Tutorial to make Full Simulation at High Pileup Scenarios

This tutorial illustrates the simulation of neutrino samples with pileup 70 considering upgraded geometry for the pixel detector (4-layers, 3-disks).

Generation of events

We are going to simulate Neutrinos with Pt 2-20 according to this configuration fragment

ssh jruizvar@lxplus5.cern.ch
cmsrel CMSSW_6_1_2_SLHC6_patch1
cd CMSSW_6_1_2_SLHC6_patch1/src
source script.sh

script.sh requires Bash terminal.

Copy the following configure file (in CMSSW_6_1_2_SLHC6_patch1/src)

• Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py

• crab.cfg

scram b -j 8
cmsRun Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py

Check edmFileUtil Neutrino_Pt2to20_gun_cff_py_GEN_SIM.root. This file is only to test and can be removed ( rm Neutrino_Pt2to20_gun_cff_py_GEN_SIM.root)

Submit to crab

crab -create
crab -submit -c Neutrino_Pt2to20_14TeV

It takes around 15 min to generate 10000 events. 0.2M events can be generated in 20 jobs. After crab getoutput, don't forget

crab -publish -c Neutrino_Pt2to20_14TeV

Create Analyzer

We are going to create an standalone analyzer (no need to checkout external libraries) that can be used with any GEN-SIM sample, not only neutrinos.

ssh jruizvar@lxplus5.cern.ch
cmsrel CMSSW_6_1_2_SLHC6_patch1
cd CMSSW_6_1_2_SLHC6_patch1/src
source script.sh

Same script.sh as above. Make a directory with a dummy name, in this case neutrinoPU70. Create an analyzer respecting the name Pxecal

mkdir neutrinoPU70
mkedanlzr Pxecal
cd Pxecal

Copy to this folder the following files

• pxecal_cfg.py
• MinBiasSource.py
• src/Pxecal.cc
• BuildFile.xml
• crab.cfg

scram b -j 8
cmsRun pxecal_cfg.py
root -l Neutrino_PU70_ntuple.root
   [0] NtupleMaker->cd()
   [1] t->Draw("sqrt(ClGx*ClGx+ClGy*ClGy):ClGz")
   [2] .q
rm Neutrino_PU70_ntuple.root

If local run succeeds, you can create and submit the jobs

crab -create
crab -submit 1-500 -c Neutrino_PU70_100K_ntuple
crab -submit 501-1000 -c Neutrino_PU70_100K_ntuple

This step is more time consuming than the generation part. No need to publish because the output are n-tuples stored at T2_BR_SPRACE for private analysis.

Example plots

• Number of pixel cluster for low pileup (neutrino sample with 5000 events)
  

• Number of pixel cluster for high pileup (neutrino sample with 5000 events)
  

Exercise

Make the full simulation for a SingleElectronFlatPt5To50 for three pileup scenarios:

• PU35 (200K events)
• PU70 (100K events)
• PU140 (50K events)

• The GEN-SIM step should be done only once.
• The python configure files should be slightly modified.
• No need to touch the Pxecal.cc source file.

Simulated events for single electron

Full simulation with CMSSW_6_1_X using steps above. Samples are @ T2_BR_SPRACE.