Fork me on GitHub

Ticket #844: delphes_card_ATLAS.tcl

File delphes_card_ATLAS.tcl, 18.0 KB (added by Cyril Becot, 9 years ago)

modified TCL card

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