PublicAnalysisDatabase: delphes_card_cms_exo_16_022.tcl

File delphes_card_cms_exo_16_022.tcl, 21.2 KB (added by Eric CONTE, 5 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 ElectronTrackingEfficiencyD0
17 MuonTrackingEfficiencyD0
18
19 TrackMerger
20
21 ECal
22 HCal
23
24 Calorimeter
25 EFlowMerger
26 EFlowFilter
27
28 PhotonEfficiency
29
30 ElectronFilter
31 ElectronEfficiency
32
33 ChargedHadronFilter
34
35 MuonEfficiency
36
37 MissingET
38
39 NeutrinoFilter
40 GenJetFinder
41 GenMissingET
42
43 FastJetFinder
44 FatJetFinder
45
46 JetEnergyScale
47
48 JetFlavorAssociation
49
50 BTagging
51 TauTagging
52
53 TreeWriter
54}
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61 set InputArray Delphes/stableParticles
62
63 set OutputArray stableParticles
64 set ChargedHadronOutputArray chargedHadrons
65 set ElectronOutputArray electrons
66 set MuonOutputArray muons
67
68 # radius of the magnetic field coverage, in m
69 set Radius 1.29
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.00
72
73 # magnetic field
74 set Bz 3.8
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82 set InputArray ParticlePropagator/chargedHadrons
83 set OutputArray chargedHadrons
84
85 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87 # tracking efficiency formula for charged hadrons
88 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
89 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
90 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
91 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
92 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
93 (abs(eta) > 2.5) * (0.00)}
94}
95
96##############################
97# Electron tracking efficiency
98##############################
99
100module Efficiency ElectronTrackingEfficiency {
101 set InputArray ParticlePropagator/electrons
102 set OutputArray electrons
103
104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
105
106 # tracking efficiency formula for electrons
107 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
108 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
109 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
110 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
112 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
113 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
114 (abs(eta) > 2.5) * (0.00)}
115}
116
117##########################
118# Muon tracking efficiency
119##########################
120
121module Efficiency MuonTrackingEfficiency {
122 set InputArray ParticlePropagator/muons
123 set OutputArray muons
124
125 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for muons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
130 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
131 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
132
133 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
134 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
135 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
136 (abs(eta) > 2.5) * (0.00)}
137}
138
139########################################
140# Momentum resolution for charged tracks
141########################################
142
143module MomentumSmearing ChargedHadronMomentumSmearing {
144 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
145 set OutputArray chargedHadrons
146
147 # set ResolutionFormula {resolution formula as a function of eta and pt}
148
149 # resolution formula for charged hadrons
150 # based on arXiv:1405.6569
151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
152 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
153 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
154}
155
156###################################
157# Momentum resolution for electrons
158###################################
159
160module MomentumSmearing ElectronMomentumSmearing {
161 set InputArray ElectronTrackingEfficiency/electrons
162 set OutputArray electrons
163
164 # set ResolutionFormula {resolution formula as a function of eta and energy}
165
166 # resolution formula for electrons
167 # based on arXiv:1405.6569
168 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
169 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
170 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
171}
172
173###############################
174# Momentum resolution for muons
175###############################
176
177module MomentumSmearing MuonMomentumSmearing {
178 set InputArray MuonTrackingEfficiency/muons
179 set OutputArray muons
180
181 # set ResolutionFormula {resolution formula as a function of eta and pt}
182
183 # resolution formula for muons
184 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
185 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
187}
188
189##############################
190# Electron tracking efficiency
191##############################
192
193module Efficiency ElectronTrackingEfficiencyD0 {
194 set InputArray ElectronMomentumSmearing/electrons
195 set OutputArray electrons
196
197 set EfficiencyFormula { (d0<=20) * (-5.06107e-7 * d0**6 + 0.0000272756 * d0**5 - 0.00049321 * d0**4 + 0.00287189 * d0**3 + 0.00522007 * d0**2 -0.0917957 * d0 + 0.924921) + (d0>20) * (0.00) }
198}
199
200
201
202##############################
203# Muon tracking efficiency
204##############################
205
206module Efficiency MuonTrackingEfficiencyD0 {
207 set InputArray MuonMomentumSmearing/muons
208 set OutputArray muons
209
210 set EfficiencyFormula { (d0<=20) * (-1.11279e-7 * d0**6 + 4.23266e-6 * d0**5 - 0.0000260292 * d0**4 - 0.000611108 * d0**3 + 0.00743217 * d0**2 - 0.0271852 * d0 + 0.99067) + (d0>20) * (0.00) }
211
212}
213
214
215
216##############
217# Track merger
218##############
219
220module Merger TrackMerger {
221# add InputArray InputArray
222 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
223 add InputArray ElectronTrackingEfficiencyD0/electrons
224 add InputArray MuonTrackingEfficiencyD0/muons
225 set OutputArray tracks
226}
227
228
229
230#############
231# ECAL
232#############
233
234module SimpleCalorimeter ECal {
235 set ParticleInputArray ParticlePropagator/stableParticles
236 set TrackInputArray TrackMerger/tracks
237
238 set TowerOutputArray ecalTowers
239 set EFlowTrackOutputArray eflowTracks
240 set EFlowTowerOutputArray eflowPhotons
241
242 set IsEcal true
243
244 set EnergyMin 0.5
245 set EnergySignificanceMin 2.0
246
247 set SmearTowerCenter true
248
249 set pi [expr {acos(-1)}]
250
251 # lists of the edges of each tower in eta and phi
252 # each list starts with the lower edge of the first tower
253 # the list ends with the higher edged of the last tower
254
255 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
256
257 set PhiBins {}
258 for {set i -180} {$i <= 180} {incr i} {
259 add PhiBins [expr {$i * $pi/180.0}]
260 }
261
262 # 0.02 unit in eta up to eta = 1.5 (barrel)
263 for {set i -85} {$i <= 86} {incr i} {
264 set eta [expr {$i * 0.0174}]
265 add EtaPhiBins $eta $PhiBins
266 }
267
268 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
269
270 set PhiBins {}
271 for {set i -180} {$i <= 180} {incr i} {
272 add PhiBins [expr {$i * $pi/180.0}]
273 }
274
275 # 0.02 unit in eta up to eta = 3
276 for {set i 1} {$i <= 84} {incr i} {
277 set eta [expr { -2.958 + $i * 0.0174}]
278 add EtaPhiBins $eta $PhiBins
279 }
280
281 for {set i 1} {$i <= 84} {incr i} {
282 set eta [expr { 1.4964 + $i * 0.0174}]
283 add EtaPhiBins $eta $PhiBins
284 }
285
286 # take present CMS granularity for HF
287
288 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
289 set PhiBins {}
290 for {set i -18} {$i <= 18} {incr i} {
291 add PhiBins [expr {$i * $pi/18.0}]
292 }
293
294 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} {
295 add EtaPhiBins $eta $PhiBins
296 }
297
298
299 add EnergyFraction {0} {0.0}
300 # energy fractions for e, gamma and pi0
301 add EnergyFraction {11} {1.0}
302 add EnergyFraction {22} {1.0}
303 add EnergyFraction {111} {1.0}
304 # energy fractions for muon, neutrinos and neutralinos
305 add EnergyFraction {12} {0.0}
306 add EnergyFraction {13} {0.0}
307 add EnergyFraction {14} {0.0}
308 add EnergyFraction {16} {0.0}
309 add EnergyFraction {1000022} {0.0}
310 add EnergyFraction {1000023} {0.0}
311 add EnergyFraction {1000025} {0.0}
312 add EnergyFraction {1000035} {0.0}
313 add EnergyFraction {1000045} {0.0}
314 # energy fractions for K0short and Lambda
315 add EnergyFraction {310} {0.3}
316 add EnergyFraction {3122} {0.3}
317
318 # set ResolutionFormula {resolution formula as a function of eta and energy}
319
320 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
321
322 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
323 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
324 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
325 (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) +
326 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
327
328}
329
330
331#############
332# HCAL
333#############
334
335module SimpleCalorimeter HCal {
336 set ParticleInputArray ParticlePropagator/stableParticles
337 set TrackInputArray ECal/eflowTracks
338
339 set TowerOutputArray hcalTowers
340 set EFlowTrackOutputArray eflowTracks
341 set EFlowTowerOutputArray eflowNeutralHadrons
342
343 set IsEcal false
344
345 set EnergyMin 1.0
346 set EnergySignificanceMin 1.0
347
348 set SmearTowerCenter true
349
350 set pi [expr {acos(-1)}]
351
352 # lists of the edges of each tower in eta and phi
353 # each list starts with the lower edge of the first tower
354 # the list ends with the higher edged of the last tower
355
356 # 5 degrees towers
357 set PhiBins {}
358 for {set i -36} {$i <= 36} {incr i} {
359 add PhiBins [expr {$i * $pi/36.0}]
360 }
361 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} {
362 add EtaPhiBins $eta $PhiBins
363 }
364
365 # 10 degrees towers
366 set PhiBins {}
367 for {set i -18} {$i <= 18} {incr i} {
368 add PhiBins [expr {$i * $pi/18.0}]
369 }
370 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} {
371 add EtaPhiBins $eta $PhiBins
372 }
373
374 # 20 degrees towers
375 set PhiBins {}
376 for {set i -9} {$i <= 9} {incr i} {
377 add PhiBins [expr {$i * $pi/9.0}]
378 }
379 foreach eta {-5 -4.7 -4.525 4.7 5} {
380 add EtaPhiBins $eta $PhiBins
381 }
382
383 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
384 add EnergyFraction {0} {1.0}
385 # energy fractions for e, gamma and pi0
386 add EnergyFraction {11} {0.0}
387 add EnergyFraction {22} {0.0}
388 add EnergyFraction {111} {0.0}
389 # energy fractions for muon, neutrinos and neutralinos
390 add EnergyFraction {12} {0.0}
391 add EnergyFraction {13} {0.0}
392 add EnergyFraction {14} {0.0}
393 add EnergyFraction {16} {0.0}
394 add EnergyFraction {1000022} {0.0}
395 add EnergyFraction {1000023} {0.0}
396 add EnergyFraction {1000025} {0.0}
397 add EnergyFraction {1000035} {0.0}
398 add EnergyFraction {1000045} {0.0}
399 # energy fractions for K0short and Lambda
400 add EnergyFraction {310} {0.7}
401 add EnergyFraction {3122} {0.7}
402
403 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
404 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
405 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
406
407}
408
409
410#################
411# Electron filter
412#################
413
414module PdgCodeFilter ElectronFilter {
415 set InputArray HCal/eflowTracks
416 set OutputArray electrons
417 set Invert true
418 add PdgCode {11}
419 add PdgCode {-11}
420}
421
422######################
423# ChargedHadronFilter
424######################
425
426module PdgCodeFilter ChargedHadronFilter {
427 set InputArray HCal/eflowTracks
428 set OutputArray chargedHadrons
429
430 add PdgCode {11}
431 add PdgCode {-11}
432 add PdgCode {13}
433 add PdgCode {-13}
434}
435
436
437###################################################
438# Tower Merger (in case not using e-flow algorithm)
439###################################################
440
441module Merger Calorimeter {
442# add InputArray InputArray
443 add InputArray ECal/ecalTowers
444 add InputArray HCal/hcalTowers
445 set OutputArray towers
446}
447
448
449
450####################
451# Energy flow merger
452####################
453
454module Merger EFlowMerger {
455# add InputArray InputArray
456 add InputArray HCal/eflowTracks
457 add InputArray ECal/eflowPhotons
458 add InputArray HCal/eflowNeutralHadrons
459 set OutputArray eflow
460}
461
462######################
463# EFlowFilter
464######################
465
466module PdgCodeFilter EFlowFilter {
467 set InputArray EFlowMerger/eflow
468 set OutputArray eflow
469
470 add PdgCode {11}
471 add PdgCode {-11}
472 add PdgCode {13}
473 add PdgCode {-13}
474}
475
476
477###################
478# Photon efficiency
479###################
480
481module Efficiency PhotonEfficiency {
482 set InputArray ECal/eflowPhotons
483 set OutputArray photons
484
485 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
486
487 # efficiency formula for photons
488 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
489 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
490 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
491 (abs(eta) > 2.5) * (0.00)}
492}
493
494
495#####################
496# Electron efficiency
497#####################
498
499module Efficiency ElectronEfficiency {
500 set InputArray ElectronFilter/electrons
501 set OutputArray electrons
502
503 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
504
505 # efficiency formula for electrons
506 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
507 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
508 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
509 (abs(eta) > 2.5) * (0.00)}
510}
511
512#################
513# Muon efficiency
514#################
515
516module Efficiency MuonEfficiency {
517 set InputArray MuonTrackingEfficiencyD0/muons
518 set OutputArray muons
519
520 # set EfficiencyFormula {efficiency as a function of eta and pt}
521
522 # efficiency formula for muons
523 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
524 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
525 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
526 (abs(eta) > 2.4) * (0.00)}
527}
528
529
530###################
531# Missing ET merger
532###################
533
534module Merger MissingET {
535# add InputArray InputArray
536 add InputArray EFlowMerger/eflow
537 set MomentumOutputArray momentum
538}
539
540
541
542#####################
543# Neutrino Filter
544#####################
545
546module PdgCodeFilter NeutrinoFilter {
547
548 set InputArray Delphes/stableParticles
549 set OutputArray filteredParticles
550
551 set PTMin 0.0
552
553 add PdgCode {12}
554 add PdgCode {14}
555 add PdgCode {16}
556 add PdgCode {-12}
557 add PdgCode {-14}
558 add PdgCode {-16}
559
560}
561
562
563#####################
564# MC truth jet finder
565#####################
566
567module FastJetFinder GenJetFinder {
568 set InputArray NeutrinoFilter/filteredParticles
569
570 set OutputArray jets
571
572 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
573 set JetAlgorithm 6
574 set ParameterR 0.5
575
576 set JetPTMin 20.0
577}
578
579#########################
580# Gen Missing ET merger
581########################
582
583module Merger GenMissingET {
584# add InputArray InputArray
585 add InputArray NeutrinoFilter/filteredParticles
586 set MomentumOutputArray momentum
587}
588
589
590
591############
592# Jet finder
593############
594
595module FastJetFinder FastJetFinder {
596# set InputArray Calorimeter/towers
597 set InputArray EFlowMerger/eflow
598
599 set OutputArray jets
600
601 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
602 set JetAlgorithm 6
603 set ParameterR 0.5
604
605 set JetPTMin 20.0
606}
607
608##################
609# Fat Jet finder
610##################
611
612module FastJetFinder FatJetFinder {
613 set InputArray EFlowMerger/eflow
614
615 set OutputArray jets
616
617 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
618 set JetAlgorithm 6
619 set ParameterR 0.8
620
621 set ComputeNsubjettiness 1
622 set Beta 1.0
623 set AxisMode 4
624
625 set ComputeTrimming 1
626 set RTrim 0.2
627 set PtFracTrim 0.05
628
629 set ComputePruning 1
630 set ZcutPrun 0.1
631 set RcutPrun 0.5
632 set RPrun 0.8
633
634 set ComputeSoftDrop 1
635 set BetaSoftDrop 0.0
636 set SymmetryCutSoftDrop 0.1
637 set R0SoftDrop 0.8
638
639 set JetPTMin 200.0
640}
641
642
643
644
645##################
646# Jet Energy Scale
647##################
648
649module EnergyScale JetEnergyScale {
650 set InputArray FastJetFinder/jets
651 set OutputArray jets
652
653 # scale formula for jets
654 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
655}
656
657########################
658# Jet Flavor Association
659########################
660
661module JetFlavorAssociation JetFlavorAssociation {
662
663 set PartonInputArray Delphes/partons
664 set ParticleInputArray Delphes/allParticles
665 set ParticleLHEFInputArray Delphes/allParticlesLHEF
666 set JetInputArray JetEnergyScale/jets
667
668 set DeltaR 0.5
669 set PartonPTMin 1.0
670 set PartonEtaMax 2.5
671
672}
673
674###########
675# b-tagging
676###########
677
678module BTagging BTagging {
679 set JetInputArray JetEnergyScale/jets
680
681 set BitNumber 0
682
683 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
684 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
685 # gluon's PDG code has the lowest priority
686
687 # based on arXiv:1211.4462
688
689 # default efficiency formula (misidentification rate)
690 add EfficiencyFormula {0} {0.01+0.000038*pt}
691
692 # efficiency formula for c-jets (misidentification rate)
693 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
694
695 # efficiency formula for b-jets
696 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
697}
698
699#############
700# tau-tagging
701#############
702
703module TauTagging TauTagging {
704 set ParticleInputArray Delphes/allParticles
705 set PartonInputArray Delphes/partons
706 set JetInputArray JetEnergyScale/jets
707
708 set DeltaR 0.5
709
710 set TauPTMin 1.0
711
712 set TauEtaMax 2.5
713
714 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
715
716 # default efficiency formula (misidentification rate)
717 add EfficiencyFormula {0} {0.01}
718 # efficiency formula for tau-jets
719 add EfficiencyFormula {15} {0.6}
720}
721
722
723##################
724# ROOT tree writer
725##################
726
727# tracks, towers and eflow objects are not stored by default in the output.
728# if needed (for jet constituent or other studies), uncomment the relevant
729# "add Branch ..." lines.
730
731module TreeWriter TreeWriter {
732# add Branch InputArray BranchName BranchClass
733 add Branch Delphes/allParticles Particle GenParticle
734
735 add Branch TrackMerger/tracks Track Track
736 add Branch Calorimeter/towers Tower Tower
737
738 add Branch HCal/eflowTracks EFlowTrack Track
739 add Branch ECal/eflowPhotons EFlowPhoton Tower
740 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
741
742 add Branch GenJetFinder/jets GenJet Jet
743 add Branch GenMissingET/momentum GenMissingET MissingET
744
745 add Branch JetEnergyScale/jets JetMA5 Jet
746 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
747 add Branch PhotonEfficiency/photons PhotonMA5 Photon
748 add Branch MuonEfficiency/muons MuonMA5 Muon
749
750 add Branch FatJetFinder/jets FatJet Jet
751
752 add Branch MissingET/momentum MissingET MissingET
753}