ALPGEN 8TeV CMS Production
Grid Files
The grid files, produced by Alex, are in his public home in afs:
~amott/public/ALPGEN/8TeV/grids/
And in my personal afs area:
~fladias/ALPGEN8TeVGrids
Unweighting and Matching Efficiency of Grid files
I have to run on the grid files to discover the efficiencies of producing:
- Unweighted events, given a number of weighted events to produce;
- Showered events, after matching procedure on PYTHIA, so we can deliver the correct amount of .lhe files.
The order of the grid files to run is:
- W + jets
- Z + jets
- ttbar + jets
- QCD
Grid Files Production
I've used Alex's script to produce all the input cards with 8TeV center of mass energy, and I am running to get the grid files according to:
Number of warm-ups = 2 + number of jets
Number of iteration/warm-up = 1.000.000
Number of events = 2 x (2+njets) x 1.000.000
Efficiencies Calculation
List of Grid Efficiencies
For calculating the grid efficiencies, have to run in all the process (through PYTHIA and so on) to get the matching efficiency. I also have to know the cross section, so I can calculate the number of events needed for a given luminosity.
I'm running 1.000.000 .wgt events to get the .unw and matching efficiency
- First, have to run on imode=1, to get the .wgt
- Then, run on imode=2, to get the .unw
Grid |
.wgt eff |
.unw eff (from .wgt) |
Matching eff (from .unw) |
Matching eff (from submitted) |
Cross section (pb) |
Events Generated |
w0j.grid2 |
0.982703 |
0.723857 |
0.859 |
calculate |
2.441e+4 +/- 1.044e+2 |
|
w0j.grid2 NEW seed |
0.504036 |
0.638775 |
0.863 |
calculate |
2.450e+4 +/- 7.005e+1 |
|
w1j_0ptw100.grid2 |
0.611141 |
0.20049 |
0.336 |
0.041180 |
3.442e+3 +/- 1.469e+1 |
|
w1j_100ptw300.grid2 |
0.424175 |
0.099343431 |
0.539 |
calculate |
1.098e+2 +/- 5.497e-1 |
|
w1j_300ptw800.grid2 |
0.128871 |
0.189328864 |
0.37 |
calculate |
9.147e-1 +/- 8.160e-3 |
|
w1j_800ptw1600.grid2 |
0.127199 |
0.132736893 |
0.228 |
0.003854 |
2.264e-06 +/- 3.581e-08 |
|
w2j_0ptw100.grid2 |
0.0284464 |
0.41629 |
0.18 |
calculate |
7.952e+2 +/- 7.498e+0 |
|
w2j_100ptw300.grid2 |
0.0158608 |
0.375 |
0.289 |
calculate |
8.401e+01 +/- 8.057e-01 |
|
w2j_300ptw800.grid2 |
0.0290226 |
0.09854 |
0.227 |
calculate |
1.168e+00 +/- 1.849e-02 |
|
w2j_800ptw1600.grid2 |
0.01706 |
0.065283 |
0.18 |
calculate |
5.196e-06 +/- 1.512e-07 |
|
w3j_0ptw100.grid2 |
0.008296 |
0.1995 |
0.539 |
calculate |
1.098e+02 +/- 5.497e-01 |
|
w3j_100ptw300.grid2 |
0.00879273 |
0.120774 |
0.15 |
0.00015933 |
3.432e+01 +/- 6.595e-01 |
|
w3j_300ptw800.grid2 |
0.008296 |
0.07119 |
0.136 |
0.0000804 |
7.600e-01 +/- 2.058e-02 |
|
w3j_800ptw1600.grid2 |
0.00372 |
0.04174 |
0.128 |
2e-5 |
4.751e-06 +/- 2.592e-07 |
|
w4j_0ptw100.grid2 |
0.00344 |
0.129 |
0.05 |
0.0000224 |
3.017e+01 +/- 1.252e+00 |
|
w4j_100ptw300.grid2 |
0.003 |
0.1144 |
0.081 |
0.00002776 |
1.023e+01 +/- 3.747e-01 |
|
w4j_300ptw800.grid2 |
0.00489 |
0.0418 |
0.086 |
0.00001768 |
3.500e-01 +/- 1.598e-02 |
|
w4j_800ptw1600.grid2 |
0.001008 |
0.012817 |
0.077 |
0.000001 |
2.159e-06 +/- 4.163e-07 |
|
w5j_0ptw100.grid2 |
0.0006288 |
0.068 |
0.035 |
0.00000148 |
6.603e+00 +/- 7.640e-01 |
|
w5j_100ptw300.grid2 |
0.001941 |
0.05814 |
0.039 |
0.00000442 |
2.203e+00 +/- 1.458e-01 |
|
w5j_300ptw800.grid2 |
0.001620 |
0.04792 |
0.047 |
0.00000368 |
1.064e-01 +/- 5.422e-03 |
|
w5j_800ptw1600.grid2 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
|
Massive Production
Initially, the aim is to produce samples equivalent to 1fb^-1 of data.
Number of matched events needed: N = sigma . L
Where sigma is the cross section for the given PROCESS, and L is the accumulated luminosity.
So, for 1fb^-1, I have to produce the number of events equal to the cross section given for the process, in
fb.
Step by step plan - Working with crab, efficiencies evaluation, etc.
I've runned the test for Alpgen on crab in my CERN afs area. Need to do the same steps in my own personal area in ACCESS from SPRACE.
- Copy the files needed (from the tarball given in this Twiki) to my local ACCESS area.
Working directory in ACCESS:
~/fladias/CMSSW_3_8_6/src/ALPGEN8TeV
- What I need to do for running with this package:
- crab_genAlpgen_grid.cfg - From the CRAB configuration file, I have to change:
-
total_number_of_events
-
events_per_job
- Output files
### The output files (comma separated list)
output_file=w1j_0ptw100.wgt, w1j_0ptw100.par, w1j_0ptw100.unw, w1j_0ptw100_unw.par
- Set the storage element and copy mode:
copy_data = 1
storage_element = T2_BR_SPRACE
user_remote_dir = FlaviaAlpgenTest
- Set the input files:
additional_input_files=alpgen.input, w1j_0ptw100.grid2
- alpgen.input - Take care to change the input card.
- runAlpgen.sh - Change the name of the executable in the line:
export ALPGEN_EXECUTABLE=wjetgen
Test
Run Alpgen with a few (60.000) events for cross-section/efficiency test.
Process chosen: w0j_0ptw100
TASK - Script to change these parameters in the script to run Alpgen for production.
Previous Alpgen Production Page
AlpgenProduction7TeV
Crab Exit Errors
Job Exit Codes
--
FlaviaDias - 14 Dec 2010