MA5SandBox: delphes_card_cms_sus_16_041.tcl

File delphes_card_cms_sus_16_041.tcl, 29.3 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 <= 10.0) * (0.00) +
498 (abs(eta) <= 0.8) * (pt > 10.0 && pt <= 15.0) * (0.320) +
499 (abs(eta) <= 0.8) * (pt > 15.0 && pt <= 20.0) * (0.450) +
500 (abs(eta) <= 0.8) * (pt > 20.0 && pt <= 25.0) * (0.555) +
501 (abs(eta) <= 0.8) * (pt > 25.0 && pt <= 30.0) * (0.606) +
502 (abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.671) +
503 (abs(eta) <= 0.8) * (pt > 40.0 && pt <= 50.0) * (0.712) +
504 (abs(eta) <= 0.8) * (pt > 50.0 && pt <= 200.0) * (0.760) +
505 (abs(eta) <= 0.8) * (pt > 200) * (0.760) +
506 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 10.0 && pt <= 15.0) * (0.268) +
507 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 15.0 && pt <= 20.0) * (0.396) +
508 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 20.0 && pt <= 25.0) * (0.500) +
509 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 25.0 && pt <= 30.0) * (0.568) +
510 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 30.0 && pt <= 40.0) * (0.653) +
511 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 40.0 && pt <= 50.0) * (0.669) +
512 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 50.0 && pt <= 200.0) * (0.738) +
513 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 200) * (0.738) +
514 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 10.0 && pt <= 15.0) * (0.083) +
515 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 15.0 && pt <= 20.0) * (0.233) +
516 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 20.0 && pt <= 25.0) * (0.193) +
517 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 25.0 && pt <= 30.0) * (0.242) +
518 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 30.0 && pt <= 40.0) * (0.274) +
519 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 40.0 && pt <= 50.0) * (0.319) +
520 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 50.0 && pt <= 200.0) * (0.376) +
521 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 200) * (0.376) +
522 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 15.0) * (0.143) +
523 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 15.0 && pt <= 20.0) * (0.250) +
524 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 25.0) * (0.380) +
525 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 25.0 && pt <= 30.0) * (0.413) +
526 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.525) +
527 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 50.0) * (0.568) +
528 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 200.0) * (0.649) +
529 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 200) * (0.649) +
530 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 15.0) * (0.130) +
531 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 15.0 && pt <= 20.0) * (0.220) +
532 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 25.0) * (0.222) +
533 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 25.0 && pt <= 30.0) * (0.370) +
534 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.455) +
535 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 50.0) * (0.517) +
536 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 200.0) * (0.559) +
537 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 200) * (0.559) +
538 (abs(eta) > 2.5) * (0.00)}
539}
540
541
542##
543## Modified by JB 20171017
544## Using CMS efficiency formula : https://twiki.cern.ch/twiki/bin/view/CMSPublic/SUSMoriond2017ObjectsEfficiency
545##
546
547####################
548# Electron isolation
549####################
550
551module Isolation ElectronIsolation {
552 set CandidateInputArray ElectronEfficiency/electrons
553 set IsolationInputArray EFlowFilter/eflow
554
555 set OutputArray electrons
556
557 set DeltaRMax 0.5
558
559 set PTMin 0.5
560
561 set PTRatioMax 0.12
562}
563
564#################
565# Muon efficiency
566#################
567
568module Efficiency MuonEfficiency {
569 set InputArray MuonMomentumSmearing/muons
570 set OutputArray muons
571
572 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
573 (abs(eta) <= 0.9) * (pt > 10.0 && pt <= 15.0) * (0.639) +
574 (abs(eta) <= 0.9) * (pt > 15.0 && pt <= 20.0) * (0.724) +
575 (abs(eta) <= 0.9) * (pt > 20.0 && pt <= 25.0) * (0.804) +
576 (abs(eta) <= 0.9) * (pt > 25.0 && pt <= 30.0) * (0.841) +
577 (abs(eta) <= 0.9) * (pt > 30.0 && pt <= 40.0) * (0.879) +
578 (abs(eta) <= 0.9) * (pt > 40.0 && pt <= 50.0) * (0.908) +
579 (abs(eta) <= 0.9) * (pt > 50.0 && pt <= 200.0) * (0.913) +
580 (abs(eta) <= 0.9) * (pt > 200) * (0.913) +
581 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 10.0 && pt <= 15.0) * (0.587) +
582 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 15.0 && pt <= 20.0) * (0.698) +
583 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 20.0 && pt <= 25.0) * (0.755) +
584 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 25.0 && pt <= 30.0) * (0.822) +
585 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 30.0 && pt <= 40.0) * (0.850) +
586 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 40.0 && pt <= 50.0) * (0.905) +
587 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 50.0 && pt <= 200.0) * (0.904) +
588 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 200) * (0.904) +
589 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 10.0 && pt <= 15.0) * (0.592) +
590 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 15.0 && pt <= 20.0) * (0.660) +
591 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 20.0 && pt <= 25.0) * (0.772) +
592 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 25.0 && pt <= 30.0) * (0.805) +
593 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 30.0 && pt <= 40.0) * (0.828) +
594 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 40.0 && pt <= 50.0) * (0.856) +
595 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 50.0 && pt <= 200.0) * (0.877) +
596 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 200) * (0.877) +
597 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 15.0) * (0.531) +
598 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 15.0 && pt <= 20.0) * (0.661) +
599 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 25.0) * (0.687) +
600 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 25.0 && pt <= 30.0) * (0.740) +
601 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40.0) * (0.723) +
602 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50.0) * (0.793) +
603 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 200.0) * (0.806) +
604 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 200) * (0.806) +
605 (abs(eta) > 2.5) * (0.00)}
606}
607
608##
609## Modified by JB 20171017
610## Using CMS efficiency formula : https://twiki.cern.ch/twiki/bin/view/CMSPublic/SUSMoriond2017ObjectsEfficiency
611##
612
613################
614# Muon isolation
615################
616
617module Isolation MuonIsolation {
618 set CandidateInputArray MuonEfficiency/muons
619 set IsolationInputArray EFlowFilter/eflow
620
621 set OutputArray muons
622
623 set DeltaRMax 0.5
624
625 set PTMin 0.5
626
627 set PTRatioMax 0.25
628}
629
630###################
631# Missing ET merger
632###################
633
634module Merger MissingET {
635# add InputArray InputArray
636 add InputArray EFlowMerger/eflow
637 set MomentumOutputArray momentum
638}
639
640##################
641# Scalar HT merger
642##################
643
644module Merger ScalarHT {
645# add InputArray InputArray
646 add InputArray UniqueObjectFinder/jets
647 add InputArray UniqueObjectFinder/electrons
648 add InputArray UniqueObjectFinder/photons
649 add InputArray UniqueObjectFinder/muons
650 set EnergyOutputArray energy
651}
652
653
654#####################
655# Neutrino Filter
656#####################
657
658module PdgCodeFilter NeutrinoFilter {
659
660 set InputArray Delphes/stableParticles
661 set OutputArray filteredParticles
662
663 set PTMin 0.0
664
665 add PdgCode {12}
666 add PdgCode {14}
667 add PdgCode {16}
668 add PdgCode {-12}
669 add PdgCode {-14}
670 add PdgCode {-16}
671 add PdgCode {1000022}
672
673}
674
675
676#####################
677# MC truth jet finder
678#####################
679
680module FastJetFinder GenJetFinder {
681 set InputArray NeutrinoFilter/filteredParticles
682
683 set OutputArray jets
684
685 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
686 set JetAlgorithm 6
687 set ParameterR 0.4
688
689 set JetPTMin 20.0
690}
691
692#########################
693# Gen Missing ET merger
694########################
695
696module Merger GenMissingET {
697# add InputArray InputArray
698 add InputArray NeutrinoFilter/filteredParticles
699 set MomentumOutputArray momentum
700}
701
702
703
704############
705# Jet finder
706############
707
708module FastJetFinder FastJetFinder {
709# set InputArray Calorimeter/towers
710 set InputArray EFlowMerger/eflow
711
712 set OutputArray jets
713
714 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
715 set JetAlgorithm 6
716 set ParameterR 0.4
717
718 set JetPTMin 20.0
719}
720
721##################
722# Fat Jet finder
723##################
724
725module FastJetFinder FatJetFinder {
726 set InputArray EFlowMerger/eflow
727
728 set OutputArray jets
729
730 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
731 set JetAlgorithm 6
732 set ParameterR 0.8
733
734 set ComputeNsubjettiness 1
735 set Beta 1.0
736 set AxisMode 4
737
738 set ComputeTrimming 1
739 set RTrim 0.2
740 set PtFracTrim 0.05
741
742 set ComputePruning 1
743 set ZcutPrun 0.1
744 set RcutPrun 0.5
745 set RPrun 0.8
746
747 set ComputeSoftDrop 1
748 set BetaSoftDrop 0.0
749 set SymmetryCutSoftDrop 0.1
750 set R0SoftDrop 0.8
751
752 set JetPTMin 200.0
753}
754
755
756
757
758##################
759# Jet Energy Scale
760##################
761
762module EnergyScale JetEnergyScale {
763 set InputArray FastJetFinder/jets
764 set OutputArray jets
765
766 # scale formula for jets
767 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
768}
769
770########################
771# Jet Flavor Association
772########################
773
774module JetFlavorAssociation JetFlavorAssociation {
775
776 set PartonInputArray Delphes/partons
777 set ParticleInputArray Delphes/allParticles
778 set ParticleLHEFInputArray Delphes/allParticlesLHEF
779 set JetInputArray JetEnergyScale/jets
780
781 set DeltaR 0.5
782 set PartonPTMin 1.0
783 set PartonEtaMax 2.5
784
785}
786
787###########
788# b-tagging
789###########
790
791module BTagging BTagging {
792 set JetInputArray JetEnergyScale/jets
793
794 set BitNumber 0
795
796 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
797 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
798 # gluon's PDG code has the lowest priority
799 #
800 # PDG code 4 = c quark
801 # 5 = b quark
802 # 0 = light flavour
803 #
804 # overall factor 0.94 used by CMS
805
806 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) +
807 0.94*(pt >= 170.0) * (5.07*10^(-3) + 6.02*10^(-5)*pt - 2.30*10^(-8)*pt^2) }
808
809 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 +
810 0.94*(pt >= 170.0) * (0.20)}
811
812 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) +
813 0.94*(pt >= 175.0) * (0.79 - 3.17*10^(-4)*pt + 1.24*10^(-7)*pt^2) }
814}
815
816#############
817# tau-tagging
818#############
819
820module TauTagging TauTagging {
821 set ParticleInputArray Delphes/allParticles
822 set PartonInputArray Delphes/partons
823 set JetInputArray JetEnergyScale/jets
824
825 set DeltaR 0.5
826
827 set TauPTMin 1.0
828
829 set TauEtaMax 2.5
830
831 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
832
833 # default efficiency formula (misidentification rate)
834 add EfficiencyFormula {0} {0.01}
835 # efficiency formula for tau-jets
836 add EfficiencyFormula {15} {0.6}
837}
838
839#####################################################
840# Find uniquely identified photons/electrons/tau/jets
841#####################################################
842
843module UniqueObjectFinder UniqueObjectFinder {
844# earlier arrays take precedence over later ones
845# add InputArray InputArray OutputArray
846 add InputArray PhotonIsolation/photons photons
847 add InputArray ElectronIsolation/electrons electrons
848 add InputArray MuonIsolation/muons muons
849 add InputArray JetEnergyScale/jets jets
850}
851
852##################
853# ROOT tree writer
854##################
855
856# tracks, towers and eflow objects are not stored by default in the output.
857# if needed (for jet constituent or other studies), uncomment the relevant
858# "add Branch ..." lines.
859
860module TreeWriter TreeWriter {
861# add Branch InputArray BranchName BranchClass
862 add Branch Delphes/allParticles Particle GenParticle
863
864 add Branch TrackMerger/tracks Track Track
865 add Branch Calorimeter/towers Tower Tower
866
867 add Branch HCal/eflowTracks EFlowTrack Track
868 add Branch ECal/eflowPhotons EFlowPhoton Tower
869 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
870
871 add Branch GenJetFinder/jets GenJet Jet
872 add Branch GenMissingET/momentum GenMissingET MissingET
873
874 add Branch UniqueObjectFinder/jets Jet Jet
875 add Branch UniqueObjectFinder/electrons Electron Electron
876 add Branch UniqueObjectFinder/photons Photon Photon
877 add Branch UniqueObjectFinder/muons Muon Muon
878
879 add Branch FatJetFinder/jets FatJet Jet
880
881 add Branch MissingET/momentum MissingET MissingET
882 add Branch ScalarHT/energy ScalarHT ScalarHT
883}