MA5SandBox: delphes_card_cms_sus_16_039.tcl

File delphes_card_cms_sus_16_039.tcl, 28.0 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 {1000022} {0.0}
287 add EnergyFraction {1000023} {0.0}
288 add EnergyFraction {1000025} {0.0}
289 add EnergyFraction {1000035} {0.0}
290 add EnergyFraction {1000045} {0.0}
291 # energy fractions for K0short and Lambda
292 add EnergyFraction {310} {0.3}
293 add EnergyFraction {3122} {0.3}
294
295 # set ResolutionFormula {resolution formula as a function of eta and energy}
296
297 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
298
299 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
300 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
301 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
302 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
303 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
304
305}
306
307
308#############
309# HCAL
310#############
311
312module SimpleCalorimeter HCal {
313 set ParticleInputArray ParticlePropagator/stableParticles
314 set TrackInputArray ECal/eflowTracks
315
316 set TowerOutputArray hcalTowers
317 set EFlowTrackOutputArray eflowTracks
318 set EFlowTowerOutputArray eflowNeutralHadrons
319
320 set IsEcal false
321
322 set EnergyMin 1.0
323 set EnergySignificanceMin 1.0
324
325 set SmearTowerCenter true
326
327 set pi [expr {acos(-1)}]
328
329 # lists of the edges of each tower in eta and phi
330 # each list starts with the lower edge of the first tower
331 # the list ends with the higher edged of the last tower
332
333 # 5 degrees towers
334 set PhiBins {}
335 for {set i -36} {$i <= 36} {incr i} {
336 add PhiBins [expr {$i * $pi/36.0}]
337 }
338 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
339 add EtaPhiBins $eta $PhiBins
340 }
341
342 # 10 degrees towers
343 set PhiBins {}
344 for {set i -18} {$i <= 18} {incr i} {
345 add PhiBins [expr {$i * $pi/18.0}]
346 }
347 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
348 add EtaPhiBins $eta $PhiBins
349 }
350
351 # 20 degrees towers
352 set PhiBins {}
353 for {set i -9} {$i <= 9} {incr i} {
354 add PhiBins [expr {$i * $pi/9.0}]
355 }
356 foreach eta {-5 -4.7 -4.525 4.7 5} {
357 add EtaPhiBins $eta $PhiBins
358 }
359
360 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
361 add EnergyFraction {0} {1.0}
362 # energy fractions for e, gamma and pi0
363 add EnergyFraction {11} {0.0}
364 add EnergyFraction {22} {0.0}
365 add EnergyFraction {111} {0.0}
366 # energy fractions for muon, neutrinos and neutralinos
367 add EnergyFraction {12} {0.0}
368 add EnergyFraction {13} {0.0}
369 add EnergyFraction {14} {0.0}
370 add EnergyFraction {16} {0.0}
371 add EnergyFraction {1000022} {0.0}
372 add EnergyFraction {1000023} {0.0}
373 add EnergyFraction {1000025} {0.0}
374 add EnergyFraction {1000035} {0.0}
375 add EnergyFraction {1000045} {0.0}
376 # energy fractions for K0short and Lambda
377 add EnergyFraction {310} {0.7}
378 add EnergyFraction {3122} {0.7}
379
380 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
381 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
382 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
383
384}
385
386
387#################
388# Electron filter
389#################
390
391module PdgCodeFilter ElectronFilter {
392 set InputArray HCal/eflowTracks
393 set OutputArray electrons
394 set Invert true
395 add PdgCode {11}
396 add PdgCode {-11}
397}
398
399######################
400# ChargedHadronFilter
401######################
402
403module PdgCodeFilter ChargedHadronFilter {
404 set InputArray HCal/eflowTracks
405 set OutputArray chargedHadrons
406
407 add PdgCode {11}
408 add PdgCode {-11}
409 add PdgCode {13}
410 add PdgCode {-13}
411}
412
413
414###################################################
415# Tower Merger (in case not using e-flow algorithm)
416###################################################
417
418module Merger Calorimeter {
419# add InputArray InputArray
420 add InputArray ECal/ecalTowers
421 add InputArray HCal/hcalTowers
422 set OutputArray towers
423}
424
425
426
427####################
428# Energy flow merger
429####################
430
431module Merger EFlowMerger {
432# add InputArray InputArray
433 add InputArray HCal/eflowTracks
434 add InputArray ECal/eflowPhotons
435 add InputArray HCal/eflowNeutralHadrons
436 set OutputArray eflow
437}
438
439######################
440# EFlowFilter
441######################
442
443module PdgCodeFilter EFlowFilter {
444 set InputArray EFlowMerger/eflow
445 set OutputArray eflow
446
447 add PdgCode {11}
448 add PdgCode {-11}
449 add PdgCode {13}
450 add PdgCode {-13}
451}
452
453
454###################
455# Photon efficiency
456###################
457
458module Efficiency PhotonEfficiency {
459 set InputArray ECal/eflowPhotons
460 set OutputArray photons
461
462 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
463
464 # efficiency formula for photons
465 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
466 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
467 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
468 (abs(eta) > 2.5) * (0.00)}
469}
470
471##################
472# Photon isolation
473##################
474
475module Isolation PhotonIsolation {
476 set CandidateInputArray PhotonEfficiency/photons
477 set IsolationInputArray EFlowFilter/eflow
478
479 set OutputArray photons
480
481 set DeltaRMax 0.5
482
483 set PTMin 0.5
484
485 set PTRatioMax 0.12
486}
487
488
489#####################
490# Electron efficiency
491#####################
492
493module Efficiency ElectronEfficiency {
494 set InputArray ElectronFilter/electrons
495 set OutputArray electrons
496
497 set EfficiencyFormula { (pt <= 15.0) * (0.00) +
498 (abs(eta) <= 0.8) * (pt > 15.0 && pt <= 20.0) * (0.507) +
499 (abs(eta) <= 0.8) * (pt > 20.0 && pt <= 25.0) * (0.619) +
500 (abs(eta) <= 0.8) * (pt > 25.0 && pt <= 30.0) * (0.682) +
501 (abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.742) +
502 (abs(eta) <= 0.8) * (pt > 40.0 && pt <= 50.0) * (0.798) +
503 (abs(eta) <= 0.8) * (pt > 50.0 && pt <= 200.0) * (0.863) +
504 (abs(eta) <= 0.8) * (pt > 200) * (0.863) +
505 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 15.0 && pt <= 20.0) * (0.429) +
506 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 20.0 && pt <= 25.0) * (0.546) +
507 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 25.0 && pt <= 30.0) * (0.619) +
508 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 30.0 && pt <= 40.0) * (0.710) +
509 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 40.0 && pt <= 50.0) * (0.734) +
510 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 50.0 && pt <= 200.0) * (0.833) +
511 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 200) * (0.833) +
512 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 15.0 && pt <= 20.0) * (0.256) +
513 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 20.0 && pt <= 25.0) * (0.221) +
514 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 25.0 && pt <= 30.0) * (0.315) +
515 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 30.0 && pt <= 40.0) * (0.351) +
516 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 40.0 && pt <= 50.0) * (0.373) +
517 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 50.0 && pt <= 200.0) * (0.437) +
518 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 200) * (0.437) +
519 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 15.0 && pt <= 20.0) * (0.249) +
520 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 25.0) * (0.404) +
521 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 25.0 && pt <= 30.0) * (0.423) +
522 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.561) +
523 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 50.0) * (0.642) +
524 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 200.0) * (0.749) +
525 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 200) * (0.749) +
526 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 15.0 && pt <= 20.0) * (0.195) +
527 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 25.0) * (0.245) +
528 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 25.0 && pt <= 30.0) * (0.380) +
529 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.441) +
530 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 50.0) * (0.533) +
531 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 200.0) * (0.644) +
532 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 200) * (0.644) +
533 (abs(eta) > 2.5) * (0.00)}
534}
535
536
537##
538## Modified by JB 20171017
539## Using CMS efficiency formula : https://twiki.cern.ch/twiki/bin/view/CMSPublic/SUSMoriond2017ObjectsEfficiency
540##
541
542####################
543# Electron isolation
544####################
545
546module Isolation ElectronIsolation {
547 set CandidateInputArray ElectronEfficiency/electrons
548 set IsolationInputArray EFlowFilter/eflow
549
550 set OutputArray electrons
551
552 set DeltaRMax 0.5
553
554 set PTMin 0.5
555
556 set PTRatioMax 0.12
557}
558
559#################
560# Muon efficiency
561#################
562
563module Efficiency MuonEfficiency {
564 set InputArray MuonMomentumSmearing/muons
565 set OutputArray muons
566
567 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
568 (abs(eta) <= 0.9) * (pt > 10.0 && pt <= 15.0) * (0.704) +
569 (abs(eta) <= 0.9) * (pt > 15.0 && pt <= 20.0) * (0.797) +
570 (abs(eta) <= 0.9) * (pt > 20.0 && pt <= 25.0) * (0.855) +
571 (abs(eta) <= 0.9) * (pt > 25.0 && pt <= 30.0) * (0.880) +
572 (abs(eta) <= 0.9) * (pt > 30.0 && pt <= 40.0) * (0.906) +
573 (abs(eta) <= 0.9) * (pt > 40.0 && pt <= 50.0) * (0.927) +
574 (abs(eta) <= 0.9) * (pt > 50.0 && pt <= 200.0) * (0.931) +
575 (abs(eta) <= 0.9) * (pt > 200) * (0.931) +
576 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 10.0 && pt <= 15.0) * (0.639) +
577 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 15.0 && pt <= 20.0) * (0.776) +
578 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 20.0 && pt <= 25.0) * (0.836) +
579 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 25.0 && pt <= 30.0) * (0.875) +
580 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 30.0 && pt <= 40.0) * (0.898) +
581 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 40.0 && pt <= 50.0) * (0.940) +
582 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 50.0 && pt <= 200.0) * (0.930) +
583 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 200) * (0.930) +
584 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 10.0 && pt <= 15.0) * (0.596) +
585 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 15.0 && pt <= 20.0) * (0.715) +
586 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 20.0 && pt <= 25.0) * (0.840) +
587 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 25.0 && pt <= 30.0) * (0.862) +
588 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 30.0 && pt <= 40.0) * (0.891) +
589 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 40.0 && pt <= 50.0) * (0.906) +
590 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 50.0 && pt <= 200.0) * (0.925) +
591 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 200) * (0.925) +
592 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 15.0) * (0.522) +
593 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 15.0 && pt <= 20.0) * (0.720) +
594 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 25.0) * (0.764) +
595 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 25.0 && pt <= 30.0) * (0.803) +
596 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40.0) * (0.807) +
597 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50.0) * (0.885) +
598 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 200.0) * (0.877) +
599 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 200) * (0.877) +
600 (abs(eta) > 2.5) * (0.00)}
601}
602
603##
604## Modified by JB 20171017
605## Using CMS efficiency formula : https://twiki.cern.ch/twiki/bin/view/CMSPublic/SUSMoriond2017ObjectsEfficiency
606##
607
608################
609# Muon isolation
610################
611
612module Isolation MuonIsolation {
613 set CandidateInputArray MuonEfficiency/muons
614 set IsolationInputArray EFlowFilter/eflow
615
616 set OutputArray muons
617
618 set DeltaRMax 0.5
619
620 set PTMin 0.5
621
622 set PTRatioMax 0.25
623}
624
625###################
626# Missing ET merger
627###################
628
629module Merger MissingET {
630# add InputArray InputArray
631 add InputArray EFlowMerger/eflow
632 set MomentumOutputArray momentum
633}
634
635##################
636# Scalar HT merger
637##################
638
639module Merger ScalarHT {
640# add InputArray InputArray
641 add InputArray UniqueObjectFinder/jets
642 add InputArray UniqueObjectFinder/electrons
643 add InputArray UniqueObjectFinder/photons
644 add InputArray UniqueObjectFinder/muons
645 set EnergyOutputArray energy
646}
647
648
649#####################
650# Neutrino Filter
651#####################
652
653module PdgCodeFilter NeutrinoFilter {
654
655 set InputArray Delphes/stableParticles
656 set OutputArray filteredParticles
657
658 set PTMin 0.0
659
660 add PdgCode {12}
661 add PdgCode {14}
662 add PdgCode {16}
663 add PdgCode {-12}
664 add PdgCode {-14}
665 add PdgCode {-16}
666 add PdgCode {1000022}
667
668}
669
670
671#####################
672# MC truth jet finder
673#####################
674
675module FastJetFinder GenJetFinder {
676 set InputArray NeutrinoFilter/filteredParticles
677
678 set OutputArray jets
679
680 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
681 set JetAlgorithm 6
682 set ParameterR 0.4
683
684 set JetPTMin 20.0
685}
686
687#########################
688# Gen Missing ET merger
689########################
690
691module Merger GenMissingET {
692# add InputArray InputArray
693 add InputArray NeutrinoFilter/filteredParticles
694 set MomentumOutputArray momentum
695}
696
697
698
699############
700# Jet finder
701############
702
703module FastJetFinder FastJetFinder {
704# set InputArray Calorimeter/towers
705 set InputArray EFlowMerger/eflow
706
707 set OutputArray jets
708
709 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
710 set JetAlgorithm 6
711 set ParameterR 0.4
712
713 set JetPTMin 20.0
714}
715
716##################
717# Fat Jet finder
718##################
719
720module FastJetFinder FatJetFinder {
721 set InputArray EFlowMerger/eflow
722
723 set OutputArray jets
724
725 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
726 set JetAlgorithm 6
727 set ParameterR 0.8
728
729 set ComputeNsubjettiness 1
730 set Beta 1.0
731 set AxisMode 4
732
733 set ComputeTrimming 1
734 set RTrim 0.2
735 set PtFracTrim 0.05
736
737 set ComputePruning 1
738 set ZcutPrun 0.1
739 set RcutPrun 0.5
740 set RPrun 0.8
741
742 set ComputeSoftDrop 1
743 set BetaSoftDrop 0.0
744 set SymmetryCutSoftDrop 0.1
745 set R0SoftDrop 0.8
746
747 set JetPTMin 200.0
748}
749
750
751
752
753##################
754# Jet Energy Scale
755##################
756
757module EnergyScale JetEnergyScale {
758 set InputArray FastJetFinder/jets
759 set OutputArray jets
760
761 # scale formula for jets
762 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
763}
764
765########################
766# Jet Flavor Association
767########################
768
769module JetFlavorAssociation JetFlavorAssociation {
770
771 set PartonInputArray Delphes/partons
772 set ParticleInputArray Delphes/allParticles
773 set ParticleLHEFInputArray Delphes/allParticlesLHEF
774 set JetInputArray JetEnergyScale/jets
775
776 set DeltaR 0.5
777 set PartonPTMin 1.0
778 set PartonEtaMax 2.5
779
780}
781
782###########
783# b-tagging
784###########
785
786module BTagging BTagging {
787 set JetInputArray JetEnergyScale/jets
788
789 set BitNumber 0
790
791 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
792 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
793 # gluon's PDG code has the lowest priority
794 #
795 # PDG code 4 = c quark
796 # 5 = b quark
797 # 0 = light flavour
798 #
799 # overall factor 0.94 used by CMS
800
801 add EfficiencyFormula {0} { 0.94*(pt >= 30.0 && pt < 170.0) * (9.59*10^(-3) - 1.96*10^(-5)*pt + 4.53*10^(-7)*pt^2 - 1.08*10^(-9)*pt^3 + 7.62*10^(-13)*pt^4) +
802 0.94*(pt >= 170.0) * (5.07*10^(-3) + 6.02*10^(-5)*pt - 2.30*10^(-8)*pt^2) }
803
804 add EfficiencyFormula {4} { 0.94*(pt >= 30.0 && pt < 170.0) * (0.13 + 1.48*10^(-3)*pt - 1.00*10^(-5)*pt^2 + 2.65*10^(-8)*pt^3) - 2.36*10^(-11)*pt^4 +
805 0.94*(pt >= 170.0) * (0.20)}
806
807 add EfficiencyFormula {5} { 0.94*(pt >= 30.0 && pt < 175.0) * (0.421 + 0.0107*pt - 1.314*10^(-4)*pt^2 + 7.628*10^(-7)*pt^3 - 1.523*10^(-9)*pt^4) +
808 0.94*(pt >= 175.0) * (0.79 - 3.17*10^(-4)*pt + 1.24*10^(-7)*pt^2) }
809}
810
811#############
812# tau-tagging
813#############
814
815module TauTagging TauTagging {
816 set ParticleInputArray Delphes/allParticles
817 set PartonInputArray Delphes/partons
818 set JetInputArray JetEnergyScale/jets
819
820 set DeltaR 0.5
821
822 set TauPTMin 1.0
823
824 set TauEtaMax 2.5
825
826 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
827
828 # default efficiency formula (misidentification rate)
829 add EfficiencyFormula {0} {0.01}
830 # efficiency formula for tau-jets
831 add EfficiencyFormula {15} {0.6}
832}
833
834#####################################################
835# Find uniquely identified photons/electrons/tau/jets
836#####################################################
837
838module UniqueObjectFinder UniqueObjectFinder {
839# earlier arrays take precedence over later ones
840# add InputArray InputArray OutputArray
841 add InputArray PhotonIsolation/photons photons
842 add InputArray ElectronIsolation/electrons electrons
843 add InputArray MuonIsolation/muons muons
844 add InputArray JetEnergyScale/jets jets
845}
846
847##################
848# ROOT tree writer
849##################
850
851# tracks, towers and eflow objects are not stored by default in the output.
852# if needed (for jet constituent or other studies), uncomment the relevant
853# "add Branch ..." lines.
854
855module TreeWriter TreeWriter {
856# add Branch InputArray BranchName BranchClass
857 add Branch Delphes/allParticles Particle GenParticle
858
859 add Branch TrackMerger/tracks Track Track
860 add Branch Calorimeter/towers Tower Tower
861
862 add Branch HCal/eflowTracks EFlowTrack Track
863 add Branch ECal/eflowPhotons EFlowPhoton Tower
864 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
865
866 add Branch GenJetFinder/jets GenJet Jet
867 add Branch GenMissingET/momentum GenMissingET MissingET
868
869 add Branch UniqueObjectFinder/jets Jet Jet
870 add Branch UniqueObjectFinder/electrons Electron Electron
871 add Branch UniqueObjectFinder/photons Photon Photon
872 add Branch UniqueObjectFinder/muons Muon Muon
873
874 add Branch FatJetFinder/jets FatJet Jet
875
876 add Branch MissingET/momentum MissingET MissingET
877 add Branch ScalarHT/energy ScalarHT ScalarHT
878}