[[snowmass2013:montecarlo| <There is a bug for b-tagging in Delphes 2.03 and it is recommended to redo jet b-tagging offline for this version. Consider using Delphes 3.0 (beta). See the bottom of this page. For the pile-up simulation, look at [[http://www.snowmass2013.org/tiki-index.php?page=Energy_Frontier_FastSimulation|Energy Frontier FastSimulation page]]. Some information for ATLAS users can be found [[https://twiki.cern.ch/twiki/bin/viewauth/Atlas/LS1PileupSamples|here]]. The output ROOT files have the truth information and the detector-level objects (jets, leptons etc). ==== How to download MC samples==== Please read [[snowmass2013:download| How to download the link]] the link ==== How to analyse MCs==== First you need to install ROOT. Then the procedure is somewhat different depending on what Delphes version you are using, "Delphes 2.03" (production) or Delphes 3.0 (beta). Go to [[https://cp3.irmp.ucl.ac.be/projects/delphes| the Delphes web page]] and download the needed version. The packages contain basic examples of how to read ROOT files generated by Delphes. One example for Delphes 3.0(beta) is given in [[snowmass2013:montecarlo#delphes_30_b_and_inclusive_ttbar_sample| this section]] Let's do some analysis using the files given on this web page. Download [[https://cp3.irmp.ucl.ac.be/projects/delphes/raw-attachment/wiki/WikiStart/Delphes_V_2.0.3.tar.gz | Delphes_V_2.0.3.tar.gz]] and compile it (should work on Scientific Linux 5 or 6): wget --no-check-certificate https://cp3.irmp.ucl.ac.be/projects/delphes/raw-attachment/wiki/WikiStart/Delphes_V_2.0.3.tar.gz tar -zvxf Delphes_V_2.0.3.tar.gz cd Delphes_V_2.0.3/ make You will see this: >> Generating tmp/Utilities/ExRootAnalysis/src/BlockClassesDict.cc >> Compiling tmp/Utilities/ExRootAnalysis/src/BlockClassesDict.cc >> Generating tmp/src/TreeClassesDict.cc >> Compiling tmp/src/TreeClassesDict.cc ........................................ many messages >> Building Frog_on_analysis_output >> Compiling Examples/Convertors_Only.cpp >> Building Convertors_Only Delphes has been compiled Ready to run The directory "Examples" contains the example source code. Let's replace it with our example program. {{:snowmass2013:snowmass_tutorial.tgz|}} mv Delphes_V_2.0.3/Examples Delphes_V_2.0.3/Examples_old cd /tmp wget --no-check-certificate https://mc.hep.anl.gov/asc/wikidoc/lib/exe/fetch.php?media=snowmass2013:snowmass_tutorial.tgz tar -zvxf snowmass_tutorial.tgz mv /tmp/snowmass_tutorial/* Delphes_V_2.0.3/ # copy the content to the Delphes installation cd Delphes_V_2.0.3/ make # recompile again The new code which is inside "Examples" calculated some basic histograms for the top production. You can run it using the script "A_RUN". Open it and define the path to the directory with ROOT files. The run it as "./A_RUN". The program runs over the list of ROOT files and creates a ROOT file with histograms inside "out/" directory. ==== TTbar production at high pT ==== The generator-level settings are given below. No pile-up was included. === PYTHIA8 samples === The events were generated with this input card: [[http://mc.hep.anl.gov/asc/snowmass2013/info/analysis.ini_top | analysis.ini]]. The luminosity info and cross section for 5000 generated events can be found here: *------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------* | | | Subprocess Code | Number of events | sigma +- delta | | | Tried Selected Accepted | (estimated) (mb) | | | | | |-----------------------------------------------------------------------------------------------------------------| | | | | | g g -> t tbar 601 | 43612 2785 2785 | 2.039e-10 1.833e-12 | | q qbar -> t tbar 602 | 21140 2215 2215 | 1.616e-10 1.798e-12 | | | | | | sum | 64752 5000 5000 | 3.655e-10 2.568e-12 | | | *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------* *------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------* | | | times message | | | | 22 Error in Pythia::next: hadronLevel failed; try again | | 22 Error in StringFragmentation::fragment: stuck in joining | | 1 Warning in MultipartonInteractions::pTnext: weight above unity | | 17 Warning in SpaceShower::pT2nextQCD: weight above unity | | 17 Warning in StringFragmentation::fragmentToJunction: bad convergence junction rest frame | | | *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------* == Cross section for this run=3.655e-01 pb == Events for this run=5000 == Luminosity for this run=1.368e+01 fb-1 == Number of top-quark events =813 sigma (fb)=5.942e+01 == Number of W events =206 sigma (fb)=1.506e+01 == Number of Z events =0 sigma (fb)=0.000e+00 Get the ROOT files from [[http://mc.hep.anl.gov/asc/snowmass2013/pythia8/ttbar650pt/ | the download link]] Example for downloading 5 files: python download.py 5 pythia8/ttbar650pt pythia8_ttbar_pt650 I total there are 80 files to download. This gives 80x5000=400,000 events. === HERWIG++ samples === The log file for 5000 generated events is here: ============================================================================== Total: 155528 155528 2.31111e+08 ============================================================================== Per matrix element breakdown: MEQCD2to2Fast 155528 155528 2.31111e+08 ============================================================================== ============================================================================== Eikonalized and soft cross sections: Model parameters: ptmin: 3.5973 GeV, mu2: 0.81 GeV2 DL mode: 2, CMenergy: 14000 GeV hard inclusive cross section (mb): 230.943 soft inclusive cross section (mb): 27.5083 total cross section (mb): 113.95 inelastic cross section (mb): 78.5474 soft inv radius (GeV2): 0.306276 slope of soft pt spectrum (1/GeV2): -0.840612 Average hard multiplicity: 5.46146 Average soft multiplicity: 0.36554 ============================================================================== >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>>>>>>>> ThePEG - Toolkit for HEP Event Generation - version 1.8.1 <<<<<<<<<< <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ============================================================================== Statistics for event handler 'LHCHandler': generated number of Cross-section events attempts (nb) ============================================================================== Total (from weighted events): including vetoed events 1.246(9)e-03 ============================================================================== Total (from unweighted events): 5000 5028 1.24(2)e-03 ============================================================================== Per matrix element breakdown: MEHeavyQuark 5000 5028 1.24(2)e-03 ============================================================================== Get the ROOT files from [[http://mc.hep.anl.gov/asc/snowmass2013/herwigpp/ttbar650pt/ | the download link]] Example for downloading 5 files: python download.py 5 herwigpp/ttbar650pt herwigpp_ttbar_pt650 In total, there are 80 files to download. This gives 80×5000=400,000 events. ==== QCD (multijet) background ==== The total luminosity for inclusive QCD jet sample is roughly 9.6 fb-1 (for each sample, PYTHIA and HERWIG). The files are grouped in 2 directories and one should combine them to get ~9.6 fb-1. === PYTHIA8 samples === The events were generated with this input card: [[http://mc.hep.anl.gov/asc/snowmass2013/info/analysis.ini_top | analysis.ini]] *------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------* | | | Subprocess Code | Number of events | sigma +- delta | | | Tried Selected Accepted | (estimated) (mb) | | | | | |-----------------------------------------------------------------------------------------------------------------| | | | | | g g -> g g 111 | 5470 831 831 | 2.753e-08 4.985e-10 | | g g -> q qbar (uds) 112 | 130 24 24 | 7.833e-10 8.864e-11 | | q g -> q g 113 | 20816 2679 2679 | 9.207e-08 8.796e-10 | | q q(bar)' -> q q(bar)' 114 | 9979 1323 1323 | 4.509e-08 6.070e-10 | | q qbar -> g g 115 | 97 17 17 | 6.241e-10 7.790e-11 | | q qbar -> q' qbar' (uds) 116 | 70 19 19 | 5.035e-10 6.732e-11 | | g g -> c cbar 121 | 52 4 4 | 1.507e-10 3.342e-11 | | q qbar -> c cbar 122 | 23 7 7 | 1.925e-10 4.638e-11 | | g g -> b bbar 123 | 33 6 6 | 2.810e-10 6.515e-11 | | q qbar -> b bbar 124 | 22 8 8 | 2.282e-10 5.441e-11 | | q g -> q gamma (udscb) 201 | 17 2 2 | 1.968e-10 5.069e-11 | | q qbar -> g gamma 202 | 8 4 4 | 6.074e-11 2.082e-11 | | g g -> g gamma 203 | 0 0 0 | 0.000e+00 0.000e+00 | | f fbar -> gamma gamma 204 | 0 0 0 | 0.000e+00 0.000e+00 | | g g -> gamma gamma 205 | 0 0 0 | 0.000e+00 0.000e+00 | | f fbar -> gamma*/Z0 221 | 59 11 11 | 4.135e-10 7.896e-11 | | f fbar' -> W+- 222 | 107 31 31 | 8.721e-10 9.339e-11 | | f fbar -> gamma*/Z0 gamma*/Z0 231 | 4 1 1 | 6.196e-12 6.196e-12 | | f fbar' -> Z0 W+- (no gamma*!) 232 | 2 0 0 | 0.000e+00 0.000e+00 | | f fbar -> W+ W- 233 | 4 0 0 | 0.000e+00 0.000e+00 | | q qbar -> gamma*/Z0 g 241 | 19 2 2 | 6.789e-11 3.454e-11 | | q g-> gamma*/Z0 q 242 | 101 11 11 | 2.479e-10 5.253e-11 | | f fbar -> gamma*/Z0 gamma 243 | 0 0 0 | 0.000e+00 0.000e+00 | | q qbar' -> W+- g 251 | 35 3 3 | 1.787e-10 3.979e-11 | | q g-> W+- q' 252 | 118 17 17 | 5.763e-10 6.943e-11 | | f fbar' -> W+- gamma 253 | 2 0 0 | 0.000e+00 0.000e+00 | | | | | | sum | 37168 5000 5000 | 1.701e-07 1.203e-09 | | | *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------* == Cross section for this run=1.701e+02 pb == Events for this run=5000 == Luminosity for this run=2.940e-02 fb-1 == Number of top-quark events =0 sigma (fb)=0.000e+00 == Number of W events =2 sigma (fb)=6.803e+01 == Number of Z events =3 sigma (fb)=1.020e+02 Get the ROOT files grouped in 2 directories: * [[http://mc.hep.anl.gov/asc/snowmass2013/pythia8/qcd650pt/run1 | the download link 1]] * [[http://mc.hep.anl.gov/asc/snowmass2013/pythia8/qcd650pt/run2 | the download link 2]] Example for downloading 5 files: python download.py 5 pythia8/qcd650pt/run1 pythia8_qcd_pt650 python download.py 5 pythia8/qcd650pt/run2 pythia8_qcd_pt650_2 === HERWIG++ samples === The MC samples were generated for inclusive jets events. The part of the log file with the cross sections is shown below (for 5000 events): Miscellaneous output from modules to the standard output: Statistics for the UE process: generated number of Cross-section events attempts (nb) ============================================================================== Total: 153340 153340 2.30923e+08 ============================================================================== Per matrix element breakdown: MEQCD2to2Fast 153340 153340 2.30923e+08 ============================================================================== ============================================================================== Eikonalized and soft cross sections: Model parameters: ptmin: 3.5973 GeV, mu2: 0.81 GeV2 DL mode: 2, CMenergy: 14000 GeV hard inclusive cross section (mb): 231.373 soft inclusive cross section (mb): 27.2365 total cross section (mb): 113.95 inelastic cross section (mb): 78.5401 soft inv radius (GeV2): 0.304806 slope of soft pt spectrum (1/GeV2): -0.884655 Average hard multiplicity: 5.47163 Average soft multiplicity: 0.360742 ============================================================================== >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>>>>>>>> ThePEG - Toolkit for HEP Event Generation - version 1.8.1 <<<<<<<<<< <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ============================================================================== Statistics for event handler 'LHCHandler': generated number of Cross-section events attempts (nb) ============================================================================== Total (from weighted events): including vetoed events 0.168(1)e+00 ============================================================================== Total (from unweighted events): 5000 5038 0.167(3)e+00 ============================================================================== Per matrix element breakdown: MEQCD2to2 4989 5027 0.166(3)e+00 MEGammaJet 4 4 0.13(7)e-03 MEZJet 1 1 30(30)e-06 MEWJet 6 6 0.20(8)e-03 ============================================================================== Get the ROOT files from these 2 locations: * [[http://mc.hep.anl.gov/asc/snowmass2013/herwigpp/qcd650pt/run1 | the download link 1]] * [[http://mc.hep.anl.gov/asc/snowmass2013/herwigpp/qcd650pt/run2 | the download link 2]] Example for downloading 5 files: python download.py 5 herwigpp/qcd650pt/run1 herwigpp_qcd_pt650 python download.py 5 herwigpp/qcd650pt/run2 herwigpp_qcd_pt650_2 --- //[[chekanov@anl.gov|Sergei Chekanov]] 2013/01/31 08:52//