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
I have to run on the grid files to discover the efficiencies of producing:
The order of the grid files to run is:
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
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 |
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.
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.
Working directory in ACCESS:
~/fladias/CMSSW_3_8_6/src/ALPGEN8TeV
total_number_of_events
events_per_job
### The output files (comma separated list)
output_file=w1j_0ptw100.wgt, w1j_0ptw100.par, w1j_0ptw100.unw, w1j_0ptw100_unw.par
copy_data = 1
storage_element = T2_BR_SPRACE
user_remote_dir = FlaviaAlpgenTest
additional_input_files=alpgen.input, w1j_0ptw100.grid2
export ALPGEN_EXECUTABLE=wjetgen
-- FlaviaDias - 14 Dec 2010
antalya escort bursa escort eskisehir escort istanbul escort izmir escort