MA5SandBox: delphes_card_cms_exo_17_015.tcl

File delphes_card_cms_exo_17_015.tcl, 21.2 KB (added by Benjamin Fuks, 6 years ago)
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
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23 EFlowFilter
24
25 PhotonEfficiency
26
27 ElectronFilter
28 ElectronEfficiency
29
30 ChargedHadronFilter
31
32 MuonEfficiency
33
34 MissingET
35
36 NeutrinoFilter
37 GenJetFinder
38 GenMissingET
39
40 FastJetFinder
41
42 JetEnergyScale
43
44 JetFlavorAssociation
45
46 BTagging
47 TauTagging
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.29
66 # half-length of the magnetic field coverage, in m
67 set HalfLength 3.00
68
69 # magnetic field
70 set Bz 3.8
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 && pt <= 1.0e3) * (0.99) +
127 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
128
129 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
130 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
132 (abs(eta) > 2.5) * (0.00)}
133}
134
135########################################
136# Momentum resolution for charged tracks
137########################################
138
139module MomentumSmearing ChargedHadronMomentumSmearing {
140 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
141 set OutputArray chargedHadrons
142
143 # set ResolutionFormula {resolution formula as a function of eta and pt}
144
145 # resolution formula for charged hadrons
146 # based on arXiv:1405.6569
147 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
148 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
150}
151
152###################################
153# Momentum resolution for electrons
154###################################
155
156module MomentumSmearing ElectronMomentumSmearing {
157 set InputArray ElectronTrackingEfficiency/electrons
158 set OutputArray electrons
159
160 # set ResolutionFormula {resolution formula as a function of eta and energy}
161
162 # resolution formula for electrons
163 # based on arXiv:1405.6569
164 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
165 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
167}
168
169###############################
170# Momentum resolution for muons
171###############################
172
173module MomentumSmearing MuonMomentumSmearing {
174 set InputArray MuonTrackingEfficiency/muons
175 set OutputArray muons
176
177 # set ResolutionFormula {resolution formula as a function of eta and pt}
178
179 # resolution formula for muons
180 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
181 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
183}
184
185
186##############
187# Track merger
188##############
189
190module Merger TrackMerger {
191# add InputArray InputArray
192 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
193 add InputArray ElectronMomentumSmearing/electrons
194 add InputArray MuonMomentumSmearing/muons
195 set OutputArray tracks
196}
197
198
199
200#############
201# ECAL
202#############
203
204module SimpleCalorimeter ECal {
205 set ParticleInputArray ParticlePropagator/stableParticles
206 set TrackInputArray TrackMerger/tracks
207
208 set TowerOutputArray ecalTowers
209 set EFlowTrackOutputArray eflowTracks
210 set EFlowTowerOutputArray eflowPhotons
211
212 set IsEcal true
213
214 set EnergyMin 0.5
215 set EnergySignificanceMin 2.0
216
217 set SmearTowerCenter true
218
219 set pi [expr {acos(-1)}]
220
221 # lists of the edges of each tower in eta and phi
222 # each list starts with the lower edge of the first tower
223 # the list ends with the higher edged of the last tower
224
225 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
226
227 set PhiBins {}
228 for {set i -180} {$i <= 180} {incr i} {
229 add PhiBins [expr {$i * $pi/180.0}]
230 }
231
232 # 0.02 unit in eta up to eta = 1.5 (barrel)
233 for {set i -85} {$i <= 86} {incr i} {
234 set eta [expr {$i * 0.0174}]
235 add EtaPhiBins $eta $PhiBins
236 }
237
238 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
239
240 set PhiBins {}
241 for {set i -180} {$i <= 180} {incr i} {
242 add PhiBins [expr {$i * $pi/180.0}]
243 }
244
245 # 0.02 unit in eta up to eta = 3
246 for {set i 1} {$i <= 84} {incr i} {
247 set eta [expr { -2.958 + $i * 0.0174}]
248 add EtaPhiBins $eta $PhiBins
249 }
250
251 for {set i 1} {$i <= 84} {incr i} {
252 set eta [expr { 1.4964 + $i * 0.0174}]
253 add EtaPhiBins $eta $PhiBins
254 }
255
256 # take present CMS granularity for HF
257
258 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
259 set PhiBins {}
260 for {set i -18} {$i <= 18} {incr i} {
261 add PhiBins [expr {$i * $pi/18.0}]
262 }
263
264 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
265 add EtaPhiBins $eta $PhiBins
266 }
267
268
269 add EnergyFraction {0} {0.0}
270 # energy fractions for e, gamma and pi0
271 add EnergyFraction {11} {1.0}
272 add EnergyFraction {22} {1.0}
273 add EnergyFraction {111} {1.0}
274 # energy fractions for muon, neutrinos and neutralinos
275 add EnergyFraction {12} {0.0}
276 add EnergyFraction {13} {0.0}
277 add EnergyFraction {14} {0.0}
278 add EnergyFraction {16} {0.0}
279 add EnergyFraction {1000022} {0.0}
280 add EnergyFraction {1000023} {0.0}
281 add EnergyFraction {1000025} {0.0}
282 add EnergyFraction {1000035} {0.0}
283 add EnergyFraction {1000045} {0.0}
284 add EnergyFraction {5000522} {0.0}
285 add EnergyFraction {9000005} {0.0}
286 # energy fractions for K0short and Lambda
287 add EnergyFraction {310} {0.3}
288 add EnergyFraction {3122} {0.3}
289
290 # set ResolutionFormula {resolution formula as a function of eta and energy}
291
292 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
293
294 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
295 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
296 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
297 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
298 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
299
300}
301
302
303#############
304# HCAL
305#############
306
307module SimpleCalorimeter HCal {
308 set ParticleInputArray ParticlePropagator/stableParticles
309 set TrackInputArray ECal/eflowTracks
310
311 set TowerOutputArray hcalTowers
312 set EFlowTrackOutputArray eflowTracks
313 set EFlowTowerOutputArray eflowNeutralHadrons
314
315 set IsEcal false
316
317 set EnergyMin 1.0
318 set EnergySignificanceMin 1.0
319
320 set SmearTowerCenter true
321
322 set pi [expr {acos(-1)}]
323
324 # lists of the edges of each tower in eta and phi
325 # each list starts with the lower edge of the first tower
326 # the list ends with the higher edged of the last tower
327
328 # 5 degrees towers
329 set PhiBins {}
330 for {set i -36} {$i <= 36} {incr i} {
331 add PhiBins [expr {$i * $pi/36.0}]
332 }
333 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
334 add EtaPhiBins $eta $PhiBins
335 }
336
337 # 10 degrees towers
338 set PhiBins {}
339 for {set i -18} {$i <= 18} {incr i} {
340 add PhiBins [expr {$i * $pi/18.0}]
341 }
342 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
343 add EtaPhiBins $eta $PhiBins
344 }
345
346 # 20 degrees towers
347 set PhiBins {}
348 for {set i -9} {$i <= 9} {incr i} {
349 add PhiBins [expr {$i * $pi/9.0}]
350 }
351 foreach eta {-5 -4.7 -4.525 4.7 5} {
352 add EtaPhiBins $eta $PhiBins
353 }
354
355 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
356 add EnergyFraction {0} {1.0}
357 # energy fractions for e, gamma and pi0
358 add EnergyFraction {11} {0.0}
359 add EnergyFraction {22} {0.0}
360 add EnergyFraction {111} {0.0}
361 # energy fractions for muon, neutrinos and neutralinos
362 add EnergyFraction {12} {0.0}
363 add EnergyFraction {13} {0.0}
364 add EnergyFraction {14} {0.0}
365 add EnergyFraction {16} {0.0}
366 add EnergyFraction {1000022} {0.0}
367 add EnergyFraction {1000023} {0.0}
368 add EnergyFraction {1000025} {0.0}
369 add EnergyFraction {1000035} {0.0}
370 add EnergyFraction {1000045} {0.0}
371 add EnergyFraction {5000522} {0.0}
372 add EnergyFraction {9000005} {0.0}
373 # energy fractions for K0short and Lambda
374 add EnergyFraction {310} {0.7}
375 add EnergyFraction {3122} {0.7}
376
377 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
378 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
379 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
380
381}
382
383
384#################
385# Electron filter
386#################
387
388module PdgCodeFilter ElectronFilter {
389 set InputArray HCal/eflowTracks
390 set OutputArray electrons
391 set Invert true
392 add PdgCode {11}
393 add PdgCode {-11}
394}
395
396######################
397# ChargedHadronFilter
398######################
399
400module PdgCodeFilter ChargedHadronFilter {
401 set InputArray HCal/eflowTracks
402 set OutputArray chargedHadrons
403
404 add PdgCode {11}
405 add PdgCode {-11}
406 add PdgCode {13}
407 add PdgCode {-13}
408}
409
410
411###################################################
412# Tower Merger (in case not using e-flow algorithm)
413###################################################
414
415module Merger Calorimeter {
416# add InputArray InputArray
417 add InputArray ECal/ecalTowers
418 add InputArray HCal/hcalTowers
419 set OutputArray towers
420}
421
422
423
424####################
425# Energy flow merger
426####################
427
428module Merger EFlowMerger {
429# add InputArray InputArray
430 add InputArray HCal/eflowTracks
431 add InputArray ECal/eflowPhotons
432 add InputArray HCal/eflowNeutralHadrons
433 set OutputArray eflow
434}
435
436######################
437# EFlowFilter
438######################
439
440module PdgCodeFilter EFlowFilter {
441 set InputArray EFlowMerger/eflow
442 set OutputArray eflow
443
444 add PdgCode {11}
445 add PdgCode {-11}
446 add PdgCode {13}
447 add PdgCode {-13}
448}
449
450
451###################
452# Photon efficiency
453###################
454
455module Efficiency PhotonEfficiency {
456 set InputArray ECal/eflowPhotons
457 set OutputArray photons
458
459 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
460
461 # efficiency formula for photons
462 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
463 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
464 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
465 (abs(eta) > 2.5) * (0.00)}
466}
467
468
469#####################
470# Electron efficiency
471#####################
472
473module Efficiency ElectronEfficiency {
474 set InputArray ElectronFilter/electrons
475 set OutputArray electrons
476
477 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
478
479 # efficiency formula for electrons
480 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
481 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
482 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
483 (abs(eta) > 2.5) * (0.00)}
484}
485
486#################
487# Muon efficiency
488#################
489
490module Efficiency MuonEfficiency {
491 set InputArray MuonMomentumSmearing/muons
492 set OutputArray muons
493
494 # set EfficiencyFormula {efficiency as a function of eta and pt}
495
496 # efficiency formula for muons
497 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
498 (abs(eta) <= 0.9) * (0.96) +
499 (abs(eta) > 0.9 && abs(eta) < 2.4) * (0.97)}
500}
501
502
503###################
504# Missing ET merger
505###################
506
507module Merger MissingET {
508# add InputArray InputArray
509 add InputArray EFlowMerger/eflow
510 set MomentumOutputArray momentum
511}
512
513
514
515#####################
516# Neutrino Filter
517#####################
518
519module PdgCodeFilter NeutrinoFilter {
520
521 set InputArray Delphes/stableParticles
522 set OutputArray filteredParticles
523
524 set PTMin 0.0
525
526 add PdgCode {12}
527 add PdgCode {14}
528 add PdgCode {16}
529 add PdgCode {-12}
530 add PdgCode {-14}
531 add PdgCode {-16}
532
533}
534
535
536#####################
537# MC truth jet finder
538#####################
539
540module FastJetFinder GenJetFinder {
541 set InputArray NeutrinoFilter/filteredParticles
542
543 set OutputArray jets
544
545 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
546 set JetAlgorithm 6
547 set ParameterR 0.4
548
549 set JetPTMin 20.0
550}
551
552#########################
553# Gen Missing ET merger
554########################
555
556module Merger GenMissingET {
557# add InputArray InputArray
558 add InputArray NeutrinoFilter/filteredParticles
559 set MomentumOutputArray momentum
560}
561
562
563
564############
565# Jet finder
566############
567
568module FastJetFinder FastJetFinder {
569# set InputArray Calorimeter/towers
570 set InputArray EFlowMerger/eflow
571
572 set OutputArray jets
573
574 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
575 set JetAlgorithm 6
576 set ParameterR 0.4
577
578 set JetPTMin 20.0
579}
580
581##################
582# Jet Energy Scale
583##################
584
585module EnergyScale JetEnergyScale {
586 set InputArray FastJetFinder/jets
587 set OutputArray jets
588
589 # scale formula for jets
590 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
591}
592
593########################
594# Jet Flavor Association
595########################
596
597module JetFlavorAssociation JetFlavorAssociation {
598
599 set PartonInputArray Delphes/partons
600 set ParticleInputArray Delphes/allParticles
601 set ParticleLHEFInputArray Delphes/allParticlesLHEF
602 set JetInputArray JetEnergyScale/jets
603
604 set DeltaR 0.5
605 set PartonPTMin 1.0
606 set PartonEtaMax 2.5
607
608}
609
610###########
611# b-tagging
612###########
613
614module BTagging BTagging {
615 set JetInputArray JetEnergyScale/jets
616
617 set BitNumber 0
618
619 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
620 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
621 # gluon's PDG code has the lowest priority
622
623 # based on arXiv:1211.4462
624
625 # default efficiency formula (misidentification rate)
626 add EfficiencyFormula {0} {(pt > 20.0 && pt <= 150.0) * (0.00284 - 8.63*(pt/10^5) + 1.38*(pt^2/10^6) - 9.69*(pt^3/10^9) + 3.19*(pt^4/10^11) - 3.97*(pt^5/10^14)) +
627 (pt > 150.0 && pt <= 1000.0) * (0.00063 + 4.51*(pt/10^6) + 2.83*(pt^2/10^9))}
628
629 # efficiency formula for c-jets (misidentification rate)
630 add EfficiencyFormula {4} {(pt > 20.0 && pt <= 65.0) * (0.0234 - 8.417*(pt/10^5) + 1.24*(pt^2/10^6) - 5.5*(pt^3/10^9) + 9.96*(pt^4/10^12) - 6.32*(pt^5/10^15)) +
631 (pt > 65.0 && pt <= 1000.0) * (0.0218 + 2.46*(pt/10^5) - 2.021*(pt^2/10^8))}
632
633 # efficiency formula for b-jets
634 add EfficiencyFormula {5} {(pt > 20.0 && pt <= 50.0) * (-0.033 + 0.0225*pt - 0.00035*pt^2 + 2.586*(pt^3/10^6) - 9.096*(pt^4/10^9) + 1.212*(pt^5/10^11)) +
635 (pt > 50.0 && pt <= 160.0) * (0.169 + 0.013*pt - 0.00019*pt^2 + 1.373*(pt^3/10^6) - 4.923*(pt^4/10^9) + 6.87*(pt^5/10^12)) +
636 (pt > 160.0 && pt <= 1000.0) * (0.62 - 0.00083*pt + 4.3078*(pt^2/10^7))}
637}
638
639#############
640# tau-tagging
641#############
642
643module TauTagging TauTagging {
644 set ParticleInputArray Delphes/allParticles
645 set PartonInputArray Delphes/partons
646 set JetInputArray JetEnergyScale/jets
647
648 set DeltaR 0.5
649
650 set TauPTMin 1.0
651
652 set TauEtaMax 2.5
653
654 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
655
656 # default efficiency formula (misidentification rate)
657 add EfficiencyFormula {0} {0.01}
658 # efficiency formula for tau-jets
659 add EfficiencyFormula {15} {0.6}
660}
661
662
663##################
664# ROOT tree writer
665##################
666
667# tracks, towers and eflow objects are not stored by default in the output.
668# if needed (for jet constituent or other studies), uncomment the relevant
669# "add Branch ..." lines.
670
671module TreeWriter TreeWriter {
672# add Branch InputArray BranchName BranchClass
673 add Branch Delphes/allParticles Particle GenParticle
674
675 add Branch TrackMerger/tracks Track Track
676 add Branch Calorimeter/towers Tower Tower
677
678 add Branch HCal/eflowTracks EFlowTrack Track
679 add Branch ECal/eflowPhotons EFlowPhoton Tower
680 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
681
682 add Branch GenJetFinder/jets GenJet Jet
683 add Branch GenMissingET/momentum GenMissingET MissingET
684
685 add Branch JetEnergyScale/jets JetMA5 Jet
686 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
687 add Branch PhotonEfficiency/photons PhotonMA5 Photon
688 add Branch MuonEfficiency/muons MuonMA5 Muon
689
690 add Branch MissingET/momentum MissingET MissingET
691}