Fork me on GitHub

Ticket #1038: delphes_card_CMS_nc.tcl

File delphes_card_CMS_nc.tcl, 22.0 KB (added by qzom190, 8 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
24 PhotonEfficiency
25 PhotonIsolation
26
27 ElectronFilter
28 ElectronEfficiency
29 ElectronIsolation
30
31 MuonEfficiency
32 MuonIsolation
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 UniqueObjectFinder
50
51 ScalarHT
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# Track merger
191##############
192
193module Merger TrackMerger {
194# add InputArray InputArray
195 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
196 add InputArray ElectronMomentumSmearing/electrons
197 add InputArray MuonMomentumSmearing/muons
198 set OutputArray tracks
199}
200
201
202
203#############
204# ECAL
205#############
206
207module SimpleCalorimeter ECal {
208 set ParticleInputArray ParticlePropagator/stableParticles
209 set TrackInputArray TrackMerger/tracks
210
211 set TowerOutputArray ecalTowers
212 set EFlowTrackOutputArray eflowTracks
213 set EFlowTowerOutputArray eflowPhotons
214
215 set IsEcal true
216
217 set EnergyMin 0.5
218 set EnergySignificanceMin 2.0
219
220 set SmearTowerCenter true
221
222 set pi [expr {acos(-1)}]
223
224 # lists of the edges of each tower in eta and phi
225 # each list starts with the lower edge of the first tower
226 # the list ends with the higher edged of the last tower
227
228 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
229
230 set PhiBins {}
231 for {set i -180} {$i <= 180} {incr i} {
232 add PhiBins [expr {$i * $pi/180.0}]
233 }
234
235 # 0.02 unit in eta up to eta = 1.5 (barrel)
236 for {set i -85} {$i <= 86} {incr i} {
237 set eta [expr {$i * 0.0174}]
238 add EtaPhiBins $eta $PhiBins
239 }
240
241 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
242
243 set PhiBins {}
244 for {set i -180} {$i <= 180} {incr i} {
245 add PhiBins [expr {$i * $pi/180.0}]
246 }
247
248 # 0.02 unit in eta up to eta = 3
249 for {set i 1} {$i <= 84} {incr i} {
250 set eta [expr { -2.958 + $i * 0.0174}]
251 add EtaPhiBins $eta $PhiBins
252 }
253
254 for {set i 1} {$i <= 84} {incr i} {
255 set eta [expr { 1.4964 + $i * 0.0174}]
256 add EtaPhiBins $eta $PhiBins
257 }
258
259 # take present CMS granularity for HF
260
261 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
262 set PhiBins {}
263 for {set i -18} {$i <= 18} {incr i} {
264 add PhiBins [expr {$i * $pi/18.0}]
265 }
266
267 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} {
268 add EtaPhiBins $eta $PhiBins
269 }
270
271
272 add EnergyFraction {0} {0.0}
273 # energy fractions for e, gamma and pi0
274 add EnergyFraction {11} {1.0}
275 add EnergyFraction {22} {1.0}
276 add EnergyFraction {111} {1.0}
277 # energy fractions for muon, neutrinos and neutralinos
278 add EnergyFraction {12} {0.0}
279 add EnergyFraction {13} {0.0}
280 add EnergyFraction {14} {0.0}
281 add EnergyFraction {16} {0.0}
282 add EnergyFraction {1000022} {0.0}
283 add EnergyFraction {1000023} {0.0}
284 add EnergyFraction {1000025} {0.0}
285 add EnergyFraction {1000035} {0.0}
286 add EnergyFraction {1000045} {0.0}
287 # energy fractions for K0short and Lambda
288 add EnergyFraction {310} {0.3}
289 add EnergyFraction {3122} {0.3}
290
291 # set ResolutionFormula {resolution formula as a function of eta and energy}
292
293 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
294
295 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
296 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
297 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
298 (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) +
299 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
300
301}
302
303
304#############
305# HCAL
306#############
307
308module SimpleCalorimeter HCal {
309 set ParticleInputArray ParticlePropagator/stableParticles
310 set TrackInputArray ECal/eflowTracks
311
312 set TowerOutputArray hcalTowers
313 set EFlowTrackOutputArray eflowTracks
314 set EFlowTowerOutputArray eflowNeutralHadrons
315
316 set IsEcal false
317
318 set EnergyMin 1.0
319 set EnergySignificanceMin 2.0
320
321 set SmearTowerCenter true
322
323 set pi [expr {acos(-1)}]
324
325 # lists of the edges of each tower in eta and phi
326 # each list starts with the lower edge of the first tower
327 # the list ends with the higher edged of the last tower
328 set PhiBins {}
329 for {set i -360} {$i <= 360} {incr i} {
330 add PhiBins [expr {$i * $pi/360.0}]
331 }
332
333 # 0.01 unit in eta up to eta = 2.5
334 for {set i -1000} {$i <= 1000} {incr i} {
335 set eta [expr {$i * 0.005}]
336 add EtaPhiBins $eta $PhiBins
337 }
338
339 # 5 degrees towers
340 set PhiBins {}
341 for {set i -36} {$i <= 36} {incr i} {
342 add PhiBins [expr {$i * $pi/36.0}]
343 }
344 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} {
345 add EtaPhiBins $eta $PhiBins
346 }
347
348 # 10 degrees towers
349 set PhiBins {}
350 for {set i -18} {$i <= 18} {incr i} {
351 add PhiBins [expr {$i * $pi/18.0}]
352 }
353 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} {
354 add EtaPhiBins $eta $PhiBins
355 }
356
357 # 20 degrees towers
358 set PhiBins {}
359 for {set i -9} {$i <= 9} {incr i} {
360 add PhiBins [expr {$i * $pi/9.0}]
361 }
362 foreach eta {-5 -4.7 -4.525 4.7 5} {
363 add EtaPhiBins $eta $PhiBins
364 }
365
366 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
367 add EnergyFraction {0} {1.0}
368 # energy fractions for e, gamma and pi0
369 add EnergyFraction {11} {0.0}
370 add EnergyFraction {22} {0.0}
371 add EnergyFraction {111} {0.0}
372 # energy fractions for muon, neutrinos and neutralinos
373 add EnergyFraction {12} {0.0}
374 add EnergyFraction {13} {0.0}
375 add EnergyFraction {14} {0.0}
376 add EnergyFraction {16} {0.0}
377 add EnergyFraction {1000022} {0.0}
378 add EnergyFraction {1000023} {0.0}
379 add EnergyFraction {1000025} {0.0}
380 add EnergyFraction {1000035} {0.0}
381 add EnergyFraction {1000045} {0.0}
382 # energy fractions for K0short and Lambda
383 add EnergyFraction {310} {0.7}
384 add EnergyFraction {3122} {0.7}
385
386 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
387 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
388 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
389
390 set ECalResolutionFormula { (abs(eta) <= 1.5) *0.05*(1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
391 (abs(eta) > 1.5 && abs(eta) <= 2.5) *0.05* (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
392 (abs(eta) > 2.5 && abs(eta) <= 5.0) *0.05* sqrt(energy^2*0.107^2 + energy*2.08^2)}
393
394 set HCalResolutionFormula { (abs(eta) <= 3.0) *0.05* sqrt(energy^2*0.050^2 + energy*1.50^2) +
395 (abs(eta) > 3.0 && abs(eta) <= 5.0) *0.05* sqrt(energy^2*0.130^2 + energy*2.70^2)}
396
397}
398
399
400#################
401# Electron filter
402#################
403
404module PdgCodeFilter ElectronFilter {
405 set InputArray HCal/eflowTracks
406 set OutputArray electrons
407 set Invert true
408 add PdgCode {11}
409 add PdgCode {-11}
410}
411
412###################################################
413# Tower Merger (in case not using e-flow algorithm)
414###################################################
415
416module Merger Calorimeter {
417# add InputArray InputArray
418 add InputArray ECal/ecalTowers
419 add InputArray HCal/hcalTowers
420 set OutputArray towers
421}
422
423
424
425####################
426# Energy flow merger
427####################
428
429module Merger EFlowMerger {
430# add InputArray InputArray
431 add InputArray HCal/eflowTracks
432 add InputArray ECal/eflowPhotons
433 add InputArray HCal/eflowNeutralHadrons
434 set OutputArray eflow
435}
436
437###################
438# Photon efficiency
439###################
440
441module Efficiency PhotonEfficiency {
442 set InputArray ECal/eflowPhotons
443 set OutputArray photons
444
445 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
446
447 # efficiency formula for photons
448 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
449 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
450 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
451 (abs(eta) > 2.5) * (0.00)}
452}
453
454##################
455# Photon isolation
456##################
457
458module Isolation PhotonIsolation {
459 set CandidateInputArray PhotonEfficiency/photons
460 set IsolationInputArray EFlowMerger/eflow
461
462 set OutputArray photons
463
464 set DeltaRMax 0.5
465
466 set PTMin 0.5
467
468 set PTRatioMax 0.12
469}
470
471
472#####################
473# Electron efficiency
474#####################
475
476module Efficiency ElectronEfficiency {
477 set InputArray ElectronFilter/electrons
478 set OutputArray electrons
479
480 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
481
482 # efficiency formula for electrons
483 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
484 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
485 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
486 (abs(eta) > 2.5) * (0.00)}
487}
488
489####################
490# Electron isolation
491####################
492
493module Isolation ElectronIsolation {
494 set CandidateInputArray ElectronEfficiency/electrons
495 set IsolationInputArray EFlowMerger/eflow
496
497 set OutputArray electrons
498
499 set DeltaRMax 0.5
500
501 set PTMin 0.5
502
503 set PTRatioMax 0.12
504}
505
506#################
507# Muon efficiency
508#################
509
510module Efficiency MuonEfficiency {
511 set InputArray MuonMomentumSmearing/muons
512 set OutputArray muons
513
514 # set EfficiencyFormula {efficiency as a function of eta and pt}
515
516 # efficiency formula for muons
517 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
518 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
519 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
520 (abs(eta) > 2.4) * (0.00)}
521}
522
523################
524# Muon isolation
525################
526
527module Isolation MuonIsolation {
528 set CandidateInputArray MuonEfficiency/muons
529 set IsolationInputArray EFlowMerger/eflow
530
531 set OutputArray muons
532
533 set DeltaRMax 0.5
534
535 set PTMin 0.5
536
537 set PTRatioMax 0.25
538}
539
540###################
541# Missing ET merger
542###################
543
544module Merger MissingET {
545# add InputArray InputArray
546 add InputArray EFlowMerger/eflow
547 set MomentumOutputArray momentum
548}
549
550##################
551# Scalar HT merger
552##################
553
554module Merger ScalarHT {
555# add InputArray InputArray
556 add InputArray UniqueObjectFinder/jets
557 add InputArray UniqueObjectFinder/electrons
558 add InputArray UniqueObjectFinder/photons
559 add InputArray UniqueObjectFinder/muons
560 set EnergyOutputArray energy
561}
562
563
564#####################
565# Neutrino Filter
566#####################
567
568module PdgCodeFilter NeutrinoFilter {
569
570 set InputArray Delphes/stableParticles
571 set OutputArray filteredParticles
572
573 set PTMin 0.0
574
575 add PdgCode {12}
576 add PdgCode {14}
577 add PdgCode {16}
578 add PdgCode {-12}
579 add PdgCode {-14}
580 add PdgCode {-16}
581
582}
583
584
585#####################
586# MC truth jet finder
587#####################
588
589module FastJetFinder GenJetFinder {
590 set InputArray NeutrinoFilter/filteredParticles
591
592 set OutputArray jets
593
594 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
595 set JetAlgorithm 6
596 set ParameterR 0.5
597
598 set JetPTMin 20.0
599}
600
601#########################
602# Gen Missing ET merger
603########################
604
605module Merger GenMissingET {
606# add InputArray InputArray
607 add InputArray NeutrinoFilter/filteredParticles
608 set MomentumOutputArray momentum
609}
610
611
612
613#############
614# Jet finder
615#############
616
617module FastJetFinder FastJetFinder {
618 set InputArray EFlowMerger/eflow
619
620 set OutputArray jets
621
622 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
623 set JetAlgorithm 5
624 set ParameterR 0.8
625
626 set ComputeNsubjettiness 1
627 set Beta 1.0
628 set AxisMode 4
629
630 set ComputeTrimming 1
631 set RTrim 0.2
632 set PtFracTrim 0.05
633
634 set ComputePruning 1
635 set ZcutPrun 0.1
636 set RcutPrun 0.5
637 set RPrun 0.8
638
639 set ComputeSoftDrop 1
640 set BetaSoftDrop 0.0
641 set SymmetryCutSoftDrop 0.1
642 set R0SoftDrop 0.8
643
644 set JetPTMin 20.0
645}
646
647##################
648# Jet Energy Scale
649##################
650
651module EnergyScale JetEnergyScale {
652 set InputArray FastJetFinder/jets
653 set OutputArray jets
654
655 # scale formula for jets
656 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
657}
658
659########################
660# Jet Flavor Association
661########################
662
663module JetFlavorAssociation JetFlavorAssociation {
664
665 set PartonInputArray Delphes/partons
666 set ParticleInputArray Delphes/allParticles
667 set ParticleLHEFInputArray Delphes/allParticlesLHEF
668 set JetInputArray JetEnergyScale/jets
669
670 set DeltaR 0.5
671 set PartonPTMin 1.0
672 set PartonEtaMax 2.5
673
674}
675
676###########
677# b-tagging
678###########
679
680module BTagging BTagging {
681 set JetInputArray JetEnergyScale/jets
682
683 set BitNumber 0
684
685 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
686 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
687 # gluon's PDG code has the lowest priority
688
689 # based on arXiv:1211.4462
690
691 # default efficiency formula (misidentification rate)
692 add EfficiencyFormula {0} {0.01+0.000038*pt}
693
694 # efficiency formula for c-jets (misidentification rate)
695 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
696
697 # efficiency formula for b-jets
698 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
699}
700
701#############
702# tau-tagging
703#############
704
705module TauTagging TauTagging {
706 set ParticleInputArray Delphes/allParticles
707 set PartonInputArray Delphes/partons
708 set JetInputArray JetEnergyScale/jets
709
710 set DeltaR 0.5
711
712 set TauPTMin 1.0
713
714 set TauEtaMax 2.5
715
716 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
717
718 # default efficiency formula (misidentification rate)
719 add EfficiencyFormula {0} {0.01}
720 # efficiency formula for tau-jets
721 add EfficiencyFormula {15} {0.6}
722}
723
724#####################################################
725# Find uniquely identified photons/electrons/tau/jets
726#####################################################
727
728module UniqueObjectFinder UniqueObjectFinder {
729# earlier arrays take precedence over later ones
730# add InputArray InputArray OutputArray
731 add InputArray PhotonIsolation/photons photons
732 add InputArray ElectronIsolation/electrons electrons
733 add InputArray MuonIsolation/muons muons
734 add InputArray JetEnergyScale/jets jets
735}
736
737##################
738# ROOT tree writer
739##################
740
741# tracks, towers and eflow objects are not stored by default in the output.
742# if needed (for jet constituent or other studies), uncomment the relevant
743# "add Branch ..." lines.
744
745module TreeWriter TreeWriter {
746# add Branch InputArray BranchName BranchClass
747 add Branch Delphes/allParticles Particle GenParticle
748
749 add Branch TrackMerger/tracks Track Track
750 add Branch Calorimeter/towers Tower Tower
751
752 add Branch HCal/eflowTracks EFlowTrack Track
753 add Branch ECal/eflowPhotons EFlowPhoton Tower
754 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
755
756 add Branch GenJetFinder/jets GenJet Jet
757 add Branch GenMissingET/momentum GenMissingET MissingET
758
759 add Branch UniqueObjectFinder/jets Jet Jet
760 add Branch UniqueObjectFinder/electrons Electron Electron
761 add Branch UniqueObjectFinder/photons Photon Photon
762 add Branch UniqueObjectFinder/muons Muon Muon
763 add Branch MissingET/momentum MissingET MissingET
764 add Branch ScalarHT/energy ScalarHT ScalarHT
765}