Alpgen Production 2008

Modus operandi

  • Alpgen ttbar + N partons, N = 0,1,2,3,4
  • ptjmin = 30, 40, 50 GeV
  • Library: 2Qwork
  • Production of grid2 files.
  • Maria's rule of thumb: 5+N iterations, where N is the number of light (udsg) partons.
  • ~ 43000 events per minute.

1          ! generation mode
ttbar4j_40 ! string labeling the output files
0          ! start with: 0=new grid, 1=previous warmup grid, 2=previous generation grid
4000000  5 ! N(events)/iteration and N(iter's) for initial grid optimiz
50000000   ! number evts to generate
*** The above 5 lines provide mandatory inputs for all processes
*** (Comment lines are introduced by the three asteriscs)
*** The lines below modify existing defaults for the hard process under study
*** For a complete list of accessible parameters and their values,
*** input 'print 1' (to display on the screen) or 'print 2' to write to file
njets 4
ebeam 5000
ih2 1
ickkw 1
ptjmin 40
etajmax 5
drjmin 0.7
mt 175
ihvy 6
itdecmode 7

  • Tried only 5 iterations for ttbar + 4 jets, trying now with 9 iterations and seeing if it makes a difference.

-- ThiagoTomei - 15 Dec 2008

  • Using the SPRACE cluster to generate 0,1,2,3 jets.
  • 4 jets still running in access.

Script to submit Alpgen generation jobs:

-- ThiagoTomei - 16 Dec 2008

Sample Events / cycle Cycles Asked events WGT events WGT Xsec (pb) Error (pb) UNW events UNW Xsec (pb) Error (pb) Unweighting efficiency Matched events Matching Xsec (pb) Error (pb) Matched efficiency
ttbar0j_30GeV 4000000 5 50000000 6462816 2.02E+03 1.06E+01 1136665 1.90E+02 9.47E-02 2.27E-002 702390 1.18E+02 0.618
ttbar0j_40GeV 4000000 5 50000000 6462816 2.02E+03 1.06E+01 1136665 1.90E+02 9.47E-02 2.27E-002 795122 1.33E+02 0.700
ttbar0j_50GeV 4000000 5 50000000 6462816 2.02E+03 1.06E+01 1136665 1.90E+02 9.47E-02 2.27E-002 860824 1.44E+02 0.757
ttbar1j_30GeV 4000000 6 50000000 2120429 1.53E+03 9.33E+00 424774 1.64E+02 1.21E-01 8.50E-003 159965 6.18E+01 0.377
ttbar1j_40GeV 4000000 6 50000000 2350711 1.17E+03 9.61E+00 453316 1.21E+02 8.66E-02 9.07E-003 194990 5.21E+01 0.430
ttbar1j_50GeV 4000000 6 50000000 2585228 9.08E+02 9.70E+00 506646 9.37E+01 6.56E-02 1.01E-002 236252 4.37E+01 0.466
ttbar2j_30GeV 4000000 7 50000000 190817 7.52E+02 7.57E+00 35558 8.73E+01 2.03E-01 7.11E-004 8410 2.06E+01 0.237
ttbar2j_40GeV 4000000 7 50000000 237619 4.35E+02 7.89E+00 46028 4.94E+01 9.63E-02 9.21E-004 12701 1.36E+01 0.276
ttbar2j_50GeV 4000000 7 50000000 298591 2.75E+02 8.07E+00 59939 3.07E+01 5.29E-02 1.20E-003 18102 9.27E+00 0.302
ttbar3j_30GeV 4000000 8 50000000 262729 4.15E+02 8.65E+00 39244 3.49E+01 6.64E-02 7.85E-004 5835 5.19E+00 0.149
ttbar3j_40GeV 4000000 8 50000000 376197 1.88E+02 9.15E+00 48051 1.54E+01 2.58E-02 9.61E-004 8515 2.74E+00 0.177
ttbar3j_50GeV 4000000 8 50000000 339983 1.00E+02 9.82E+00 56262 7.87E+00 1.26E-02 1.13E-003 11006 1.54E+00 0.196
ttbar4j_30GeV 4000000 9 50000000 208188 2.02E+02 1.03E+01 20263 1.15E+01 2.53E-02 4.05E-004 2760 1.56E+00 0.136
ttbar4j_40GeV 4000000 9 50000000 155193 7.33E+01 1.13E+01 25736 4.04E+00 8.59E-03 5.15E-004 3854 6.05E-01 0.150
ttbar4j_50GeV 4000000 9 50000000 180885 3.09E+01 1.14E+01 30121 1.71E+00 3.56E-03 6.02E-004 4854 2.75E-01 0.161

-- ThiagoTomei - 17 Dec 2008

Grids done. Running showering and matching now. Efficiencies seem to come down a little bit too low. The matching is done with:

  • ETCLUS = 1.2*ptjmin, i.e., for ptjmin = (30,40,50) GeV we have ETCLUS = (36,48,60) GeV, as per the ALPGEN defaults.
  • DR = 0.7 (fixed).
  • ETACLUS = 5.0 (fixed).

Scripts for showering and matching:

-- ThiagoTomei - 17 Dec 2008

All finished. The AlpgenProduction2008_ttbarjets.tar.gz tarball contains a directory for each sample, and in each directory we find the .input file, used to generate the grids, and the .grid2 file.

-- ThiagoTomei - 18 Dec 2008

Final Results

Grids, input files, customized functions

In the attachments below you can find the grids for ttbar, W and Z - ttbar was done by Thiago, W and Z were done by Maurizio.

Name Description
AlpgenProduction2008_ttbarjets.tar.gz Grids and input files for ttbar + jets, with (30,40,50) GeV parton thresholds
AlpgenProduction2008_Wjets.tar.gz Grids, input files and custom functions for W + jets, in W pT bins
AlpgenProduction2008_Zjets.tar.gz Grids, input files and custom functions for Z + jets, in Z pT bins

Plots - cross section as a function of parton multiplicity

NOTICE : slopes were fitted in the bins 1 to 4 partons only.

ttbar - Xsec vs number of jets - exclusive:

ttbar_exclusive.png

ttbar - Xsec vs number of jets - inclusive:

ttbar_inclusive.png

-- ThiagoTomei - 19 Dec 2008

Tables in consolidated format

ttbar + jets

Sample name Exc xsec (pb) Error Inc Xsec (pb) Error
ttbar0j_30GeV 117.529 0.059 206.693 0.407
ttbar1j_30GeV 61.770 0.046 89.165 0.327
ttbar2j_30GeV 20.639 0.048 27.395 0.248
ttbar3j_30GeV 5.191 0.010 6.755 0.115
ttbar4j_30GeV 1.564 0.003 1.564 0.059
ttbar0j_40GeV 133.045 0.066 202.070 0.369
ttbar1j_40GeV 52.052 0.037 69.025 0.264
ttbar2j_40GeV 13.631 0.027 16.973 0.180
ttbar3j_40GeV 2.737 0.005 3.342 0.077
ttbar4j_40GeV 0.605 0.001 0.605 0.036
ttbar0j_50GeV 144.039 0.072 198.817 0.348
ttbar1j_50GeV 43.692 0.031 54.778 0.223
ttbar2j_50GeV 9.270 0.016 11.086 0.138
ttbar3j_50GeV 1.541 0.003 1.815 0.055
ttbar4j_50GeV 0.275 0.001 0.275 0.024

W and Z + jets

W+jets xsec excl error xsec incl error
1j 30538.000 3.058 39321.200 590.948
2j 6373.580 20.461 8783.190 131.997
3j 1796.860 12.516 2409.610 35.377
4j 469.031 3.463 612.751 7.919
5j 110.275 0.809 143.720 1.525
6j 33.445 0.256 33.445 0.256
Z+jets xsec excl error xsec incl error
1j 2950.500 1.353 3872.020 53.125
2j 656.510 1.494 921.517 12.636
3j 192.907 1.205 265.008 3.584
4j 51.211 0.305 72.101 0.865
5j 17.240 0.107 20.890 0.217
6j 3.650 0.031 3.650 0.031

-- ThiagoTomei - 22 Dec 2008

Scale studies

  • Started today with W+jets, iqopt=4, based on Maurizio's code.
  • Takes a little over one day to run for W+1,2,3,4 jets.
  • Takes a little over two days and a half to run for W+5 jets in the (0,100),(100,300),(300,800) bins.
  • Takes expected TWO MONTHS (55 days!!!) to run for W+5 jets in the (800,1600),(1600,infinity) bins. Something must be done.
  • Absolutely NO OUTPUT from W+0 jets...

-- ThiagoTomei - 23 Dec 2008

  • W+5 jets in the (0,100),(100,300),(300,800) bins finished today
  • W+5 jets in the (800,1600),(1600,infinity) new estimation: eleven days (two already passed). Still not good.
  • NO OUTPUT from W+0 jets...

-- ThiagoTomei - 24 Dec 2008

  • W+5 jets in the (800,1600),(1600,infinity) new estimation:
    • 73.2 MEvents already processed, of a total of 200 MEvents.
    • Average rate of events: 1.18 MEvents / hour.
      • Estimate to completion: 107.5 hours (four days and a half).
      • Estimate of total time @ average rate: seven days.
      • Estimate changing over time: 55 days -> 11 days -> 7 days.
      • GUESS: I believe the estimate of total time would freeze at 6.6 days = 6 days, 14 hours.
        • This puts the completion time @ Mon Dec 29, 2 PM
        • BUT, there is a power cut expected @ Sun, Dec 28 ... will we make it?
  • NO OUTPUT from W+0 jets...

-- ThiagoTomei - 25 Dec 2008

  • W+5 jets in the (800,1600),(1600,infinity) new estimation:
    • 139.1 MEvents already processed, of a total of 200 MEvents.
    • Average rate of events: 1.53 MEvents / hour.
      • Estimate to completion: 40 hours. Completion time @ Sun, Dec 28, 10 AM
      • Estimate of total time @ average rate: five days and a half.
      • Estimate changing over time: 55 days -> 11 days -> 7 days -> 5.5 days.
    • Current rate of events: 2.32 MEvents / hour.
      • Estimate to completion: 26 hours. Completion time @ Sat, Dec 27, 8 PM
      • This puts the completion time BEFORE the power cut. Hooray!!!
  • NO OUTPUT from W+0 jets...

-- ThiagoTomei - 26 Dec 2008

  • W+5 jets in the (800,1600),(1600,infinity) new estimation:
    • 177.4 MEvents already processed, of a total of 200 MEvents.
    • Current rate of events: 2.25 MEvents / hour.
      • Estimate to completion: 10 hours. Completion time @ Sat, Dec 27, 8 PM
  • NO OUTPUT from W+0 jets...

-- ThiagoTomei - 27 Dec 2008

Continues in AlpgenProduction2009.

-- ThiagoTomei - 04 Jan 2009

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng ttbar_exclusive.png r2 r1 manage 16.5 K 2008-12-22 - 11:47 ThiagoTomei ttbar - Xsec vs number of jets - exclusive
PNGpng ttbar_inclusive.png r2 r1 manage 17.1 K 2008-12-22 - 11:48 ThiagoTomei ttbar - Xsec vs number of jets - inclusive
Topic revision: r24 - 2009-01-12 - ThiagoTomei
 

This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback

antalya escort bursa escort eskisehir escort istanbul escort izmir escort