Fork me on GitHub

Ticket #1333: delphes_card.tcl

File delphes_card.tcl, 18.0 KB (added by Michael Albergo, 6 years ago)

My Delphes 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.01^2 + pt^2*1.5e-4^2) +
144 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^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 1.0
484
485 set ComputeTrimming 1
486 set RTrim 0.3
487 set PtFracTrim 0.05
488
489 set JetPTMin 250.0
490 set JetPTMax 300.0
491}
492
493##################
494# Jet Energy Scale
495##################
496
497module EnergyScale JetEnergyScale {
498 set InputArray FastJetFinder/jets
499 set OutputArray jets
500
501 # scale formula for jets
502 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
503}
504
505########################
506# Jet Flavor Association
507########################
508
509module JetFlavorAssociation JetFlavorAssociation {
510
511 set PartonInputArray Delphes/partons
512 set ParticleInputArray Delphes/allParticles
513 set ParticleLHEFInputArray Delphes/allParticlesLHEF
514 set JetInputArray JetEnergyScale/jets
515
516 set DeltaR 0.5
517 set PartonPTMin 1.0
518 set PartonEtaMax 2.5
519
520}
521
522###########
523# b-tagging
524###########
525
526module BTagging BTagging {
527 set JetInputArray JetEnergyScale/jets
528
529 set BitNumber 0
530
531 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
532 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
533 # gluon's PDG code has the lowest priority
534
535 # based on ATL-PHYS-PUB-2015-022
536
537 # default efficiency formula (misidentification rate)
538 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
539
540 # efficiency formula for c-jets (misidentification rate)
541 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
542
543 # efficiency formula for b-jets
544 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
545}
546
547#############
548# tau-tagging
549#############
550
551module TrackCountingTauTagging TauTagging {
552
553 set ParticleInputArray Delphes/allParticles
554 set PartonInputArray Delphes/partons
555 set TrackInputArray TrackMerger/tracks
556 set JetInputArray JetEnergyScale/jets
557
558 set DeltaR 0.2
559 set DeltaRTrack 0.2
560
561 set TrackPTMin 1.0
562
563 set TauPTMin 1.0
564 set TauEtaMax 2.5
565
566 # instructions: {n-prongs} {eff}
567
568 # 1 - one prong efficiency
569 # 2 - two or more efficiency
570 # -1 - one prong mistag rate
571 # -2 - two or more mistag rate
572
573 set BitNumber 0
574
575 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
576 add EfficiencyFormula {1} {0.70}
577 add EfficiencyFormula {2} {0.60}
578 add EfficiencyFormula {-1} {0.02}
579 add EfficiencyFormula {-2} {0.01}
580
581}
582
583#####################################################
584# Find uniquely identified photons/electrons/tau/jets
585#####################################################
586
587module UniqueObjectFinder UniqueObjectFinder {
588# earlier arrays take precedence over later ones
589# add InputArray InputArray OutputArray
590 add InputArray PhotonIsolation/photons photons
591 add InputArray ElectronIsolation/electrons electrons
592 add InputArray MuonIsolation/muons muons
593 add InputArray JetEnergyScale/jets jets
594}
595
596##################
597# ROOT tree writer
598##################
599
600# tracks, towers and eflow objects are not stored by default in the output.
601# if needed (for jet constituent or other studies), uncomment the relevant
602# "add Branch ..." lines.
603
604module TreeWriter TreeWriter {
605# add Branch InputArray BranchName BranchClass
606 add Branch Delphes/allParticles Particle GenParticle
607
608 add Branch TrackMerger/tracks Track Track
609 add Branch Calorimeter/towers Tower Tower
610
611 add Branch Calorimeter/eflowTracks EFlowTrack Track
612 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
613 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
614
615 add Branch GenJetFinder/jets GenJet Jet
616# add Branch GenMissingET/momentum GenMissingET MissingET
617
618 add Branch UniqueObjectFinder/jets Jet Jet
619 add Branch UniqueObjectFinder/electrons Electron Electron
620 add Branch UniqueObjectFinder/photons Photon Photon
621 add Branch UniqueObjectFinder/muons Muon Muon
622 add Branch MissingET/momentum MissingET MissingET
623# add Branch ScalarHT/energy ScalarHT ScalarHT
624}