MA5SandBox: delphes_card_cms_top_18_003.tcl

File delphes_card_cms_top_18_003.tcl, 29.2 KB (added by Benjamin Fuks, 5 years ago)
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6
7 PileUpMerger
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19
20 ECal
21 HCal
22
23 ElectronFilter
24 TrackPileUpSubtractor
25 NeutralTowerMerger
26 EFlowMergerAllTracks
27 EFlowMerger
28 EFlowFilter
29
30 NeutrinoFilter
31 GenJetFinder
32 GenMissingET
33
34 Rho
35 FastJetFinder
36 PileUpJetID
37 JetPileUpSubtractor
38
39 JetEnergyScale
40
41 PhotonEfficiency
42 PhotonIsolation
43
44 ElectronEfficiency
45 ElectronIsolation
46
47 MuonEfficiency
48 MuonIsolation
49
50 MissingET
51
52 JetFlavorAssociation
53
54 BTagging
55 TauTagging
56
57 UniqueObjectFinder
58
59 ScalarHT
60
61 TreeWriter
62}
63
64###############
65# PileUp Merger
66###############
67
68module PileUpMerger PileUpMerger {
69 set InputArray Delphes/stableParticles
70
71 set ParticleOutputArray stableParticles
72 set VertexOutputArray vertices
73
74 # pre-generated minbias input file
75 set PileUpFile MinBias.pileup
76
77 # average expected pile up 12 for 2015 run, 32 afterwards
78 set MeanPileUp 32
79
80 # maximum spread in the beam direction in m
81 set ZVertexSpread 0.25
82
83 # maximum spread in time in s
84 set TVertexSpread 800E-12
85
86 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
87 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
88
89
90}
91
92#################################
93# Propagate particles in cylinder
94#################################
95
96module ParticlePropagator ParticlePropagator {
97 set InputArray PileUpMerger/stableParticles
98
99 set OutputArray stableParticles
100 set ChargedHadronOutputArray chargedHadrons
101 set ElectronOutputArray electrons
102 set MuonOutputArray muons
103
104 # radius of the magnetic field coverage, in m
105 set Radius 1.29
106 # half-length of the magnetic field coverage, in m
107 set HalfLength 3.00
108
109 # magnetic field
110 set Bz 3.8
111}
112
113####################################
114# Charged hadron tracking efficiency
115####################################
116
117module Efficiency ChargedHadronTrackingEfficiency {
118 set InputArray ParticlePropagator/chargedHadrons
119 set OutputArray chargedHadrons
120
121 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
122
123 # tracking efficiency formula for charged hadrons
124 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
125 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
126 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
128 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
129 (abs(eta) > 2.5) * (0.00)}
130}
131
132##############################
133# Electron tracking efficiency
134##############################
135
136module Efficiency ElectronTrackingEfficiency {
137 set InputArray ParticlePropagator/electrons
138 set OutputArray electrons
139
140 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
141
142 # tracking efficiency formula for electrons
143 set EfficiencyFormula {
144 (pt <= 0.1) * (0.00) +
145 (abs(eta) <= 1.) * (pt > 0.1 && pt <= 20.0) * (0.73) +
146 (abs(eta) <= 1.) * (pt > 20.0 && pt <= 75) * (0.96) +
147 (abs(eta) <= 1.) * (pt > 75) * (0.98) +
148 (abs(eta) > 1. && abs(eta) <= 2.) * (pt > 0.1 && pt <=20.0) * (0.50) +
149 (abs(eta) > 1. && abs(eta) <= 2.) * (pt > 20.0 && pt <= 75) * (0.95) +
150 (abs(eta) > 1. && abs(eta) <= 2.) * (pt > 75) * (0.97) +
151 (abs(eta) > 2. && abs(eta) <= 2.5) * (pt > 0.1 && pt <=20.0) * (0.50) +
152 (abs(eta) > 2. && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 75) * (0.96) +
153 (abs(eta) > 2. && abs(eta) <= 2.5) * (pt > 75) * (0.98) +
154 (abs(eta) > 2.5) * (0.00)}
155}
156
157##########################
158# Muon tracking efficiency
159##########################
160
161module Efficiency MuonTrackingEfficiency {
162 set InputArray ParticlePropagator/muons
163 set OutputArray muons
164
165 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
166
167 # tracking efficiency formula for muons
168 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
169 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
170 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
171 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
172
173 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
174 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
175 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
176 (abs(eta) > 2.5) * (0.00)}
177}
178
179########################################
180# Momentum resolution for charged tracks
181########################################
182
183module MomentumSmearing ChargedHadronMomentumSmearing {
184 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
185 set OutputArray chargedHadrons
186
187 # set ResolutionFormula {resolution formula as a function of eta and pt}
188
189 # resolution formula for charged hadrons
190 # based on arXiv:1405.6569
191 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
192 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
193 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
194}
195
196###################################
197# Momentum resolution for electrons
198###################################
199
200module MomentumSmearing ElectronMomentumSmearing {
201 set InputArray ElectronTrackingEfficiency/electrons
202 set OutputArray electrons
203
204 # set ResolutionFormula {resolution formula as a function of eta and energy}
205
206 # resolution formula for electrons
207 # based on arXiv:1405.6569
208 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
209 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
210 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
211}
212
213###############################
214# Momentum resolution for muons
215###############################
216
217module MomentumSmearing MuonMomentumSmearing {
218 set InputArray MuonTrackingEfficiency/muons
219 set OutputArray muons
220
221 # set ResolutionFormula {resolution formula as a function of eta and pt}
222
223 # resolution formula for muons
224 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
225 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
226 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
227}
228
229##############
230# Track merger
231##############
232
233module Merger TrackMerger {
234# add InputArray InputArray
235 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
236 add InputArray ElectronMomentumSmearing/electrons
237 add InputArray MuonMomentumSmearing/muons
238 set OutputArray tracks
239}
240
241#############
242# ECAL
243#############
244
245module SimpleCalorimeter ECal {
246 set ParticleInputArray ParticlePropagator/stableParticles
247 set TrackInputArray TrackMerger/tracks
248
249 set TowerOutputArray ecalTowers
250 set EFlowTrackOutputArray eflowTracks
251 set EFlowTowerOutputArray eflowPhotons
252
253 set IsEcal true
254
255 set EnergyMin 0.5
256 set EnergySignificanceMin 2.0
257
258 set SmearTowerCenter true
259
260 set pi [expr {acos(-1)}]
261
262 # lists of the edges of each tower in eta and phi
263 # each list starts with the lower edge of the first tower
264 # the list ends with the higher edged of the last tower
265
266 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
267
268 set PhiBins {}
269 for {set i -180} {$i <= 180} {incr i} {
270 add PhiBins [expr {$i * $pi/180.0}]
271 }
272
273 # 0.02 unit in eta up to eta = 1.5 (barrel)
274 for {set i -85} {$i <= 86} {incr i} {
275 set eta [expr {$i * 0.0174}]
276 add EtaPhiBins $eta $PhiBins
277 }
278
279 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
280
281 set PhiBins {}
282 for {set i -180} {$i <= 180} {incr i} {
283 add PhiBins [expr {$i * $pi/180.0}]
284 }
285
286 # 0.02 unit in eta up to eta = 3
287 for {set i 1} {$i <= 84} {incr i} {
288 set eta [expr { -2.958 + $i * 0.0174}]
289 add EtaPhiBins $eta $PhiBins
290 }
291
292 for {set i 1} {$i <= 84} {incr i} {
293 set eta [expr { 1.4964 + $i * 0.0174}]
294 add EtaPhiBins $eta $PhiBins
295 }
296
297 # take present CMS granularity for HF
298
299 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
300 set PhiBins {}
301 for {set i -18} {$i <= 18} {incr i} {
302 add PhiBins [expr {$i * $pi/18.0}]
303 }
304
305 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} {
306 add EtaPhiBins $eta $PhiBins
307 }
308
309
310 add EnergyFraction {0} {0.0}
311 # energy fractions for e, gamma and pi0
312 add EnergyFraction {11} {1.0}
313 add EnergyFraction {22} {1.0}
314 add EnergyFraction {111} {1.0}
315 # energy fractions for muon, neutrinos and neutralinos
316 add EnergyFraction {12} {0.0}
317 add EnergyFraction {13} {0.0}
318 add EnergyFraction {14} {0.0}
319 add EnergyFraction {16} {0.0}
320 add EnergyFraction {1000022} {0.0}
321 add EnergyFraction {1000023} {0.0}
322 add EnergyFraction {1000025} {0.0}
323 add EnergyFraction {1000035} {0.0}
324 add EnergyFraction {1000045} {0.0}
325 # energy fractions for K0short and Lambda
326 add EnergyFraction {310} {0.3}
327 add EnergyFraction {3122} {0.3}
328
329 # set ResolutionFormula {resolution formula as a function of eta and energy}
330
331 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
332
333 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
334 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
335 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
336 (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) +
337 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
338
339}
340
341
342#############
343# HCAL
344#############
345
346module SimpleCalorimeter HCal {
347 set ParticleInputArray ParticlePropagator/stableParticles
348 set TrackInputArray ECal/eflowTracks
349
350 set TowerOutputArray hcalTowers
351 set EFlowTrackOutputArray eflowTracks
352 set EFlowTowerOutputArray eflowNeutralHadrons
353
354 set IsEcal false
355
356 set EnergyMin 1.0
357 set EnergySignificanceMin 1.0
358
359 set SmearTowerCenter true
360
361 set pi [expr {acos(-1)}]
362
363 # lists of the edges of each tower in eta and phi
364 # each list starts with the lower edge of the first tower
365 # the list ends with the higher edged of the last tower
366
367 # 5 degrees towers
368 set PhiBins {}
369 for {set i -36} {$i <= 36} {incr i} {
370 add PhiBins [expr {$i * $pi/36.0}]
371 }
372 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} {
373 add EtaPhiBins $eta $PhiBins
374 }
375
376 # 10 degrees towers
377 set PhiBins {}
378 for {set i -18} {$i <= 18} {incr i} {
379 add PhiBins [expr {$i * $pi/18.0}]
380 }
381 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} {
382 add EtaPhiBins $eta $PhiBins
383 }
384
385 # 20 degrees towers
386 set PhiBins {}
387 for {set i -9} {$i <= 9} {incr i} {
388 add PhiBins [expr {$i * $pi/9.0}]
389 }
390 foreach eta {-5 -4.7 -4.525 4.7 5} {
391 add EtaPhiBins $eta $PhiBins
392 }
393
394 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
395 add EnergyFraction {0} {1.0}
396 # energy fractions for e, gamma and pi0
397 add EnergyFraction {11} {0.0}
398 add EnergyFraction {22} {0.0}
399 add EnergyFraction {111} {0.0}
400 # energy fractions for muon, neutrinos and neutralinos
401 add EnergyFraction {12} {0.0}
402 add EnergyFraction {13} {0.0}
403 add EnergyFraction {14} {0.0}
404 add EnergyFraction {16} {0.0}
405 add EnergyFraction {1000022} {0.0}
406 add EnergyFraction {1000023} {0.0}
407 add EnergyFraction {1000025} {0.0}
408 add EnergyFraction {1000035} {0.0}
409 add EnergyFraction {1000045} {0.0}
410 # energy fractions for K0short and Lambda
411 add EnergyFraction {310} {0.7}
412 add EnergyFraction {3122} {0.7}
413
414 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
415 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
416 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
417
418}
419
420#################
421# Electron filter
422#################
423
424module PdgCodeFilter ElectronFilter {
425 set InputArray HCal/eflowTracks
426 set OutputArray electrons
427 set Invert true
428 add PdgCode {11}
429 add PdgCode {-11}
430}
431
432###################################################
433# Tower Merger (in case not using e-flow algorithm)
434###################################################
435
436module Merger Calorimeter {
437# add InputArray InputArray
438 add InputArray ECal/ecalTowers
439 add InputArray HCal/hcalTowers
440 set OutputArray towers
441}
442
443######################
444# EFlowFilter
445######################
446
447module PdgCodeFilter EFlowFilter {
448 set InputArray EFlowMergerAllTracks/eflow
449 set OutputArray eflow
450
451 add PdgCode {11}
452 add PdgCode {-11}
453 add PdgCode {13}
454 add PdgCode {-13}
455}
456
457##########################
458# Track pile-up subtractor
459##########################
460
461module TrackPileUpSubtractor TrackPileUpSubtractor {
462# add InputArray InputArray OutputArray
463 add InputArray HCal/eflowTracks eflowTracks
464 add InputArray ElectronFilter/electrons electrons
465 add InputArray MuonMomentumSmearing/muons muons
466
467 set VertexInputArray PileUpMerger/vertices
468 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
469 # Z vertex resolution in m
470 set ZVertexResolution {0.0001}
471}
472
473
474####################
475# Neutral Tower merger
476####################
477
478module Merger NeutralTowerMerger {
479# add InputArray InputArray
480 add InputArray ECal/eflowPhotons
481 add InputArray HCal/eflowNeutralHadrons
482 set OutputArray towers
483}
484
485
486####################
487# Energy flow merger
488####################
489
490module Merger EFlowMergerAllTracks {
491# add InputArray InputArray
492 add InputArray HCal/eflowTracks
493 add InputArray ECal/eflowPhotons
494 add InputArray HCal/eflowNeutralHadrons
495 set OutputArray eflow
496}
497
498
499
500
501####################
502# Energy flow merger
503####################
504
505module Merger EFlowMerger {
506# add InputArray InputArray
507 add InputArray TrackPileUpSubtractor/eflowTracks
508 add InputArray ECal/eflowPhotons
509 add InputArray HCal/eflowNeutralHadrons
510 set OutputArray eflow
511}
512
513#############
514# Rho pile-up
515#############
516
517module FastJetGridMedianEstimator Rho {
518
519 set InputArray EFlowMerger/eflow
520 set RhoOutputArray rho
521
522 # add GridRange rapmin rapmax drap dphi
523 # rapmin - the minimum rapidity extent of the grid
524 # rapmax - the maximum rapidity extent of the grid
525 # drap - the grid spacing in rapidity
526 # dphi - the grid spacing in azimuth
527
528 add GridRange -5.0 -2.5 1.0 1.0
529 add GridRange -2.5 2.5 1.0 1.0
530 add GridRange 2.5 5.0 1.0 1.0
531
532}
533
534#####################
535# Neutrino Filter
536#####################
537
538module PdgCodeFilter NeutrinoFilter {
539
540 set InputArray Delphes/stableParticles
541 set OutputArray filteredParticles
542
543 set PTMin 0.0
544
545 add PdgCode {12}
546 add PdgCode {14}
547 add PdgCode {16}
548 add PdgCode {-12}
549 add PdgCode {-14}
550 add PdgCode {-16}
551
552}
553
554
555
556#####################
557# MC truth jet finder
558#####################
559
560module FastJetFinder GenJetFinder {
561 set InputArray NeutrinoFilter/filteredParticles
562
563 set OutputArray jets
564
565 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
566 set JetAlgorithm 6
567 set ParameterR 0.4
568
569 set JetPTMin 20.0
570}
571
572#########################
573# Gen Missing ET merger
574########################
575
576module Merger GenMissingET {
577# add InputArray InputArray
578 add InputArray NeutrinoFilter/filteredParticles
579 set MomentumOutputArray momentum
580}
581
582############
583# Jet finder
584############
585
586module FastJetFinder FastJetFinder {
587# set InputArray Calorimeter/towers
588 set InputArray EFlowMerger/eflow
589
590 set OutputArray jets
591
592 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
593 set AreaAlgorithm 5
594
595 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
596 set JetAlgorithm 6
597 set ParameterR 0.4
598
599 set JetPTMin 20.0
600}
601
602###########################
603# Jet Pile-Up ID
604###########################
605
606module PileUpJetID PileUpJetID {
607 set JetInputArray FastJetFinder/jets
608 set TrackInputArray HCal/eflowTracks
609 set NeutralInputArray NeutralTowerMerger/towers
610
611 set VertexInputArray PileUpMerger/vertices
612 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
613 # Z vertex resolution in m
614 set ZVertexResolution 0.0001
615
616 set OutputArray jets
617
618 set UseConstituents 0
619 set ParameterR 0.4
620
621 set JetPTMin 20.0
622}
623
624###########################
625# Jet Pile-Up Subtraction
626###########################
627
628module JetPileUpSubtractor JetPileUpSubtractor {
629 set JetInputArray PileUpJetID/jets
630 set RhoInputArray Rho/rho
631
632 set OutputArray jets
633
634 set JetPTMin 20.0
635}
636
637##################
638# Jet Energy Scale
639##################
640
641module EnergyScale JetEnergyScale {
642 set InputArray JetPileUpSubtractor/jets
643 set OutputArray jets
644
645 # scale formula for jets
646 set ScaleFormula {1.0}
647}
648
649###################
650# Photon efficiency
651###################
652
653module Efficiency PhotonEfficiency {
654 set InputArray ECal/eflowPhotons
655 set OutputArray photons
656
657 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
658
659 # efficiency formula for photons
660 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
661 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
662 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
663 (abs(eta) > 2.5) * (0.00)}
664}
665
666
667##################
668# Photon isolation
669##################
670
671module Isolation PhotonIsolation {
672 set CandidateInputArray PhotonEfficiency/photons
673 set IsolationInputArray EFlowFilter/eflow
674 set RhoInputArray Rho/rho
675
676 set OutputArray photons
677
678 set DeltaRMax 0.4
679
680 set PTMin 0.5
681
682 set PTRatioMax 0.12
683}
684
685#####################
686# Electron efficiency
687#####################
688
689module Efficiency ElectronEfficiency {
690 set InputArray TrackPileUpSubtractor/electrons
691 set OutputArray electrons
692
693 # efficiency formula for electrons, previously 0.85/95 for eta>1.5/<1.5
694 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
695 (abs(eta) <= 1.5) * (pt > 10.0 ) * (pt < 20.0 ) * (0.73) +
696 (abs(eta) <= 1.5) * (pt > 20.0 ) * (pt < 35.0 ) * (0.80) +
697 (abs(eta) <= 1.5) * (pt > 35.0 ) * (pt < 50.0 ) * (0.85) +
698 (abs(eta) <= 1.5) * (pt > 50.0 ) * (0.87) +
699 (abs(eta) > 1.5 && abs(eta) <= 2.) * (pt > 10.0 )* (pt < 35.0 ) * (0.80) +
700 (abs(eta) > 1.5 && abs(eta) <= 2.) * (pt > 35.0 )* (pt < 100.0 ) * (0.88) +
701 (abs(eta) > 1.5 && abs(eta) <= 2.) * (pt > 100.0 ) * (0.93) +
702 (abs(eta) > 2 && abs(eta) <= 2.5) * (pt > 10.0 )* (pt < 20.0 ) * (0.70) +
703 (abs(eta) > 2 && abs(eta) <= 2.5) * (pt > 20.0 )* (pt < 35.0 ) * (0.75) +
704 (abs(eta) > 2 && abs(eta) <= 2.5) * (pt > 35.0 )* (pt < 100.0 ) * (0.83) +
705 (abs(eta) > 2 && abs(eta) <= 2.5) * (pt > 100.0 ) * (0.87) +
706 (abs(eta) > 2.5) * (0.00)}
707}
708
709####################
710# Electron isolation
711####################
712
713module Isolation ElectronIsolation {
714 set CandidateInputArray ElectronEfficiency/electrons
715 set IsolationInputArray EFlowFilter/eflow
716 set RhoInputArray Rho/rho
717
718 set OutputArray electrons
719
720 set DeltaRMax 0.4
721
722 set PTMin 0.5
723
724 set PTRatioMax 0.12
725}
726
727#################
728# Muon efficiency
729#################
730
731module Efficiency MuonEfficiency {
732 set InputArray TrackPileUpSubtractor/muons
733 set OutputArray muons
734
735 # set EfficiencyFormula {efficiency as a function of eta and pt}
736
737#
738
739#set EfficiencyFormula { (pt <= 5.0) * (0.00) +
740# (abs(eta) <= 2.4) * (pt > 5.0 && pt < 60) * (0.70 + 0.2*(pt-5)/55.) +
741# (abs(eta) <= 2.4) * (pt > 60.0) * (0.90) +
742# (abs(eta) > 2.4) * (0.00)}
743
744
745 set EfficiencyFormula { (pt <= 5.0) * (0.00) +
746 (abs(eta) <= 1.5) * (pt > 5.0) * (0.95) +
747 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 5.0) * (0.95) +
748 (abs(eta) > 2.4) * (0.00)}
749
750# set EfficiencyFormula { (pt <= 5.0) * (0.00) +
751# (abs(eta) <= 1.5) * (pt > 5.0) * (0.90) +
752# (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 5.0) * (0.90) +
753# (abs(eta) > 2.4) * (0.00)}
754
755
756
757 # efficiency formula for muons # Previously 95
758# set EfficiencyFormula { (pt <= 5.0) * (0.00) +
759# (abs(eta) <= 1.5) * (pt > 5.0) * (0.98) +
760# (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 5.0) * (0.98) +
761# (abs(eta) > 2.4) * (0.00)}
762
763}
764
765################
766#a Muon isolation
767################
768
769module Isolation MuonIsolation {
770 set CandidateInputArray MuonEfficiency/muons
771 set IsolationInputArray EFlowFilter/eflow
772 set RhoInputArray Rho/rho
773
774 set OutputArray muons
775
776 set DeltaRMax 0.4
777
778 set PTMin 0.5
779
780 set PTRatioMax 0.25
781}
782
783###################
784# Missing ET merger
785###################
786
787module Merger MissingET {
788# add InputArray InputArray
789 add InputArray EFlowMergerAllTracks/eflow
790 set MomentumOutputArray momentum
791}
792
793
794
795##################
796# Scalar HT merger
797##################
798
799module Merger ScalarHT {
800# add InputArray InputArray
801 add InputArray UniqueObjectFinder/jets
802 add InputArray UniqueObjectFinder/electrons
803 add InputArray UniqueObjectFinder/photons
804 add InputArray UniqueObjectFinder/muons
805 set EnergyOutputArray energy
806}
807
808########################
809# Jet Flavor Association
810########################
811
812module JetFlavorAssociation JetFlavorAssociation {
813
814 set PartonInputArray Delphes/partons
815 set ParticleInputArray Delphes/allParticles
816 set ParticleLHEFInputArray Delphes/allParticlesLHEF
817 set JetInputArray JetEnergyScale/jets
818
819 set DeltaR 0.4
820 set PartonPTMin 1.0
821 set PartonEtaMax 2.5
822
823}
824
825###########
826# b-tagging
827###########
828
829module BTagging BTagging {
830 set JetInputArray JetEnergyScale/jets
831
832 set BitNumber 0
833
834 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
835 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
836 # gluon's PDG code has the lowest priority
837
838#CSVv2 efficiencies
839
840# add EfficiencyFormula {0} { (pt > 20.0 && pt < 105.0) * (-1.770*10^-9*pt^4 + 4.492*10^-7*pt^3-0.00003935*pt^2+0.00136*pt-0.006) +
841# (pt >= 105.0 && pt < 235.0) * (-7.638*10^-12*pt^5 + 6.486*10^-9*pt^4 -2.15015*10^-6*pt^3 + 0.000347053*pt^2 - 0.0271883*pt + 0.8345)+
842# (pt >= 235.0) * (3.9855*10^-17*pt^6 - 1.3094*10^-13*pt^5 + 1.7179*10^-10*pt^4 -
843# 1.14705*10^-7*pt^3 + 0.0000409144*pt^2 - 0.00733112*pt + 0.52874) }
844
845# add EfficiencyFormula {4} { (pt > 20.0 && pt <45.0) * (-0.000159*pt^2+0.0117*pt-0.087) +
846# (pt >= 45.0 && pt <145.0) * (2.76*10^-8 * pt^3 - 6.91*10^-6*pt^2 + 0.000723*pt+0.09675) +
847# (pt >= 145.0 ) * (5.21785*10^-17*pt^6 - 1.60294*10^-13*pt^5 + 1.94745*10^-10*pt^4 -
848# 1.19154*10^-7*pt^3 + 0.0000386176*pt^2 - 0.0063131*pt + 0.532) }
849
850# add EfficiencyFormula {5} { (pt > 20.0 && pt < 80.0)*(-7.649*10^-8*pt^4 + 0.00001932*pt^3 - 0.001797*pt^2 + 0.07384*pt - 0.5) +
851# (pt >= 80.0 && pt < 143.0) * (0.665) +
852# (pt >= 143.0) * (-3.628*10^-15*pt^5 + 8.607*10^-12*pt^4 - 7.86*10^-9*pt^3 + 3.92*10^-6*pt^2 - 0.001468*pt + 0.8147)}
853
854# DeepCSV efficiencies
855
856 add EfficiencyFormula {0} { (pt >= 25.0 && pt < 360.0) * (-8.56347 * 10^-17 * pt^6 + 4.81506 * 10^-14 * pt^5 + 1.26688 * 10^-11 * pt^4 - 1.43111 * 10^-8 * pt^3 + 3.66113 * 10^-6 * pt^2 - 0.000343664 * pt + 0.0194182) +
857 (pt >= 360.0 && pt < 950 ) * (2.47657 * 10^-17 * pt^6 - 9.79073 * 10^-14 * pt^5 + 1.56687 * 10^-10 * pt^4 - 1.29616 * 10^-7 * pt^3 + 0.0000582789 * pt^2 - 0.0134288 * pt + 1.25506) + (pt >= 950 ) * (0.035)}
858
859 add EfficiencyFormula {4} { (pt >= 25.0 && pt <155.0) * (2.49952*10^-11*pt^5 - 1.19968 * 10^-8 * pt^4 + 2.15425 * 10^-6 * pt^3 - 0.000174775 * pt^2 + 0.00603344 * pt + 0.0571752) +
860 (pt >= 155.0 && pt <318.0) * (1.14416*10^-13 * pt^6 - 1.62168 * 10^-10 * pt^5 + 9.434 * 10^-8 * pt^4 - 0.000028816 * pt^3 + 0.00487191 * pt^2 - 0.432125 * pt + 15.8347) +
861 (pt >= 318.0 && pt < 950) * (1.10863*10^-16 * pt^5 + 7.37806 * 10^-13 * pt^4 - 1.96006 * 10^-9 * pt^3 + 1.36171 * 10^-6 * pt^2 - 0.000225177 * pt + 0.118996) +
862 (pt >= 950 ) * (0.14)}
863
864 add EfficiencyFormula {5} { (pt > 25.0 && pt < 115.0)*(-1.20579*10^-8 * pt^4+4.07131*10^-6 * pt^3-0.00050732 * pt^2+0.0279171 * pt + 0.1297235) +
865 (pt >= 115.0 && pt < 950 ) * (1.11352*10^-14 * pt^5 - 2.78119 * 10^-11 * pt^4 + 2.55081 * 10^-8 * pt^3 - 0.000010315 * pt^2 + 0.00142765 * pt + 0.651055)+
866 (pt >= 950 ) * (0.50)}
867
868
869
870# 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) +
871# (pt >= 130.0) * (0.055 + 4.53*10^-4*pt - 1.60*10^-7*pt^2) }
872
873# 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) +
874# (pt >= 205.0) * (0.478 + 1.573*10^-4*pt)}
875
876# 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) +
877# (pt >= 150.0) * (0.906 - 6.39*10^-5*pt + 4.11*10^-8*pt^2) }
878}
879# tau-tagging
880#############
881
882module TauTagging TauTagging {
883 set ParticleInputArray Delphes/allParticles
884 set PartonInputArray Delphes/partons
885 set JetInputArray JetEnergyScale/jets
886
887 set DeltaR 0.4
888
889 set TauPTMin 1.0
890
891 set TauEtaMax 2.5
892
893 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
894
895 # default efficiency formula (misidentification rate)
896 add EfficiencyFormula {0} {0.01}
897 # efficiency formula for tau-jets
898 add EfficiencyFormula {15} {0.6}
899}
900
901#####################################################
902# Find uniquely identified photons/electrons/tau/jets
903#####################################################
904
905module UniqueObjectFinder UniqueObjectFinder {
906# earlier arrays take precedence over later ones
907# add InputArray InputArray OutputArray
908 add InputArray PhotonIsolation/photons photons
909 add InputArray ElectronIsolation/electrons electrons
910 add InputArray MuonIsolation/muons muons
911 add InputArray JetEnergyScale/jets jets
912}
913
914##################
915# ROOT tree writer
916##################
917
918# tracks, towers and eflow objects are not stored by default in the output.
919# if needed (for jet constituent or other studies), uncomment the relevant
920# "add Branch ..." lines.
921
922module TreeWriter TreeWriter {
923# add Branch InputArray BranchName BranchClass
924 add Branch Delphes/allParticles Particle GenParticle
925
926# add Branch TrackMerger/tracks Track Track
927# add Branch Calorimeter/towers Tower Tower
928
929# add Branch Calorimeter/eflowTracks EFlowTrack Track
930# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
931# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
932
933 add Branch GenJetFinder/jets GenJet Jet
934 add Branch GenMissingET/momentum GenMissingET MissingET
935
936 add Branch UniqueObjectFinder/jets Jet Jet
937 add Branch UniqueObjectFinder/electrons Electron Electron
938 add Branch UniqueObjectFinder/photons Photon Photon
939 add Branch UniqueObjectFinder/muons Muon Muon
940 add Branch MissingET/momentum MissingET MissingET
941 add Branch ScalarHT/energy ScalarHT ScalarHT
942 add Branch Rho/rho Rho Rho
943 add Branch PileUpMerger/vertices Vertex Vertex
944}