(r20) AlpgenWorkLog2 < Main

Tags: , view all tags
---+ 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)* |
|  =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 [[https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideCMSAlpgen#Using_CRAB_to_produce_ALPGEN_fil][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 <br>
         =### The output files (comma separated list)= <br>
         =output_file=w1j_0ptw100.wgt, w1j_0ptw100.par, w1j_0ptw100.unw, w1j_0ptw100_unw.par= <br>
      * Set the storage element and copy mode: <br>
         =copy_data = 1= <br>
         =storage_element = T2_BR_SPRACE= <br>
         =user_remote_dir = FlaviaAlpgenTest= <br>
      * Set the input files: <br>
         =additional_input_files=alpgen.input, w1j_0ptw100.grid2= <br>

   * *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
[[http://www.sprace.org.br/Twiki/bin/view/Main/AlpgenEventProductionPage][AlpgenProduction7TeV]]

---+++Crab Exit Errors
[[https://twiki.cern.ch/twiki/bin/view/CMS/JobExitCodes][Job Exit Codes]]

 -- Main.FlaviaDias - 14 Dec 2010
Topic revision: r20 - 2010-12-20 - FlaviaDias
antalya escort bursa escort eskisehir escort istanbul escort izmir escort