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.

Template file Hide template file

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:

• sprace_ALPGEN_submit.py.txt

-- ThiagoTomei - 16 Dec 2008

-- 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:

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

-- 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

Plots Hide plots

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:

-- ThiagoTomei - 19 Dec 2008

Tables in consolidated format

ttbar + jets

W and Z + jets

-- 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