Fork me on GitHub

Ticket #1349: delphes_card_CMS_smear.tcl

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