Fork me on GitHub

source: svn/trunk/examples/delphes_card_CMS_PileUp_Timing.tcl@ 1348

Last change on this file since 1348 was 1345, checked in by Michele Selvaggi, 11 years ago

timing implemented
File size: 9.0 KB
Line 
1set MaxEvents 1000
2
3#######################################
4# Order of execution of various modules
5#######################################
6
7set ExecutionPath {
8
9 PileUpMerger
10 ParticlePropagator
11
12 ChargedHadronTrackingEfficiency
13 ElectronTrackingEfficiency
14 MuonTrackingEfficiency
15
16 ChargedHadronMomentumSmearing
17 ElectronEnergySmearing
18 MuonMomentumSmearing
19
20 MuonTimeSmearing
21
22 TreeWriter
23}
24
25
26###############
27# PileUp Merger
28###############
29
30module PileUpMerger PileUpMerger {
31 set InputArray Delphes/stableParticles
32
33 set ParticleOutputArray stableParticles
34 set VertexOutputArray vertices
35
36 # pre-generated minbias input file
37 set PileUpFile MinBias.pileup
38
39 # average expected pile up
40 set MeanPileUp 20
41
42 # maximum spread in the beam direction in m
43 set ZVertexSpread 0.05
44
45 # maximum spread in time in s
46 set TVertexSpread 0.5E-09
47
48 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
49
50 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
51
52 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
53 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
54 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \
55 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
56
57}
58
59#################################
60# Propagate particles in cylinder
61#################################
62
63module ParticlePropagator ParticlePropagator {
64 set InputArray PileUpMerger/stableParticles
65
66 set OutputArray stableParticles
67 set ChargedHadronOutputArray chargedHadrons
68 set ElectronOutputArray electrons
69 set MuonOutputArray muons
70
71 # radius of the magnetic field coverage, in m
72 set Radius 1.29
73 # half-length of the magnetic field coverage, in m
74 set HalfLength 3.00
75
76 # magnetic field
77 set Bz 3.8
78}
79
80####################################
81# Charged hadron tracking efficiency
82####################################
83
84module Efficiency ChargedHadronTrackingEfficiency {
85 set InputArray ParticlePropagator/chargedHadrons
86 set OutputArray chargedHadrons
87
88 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
89
90 # tracking efficiency formula for charged hadrons
91 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
92 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
93 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
94 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
95 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
96 (abs(eta) > 2.5) * (0.00)}
97}
98
99##############################
100# Electron tracking efficiency
101##############################
102
103module Efficiency ElectronTrackingEfficiency {
104 set InputArray ParticlePropagator/electrons
105 set OutputArray electrons
106
107 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
108
109 # tracking efficiency formula for electrons
110 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
111 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
112 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
113 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
114 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
115 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
117 (abs(eta) > 2.5) * (0.00)}
118}
119
120##########################
121# Muon tracking efficiency
122##########################
123
124module Efficiency MuonTrackingEfficiency {
125 set InputArray ParticlePropagator/muons
126 set OutputArray muons
127
128 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
129
130 # tracking efficiency formula for muons
131 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
132 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
133 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
134 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
135 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
136 (abs(eta) > 2.5) * (0.00)}
137}
138
139########################################
140# Momentum resolution for charged tracks
141########################################
142
143module MomentumSmearing ChargedHadronMomentumSmearing {
144 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
145 set OutputArray chargedHadrons
146
147 # set ResolutionFormula {resolution formula as a function of eta and pt}
148
149 # resolution formula for charged hadrons
150 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
151 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
152 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
153 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
154 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
155 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
156 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
157 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
158}
159
160#################################
161# Energy resolution for electrons
162#################################
163
164module EnergySmearing ElectronEnergySmearing {
165 set InputArray ElectronTrackingEfficiency/electrons
166 set OutputArray electrons
167
168 # set ResolutionFormula {resolution formula as a function of eta and energy}
169
170 # resolution formula for electrons
171 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
172 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
173 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
174 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
175
176}
177
178###############################
179# Momentum resolution for muons
180###############################
181
182module MomentumSmearing MuonMomentumSmearing {
183 set InputArray MuonTrackingEfficiency/muons
184 set OutputArray muons
185
186 # set ResolutionFormula {resolution formula as a function of eta and pt}
187
188 # resolution formula for muons
189 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
190 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
191 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
192 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
193 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
194 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
195 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
196 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
197 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
198 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
199 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
200 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
201}
202
203
204##############
205# Muon Timing
206##############
207
208module TimeSmearing MuonTimeSmearing {
209 set InputArray ParticlePropagator/muons
210 set OutputArray muons
211
212 set TimeResolution 1.0e-10
213
214}
215
216##################
217# ROOT tree writer
218##################
219
220module TreeWriter TreeWriter {
221# add Branch InputArray BranchName BranchClass
222 add Branch Delphes/allParticles Particle GenParticle
223 add Branch MuonTimeSmearing/muons MuonTimeSmeared Muon
224 add Branch PileUpMerger/vertices Vertex Vertex
225}
226
Note: See TracBrowser for help on using the repository browser.