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
)
Before submitting the job to crab using
do a local test
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
Run a local test and check the resulting n-tuple
Neutrino_PU70_ntuple.root
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)
Hint:
- 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.