MA5SandBox: delphes_card_cms_exo_16_012.tcl

File delphes_card_cms_exo_16_012.tcl, 22.4 KB (added by Benjamin Fuks, 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 {18} {0.0}
287 add EnergyFraction {1000022} {0.0}
288 add EnergyFraction {1000023} {0.0}
289 add EnergyFraction {1000025} {0.0}
290 add EnergyFraction {1000035} {0.0}
291 add EnergyFraction {1000045} {0.0}
292 # energy fractions for K0short and Lambda
293 add EnergyFraction {310} {0.3}
294 add EnergyFraction {3122} {0.3}
295
296 # set ResolutionFormula {resolution formula as a function of eta and energy}
297
298 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
299
300 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
301 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
302 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
303 (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) +
304 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
305
306}
307
308
309#############
310# HCAL
311#############
312
313module SimpleCalorimeter HCal {
314 set ParticleInputArray ParticlePropagator/stableParticles
315 set TrackInputArray ECal/eflowTracks
316
317 set TowerOutputArray hcalTowers
318 set EFlowTrackOutputArray eflowTracks
319 set EFlowTowerOutputArray eflowNeutralHadrons
320
321 set IsEcal false
322
323 set EnergyMin 1.0
324 set EnergySignificanceMin 1.0
325
326 set SmearTowerCenter true
327
328 set pi [expr {acos(-1)}]
329
330 # lists of the edges of each tower in eta and phi
331 # each list starts with the lower edge of the first tower
332 # the list ends with the higher edged of the last tower
333
334 # 5 degrees towers
335 set PhiBins {}
336 for {set i -36} {$i <= 36} {incr i} {
337 add PhiBins [expr {$i * $pi/36.0}]
338 }
339 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} {
340 add EtaPhiBins $eta $PhiBins
341 }
342
343 # 10 degrees towers
344 set PhiBins {}
345 for {set i -18} {$i <= 18} {incr i} {
346 add PhiBins [expr {$i * $pi/18.0}]
347 }
348 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} {
349 add EtaPhiBins $eta $PhiBins
350 }
351
352 # 20 degrees towers
353 set PhiBins {}
354 for {set i -9} {$i <= 9} {incr i} {
355 add PhiBins [expr {$i * $pi/9.0}]
356 }
357 foreach eta {-5 -4.7 -4.525 4.7 5} {
358 add EtaPhiBins $eta $PhiBins
359 }
360
361 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
362 add EnergyFraction {0} {1.0}
363 # energy fractions for e, gamma and pi0
364 add EnergyFraction {11} {0.0}
365 add EnergyFraction {22} {0.0}
366 add EnergyFraction {111} {0.0}
367 # energy fractions for muon, neutrinos and neutralinos
368 add EnergyFraction {12} {0.0}
369 add EnergyFraction {13} {0.0}
370 add EnergyFraction {14} {0.0}
371 add EnergyFraction {16} {0.0}
372 add EnergyFraction {18} {0.0}
373 add EnergyFraction {1000022} {0.0}
374 add EnergyFraction {1000023} {0.0}
375 add EnergyFraction {1000025} {0.0}
376 add EnergyFraction {1000035} {0.0}
377 add EnergyFraction {1000045} {0.0}
378 # energy fractions for K0short and Lambda
379 add EnergyFraction {310} {0.7}
380 add EnergyFraction {3122} {0.7}
381
382 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
383 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
384 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
385
386}
387
388
389#################
390# Electron filter
391#################
392
393module PdgCodeFilter ElectronFilter {
394 set InputArray HCal/eflowTracks
395 set OutputArray electrons
396 set Invert true
397 add PdgCode {11}
398 add PdgCode {-11}
399}
400
401######################
402# ChargedHadronFilter
403######################
404
405module PdgCodeFilter ChargedHadronFilter {
406 set InputArray HCal/eflowTracks
407 set OutputArray chargedHadrons
408
409 add PdgCode {11}
410 add PdgCode {-11}
411 add PdgCode {13}
412 add PdgCode {-13}
413}
414
415
416###################################################
417# Tower Merger (in case not using e-flow algorithm)
418###################################################
419
420module Merger Calorimeter {
421# add InputArray InputArray
422 add InputArray ECal/ecalTowers
423 add InputArray HCal/hcalTowers
424 set OutputArray towers
425}
426
427
428
429####################
430# Energy flow merger
431####################
432
433module Merger EFlowMerger {
434# add InputArray InputArray
435 add InputArray HCal/eflowTracks
436 add InputArray ECal/eflowPhotons
437 add InputArray HCal/eflowNeutralHadrons
438 set OutputArray eflow
439}
440
441######################
442# EFlowFilter
443######################
444
445module PdgCodeFilter EFlowFilter {
446 set InputArray EFlowMerger/eflow
447 set OutputArray eflow
448
449 add PdgCode {11}
450 add PdgCode {-11}
451 add PdgCode {13}
452 add PdgCode {-13}
453}
454
455
456###################
457# Photon efficiency
458###################
459
460module Efficiency PhotonEfficiency {
461 set InputArray ECal/eflowPhotons
462 set OutputArray photons
463
464 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
465
466 # efficiency formula for photons
467 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
468 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
469 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
470 (abs(eta) > 2.5) * (0.00)}
471}
472
473##################
474# Photon isolation
475##################
476
477module Isolation PhotonIsolation {
478 set CandidateInputArray PhotonEfficiency/photons
479 set IsolationInputArray EFlowFilter/eflow
480
481 set OutputArray photons
482
483 set DeltaRMax 0.3
484
485 set PTMin 0.5
486
487 set PTRatioMax 0.12
488}
489
490
491#####################
492# Electron efficiency
493#####################
494
495module Efficiency ElectronEfficiency {
496 set InputArray ElectronFilter/electrons
497 set OutputArray electrons
498
499 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
500
501 # efficiency formula for electrons
502 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
503 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
504 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
505 (abs(eta) > 2.5) * (0.00)}
506}
507
508####################
509# Electron isolation
510####################
511
512module Isolation ElectronIsolation {
513 set CandidateInputArray ElectronEfficiency/electrons
514 set IsolationInputArray EFlowFilter/eflow
515
516 set OutputArray electrons
517
518 set DeltaRMax 0.5
519
520 set PTMin 0.5
521
522 set PTRatioMax 0.12
523}
524
525#################
526# Muon efficiency
527#################
528
529module Efficiency MuonEfficiency {
530 set InputArray MuonMomentumSmearing/muons
531 set OutputArray muons
532
533 # set EfficiencyFormula {efficiency as a function of eta and pt}
534
535 # efficiency formula for muons
536 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
537 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
538 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
539 (abs(eta) > 2.4) * (0.00)}
540}
541
542################
543# Muon isolation
544################
545
546module Isolation MuonIsolation {
547 set CandidateInputArray MuonEfficiency/muons
548 set IsolationInputArray EFlowFilter/eflow
549
550 set OutputArray muons
551
552 set DeltaRMax 0.5
553
554 set PTMin 0.5
555
556 set PTRatioMax 0.25
557}
558
559###################
560# Missing ET merger
561###################
562
563module Merger MissingET {
564# add InputArray InputArray
565 add InputArray EFlowMerger/eflow
566 set MomentumOutputArray momentum
567}
568
569##################
570# Scalar HT merger
571##################
572
573module Merger ScalarHT {
574# add InputArray InputArray
575 add InputArray UniqueObjectFinder/jets
576 add InputArray UniqueObjectFinder/electrons
577 add InputArray UniqueObjectFinder/photons
578 add InputArray UniqueObjectFinder/muons
579 set EnergyOutputArray energy
580}
581
582
583#####################
584# Neutrino Filter
585#####################
586
587module PdgCodeFilter NeutrinoFilter {
588
589 set InputArray Delphes/stableParticles
590 set OutputArray filteredParticles
591
592 set PTMin 0.0
593
594 add PdgCode {12}
595 add PdgCode {14}
596 add PdgCode {16}
597 add PdgCode {-12}
598 add PdgCode {-14}
599 add PdgCode {-16}
600
601}
602
603
604#####################
605# MC truth jet finder
606#####################
607
608module FastJetFinder GenJetFinder {
609 set InputArray NeutrinoFilter/filteredParticles
610
611 set OutputArray jets
612
613 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
614 set JetAlgorithm 6
615 set ParameterR 0.4
616
617 set JetPTMin 20.0
618}
619
620#########################
621# Gen Missing ET merger
622########################
623
624module Merger GenMissingET {
625# add InputArray InputArray
626 add InputArray NeutrinoFilter/filteredParticles
627 set MomentumOutputArray momentum
628}
629
630
631
632############
633# Jet finder
634############
635
636module FastJetFinder FastJetFinder {
637# set InputArray Calorimeter/towers
638 set InputArray EFlowMerger/eflow
639
640 set OutputArray jets
641
642 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
643 set JetAlgorithm 6
644 set ParameterR 0.4
645
646 set JetPTMin 20.0
647}
648
649##################
650# Fat Jet finder
651##################
652
653module FastJetFinder FatJetFinder {
654 set InputArray EFlowMerger/eflow
655
656 set OutputArray jets
657
658 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
659 set JetAlgorithm 6
660 set ParameterR 0.8
661
662 set ComputeNsubjettiness 1
663 set Beta 1.0
664 set AxisMode 4
665
666 set ComputeTrimming 1
667 set RTrim 0.2
668 set PtFracTrim 0.05
669
670 set ComputePruning 1
671 set ZcutPrun 0.1
672 set RcutPrun 0.5
673 set RPrun 0.8
674
675 set ComputeSoftDrop 1
676 set BetaSoftDrop 0.0
677 set SymmetryCutSoftDrop 0.1
678 set R0SoftDrop 0.8
679
680 set JetPTMin 200.0
681}
682
683
684
685
686##################
687# Jet Energy Scale
688##################
689
690module EnergyScale JetEnergyScale {
691 set InputArray FastJetFinder/jets
692 set OutputArray jets
693
694 # scale formula for jets
695 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
696}
697
698########################
699# Jet Flavor Association
700########################
701
702module JetFlavorAssociation JetFlavorAssociation {
703
704 set PartonInputArray Delphes/partons
705 set ParticleInputArray Delphes/allParticles
706 set ParticleLHEFInputArray Delphes/allParticlesLHEF
707 set JetInputArray JetEnergyScale/jets
708
709 set DeltaR 0.5
710 set PartonPTMin 1.0
711 set PartonEtaMax 2.5
712
713}
714
715###########
716# b-tagging
717###########
718
719module BTagging BTagging {
720 set JetInputArray JetEnergyScale/jets
721
722 set BitNumber 0
723
724 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
725 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
726 # gluon's PDG code has the lowest priority
727
728 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) +
729 (pt >= 130.0) * (0.055 + 4.53*10^-4*pt - 1.60*10^-7*pt^2) }
730
731 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) +
732 (pt >= 205.0) * (0.478 + 1.573*10^-4*pt)}
733
734 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) +
735 (pt >= 150.0) * (0.906 - 6.39*10^-5*pt + 4.11*10^-8*pt^2) }
736}
737
738#############
739# tau-tagging
740#############
741
742module TauTagging TauTagging {
743 set ParticleInputArray Delphes/allParticles
744 set PartonInputArray Delphes/partons
745 set JetInputArray JetEnergyScale/jets
746
747 set DeltaR 0.5
748
749 set TauPTMin 1.0
750
751 set TauEtaMax 2.5
752
753 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
754
755 # default efficiency formula (misidentification rate)
756 add EfficiencyFormula {0} {0.01}
757 # efficiency formula for tau-jets
758 add EfficiencyFormula {15} {0.6}
759}
760
761#####################################################
762# Find uniquely identified photons/electrons/tau/jets
763#####################################################
764
765module UniqueObjectFinder UniqueObjectFinder {
766# earlier arrays take precedence over later ones
767# add InputArray InputArray OutputArray
768 add InputArray PhotonIsolation/photons photons
769 add InputArray ElectronIsolation/electrons electrons
770 add InputArray MuonIsolation/muons muons
771 add InputArray JetEnergyScale/jets jets
772}
773
774##################
775# ROOT tree writer
776##################
777
778# tracks, towers and eflow objects are not stored by default in the output.
779# if needed (for jet constituent or other studies), uncomment the relevant
780# "add Branch ..." lines.
781
782module TreeWriter TreeWriter {
783# add Branch InputArray BranchName BranchClass
784 add Branch Delphes/allParticles Particle GenParticle
785
786 add Branch TrackMerger/tracks Track Track
787 add Branch Calorimeter/towers Tower Tower
788
789 add Branch HCal/eflowTracks EFlowTrack Track
790 add Branch ECal/eflowPhotons EFlowPhoton Tower
791 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
792
793 add Branch GenJetFinder/jets GenJet Jet
794 add Branch GenMissingET/momentum GenMissingET MissingET
795
796 add Branch UniqueObjectFinder/jets Jet Jet
797 add Branch UniqueObjectFinder/electrons Electron Electron
798 add Branch UniqueObjectFinder/photons Photon Photon
799 add Branch UniqueObjectFinder/muons Muon Muon
800
801 add Branch FatJetFinder/jets FatJet Jet
802
803 add Branch MissingET/momentum MissingET MissingET
804 add Branch ScalarHT/energy ScalarHT ScalarHT
805}