(r3) PixelEcal < Main

Tags: , view all tags

---+++ 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 [[http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/Configuration/GenProduction/python/FourteenTeV/Neutrino_Pt2to20_gun_cff.py?view=markup][configuration fragment]]

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

[[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/script.sh][script.sh]] requires *Bash* terminal. 

Copy the following configure file (in =CMSSW_6_1_2_SLHC6_patch1/src=)

   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py.txt][Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py]]

Before submitting the job to crab using 

   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/crabGENSIM.cfg][crab.cfg]]

do a local test

<verbatim>
scram b -j 8
cmsRun Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py
</verbatim>

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

<verbatim>
crab -create
crab -submit -c Neutrino_Pt2to20_14TeV
</verbatim>

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

<verbatim>
crab -publish -c Neutrino_Pt2to20_14TeV
</verbatim>

---++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.

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

Same [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/script.sh][script.sh]] as above. Make a directory with a dummy name, in this case =neutrinoPU70=. Create an analyzer respecting the name ==Pxecal==
 
<verbatim>
mkdir neutrinoPU70
mkedanlzr Pxecal
cd Pxecal
</verbatim> 

Copy to this folder the following files 

   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/pxecal_cfg.py.txt][pxecal_cfg.py]]
   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/MinBiasSource.py.txt][MinBiasSource.py]]
   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/Pxecal.cc][src/Pxecal.cc]]
   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/BuildFile.xml.txt][BuildFile.xml]]
   * [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/crab.cfg][crab.cfg]]

Run a local test and check the resulting n-tuple *Neutrino_PU70_ntuple.root*

<verbatim>
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
</verbatim>

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

<verbatim>
crab -create
crab -submit 1-500 -c Neutrino_PU70_100K_ntuple
crab -submit 501-1000 -c Neutrino_PU70_100K_ntuple
</verbatim>

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

   * Distribution of the number of clusters (5000 minbias events): <br />
     <img src="%PUBURLPATH%/%WEB%/%TOPIC%/dist.gif" alt="dist.gif" width="1000" height="600" />

---++Exercise

Make the full simulation for a [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/Screen_Shot_2013-08-28_at_10.02.51_AM.png][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.

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
txt	BuildFile.xml.txt	r1	manage	0.4 K	2013-08-28 - 07:26	JoseRuiz	BuildFile
txt	MinBiasSource.py.txt	r1	manage	108.2 K	2013-08-28 - 05:29	JoseRuiz	Official MinBias
txt	Neutrino_Pt2to20_gun_cff_py_GEN_SIM.py.txt	r1	manage	4.8 K	2013-08-28 - 06:22	JoseRuiz	Configure file to make GEN-SIM
cc	Pxecal.cc	r1	manage	14.1 K	2013-08-28 - 05:31	JoseRuiz	Source code
png	Screen_Shot_2013-08-28_at_10.02.51_AM.png	r1	manage	84.1 K	2013-08-28 - 08:05	JoseRuiz	Generator fragment
cfg	crab.cfg	r1	manage	0.8 K	2013-08-28 - 05:28	JoseRuiz	Crab configure file
cfg	crabGENSIM.cfg	r1	manage	0.7 K	2013-08-28 - 06:30	JoseRuiz	Crab configure file to generate and publish GEN-SIM
gif	dist.gif	r1	manage	389.6 K	2013-09-06 - 04:15	JoseRuiz	Distribution of the number of clusters (5000 events)
txt	pxecal_cfg.py.txt	r1	manage	6.4 K	2013-08-28 - 05:29	JoseRuiz	Configure file
sh	script.sh	r1	manage	0.1 K	2013-08-28 - 05:29	JoseRuiz	Setup crab at lxplus

Topic revision: r3 - 2013-09-06 - JoseRuiz

antalya escort bursa escort eskisehir escort istanbul escort izmir escort