Fork me on GitHub

source: git/cards/delphes_card_LHCb.tcl@ d870fc5

ImprovedOutputFile Timing dual_readout llp
Last change on this file since d870fc5 was 8624f58, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

replace Ecal with ECal and Hcal with HCal
  • Property mode set to 100644
File size: 10.9 KB
Line 
1#set MaxEvents 1000
2#set RandomSeed 123
3
4
5#######################################
6# Order of execution of various modules
7#######################################
8
9set ExecutionPath {
10
11 ParticlePropagator
12
13 ChargedHadronMomentumSmearing
14 ElectronEnergySmearing
15 MuonMomentumSmearing
16
17 TrackMerger
18 ImpactParameterSmearing
19
20 IdentificationMap
21
22 ECal
23 HCal
24
25 TreeWriter
26}
27
28
29
30#################################
31# Propagate particles in cylinder
32#################################
33
34module ParticlePropagator ParticlePropagator {
35 set InputArray Delphes/stableParticles
36
37 set OutputArray stableParticles
38 set ChargedHadronOutputArray chargedHadrons
39 set ElectronOutputArray electrons
40 set MuonOutputArray muons
41
42 # radius of the magnetic field coverage, in m
43
44 set Radius 3.31
45
46
47 # half-length of the magnetic field coverage, in m
48 set HalfLength 12.0
49
50 # magnetic field
51 set Bz 1.1
52
53 # Need to veto anything with theta > 0.269 rad -> eta = 2
54 # theta < 0.0135 rad -> eta = 5
55
56 # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269)
57
58}
59
60
61########################################
62# Momentum resolution for charged tracks
63########################################
64
65module MomentumSmearing ChargedHadronMomentumSmearing {
66 set InputArray ParticlePropagator/chargedHadrons
67 set OutputArray chargedHadrons
68
69 # set ResolutionFormula {resolution formula as a function of eta and pt}
70
71 # resolution formula for charged hadrons
72 set ResolutionFormula {(eta > 2.0 && eta <= 5.0) * (pt > 0.5) * (0.005)}
73}
74
75#################################
76# Energy resolution for electrons
77#################################
78
79module EnergySmearing ElectronEnergySmearing {
80 set InputArray ParticlePropagator/electrons
81 set OutputArray electrons
82
83 # set ResolutionFormula {resolution formula as a function of eta and energy}
84
85 # resolution formula for electrons
86 set ResolutionFormula { (eta > 2.0 && eta <= 5.0) * (energy > 0.1 && energy <= 8.0) * (energy*0.05) + \
87 (eta > 2.0 && eta <= 5.0) * (energy > 8.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
88 }
89
90###############################
91# Momentum resolution for muons
92###############################
93
94module MomentumSmearing MuonMomentumSmearing {
95 set InputArray ParticlePropagator/muons
96 set OutputArray muons
97
98 # resolution formula for muons
99 set ResolutionFormula {(eta > 2.0 && eta <= 5.0) * (pt > 0.5)* (0.005)}
100}
101
102
103##############
104# Track merger
105##############
106
107module Merger TrackMerger {
108# add InputArray InputArray
109 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
110 add InputArray ElectronEnergySmearing/electrons
111 add InputArray MuonMomentumSmearing/muons
112 set OutputArray tracks
113}
114
115
116
117################################
118# Track impact parameter smearing
119################################
120
121module ImpactParameterSmearing ImpactParameterSmearing {
122 set InputArray TrackMerger/tracks
123 set OutputArray tracks
124
125
126# absolute impact parameter smearing formula (in mm) as a function of pt and eta
127set ResolutionFormula {0.0116 + 0.0234/pt}
128
129}
130
131###### ref. JINST 9 C01065 #############
132
133
134####################################
135# Charged hadron PID
136####################################
137
138module IdentificationMap IdentificationMap {
139 set InputArray ImpactParameterSmearing/tracks
140 set OutputArray tracks
141
142 # {PID in} {PID out} {formula}
143 # make sure "PID in" and "PID out" have the same charge (e.g {-13} {211} or {-321} {211})
144 # {211} {-13} is equivalent to {-211} {13} (and needs to be written once only...)
145
146
147
148
149
150
151 # --- pions ---
152
153 add EfficiencyFormula {211} {211} { (eta <= 2.0) * (0.00) + \
154 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
155 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \
156 (eta > 5.0) * (0.00)}
157
158 add EfficiencyFormula {211} {-13} { (eta <= 2.0) * (0.00) + \
159 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
160 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) + \
161 (eta > 5.0) * (0.00)}
162
163
164 # --- kaons ---
165
166
167 add EfficiencyFormula {321} {321} { (eta <= 2.0) * (0.00) + \
168 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
169 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \
170 (eta > 5.0) * (0.00)}
171
172 add EfficiencyFormula {321} {-13} { (eta <= 2.0) * (0.00) + \
173 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
174 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) + \
175 (eta > 5.0) * (0.00)}
176
177
178 # --- protons ---
179
180
181 add EfficiencyFormula {2212} {2212} { (eta <= 2.0) * (0.00) + \
182 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
183 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.95) + \
184 (eta > 5.0) * (0.00)}
185
186 add EfficiencyFormula {2212} {-13} { (eta <= 2.0) * (0.00) + \
187 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
188 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.002) + \
189 (eta > 5.0) * (0.00)}
190
191
192
193 # --- muons ---
194
195
196
197 add EfficiencyFormula {-13} {-13} { (eta <= 2.0) * (0.00) + \
198 (eta > 2.0 && eta <= 5.0) * (pt < 0.8)* (0.00) + \
199 (eta > 2.0 && eta <= 5.0) * (pt >= 0.8)* (0.97) + \
200 (eta > 5.0) * (0.00)}
201
202
203 # efficiency for other charged particles (should be always {0} {0} {formula})
204
205 add EfficiencyFormula {0} {0} { (eta <= 2.0) * (0.00) + \
206 (eta > 2.0 && eta <= 5.0) * (pt < 0.8) * (0.00) + \
207 (eta > 2.0 && eta <= 5.0) * (pt > 0.8) * (0.95) + \
208 (eta > 5.0) * (0.00)}
209
210}
211
212###### ref. JINST 8 P10020 #############
213
214
215#############
216# ECAL
217#############
218
219module SimpleCalorimeter ECal {
220 set ParticleInputArray ParticlePropagator/stableParticles
221 set TrackInputArray IdentificationMap/tracks
222
223 set TowerOutputArray ecalTowers
224 set EFlowTowerOutputArray eflowPhotons
225
226 set pi [expr {acos(-1)}]
227
228 # lists of the edges of each tower in eta and phi
229 # each list starts with the lower edge of the first tower
230 # the list ends with the higher edged of the last tower
231
232 # 1 degree towers
233 set PhiBins {}
234 for {set i -180} {$i <= 180} {incr i} {
235 add PhiBins [expr {$i * $pi/180.0}]
236 }
237
238 # 0.02 unit in eta from eta = 3.2 to eta = 5.0
239 for {set i 1} {$i <= 90} {incr i} {
240 set eta [expr {3.2 + $i * 0.02}]
241 add EtaPhiBins $eta $PhiBins
242 }
243
244 # 1.25 degree towers
245 set PhiBins {}
246 for {set i -135} {$i <= 135} {incr i} {
247 add PhiBins [expr {$i * $pi/135.0}]
248 }
249
250 # 0.025 unit in eta from eta = 2.6 to eta = 3.2
251 for {set i 1} {$i <= 24} {incr i} {
252 set eta [expr {2.6 + $i * 0.025}]
253 add EtaPhiBins $eta $PhiBins
254 }
255
256 # 1.25 degree towers
257 set PhiBins {}
258 for {set i -100} {$i <= 100} {incr i} {
259 add PhiBins [expr {$i * $pi/100.0}]
260 }
261
262 # 0.04 unit in eta from eta = 2.0 to eta = 2.6
263 for {set i 0} {$i <= 24} {incr i} {
264 set eta [expr {2.0 + $i * 0.04}]
265 add EtaPhiBins $eta $PhiBins
266 }
267
268 add EnergyFraction {0} {0.0}
269 # energy fractions for e, gamma and pi0
270 add EnergyFraction {11} {1.0}
271 add EnergyFraction {22} {1.0}
272 add EnergyFraction {111} {1.0}
273 # energy fractions for muon, neutrinos and neutralinos
274 add EnergyFraction {12} {0.0}
275 add EnergyFraction {13} {0.0}
276 add EnergyFraction {14} {0.0}
277 add EnergyFraction {16} {0.0}
278 add EnergyFraction {1000022} {0.0}
279 add EnergyFraction {1000023} {0.0}
280 add EnergyFraction {1000025} {0.0}
281 add EnergyFraction {1000035} {0.0}
282 add EnergyFraction {1000045} {0.0}
283 # energy fractions for K0short and Lambda
284 add EnergyFraction {310} {0.3}
285 add EnergyFraction {3122} {0.3}
286
287 set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
288}
289
290#############
291# HCAL
292#############
293
294module SimpleCalorimeter HCal {
295 set ParticleInputArray ParticlePropagator/stableParticles
296 set TrackInputArray IdentificationMap/tracks
297
298 set TowerOutputArray hcalTowers
299 set EFlowTowerOutputArray eflowNeutralHadrons
300
301 set pi [expr {acos(-1)}]
302
303 # lists of the edges of each tower in eta and phi
304 # each list starts with the lower edge of the first tower
305 # the list ends with the higher edged of the last tower
306
307 # 1 degree towers
308 set PhiBins {}
309 for {set i -16} {$i <= 16} {incr i} {
310 add PhiBins [expr {$i * $pi/16.0}]
311 }
312
313 # 0.20 unit in eta from eta = 2.6 to eta = 5.0
314 for {set i 1} {$i <= 12} {incr i} {
315 set eta [expr {2.6 + $i * 0.2}]
316 add EtaPhiBins $eta $PhiBins
317 }
318
319 # 1 degree towers
320 set PhiBins {}
321 for {set i -32} {$i <= 32} {incr i} {
322 add PhiBins [expr {$i * $pi/32.0}]
323 }
324
325 # 0.1 unit in eta from eta = 2 to eta = 2.6
326 for {set i 0} {$i <= 6} {incr i} {
327 set eta [expr {2.0 + $i * 0.1}]
328 add EtaPhiBins $eta $PhiBins
329 }
330
331
332
333 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
334 add EnergyFraction {0} {1.0}
335 # energy fractions for e, gamma and pi0
336 add EnergyFraction {11} {0.0}
337 add EnergyFraction {22} {0.0}
338 add EnergyFraction {111} {0.0}
339 # energy fractions for muon, neutrinos and neutralinos
340 add EnergyFraction {12} {0.0}
341 add EnergyFraction {13} {0.0}
342 add EnergyFraction {14} {0.0}
343 add EnergyFraction {16} {0.0}
344 add EnergyFraction {1000022} {0.0}
345 add EnergyFraction {1000023} {0.0}
346 add EnergyFraction {1000025} {0.0}
347 add EnergyFraction {1000035} {0.0}
348 add EnergyFraction {1000045} {0.0}
349 # energy fractions for K0short and Lambda
350 add EnergyFraction {310} {0.7}
351 add EnergyFraction {3122} {0.7}
352
353 set ResolutionFormula { (eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.05^2 + energy*0.80^2)}
354 }
355
356
357##################
358# ROOT tree writer
359##################
360
361# tracks, towers and eflow objects are not stored by default in the output.
362# if needed (for jet constituent or other studies), uncomment the relevant
363# "add Branch ..." lines.
364
365module TreeWriter TreeWriter {
366# add Branch InputArray BranchName BranchClass
367
368 add Branch Delphes/allParticles Particle GenParticle
369
370 add Branch IdentificationMap/tracks Track Track
371 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
372 add Branch ECal/eflowPhotons Photon Photon
373
374}
375
Note: See TracBrowser for help on using the repository browser.