Fork me on GitHub

source: git/cards/delphes_card_ATLAS.tcl@ a933829

ImprovedOutputFile Timing dual_readout llp
Last change on this file since a933829 was a0c065d, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

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