Fork me on GitHub

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

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

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