Fork me on GitHub

source: git/cards/delphes_card_ATLAS.tcl@ 122e1e5

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 122e1e5 was 6153fb0, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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