Fork me on GitHub

source: svn/trunk/examples/delphes_card_ATLAS.tcl@ 1307

Last change on this file since 1307 was 1279, checked in by Pavel Demin, 11 years ago

revert to pre 3.0.10 charged hadron resolution for ATLAS
File size: 17.9 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 ElectronEnergySmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 Calorimeter
18 EFlowMerger
19
20 PhotonEfficiency
21 PhotonIsolation
22
23 ElectronEfficiency
24 ElectronIsolation
25
26 MuonEfficiency
27 MuonIsolation
28
29 MissingET
30
31 GenJetFinder
32 FastJetFinder
33
34 JetEnergyScale
35
36 BTagging
37 TauTagging
38
39 UniqueObjectFinder
40
41 ScalarHT
42
43 TreeWriter
44}
45
46#################################
47# Propagate particles in cylinder
48#################################
49
50module ParticlePropagator ParticlePropagator {
51 set InputArray Delphes/stableParticles
52
53 set OutputArray stableParticles
54 set ChargedHadronOutputArray chargedHadrons
55 set ElectronOutputArray electrons
56 set MuonOutputArray muons
57
58 # radius of the magnetic field coverage, in m
59 set Radius 1.15
60 # half-length of the magnetic field coverage, in m
61 set HalfLength 3.51
62
63 # magnetic field
64 set Bz 2.0
65}
66
67####################################
68# Charged hadron tracking efficiency
69####################################
70
71module Efficiency ChargedHadronTrackingEfficiency {
72 set InputArray ParticlePropagator/chargedHadrons
73 set OutputArray chargedHadrons
74
75 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
76
77 # tracking efficiency formula for charged hadrons
78 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
79 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
80 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
81 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
83 (abs(eta) > 2.5) * (0.00)}
84}
85
86##############################
87# Electron tracking efficiency
88##############################
89
90module Efficiency ElectronTrackingEfficiency {
91 set InputArray ParticlePropagator/electrons
92 set OutputArray electrons
93
94 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
95
96 # tracking efficiency formula for electrons
97 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
98 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
99 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
100 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
101 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
104 (abs(eta) > 2.5) * (0.00)}
105}
106
107##########################
108# Muon tracking efficiency
109##########################
110
111module Efficiency MuonTrackingEfficiency {
112 set InputArray ParticlePropagator/muons
113 set OutputArray muons
114
115 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
116
117 # tracking efficiency formula for muons
118 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
119 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
120 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
123 (abs(eta) > 2.5) * (0.00)}
124}
125
126########################################
127# Momentum resolution for charged tracks
128########################################
129
130module MomentumSmearing ChargedHadronMomentumSmearing {
131 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
132 set OutputArray chargedHadrons
133
134 # set ResolutionFormula {resolution formula as a function of eta and pt}
135
136 # resolution formula for charged hadrons
137 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
138 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
139 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
140 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
145}
146
147#################################
148# Energy resolution for electrons
149#################################
150
151module EnergySmearing ElectronEnergySmearing {
152 set InputArray ElectronTrackingEfficiency/electrons
153 set OutputArray electrons
154
155 # set ResolutionFormula {resolution formula as a function of eta and energy}
156
157 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
158 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
159 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
160 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
161
162}
163
164###############################
165# Momentum resolution for muons
166###############################
167
168module MomentumSmearing MuonMomentumSmearing {
169 set InputArray MuonTrackingEfficiency/muons
170 set OutputArray muons
171
172 # set ResolutionFormula {resolution formula as a function of eta and pt}
173
174 # resolution formula for muons
175 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
176 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
177 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
178 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
179 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
183}
184
185##############
186# Track merger
187##############
188
189module Merger TrackMerger {
190# add InputArray InputArray
191 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
192 add InputArray ElectronEnergySmearing/electrons
193 set OutputArray tracks
194}
195
196#############
197# Calorimeter
198#############
199
200module Calorimeter Calorimeter {
201 set ParticleInputArray ParticlePropagator/stableParticles
202 set TrackInputArray TrackMerger/tracks
203
204 set TowerOutputArray towers
205 set PhotonOutputArray photons
206
207 set EFlowTrackOutputArray eflowTracks
208 set EFlowTowerOutputArray eflowTowers
209
210 set pi [expr {acos(-1)}]
211
212 # lists of the edges of each tower in eta and phi
213 # each list starts with the lower edge of the first tower
214 # the list ends with the higher edged of the last tower
215
216 # 10 degrees towers
217 set PhiBins {}
218 for {set i -18} {$i <= 18} {incr i} {
219 add PhiBins [expr {$i * $pi/18.0}]
220 }
221 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
222 add EtaPhiBins $eta $PhiBins
223 }
224
225 # 20 degrees towers
226 set PhiBins {}
227 for {set i -9} {$i <= 9} {incr i} {
228 add PhiBins [expr {$i * $pi/9.0}]
229 }
230 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
231 add EtaPhiBins $eta $PhiBins
232 }
233
234 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
235 add EnergyFraction {0} {0.0 1.0}
236 # energy fractions for e, gamma and pi0
237 add EnergyFraction {11} {1.0 0.0}
238 add EnergyFraction {22} {1.0 0.0}
239 add EnergyFraction {111} {1.0 0.0}
240 # energy fractions for muon, neutrinos and neutralinos
241 add EnergyFraction {12} {0.0 0.0}
242 add EnergyFraction {13} {0.0 0.0}
243 add EnergyFraction {14} {0.0 0.0}
244 add EnergyFraction {16} {0.0 0.0}
245 add EnergyFraction {1000022} {0.0 0.0}
246 add EnergyFraction {1000023} {0.0 0.0}
247 add EnergyFraction {1000025} {0.0 0.0}
248 add EnergyFraction {1000035} {0.0 0.0}
249 add EnergyFraction {1000045} {0.0 0.0}
250 # energy fractions for K0short and Lambda
251 add EnergyFraction {310} {0.3 0.7}
252 add EnergyFraction {3122} {0.3 0.7}
253
254 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
255 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
256 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
257 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
258 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
259 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
260
261 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
262 # http://arxiv.org/pdf/hep-ex/0004009v1
263 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
264 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
265 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
266 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
267}
268
269####################
270# Energy flow merger
271####################
272
273module Merger EFlowMerger {
274# add InputArray InputArray
275 add InputArray Calorimeter/eflowTracks
276 add InputArray Calorimeter/eflowTowers
277 add InputArray MuonMomentumSmearing/muons
278 set OutputArray eflow
279}
280
281###################
282# Photon efficiency
283###################
284
285module Efficiency PhotonEfficiency {
286 set InputArray Calorimeter/photons
287 set OutputArray photons
288
289 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
290
291 # efficiency formula for photons
292 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
293 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
294 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
295 (abs(eta) > 2.5) * (0.00)}
296}
297
298##################
299# Photon isolation
300##################
301
302module Isolation PhotonIsolation {
303 set CandidateInputArray PhotonEfficiency/photons
304 set IsolationInputArray EFlowMerger/eflow
305
306 set OutputArray photons
307
308 set DeltaRMax 0.5
309
310 set PTMin 0.5
311
312 set PTRatioMax 0.1
313}
314
315#####################
316# Electron efficiency
317#####################
318
319module Efficiency ElectronEfficiency {
320 set InputArray ElectronEnergySmearing/electrons
321 set OutputArray electrons
322
323 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
324
325 # efficiency formula for electrons
326 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
327 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
328 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
329 (abs(eta) > 2.5) * (0.00)}
330}
331
332####################
333# Electron isolation
334####################
335
336module Isolation ElectronIsolation {
337 set CandidateInputArray ElectronEfficiency/electrons
338 set IsolationInputArray EFlowMerger/eflow
339
340 set OutputArray electrons
341
342 set DeltaRMax 0.5
343
344 set PTMin 0.5
345
346 set PTRatioMax 0.1
347}
348
349#################
350# Muon efficiency
351#################
352
353module Efficiency MuonEfficiency {
354 set InputArray MuonMomentumSmearing/muons
355 set OutputArray muons
356
357 # set EfficiencyFormula {efficiency as a function of eta and pt}
358
359 # efficiency formula for muons
360 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
361 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
362 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
363 (abs(eta) > 2.7) * (0.00)}
364}
365
366################
367# Muon isolation
368################
369
370module Isolation MuonIsolation {
371 set CandidateInputArray MuonEfficiency/muons
372 set IsolationInputArray EFlowMerger/eflow
373
374 set OutputArray muons
375
376 set DeltaRMax 0.5
377
378 set PTMin 0.5
379
380 set PTRatioMax 0.1
381}
382
383###################
384# Missing ET merger
385###################
386
387module Merger MissingET {
388# add InputArray InputArray
389 add InputArray EFlowMerger/eflow
390 set MomentumOutputArray momentum
391}
392
393##################
394# Scalar HT merger
395##################
396
397module Merger ScalarHT {
398# add InputArray InputArray
399 add InputArray UniqueObjectFinder/jets
400 add InputArray UniqueObjectFinder/electrons
401 add InputArray UniqueObjectFinder/photons
402 add InputArray UniqueObjectFinder/muons
403 set EnergyOutputArray energy
404}
405
406#####################
407# MC truth jet finder
408#####################
409
410module FastJetFinder GenJetFinder {
411 set InputArray Delphes/stableParticles
412
413 set OutputArray jets
414
415 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
416 set JetAlgorithm 6
417 set ParameterR 0.6
418
419 set JetPTMin 20.0
420}
421
422############
423# Jet finder
424############
425
426module FastJetFinder FastJetFinder {
427 set InputArray Calorimeter/towers
428
429 set OutputArray jets
430
431 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
432 set JetAlgorithm 6
433 set ParameterR 0.6
434
435 set JetPTMin 20.0
436}
437
438##################
439# Jet Energy Scale
440##################
441
442module EnergyScale JetEnergyScale {
443 set InputArray FastJetFinder/jets
444 set OutputArray jets
445
446 # scale formula for jets
447 set ScaleFormula {1.08}
448}
449
450###########
451# b-tagging
452###########
453
454module BTagging BTagging {
455 set PartonInputArray Delphes/partons
456 set JetInputArray JetEnergyScale/jets
457
458 set BitNumber 0
459
460 set DeltaR 0.5
461
462 set PartonPTMin 1.0
463
464 set PartonEtaMax 2.5
465
466 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
467 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
468 # gluon's PDG code has the lowest priority
469
470 # default efficiency formula (misidentification rate)
471 add EfficiencyFormula {0} {0.001}
472
473 # efficiency formula for c-jets (misidentification rate)
474 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
475 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
476 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
477 (abs(eta) > 2.5) * (0.000)}
478
479 # efficiency formula for b-jets
480 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
481 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
482 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
483 (abs(eta) > 2.5) * (0.000)}
484}
485
486module TauTagging TauTagging {
487 set ParticleInputArray Delphes/allParticles
488 set PartonInputArray Delphes/partons
489 set JetInputArray JetEnergyScale/jets
490
491 set DeltaR 0.5
492
493 set TauPTMin 1.0
494
495 set TauEtaMax 2.5
496
497 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
498
499 # default efficiency formula (misidentification rate)
500 add EfficiencyFormula {0} {0.001}
501 # efficiency formula for tau-jets
502 add EfficiencyFormula {15} {0.4}
503}
504
505#####################################################
506# Find uniquely identified photons/electrons/tau/jets
507#####################################################
508
509module UniqueObjectFinder UniqueObjectFinder {
510# earlier arrays take precedence over later ones
511# add InputArray InputArray OutputArray
512 add InputArray PhotonIsolation/photons photons
513 add InputArray ElectronIsolation/electrons electrons
514 add InputArray MuonIsolation/muons muons
515 add InputArray JetEnergyScale/jets jets
516}
517
518##################
519# ROOT tree writer
520##################
521
522module TreeWriter TreeWriter {
523# add Branch InputArray BranchName BranchClass
524 add Branch Delphes/allParticles Particle GenParticle
525 add Branch TrackMerger/tracks Track Track
526 add Branch Calorimeter/towers Tower Tower
527# add Branch Calorimeter/eflowTracks EFlowTrack Track
528# add Branch Calorimeter/eflowTowers EFlowTower Tower
529# add Branch MuonMomentumSmearing/muons EFlowMuon Muon
530 add Branch GenJetFinder/jets GenJet Jet
531 add Branch UniqueObjectFinder/jets Jet Jet
532 add Branch UniqueObjectFinder/electrons Electron Electron
533 add Branch UniqueObjectFinder/photons Photon Photon
534 add Branch UniqueObjectFinder/muons Muon Muon
535 add Branch MissingET/momentum MissingET MissingET
536 add Branch ScalarHT/energy ScalarHT ScalarHT
537}
538
Note: See TracBrowser for help on using the repository browser.