Fork me on GitHub

Ticket #1521: delphes_card_CMS_AK4_AK8_isolation.tcl

File delphes_card_CMS_AK4_AK8_isolation.tcl, 28.2 KB (added by bickendorf, 3 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 ElectronIsolation1
32 ElectronIsolation2
33 ElectronIsolation3
34 ElectronIsolation4
35 ElectronIsolation5
36 ElectronIsolation6
37 ElectronIsolation7
38
39 ChargedHadronFilter
40
41 MuonEfficiency
42 MuonIsolation
43 MuonIsolation1
44 MuonIsolation2
45 MuonIsolation3
46 MuonIsolation4
47 MuonIsolation5
48 MuonIsolation6
49 MuonIsolation7
50
51 MissingET
52
53 NeutrinoFilter
54 GenJetFinder
55 GenMissingET
56
57 FastJetFinderAK4
58 FastJetFinderAK8
59 FatJetFinder
60
61 JetEnergyScaleAK4
62 JetEnergyScaleAK8
63
64 JetFlavorAssociation
65
66 BTagging
67 TauTagging
68
69 UniqueObjectFinderLeptonIsolation
70 UniqueObjectFinder
71
72 ElectronMerger
73 MuonMerger
74
75 ScalarHT
76
77 TreeWriter
78}
79
80#################################
81# Propagate particles in cylinder
82#################################
83
84module ParticlePropagator ParticlePropagator {
85 set InputArray Delphes/stableParticles
86
87 set OutputArray stableParticles
88 set ChargedHadronOutputArray chargedHadrons
89 set ElectronOutputArray electrons
90 set MuonOutputArray muons
91
92 # radius of the magnetic field coverage, in m
93 set Radius 1.29
94 # half-length of the magnetic field coverage, in m
95 set HalfLength 3.00
96
97 # magnetic field
98 set Bz 3.8
99}
100
101####################################
102# Charged hadron tracking efficiency
103####################################
104
105module Efficiency ChargedHadronTrackingEfficiency {
106 set InputArray ParticlePropagator/chargedHadrons
107 set OutputArray chargedHadrons
108
109 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
110
111 # tracking efficiency formula for charged hadrons
112 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
113 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
114 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
115 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
117 (abs(eta) > 2.5) * (0.00)}
118}
119
120##############################
121# Electron tracking efficiency
122##############################
123
124module Efficiency ElectronTrackingEfficiency {
125 set InputArray ParticlePropagator/electrons
126 set OutputArray electrons
127
128 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
129
130 # tracking efficiency formula for electrons
131 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
132 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
133 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
134 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
135 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
136 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
137 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
138 (abs(eta) > 2.5) * (0.00)}
139}
140
141##########################
142# Muon tracking efficiency
143##########################
144
145module Efficiency MuonTrackingEfficiency {
146 set InputArray ParticlePropagator/muons
147 set OutputArray muons
148
149 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
150
151 # tracking efficiency formula for muons
152 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
153 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
154 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
155 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
156
157 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
158 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
159 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
160 (abs(eta) > 2.5) * (0.00)}
161}
162
163########################################
164# Momentum resolution for charged tracks
165########################################
166
167module MomentumSmearing ChargedHadronMomentumSmearing {
168 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
169 set OutputArray chargedHadrons
170
171 # set ResolutionFormula {resolution formula as a function of eta and pt}
172
173 # resolution formula for charged hadrons
174 # based on arXiv:1405.6569
175 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
176 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
177 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
178}
179
180###################################
181# Momentum resolution for electrons
182###################################
183
184module MomentumSmearing ElectronMomentumSmearing {
185 set InputArray ElectronTrackingEfficiency/electrons
186 set OutputArray electrons
187
188 # set ResolutionFormula {resolution formula as a function of eta and energy}
189
190 # resolution formula for electrons
191 # based on arXiv:1502.02701
192 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
193 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
194 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
195}
196
197###############################
198# Momentum resolution for muons
199###############################
200
201module MomentumSmearing MuonMomentumSmearing {
202 set InputArray MuonTrackingEfficiency/muons
203 set OutputArray muons
204
205 # set ResolutionFormula {resolution formula as a function of eta and pt}
206
207 # resolution formula for muons
208 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
209 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
210 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
211}
212
213##############
214# Track merger
215##############
216
217module Merger TrackMerger {
218# add InputArray InputArray
219 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
220 add InputArray ElectronMomentumSmearing/electrons
221 add InputArray MuonMomentumSmearing/muons
222 set OutputArray tracks
223}
224
225
226
227#############
228# ECAL
229#############
230
231module SimpleCalorimeter ECal {
232 set ParticleInputArray ParticlePropagator/stableParticles
233 set TrackInputArray TrackMerger/tracks
234
235 set TowerOutputArray ecalTowers
236 set EFlowTrackOutputArray eflowTracks
237 set EFlowTowerOutputArray eflowPhotons
238
239 set IsEcal true
240
241 set EnergyMin 0.5
242 set EnergySignificanceMin 2.0
243
244 set SmearTowerCenter true
245
246 set pi [expr {acos(-1)}]
247
248 # lists of the edges of each tower in eta and phi
249 # each list starts with the lower edge of the first tower
250 # the list ends with the higher edged of the last tower
251
252 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
253
254 set PhiBins {}
255 for {set i -180} {$i <= 180} {incr i} {
256 add PhiBins [expr {$i * $pi/180.0}]
257 }
258
259 # 0.02 unit in eta up to eta = 1.5 (barrel)
260 for {set i -85} {$i <= 86} {incr i} {
261 set eta [expr {$i * 0.0174}]
262 add EtaPhiBins $eta $PhiBins
263 }
264
265 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
266
267 set PhiBins {}
268 for {set i -180} {$i <= 180} {incr i} {
269 add PhiBins [expr {$i * $pi/180.0}]
270 }
271
272 # 0.02 unit in eta up to eta = 3
273 for {set i 1} {$i <= 84} {incr i} {
274 set eta [expr { -2.958 + $i * 0.0174}]
275 add EtaPhiBins $eta $PhiBins
276 }
277
278 for {set i 1} {$i <= 84} {incr i} {
279 set eta [expr { 1.4964 + $i * 0.0174}]
280 add EtaPhiBins $eta $PhiBins
281 }
282
283 # take present CMS granularity for HF
284
285 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
286 set PhiBins {}
287 for {set i -18} {$i <= 18} {incr i} {
288 add PhiBins [expr {$i * $pi/18.0}]
289 }
290
291 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} {
292 add EtaPhiBins $eta $PhiBins
293 }
294
295
296 add EnergyFraction {0} {0.0}
297 # energy fractions for e, gamma and pi0
298 add EnergyFraction {11} {1.0}
299 add EnergyFraction {22} {1.0}
300 add EnergyFraction {111} {1.0}
301 # energy fractions for muon, neutrinos and neutralinos
302 add EnergyFraction {12} {0.0}
303 add EnergyFraction {13} {0.0}
304 add EnergyFraction {14} {0.0}
305 add EnergyFraction {16} {0.0}
306 add EnergyFraction {1000022} {0.0}
307 add EnergyFraction {1000023} {0.0}
308 add EnergyFraction {1000025} {0.0}
309 add EnergyFraction {1000035} {0.0}
310 add EnergyFraction {1000045} {0.0}
311 # energy fractions for K0short and Lambda
312 add EnergyFraction {310} {0.3}
313 add EnergyFraction {3122} {0.3}
314
315 # set ResolutionFormula {resolution formula as a function of eta and energy}
316
317 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
318
319 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
320 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
321 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
322 (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) +
323 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
324
325}
326
327
328#############
329# HCAL
330#############
331
332module SimpleCalorimeter HCal {
333 set ParticleInputArray ParticlePropagator/stableParticles
334 set TrackInputArray ECal/eflowTracks
335
336 set TowerOutputArray hcalTowers
337 set EFlowTrackOutputArray eflowTracks
338 set EFlowTowerOutputArray eflowNeutralHadrons
339
340 set IsEcal false
341
342 set EnergyMin 1.0
343 set EnergySignificanceMin 1.0
344
345 set SmearTowerCenter true
346
347 set pi [expr {acos(-1)}]
348
349 # lists of the edges of each tower in eta and phi
350 # each list starts with the lower edge of the first tower
351 # the list ends with the higher edged of the last tower
352
353 # 5 degrees towers
354 set PhiBins {}
355 for {set i -36} {$i <= 36} {incr i} {
356 add PhiBins [expr {$i * $pi/36.0}]
357 }
358 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} {
359 add EtaPhiBins $eta $PhiBins
360 }
361
362 # 10 degrees towers
363 set PhiBins {}
364 for {set i -18} {$i <= 18} {incr i} {
365 add PhiBins [expr {$i * $pi/18.0}]
366 }
367 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} {
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 # 20 degrees towers
372 set PhiBins {}
373 for {set i -9} {$i <= 9} {incr i} {
374 add PhiBins [expr {$i * $pi/9.0}]
375 }
376 foreach eta {-5 -4.7 -4.525 4.7 5} {
377 add EtaPhiBins $eta $PhiBins
378 }
379
380 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
381 add EnergyFraction {0} {1.0}
382 # energy fractions for e, gamma and pi0
383 add EnergyFraction {11} {0.0}
384 add EnergyFraction {22} {0.0}
385 add EnergyFraction {111} {0.0}
386 # energy fractions for muon, neutrinos and neutralinos
387 add EnergyFraction {12} {0.0}
388 add EnergyFraction {13} {0.0}
389 add EnergyFraction {14} {0.0}
390 add EnergyFraction {16} {0.0}
391 add EnergyFraction {1000022} {0.0}
392 add EnergyFraction {1000023} {0.0}
393 add EnergyFraction {1000025} {0.0}
394 add EnergyFraction {1000035} {0.0}
395 add EnergyFraction {1000045} {0.0}
396 # energy fractions for K0short and Lambda
397 add EnergyFraction {310} {0.7}
398 add EnergyFraction {3122} {0.7}
399
400 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
401 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
402 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
403
404}
405
406
407#################
408# Electron filter
409#################
410
411module PdgCodeFilter ElectronFilter {
412 set InputArray HCal/eflowTracks
413 set OutputArray electrons
414 set Invert true
415 add PdgCode {11}
416 add PdgCode {-11}
417}
418
419######################
420# ChargedHadronFilter
421######################
422
423module PdgCodeFilter ChargedHadronFilter {
424 set InputArray HCal/eflowTracks
425 set OutputArray chargedHadrons
426
427 add PdgCode {11}
428 add PdgCode {-11}
429 add PdgCode {13}
430 add PdgCode {-13}
431}
432
433
434###################################################
435# Tower Merger (in case not using e-flow algorithm)
436###################################################
437
438module Merger Calorimeter {
439# add InputArray InputArray
440 add InputArray ECal/ecalTowers
441 add InputArray HCal/hcalTowers
442 set OutputArray towers
443}
444
445
446
447####################
448# Energy flow merger
449####################
450
451module Merger EFlowMerger {
452# add InputArray InputArray
453 add InputArray HCal/eflowTracks
454 add InputArray ECal/eflowPhotons
455 add InputArray HCal/eflowNeutralHadrons
456 set OutputArray eflow
457}
458
459######################
460# EFlowFilter
461######################
462
463module PdgCodeFilter EFlowFilter {
464 set InputArray EFlowMerger/eflow
465 set OutputArray eflow
466
467 add PdgCode {11}
468 add PdgCode {-11}
469 add PdgCode {13}
470 add PdgCode {-13}
471}
472
473
474###################
475# Photon efficiency
476###################
477
478module Efficiency PhotonEfficiency {
479 set InputArray ECal/eflowPhotons
480 set OutputArray photons
481
482 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
483
484 # efficiency formula for photons
485 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
486 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
487 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
488 (abs(eta) > 2.5) * (0.00)}
489}
490
491##################
492# Photon isolation
493##################
494
495module Isolation PhotonIsolation {
496 set CandidateInputArray PhotonEfficiency/photons
497 set IsolationInputArray EFlowFilter/eflow
498
499 set OutputArray photons
500
501 set DeltaRMax 0.5
502
503 set PTMin 0.5
504
505 set PTRatioMax 0.12
506}
507
508
509#####################
510# Electron efficiency
511#####################
512
513module Efficiency ElectronEfficiency {
514 set InputArray ElectronFilter/electrons
515 set OutputArray electrons
516
517 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
518
519 # efficiency formula for electrons
520 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
521 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
522 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
523 (abs(eta) > 2.5) * (0.00)}
524}
525
526####################
527# Electron isolation
528####################
529
530module Isolation ElectronIsolation {
531 set CandidateInputArray ElectronEfficiency/electrons
532 set IsolationInputArray EFlowFilter/eflow
533
534 set OutputArray electrons
535
536 set DeltaRMax 0.2
537
538 set PTMin 0.5
539
540 set PTRatioMax 0.12
541}
542
543module Isolation ElectronIsolation1 {
544 set CandidateInputArray ElectronEfficiency/electrons
545 set IsolationInputArray EFlowFilter/eflow
546 set OutputArray electrons
547 set DeltaRMax 0.05
548 set PTMin 200
549 set PTRatioMax 0.1
550}
551
552module Isolation ElectronIsolation2 {
553 set CandidateInputArray ElectronEfficiency/electrons
554 set IsolationInputArray EFlowFilter/eflow
555 set OutputArray electrons
556 set DeltaRMax 0.075
557 set PTMin 133
558 set PTRatioMax 0.1
559}
560
561module Isolation ElectronIsolation3 {
562 set CandidateInputArray ElectronEfficiency/electrons
563 set IsolationInputArray EFlowFilter/eflow
564 set OutputArray electrons
565 set DeltaRMax 0.1
566 set PTMin 100
567 set PTRatioMax 0.1
568}
569
570module Isolation ElectronIsolation4 {
571 set CandidateInputArray ElectronEfficiency/electrons
572 set IsolationInputArray EFlowFilter/eflow
573 set OutputArray electrons
574 set DeltaRMax 0.125
575 set PTMin 80
576 set PTRatioMax 0.1
577}
578
579module Isolation ElectronIsolation5 {
580 set CandidateInputArray ElectronEfficiency/electrons
581 set IsolationInputArray EFlowFilter/eflow
582 set OutputArray electrons
583 set DeltaRMax 0.15
584 set PTMin 66
585 set PTRatioMax 0.1
586}
587
588module Isolation ElectronIsolation6 {
589 set CandidateInputArray ElectronEfficiency/electrons
590 set IsolationInputArray EFlowFilter/eflow
591 set OutputArray electrons
592 set DeltaRMax 0.175
593 set PTMin 57
594 set PTRatioMax 0.1
595}
596
597module Isolation ElectronIsolation7 {
598 set CandidateInputArray ElectronEfficiency/electrons
599 set IsolationInputArray EFlowFilter/eflow
600 set OutputArray electrons
601 set DeltaRMax 0.2
602 set PTMin 10
603 set PTRatioMax 0.1
604}
605
606
607#################
608# Muon efficiency
609#################
610
611module Efficiency MuonEfficiency {
612 set InputArray MuonMomentumSmearing/muons
613 set OutputArray muons
614
615 # set EfficiencyFormula {efficiency as a function of eta and pt}
616
617 # efficiency formula for muons
618 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
619 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
620 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
621 (abs(eta) > 2.4) * (0.00)}
622}
623
624################
625# Muon isolation
626################
627
628module Isolation MuonIsolation {
629 set CandidateInputArray MuonEfficiency/muons
630 set IsolationInputArray EFlowFilter/eflow
631 set OutputArray muons
632
633 set DeltaRMax 0.2
634
635 set PTMin 0.5
636
637 set PTRatioMax 0.25
638}
639
640module Isolation MuonIsolation1 {
641 set CandidateInputArray MuonEfficiency/muons
642 set IsolationInputArray EFlowFilter/eflow
643 set OutputArray muons
644 set DeltaRMax 0.05
645 set PTMin 200
646 set PTRatioMax 0.2
647}
648
649module Isolation MuonIsolation2 {
650 set CandidateInputArray MuonEfficiency/muons
651 set IsolationInputArray EFlowFilter/eflow
652 set OutputArray muons
653 set DeltaRMax 0.075
654 set PTMin 133
655 set PTRatioMax 0.2
656}
657
658module Isolation MuonIsolation3 {
659 set CandidateInputArray MuonEfficiency/muons
660 set IsolationInputArray EFlowFilter/eflow
661 set OutputArray muons
662 set DeltaRMax 0.1
663 set PTMin 100
664 set PTRatioMax 0.2
665}
666
667module Isolation MuonIsolation4 {
668 set CandidateInputArray MuonEfficiency/muons
669 set IsolationInputArray EFlowFilter/eflow
670 set OutputArray muons
671 set DeltaRMax 0.125
672 set PTMin 80
673 set PTRatioMax 0.2
674}
675
676module Isolation MuonIsolation5 {
677 set CandidateInputArray MuonEfficiency/muons
678 set IsolationInputArray EFlowFilter/eflow
679 set OutputArray muons
680 set DeltaRMax 0.15
681 set PTMin 66
682 set PTRatioMax 0.2
683}
684
685module Isolation MuonIsolation6 {
686 set CandidateInputArray MuonEfficiency/muons
687 set IsolationInputArray EFlowFilter/eflow
688 set OutputArray muons
689 set DeltaRMax 0.175
690 set PTMin 57
691 set PTRatioMax 0.2
692}
693
694module Isolation MuonIsolation7 {
695 set CandidateInputArray MuonEfficiency/muons
696 set IsolationInputArray EFlowFilter/eflow
697 set OutputArray muons
698 set DeltaRMax 0.2
699 set PTMin 10
700 set PTRatioMax 0.2
701}
702
703###################
704# Missing ET merger
705###################
706
707module Merger MissingET {
708# add InputArray InputArray
709 add InputArray EFlowMerger/eflow
710 set MomentumOutputArray momentum
711}
712
713##################
714# Scalar HT merger
715##################
716
717module Merger ScalarHT {
718# add InputArray InputArray
719 add InputArray UniqueObjectFinder/jets
720 add InputArray UniqueObjectFinder/electrons
721 add InputArray UniqueObjectFinder/photons
722 add InputArray UniqueObjectFinder/muons
723 set EnergyOutputArray energy
724}
725
726
727#####################
728# Neutrino Filter
729#####################
730
731module PdgCodeFilter NeutrinoFilter {
732
733 set InputArray Delphes/stableParticles
734 set OutputArray filteredParticles
735
736 set PTMin 0.0
737
738 add PdgCode {12}
739 add PdgCode {14}
740 add PdgCode {16}
741 add PdgCode {-12}
742 add PdgCode {-14}
743 add PdgCode {-16}
744
745}
746
747
748#####################
749# MC truth jet finder
750#####################
751
752module FastJetFinder GenJetFinder {
753 set InputArray NeutrinoFilter/filteredParticles
754
755 set OutputArray jets
756
757 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
758 set JetAlgorithm 6
759 set ParameterR 0.5
760
761 set JetPTMin 20.0
762}
763
764#########################
765# Gen Missing ET merger
766########################
767
768module Merger GenMissingET {
769# add InputArray InputArray
770 add InputArray NeutrinoFilter/filteredParticles
771 set MomentumOutputArray momentum
772}
773
774
775
776############
777# Jet finder
778############
779
780module FastJetFinder FastJetFinderAK4 {
781# set InputArray Calorimeter/towers
782 set InputArray EFlowMerger/eflow
783
784 set OutputArray jets
785
786 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
787 set JetAlgorithm 6
788 set ParameterR 0.4
789
790 set JetPTMin 30.0
791}
792
793module FastJetFinder FastJetFinderAK8 {
794# set InputArray Calorimeter/towers
795 set InputArray EFlowMerger/eflow
796
797 set OutputArray jets
798
799 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
800 set JetAlgorithm 6
801 set ParameterR 0.8
802
803 set ComputeNsubjettiness 1
804 set Beta 1.0
805 set AxisMode 4
806
807 set ComputeTrimming 1
808 set RTrim 0.2
809 set PtFracTrim 0.05
810
811 set ComputePruning 1
812 set ZcutPrun 0.1
813 set RcutPrun 0.5
814 set RPrun 0.8
815
816 set ComputeSoftDrop 1
817 set BetaSoftDrop 0.0
818 set SymmetryCutSoftDrop 0.1
819 set R0SoftDrop 0.8
820
821 set JetPTMin 30.0
822}
823
824##################
825# Fat Jet finder
826##################
827
828module FastJetFinder FatJetFinder {
829 set InputArray EFlowMerger/eflow
830
831 set OutputArray jets
832
833 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
834 set JetAlgorithm 6
835 set ParameterR 0.8
836
837 set ComputeNsubjettiness 1
838 set Beta 1.0
839 set AxisMode 4
840
841 set ComputeTrimming 1
842 set RTrim 0.2
843 set PtFracTrim 0.05
844
845 set ComputePruning 1
846 set ZcutPrun 0.1
847 set RcutPrun 0.5
848 set RPrun 0.8
849
850 set ComputeSoftDrop 1
851 set BetaSoftDrop 0.0
852 set SymmetryCutSoftDrop 0.1
853 set R0SoftDrop 0.8
854
855 set JetPTMin 200.0
856}
857
858
859
860
861##################
862# Jet Energy Scale
863##################
864
865module EnergyScale JetEnergyScaleAK4 {
866 set InputArray FastJetFinderAK4/jets
867 set OutputArray jets
868
869 # scale formula for jets
870 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
871}
872
873module EnergyScale JetEnergyScaleAK8 {
874 set InputArray FastJetFinderAK8/jets
875 set OutputArray jets
876
877 # scale formula for jets
878 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
879}
880
881########################
882# Jet Flavor Association
883########################
884
885module JetFlavorAssociation JetFlavorAssociation {
886
887 set PartonInputArray Delphes/partons
888 set ParticleInputArray Delphes/allParticles
889 set ParticleLHEFInputArray Delphes/allParticlesLHEF
890 set JetInputArray JetEnergyScaleAK4/jets
891
892 set DeltaR 0.5
893 set PartonPTMin 1.0
894 set PartonEtaMax 2.5
895
896}
897
898###########
899# b-tagging
900###########
901
902module BTagging BTagging {
903 set JetInputArray JetEnergyScaleAK4/jets
904
905 set BitNumber 0
906
907 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
908 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
909 # gluon's PDG code has the lowest priority
910
911 # based on arXiv:1211.4462
912
913 # default efficiency formula (misidentification rate)
914 add EfficiencyFormula {0} {0.01+0.000038*pt}
915
916 # efficiency formula for c-jets (misidentification rate)
917 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
918
919 # efficiency formula for b-jets
920 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
921}
922
923#############
924# tau-tagging
925#############
926
927module TauTagging TauTagging {
928 set ParticleInputArray Delphes/allParticles
929 set PartonInputArray Delphes/partons
930 set JetInputArray JetEnergyScaleAK4/jets
931
932 set DeltaR 0.5
933
934 set TauPTMin 1.0
935
936 set TauEtaMax 2.5
937
938 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
939
940 # default efficiency formula (misidentification rate)
941 add EfficiencyFormula {0} {0.01}
942 # efficiency formula for tau-jets
943 add EfficiencyFormula {15} {0.6}
944}
945
946#####################################################
947# Find uniquely identified photons/electrons/tau/jets
948#####################################################
949
950module UniqueObjectFinder UniqueObjectFinderLeptonIsolation {
951# earlier arrays take precedence over later ones
952# add InputArray InputArray OutputArray
953 add InputArray ElectronIsolation1/electrons electrons1
954 add InputArray ElectronIsolation2/electrons electrons2
955 add InputArray ElectronIsolation3/electrons electrons3
956 add InputArray ElectronIsolation4/electrons electrons4
957 add InputArray ElectronIsolation5/electrons electrons5
958 add InputArray ElectronIsolation6/electrons electrons6
959 add InputArray ElectronIsolation7/electrons electrons7
960
961 add InputArray MuonIsolation1/muons muons1
962 add InputArray MuonIsolation2/muons muons2
963 add InputArray MuonIsolation3/muons muons3
964 add InputArray MuonIsolation4/muons muons4
965 add InputArray MuonIsolation5/muons muons5
966 add InputArray MuonIsolation6/muons muons6
967 add InputArray MuonIsolation7/muons muons7
968}
969
970module UniqueObjectFinder UniqueObjectFinder {
971# earlier arrays take precedence over later ones
972# add InputArray InputArray OutputArray
973 add InputArray PhotonIsolation/photons photons
974 add InputArray ElectronIsolation/electrons electrons
975 add InputArray MuonIsolation/muons muons
976 add InputArray JetEnergyScaleAK4/jets jets
977}
978
979module Merger ElectronMerger {
980# add InputArray InputArray
981 add InputArray UniqueObjectFinderLeptonIsolation/electrons1
982 add InputArray UniqueObjectFinderLeptonIsolation/electrons2
983 add InputArray UniqueObjectFinderLeptonIsolation/electrons3
984 add InputArray UniqueObjectFinderLeptonIsolation/electrons4
985 add InputArray UniqueObjectFinderLeptonIsolation/electrons5
986 add InputArray UniqueObjectFinderLeptonIsolation/electrons6
987 add InputArray UniqueObjectFinderLeptonIsolation/electrons7
988 set OutputArray Electron
989}
990
991module Merger MuonMerger {
992# add InputArray InputArray
993 add InputArray UniqueObjectFinderLeptonIsolation/muons1
994 add InputArray UniqueObjectFinderLeptonIsolation/muons2
995 add InputArray UniqueObjectFinderLeptonIsolation/muons3
996 add InputArray UniqueObjectFinderLeptonIsolation/muons4
997 add InputArray UniqueObjectFinderLeptonIsolation/muons5
998 add InputArray UniqueObjectFinderLeptonIsolation/muons6
999 add InputArray UniqueObjectFinderLeptonIsolation/muons7
1000 set OutputArray Muon
1001}
1002
1003##################
1004# ROOT tree writer
1005##################
1006
1007# tracks, towers and eflow objects are not stored by default in the output.
1008# if needed (for jet constituent or other studies), uncomment the relevant
1009# "add Branch ..." lines.
1010
1011module TreeWriter TreeWriter {
1012# add Branch InputArray BranchName BranchClass
1013 #add Branch Delphes/allParticles Particle GenParticle
1014
1015 add Branch TrackMerger/tracks Track Track
1016 #add Branch Calorimeter/towers Tower Tower
1017
1018 #add Branch HCal/eflowTracks EFlowTrack Track
1019 #add Branch ECal/eflowPhotons EFlowPhoton Tower
1020 #add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1021
1022 #add Branch GenJetFinder/jets GenJet Jet
1023 #add Branch GenMissingET/momentum GenMissingET MissingET
1024 add Branch ElectronMerger/Electron IsoElectrons Electron
1025 add Branch MuonMerger/Muon IsoMuon Muon
1026
1027 add Branch UniqueObjectFinder/jets Jet Jet
1028 add Branch UniqueObjectFinder/electrons Electron Electron
1029 add Branch UniqueObjectFinder/photons Photon Photon
1030 add Branch UniqueObjectFinder/muons Muon Muon
1031
1032 #add Branch FatJetFinder/jets FatJet Jet
1033
1034 add Branch JetEnergyScaleAK4/jets AK4Jets Jet
1035 add Branch JetEnergyScaleAK8/jets AK8Jets Jet
1036
1037 add Branch MissingET/momentum MissingET MissingET
1038 #add Branch ScalarHT/energy ScalarHT ScalarHT
1039}