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)Sorted ascending Events Generated
w1j_100ptw300.grid2 0.424175 0.099343431 0.539 calculate 1.098e+2 +/- 5.497e-1  
w3j_0ptw100.grid2 0.008296 0.1995 0.539 calculate 1.098e+02 +/- 5.497e-01  
w4j_100ptw300.grid2 0.003 0.1144 0.081 0.00002776 1.023e+01 +/- 3.747e-01  
w5j_300ptw800.grid2 0.001620 0.04792 0.047 0.00000368 1.064e-01 +/- 5.422e-03  
w2j_300ptw800.grid2 0.0290226 0.09854 0.227 calculate 1.168e+00 +/- 1.849e-02  
w4j_800ptw1600.grid2 0.001008 0.012817 0.077 0.000001 2.159e-06 +/- 4.163e-07  
w1j_800ptw1600.grid2 0.127199 0.132736893 0.228 0.003854 2.264e-06 +/- 3.581e-08  
w5j_100ptw300.grid2 0.001941 0.05814 0.039 0.00000442 2.203e+00 +/- 1.458e-01  
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  
w4j_0ptw100.grid2 0.00344 0.129 0.05 0.0000224 3.017e+01 +/- 1.252e+00  
w1j_0ptw100.grid2 0.611141 0.20049 0.336 0.041180 3.442e+3 +/- 1.469e+1  
w3j_100ptw300.grid2 0.00879273 0.120774 0.15 0.00015933 3.432e+01 +/- 6.595e-01  
w4j_300ptw800.grid2 0.00489 0.0418 0.086 0.00001768 3.500e-01 +/- 1.598e-02  
w3j_800ptw1600.grid2 0.00372 0.04174 0.128 2e-5 4.751e-06 +/- 2.592e-07  
w2j_800ptw1600.grid2 0.01706 0.065283 0.18 calculate 5.196e-06 +/- 1.512e-07  
w5j_0ptw100.grid2 0.0006288 0.068 0.035 0.00000148 6.603e+00 +/- 7.640e-01  
w3j_300ptw800.grid2 0.008296 0.07119 0.136 0.0000804 7.600e-01 +/- 2.058e-02  
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  
w1j_300ptw800.grid2 0.128871 0.189328864 0.37 calculate 9.147e-1 +/- 8.160e-3  

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

Topic revision: r23 - 2014-06-17 - trtomei
 

This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback

antalya escort bursa escort eskisehir escort istanbul escort izmir escort