Fork me on GitHub

source: git/cards/delphes_card_CMS_NoFastJet.tcl

Last change on this file was d091310, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

replace calo by simple calo in CMS cards
  • Property mode set to 100644
File size: 12.7 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23
24}
25
26#################################
27# Propagate particles in cylinder
28#################################
29
30module ParticlePropagator ParticlePropagator {
31 set InputArray Delphes/stableParticles
32
33 set OutputArray stableParticles
34 set ChargedHadronOutputArray chargedHadrons
35 set ElectronOutputArray electrons
36 set MuonOutputArray muons
37
38 # radius of the magnetic field coverage, in m
39 set Radius 1.29
40 # half-length of the magnetic field coverage, in m
41 set HalfLength 3.00
42
43 # magnetic field
44 set Bz 3.8
45}
46
47####################################
48# Charged hadron tracking efficiency
49####################################
50
51module Efficiency ChargedHadronTrackingEfficiency {
52 set InputArray ParticlePropagator/chargedHadrons
53 set OutputArray chargedHadrons
54
55 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
56
57 # tracking efficiency formula for charged hadrons
58 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
59 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
60 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
61 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
62 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
63 (abs(eta) > 2.5) * (0.00)}
64}
65
66##############################
67# Electron tracking efficiency
68##############################
69
70module Efficiency ElectronTrackingEfficiency {
71 set InputArray ParticlePropagator/electrons
72 set OutputArray electrons
73
74 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
75
76 # tracking efficiency formula for electrons
77 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
78 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
79 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
80 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
81 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
83 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
84 (abs(eta) > 2.5) * (0.00)}
85}
86
87##########################
88# Muon tracking efficiency
89##########################
90
91module Efficiency MuonTrackingEfficiency {
92 set InputArray ParticlePropagator/muons
93 set OutputArray muons
94
95 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
96
97 # tracking efficiency formula for muons
98 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
100 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
101 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
102
103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
105 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
106 (abs(eta) > 2.5) * (0.00)}
107
108}
109
110########################################
111# Momentum resolution for charged tracks
112########################################
113
114module MomentumSmearing ChargedHadronMomentumSmearing {
115 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
116 set OutputArray chargedHadrons
117
118 # set ResolutionFormula {resolution formula as a function of eta and pt}
119
120 # resolution formula for charged hadrons
121 # based on arXiv:1405.6569
122 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
123 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
124 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
125}
126
127###################################
128# Momentum resolution for electrons
129###################################
130
131module MomentumSmearing ElectronMomentumSmearing {
132 set InputArray ElectronTrackingEfficiency/electrons
133 set OutputArray electrons
134
135 # set ResolutionFormula {resolution formula as a function of eta and energy}
136
137 # resolution formula for electrons
138 # based on arXiv:1405.6569
139 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
140 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
142}
143
144###############################
145# Momentum resolution for muons
146###############################
147
148module MomentumSmearing MuonMomentumSmearing {
149 set InputArray MuonTrackingEfficiency/muons
150 set OutputArray muons
151
152 # set ResolutionFormula {resolution formula as a function of eta and pt}
153
154 # resolution formula for muons
155 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
156 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
157 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
158}
159
160
161##############
162# Track merger
163##############
164
165module Merger TrackMerger {
166# add InputArray InputArray
167 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
168 add InputArray ElectronMomentumSmearing/electrons
169 add InputArray MuonMomentumSmearing/muons
170 set OutputArray tracks
171}
172
173#############
174# ECAL
175#############
176
177module SimpleCalorimeter ECal {
178 set ParticleInputArray ParticlePropagator/stableParticles
179 set TrackInputArray TrackMerger/tracks
180
181 set TowerOutputArray ecalTowers
182 set EFlowTrackOutputArray eflowTracks
183 set EFlowTowerOutputArray eflowPhotons
184
185 set IsEcal true
186
187 set EnergyMin 0.5
188 set EnergySignificanceMin 2.0
189
190 set SmearTowerCenter true
191
192 set pi [expr {acos(-1)}]
193
194 # lists of the edges of each tower in eta and phi
195 # each list starts with the lower edge of the first tower
196 # the list ends with the higher edged of the last tower
197
198 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
199
200 set PhiBins {}
201 for {set i -180} {$i <= 180} {incr i} {
202 add PhiBins [expr {$i * $pi/180.0}]
203 }
204
205 # 0.02 unit in eta up to eta = 1.5 (barrel)
206 for {set i -85} {$i <= 86} {incr i} {
207 set eta [expr {$i * 0.0174}]
208 add EtaPhiBins $eta $PhiBins
209 }
210
211 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
212
213 set PhiBins {}
214 for {set i -180} {$i <= 180} {incr i} {
215 add PhiBins [expr {$i * $pi/180.0}]
216 }
217
218 # 0.02 unit in eta up to eta = 3
219 for {set i 1} {$i <= 84} {incr i} {
220 set eta [expr { -2.958 + $i * 0.0174}]
221 add EtaPhiBins $eta $PhiBins
222 }
223
224 for {set i 1} {$i <= 84} {incr i} {
225 set eta [expr { 1.4964 + $i * 0.0174}]
226 add EtaPhiBins $eta $PhiBins
227 }
228
229 # take present CMS granularity for HF
230
231 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
232 set PhiBins {}
233 for {set i -18} {$i <= 18} {incr i} {
234 add PhiBins [expr {$i * $pi/18.0}]
235 }
236
237 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
238 add EtaPhiBins $eta $PhiBins
239 }
240
241
242 add EnergyFraction {0} {0.0}
243 # energy fractions for e, gamma and pi0
244 add EnergyFraction {11} {1.0}
245 add EnergyFraction {22} {1.0}
246 add EnergyFraction {111} {1.0}
247 # energy fractions for muon, neutrinos and neutralinos
248 add EnergyFraction {12} {0.0}
249 add EnergyFraction {13} {0.0}
250 add EnergyFraction {14} {0.0}
251 add EnergyFraction {16} {0.0}
252 add EnergyFraction {1000022} {0.0}
253 add EnergyFraction {1000023} {0.0}
254 add EnergyFraction {1000025} {0.0}
255 add EnergyFraction {1000035} {0.0}
256 add EnergyFraction {1000045} {0.0}
257 # energy fractions for K0short and Lambda
258 add EnergyFraction {310} {0.3}
259 add EnergyFraction {3122} {0.3}
260
261 # set ResolutionFormula {resolution formula as a function of eta and energy}
262
263 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
264
265 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
266 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
267 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
268 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
269 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
270
271}
272
273
274#############
275# HCAL
276#############
277
278module SimpleCalorimeter HCal {
279 set ParticleInputArray ParticlePropagator/stableParticles
280 set TrackInputArray ECal/eflowTracks
281
282 set TowerOutputArray hcalTowers
283 set EFlowTrackOutputArray eflowTracks
284 set EFlowTowerOutputArray eflowNeutralHadrons
285
286 set IsEcal false
287
288 set EnergyMin 1.0
289 set EnergySignificanceMin 2.0
290
291 set SmearTowerCenter true
292
293 set pi [expr {acos(-1)}]
294
295 # lists of the edges of each tower in eta and phi
296 # each list starts with the lower edge of the first tower
297 # the list ends with the higher edged of the last tower
298
299 # 5 degrees towers
300 set PhiBins {}
301 for {set i -36} {$i <= 36} {incr i} {
302 add PhiBins [expr {$i * $pi/36.0}]
303 }
304 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
305 add EtaPhiBins $eta $PhiBins
306 }
307
308 # 10 degrees towers
309 set PhiBins {}
310 for {set i -18} {$i <= 18} {incr i} {
311 add PhiBins [expr {$i * $pi/18.0}]
312 }
313 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
314 add EtaPhiBins $eta $PhiBins
315 }
316
317 # 20 degrees towers
318 set PhiBins {}
319 for {set i -9} {$i <= 9} {incr i} {
320 add PhiBins [expr {$i * $pi/9.0}]
321 }
322 foreach eta {-5 -4.7 -4.525 4.7 5} {
323 add EtaPhiBins $eta $PhiBins
324 }
325
326 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
327 add EnergyFraction {0} {1.0}
328 # energy fractions for e, gamma and pi0
329 add EnergyFraction {11} {0.0}
330 add EnergyFraction {22} {0.0}
331 add EnergyFraction {111} {0.0}
332 # energy fractions for muon, neutrinos and neutralinos
333 add EnergyFraction {12} {0.0}
334 add EnergyFraction {13} {0.0}
335 add EnergyFraction {14} {0.0}
336 add EnergyFraction {16} {0.0}
337 add EnergyFraction {1000022} {0.0}
338 add EnergyFraction {1000023} {0.0}
339 add EnergyFraction {1000025} {0.0}
340 add EnergyFraction {1000035} {0.0}
341 add EnergyFraction {1000045} {0.0}
342 # energy fractions for K0short and Lambda
343 add EnergyFraction {310} {0.7}
344 add EnergyFraction {3122} {0.7}
345
346 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
347 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
348 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
349
350}
351
352###################################################
353# Tower Merger (in case not using e-flow algorithm)
354###################################################
355
356module Merger Calorimeter {
357# add InputArray InputArray
358 add InputArray ECal/ecalTowers
359 add InputArray HCal/hcalTowers
360 set OutputArray towers
361}
362
363
364####################
365# Energy flow merger
366####################
367
368module Merger EFlowMerger {
369# add InputArray InputArray
370 add InputArray HCal/eflowTracks
371 add InputArray ECal/eflowPhotons
372 add InputArray HCal/eflowNeutralHadrons
373 set OutputArray eflow
374}
375
Note: See TracBrowser for help on using the repository browser.