PixelEcal < Main

Tags: , view all tags
---+!! Tutorial to make Full Simulation at High Pileup Scenarios

%TOC%

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

   * Number of pixel cluster for low pileup (neutrino sample with 5000 events) <br /> <img width="996" alt="lowPU.gif" src="%ATTACHURLPATH%/lowPU.gif" height="572" />

   * Number of pixel cluster for high pileup (neutrino sample with 5000 events) <br /> <img width="996" alt="highPU.gif" src="%ATTACHURLPATH%/highPU.gif" height="572" />

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

---++ Simulated events for single electron

Full simulation with CMSSW_6_1_X using steps above. <span style="background-color: transparent;">Samples are @ T2_BR_SPRACE.</span>

| Step | Path to dataset | Number of Events |
| GEN-SIM | /pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectron_Pt5to50_14TeV_v2/ |  200K  |
| PileUP35 | /pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectron_PU35_20K_ntuple_v4/ |  20K  |
| PileUP70 | /pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectron_PU70_10K_ntuple_v4/ |  10K  |
| PileUP140 | /pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectron_PU140_5K_ntuple_v4/ |  5K  |

---++ New Simulated events for single electron

Files in the table bellow were produced with CMSSW_6_1_2_SLHC6_patch1 from these official GEN-SIM samples (50K events)
<pre>/SingleElectronFlatPt0p2To50/UpgFall13-POSTLS261_V2-v1/GEN-SIM</pre>

and using the official MinBias samples (30M events)
<pre>/MinBias_TuneZ2star_14TeV-pythia6/UpgFall13-POSTLS261_V2-v1/GEN-SIM</pre>

To open a merged file in lxplus , inside root mode, one can use
<pre><pre>TFile *file = TFile::Open("dcap://osg-se.sprace.org.br:/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/MergedNtuples/SingleElectron_Pt2to50_noPU_ntuple.root");</pre>
</pre>

Ntuples listed here are located at T2_SPRACE and contain information about both 5x5 and 1x1 EM crystals.
<table align="center" cellspacing="1" cellpadding="0" border="1"> <thead> <tr style="background-color: #1e90ff;"> <td align="center" valign="middle" scope="rowgroup"> *Scenario* </td> <td align="center" valign="middle" scope="rowgroup"> *Format* </td> <td align="center" valign="middle" scope="rowgroup"> *Size* </td> <td align="center" valign="middle" scope="rowgroup"> *Events* </td> <td align="center" valign="middle" scope="rowgroup"> *Details* </td> <td align="left" valign="middle" scope="col"> *Paths* </td> </tr> </thead> <tbody> <tr style="background-color: #66cdaa;"> <td rowspan="2" align="center" valign="middle" scope="rowgroup">No PileUP</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">DIGIRAW</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">12GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K</td> <td align="center" valign="middle" scope="rowgroup">Location</td> <td align="left" valign="middle" scope="col">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_noPileUP_50K_DIGI_RAW_v6/a49c986a608faae0f4ff11329e2bb83f/</td> </tr> <tr style="background-color: #66cdaa;"> <td align="center" valign="middle" scope="rowgroup">datasetpath</td> <td align="left" valign="middle" scope="rowgroup">/SingleElectronFlatPt0p2To50/adesouza-SingleElectron_noPileUP_50K_DIGI_RAW_v6-a49c986a608faae0f4ff11329e2bb83f/USER</td> </tr> <tr> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP35</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">DIGIRAW<br /></td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">110GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">

<span style="background-color: transparent;">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU35_50K_DIGI_RAW_v6/61196276663e836a1b830c8b84a619d4/</span>
</td> </tr> <tr> <td align="center" valign="middle" scope="rowgroup">datasetpath<br /></td> <td align="left" valign="middle" scope="rowgroup">

<span style="background-color: transparent;">/SingleElectronFlatPt0p2To50/adesouza-SingleElectron_PU35_50K_DIGI_RAW_v6-61196276663e836a1b830c8b84a619d4/USER</span>
</td> </tr> <tr style="background-color: #66cdaa;"> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP70</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">DIGIRAW</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">210GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location</td> <td align="left" valign="middle" scope="rowgroup">

<span style="background-color: transparent;">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU35_50K_DIGI_RAW_v6/61196276663e836a1b830c8b84a619d4/</span>
</td> </tr> <tr style="background-color: #66cdaa;"> <td align="center" valign="middle" scope="rowgroup">datasetpath<br /></td> <td align="left" valign="middle" scope="rowgroup">/SingleElectronFlatPt0p2To50/adesouza-SingleElectron_PU70_50K_DIGI_RAW_v6-3752475e524a60330e4a83a2952b3625/USER<br /></td> </tr> <tr> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP140</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">DIGIRAW</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">400GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">

<span style="background-color: transparent;">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU140_50K_DIGI_RAW_v6/f74d2f3945f4866a2ea25695b38d79f0/</span>
</td> </tr> <tr> <td align="center" valign="middle" scope="rowgroup">datasetpath<br /></td> <td align="left" valign="middle" scope="rowgroup">/SingleElectronFlatPt0p2To50/adesouza-SingleElectron_PU140_50K_DIGI_RAW_v6-f74d2f3945f4866a2ea25695b38d79f0/USER<br /></td> </tr> <tr style="background-color: #778899;"> <td align="center" valign="middle" scope="rowgroup"><br /></td> <td align="center" valign="middle" scope="rowgroup"><br /></td> <td align="center" valign="middle" scope="rowgroup"><br /></td> <td align="center" valign="middle" scope="rowgroup"><br /></td> <td align="center" valign="middle" scope="rowgroup"><br /></td> <td align="center" valign="middle" scope="rowgroup"><br /></td> </tr> <tr style="background-color: #87ceeb;"> <td rowspan="2" align="center" valign="middle" scope="rowgroup">No PileUP</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">Ntuple</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">24MB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_noPileUP_50K_ntuple_v6/7fef2304deb2968f78764a0244f910ee/<br /></td> </tr> <tr style="background-color: #87ceeb;"> <td align="center" valign="middle" scope="rowgroup">Merged file</td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/MergedNtuples/SingleElectron_Pt2to50_noPU_ntuple.root<br /></td> </tr> <tr> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP35</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">Ntuple</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">5.6GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU35_50K_ntuple_v6/7fef2304deb2968f78764a0244f910ee/<br /></td> </tr> <tr> <td align="center" valign="middle" scope="rowgroup">Merged file<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/MergedNtuples/SingleElectron_Pt2to50_PU35_ntuple_50K.root<br /></td> </tr> <tr style="background-color: #87ceeb;"> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP70</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">Ntuple</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">11GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU70_50K_ntuple_v6/7fef2304deb2968f78764a0244f910ee/<br /></td> </tr> <tr style="background-color: #87ceeb;"> <td align="center" valign="middle" scope="rowgroup">Merged file<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/MergedNtuples/SingleElectron_Pt2to50_PU70_ntuple_50K.root<br /></td> </tr> <tr> <td rowspan="2" align="center" valign="middle" scope="rowgroup">PileUP140</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">Ntuple</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">22GB</td> <td rowspan="2" align="center" valign="middle" scope="rowgroup">50K<br /></td> <td align="center" valign="middle" scope="rowgroup">Location<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/SingleElectron_PU140_50K_ntuple_v6/7fef2304deb2968f78764a0244f910ee/<br /></td> </tr> <tr> <td align="center" valign="middle" scope="rowgroup">Merged file<br /></td> <td align="left" valign="middle" scope="rowgroup">/pnfs/sprace.org.br/data/cms/store/user/adesouza/SingleElectronFlatPt0p2To50/MergedNtuples/SingleElectron_Pt2to50_PU140_ntuple_50K.root<br /></td> </tr> </tbody> </table>

For rejection factor
| Scenario | Format | Size | Events | Details | Paths |
| No PileUp | DIGIRAW | 560GB | 1M    | 400 files   | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/MinBias_TuneZ2star_14TeV-pythia6/MinBias_NoPU_1M_DIGI_RAW/ |
| No PileUp | DIGIRAW | 5.60TB | 9.9M | 4000 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/MinBias_TuneZ2star_14TeV-pythia6/MinBias_NoPU_10M_DIGI_RAW/ |
| PileUp35 | DIGIRAW | 440GB | 0.2M | 2000 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/Neutrino_Pt2to20_gun/Neutrino_PU35_200K_DIGI_RAW/ |
| PileUp70 | DIGIRAW | 880GB | 0.2M | 2000 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/Neutrino_Pt2to20_gun/Neutrino_PU70_200K_DIGI_RAW/ |
| No PileUp| Ntuple | 19GB | 1M | 100 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/MinBias_NoPU_1M_ntuple/ |
| No PileUp |Ntuple | 280GB | 9.9M | 500 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/MinBias_NoPU_10M_ntuple/ |
| PileUp35 | Ntuple | 30GB | 0.2M | 200 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/Neutrino_PU35_200K_ntuple/ |
| PileUp70 | Ntuple | 60GB | 0.2M | 200 files | /pnfs/sprace.org.br/data/cms/store/user/jruizvar/Neutrino_PU70_200K_ntuple/ |

Updated code available in [[https://github.com/jruizvar/pixel-analysis][GIT]].

---++ Raw tables with pixel and ecal information

To produce a raw table use the following

<verbatim>
root -b SingleElectron_Pt2to50_noPU_ntuple.root
   [0] NtupleMaker->cd() 
   [1] t->Scan("bHitLayer:bHitGx:bHitGy:bHitGz:egEt:egGx:egGy:egGz","event<=10 && bHitLayer==1"); > singleElectron_withoutPileup.log
</verbatim>

The resulting table is

* [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/singleElectron_withoutPileup.log][singleElectron_withoutPileup.log]]

| Name | Description |
| Row | Event identifier |
| Instance | Value 0 indicates no pileup |
| bHitLayer | Value 1 indicates first layer of the pixel detector |
| bHitGx | Global x position (in cm) of pixel cluster |
| bHitGy | Global y position of pixel cluster |
| bHitGz | Global z position of pixel cluster |
| egEt | Transverse energy (in GeV) of the cluster in the ECAL |
| egGx | Global x position of ecal cluster |
| egGy | Global y position of ecal cluster |
| egGz | Global z position of ecal cluster |

Adding pileup to the event, the resulting table is

* [[https://www.sprace.org.br/twiki/pub/Main/PixelEcal/singleElectron_35Pileup.log][singleElectron_35Pileup.log]]
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	16.6 K	2013-09-12 - 02:43	JoseRuiz	Analyzer
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	highPU.gif	r1	manage	130.5 K	2013-09-10 - 21:55	JoseRuiz	Number of pixel cluster for high pileup scenarios
gif	lowPU.gif	r1	manage	129.4 K	2013-09-10 - 21:56	JoseRuiz	Number of pixel cluster for low pileup scenarios:
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
log	singleElectron_35Pileup.log	r1	manage	9.6 K	2013-11-09 - 20:37	JoseRuiz	raw table with pixel and ecal information considering pileup
log	singleElectron_withoutPileup.log	r1	manage	1.5 K	2013-11-09 - 20:28	JoseRuiz	Raw table with pixel and ecal information
Topic revision: r16 - 2014-01-18 - JoseRuiz
antalya escort bursa escort eskisehir escort istanbul escort izmir escort