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

ttbar - Xsec vs number of jets - exclusive:

ttbar - Xsec vs number of jets - inclusive:

Plots - cross section as a function of parton multiplicity

Plots Hide plots

<--/twistyPlugin twikiMakeVisibleInline-->

-- ThiagoTomei - 04 Jan 2009

• 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

• • Average rate of events: 1.18 MEvents / hour.

• • • Estimate of total time @ average rate: seven days.

• 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

• Takes expected TWO MONTHS (55 days!!!) to run for W+5 jets in the (800,1600),(1600,infinity) bins. Something must be done.

• W+5 jets in the (800,1600),(1600,infinity) new estimation: eleven days (two already passed). Still not good.

• W+5 jets in the (800,1600),(1600,infinity) new estimation:
  • 73.2 MEvents already processed, of a total of 200 MEvents.
  • Current rate of events: 1.18 MEvents / hour.
  • Estimate to completion: 107.5 hours (four days and a half).
  • Estimate of total time @ current 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

Modus operandi

Final Results

Grids, input files, customized functions

Scale studies

• Absolutely NO OUTPUT from W+0 jets...

Scale studies

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

-- ThiagoTomei - 24 Dec 2008

• 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 (!!!) 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

Template file Hide template file

<--/twistyPlugin twikiMakeVisibleInline-->

<--/twistyPlugin-->

Plots - cross section as a function of parton multiplicity Hide plots

<--/twistyPlugin twikiMakeVisibleInline-->

ttbar - Xsec vs number of jets - inclusive:

<--/twistyPlugin-->

Tables in consolidated format

ttbar + jets

W and Z + jets

-- ThiagoTomei - 22 Dec 2008

• ttbar - Xsec vs number of jets - exclusive:
  

• ttbar - Xsec vs number of jets - inclusive:
  

-- ThiagoTomei - 19 Dec 2008

• 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).

-- ThiagoTomei - 18 Dec 2008

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

• sprace_ALPGEN_submit.py.txt

• sprace_CMSSW_alpgmatch.py.txt
• Alpgen_cfg.py.txt

-- ThiagoTomei - 17 Dec 2008

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

-- ThiagoTomei - 17 Dec 2008

-- ThiagoTomei - 16 Dec 2008

• ptjmin = 30, 40, 50 GeV

• Using the SPRACE cluster to generate 0,1,2,3 jets.
• 4 jets still running in access.
• sprace_ALPGEN_submit.py.txt: Script to submit Alpgen generation jobs.

-- ThiagoTomei - 16 Dec 2008

• ptjmin = 30, 40, 50 GeV

• Alpgen ttbar + N partons, N = 0,1,2,3,4
• 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.

Alpgen Production 2008

-- ThiagoTomei - 13 Dec 2008