Korea2017-04-CMS-SUS-16-041: delphes_card_cms_exo_16_010.tcl

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