MA5SandBox: delphes_card_cms_sus_16_048.tcl

File delphes_card_cms_sus_16_048.tcl, 19.1 KB (added by Benjamin Fuks, 5 years ago)
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9
10 ChargedHadronMomentumSmearing
11 ElectronEnergySmearing
12 MuonMomentumSmearing
13
14 TrackMerger
15
16 ECal
17 HCal
18 Calorimeter
19
20 EFlowMerger
21 EFlowFilter
22
23 PhotonEfficiency
24
25 ElectronFilter
26 ElectronEfficiency
27
28 ChargedHadronFilter
29
30 MuonEfficiency
31
32 MissingET
33
34 NeutrinoFilter
35 GenJetFinder
36 GenMissingET
37
38 FastJetFinder
39 JetEnergyScale
40 JetFlavorAssociation
41
42 BTagging
43 TauTagging
44
45 TreeWriter
46}
47
48#################################
49# Propagate particles in cylinder
50#################################
51
52module ParticlePropagator ParticlePropagator {
53 set InputArray Delphes/stableParticles
54
55 set OutputArray stableParticles
56 set ChargedHadronOutputArray chargedHadrons
57 set ElectronOutputArray electrons
58 set MuonOutputArray muons
59
60 # radius of the magnetic field coverage, in m
61 set Radius 1.29
62 # half-length of the magnetic field coverage, in m
63 set HalfLength 3.00
64
65 # magnetic field
66 set Bz 3.8
67}
68
69####################################
70# Charged hadron tracking efficiency
71####################################
72
73module Efficiency ChargedHadronTrackingEfficiency {
74 set InputArray ParticlePropagator/chargedHadrons
75 set OutputArray chargedHadrons
76
77 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
78 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
79 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
80 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
81 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
82 (abs(eta) > 2.5) * (0.00)}
83}
84
85########################################
86# Momentum resolution for charged tracks
87########################################
88
89module MomentumSmearing ChargedHadronMomentumSmearing {
90 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
91 set OutputArray chargedHadrons
92
93 # resolution formula for charged hadrons
94 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
95 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
96 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
97}
98
99###################################
100# Momentum resolution for electrons
101###################################
102module EnergySmearing ElectronEnergySmearing {
103 set InputArray ParticlePropagator/electrons
104 set OutputArray electrons
105
106 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
107 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
108 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
109}
110
111
112###############################
113# Momentum resolution for muons
114###############################
115module MomentumSmearing MuonMomentumSmearing {
116 set InputArray ParticlePropagator/muons
117 set OutputArray muons
118
119 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
120 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
122}
123
124##############
125# Track merger
126##############
127
128module Merger TrackMerger {
129 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
130 add InputArray ElectronEnergySmearing/electrons
131 add InputArray MuonMomentumSmearing/muons
132 set OutputArray tracks
133}
134
135
136
137#############
138# ECAL
139#############
140
141module SimpleCalorimeter ECal {
142 set ParticleInputArray ParticlePropagator/stableParticles
143 set TrackInputArray TrackMerger/tracks
144
145 set TowerOutputArray ecalTowers
146 set EFlowTrackOutputArray eflowTracks
147 set EFlowTowerOutputArray eflowPhotons
148
149 set IsEcal true
150
151 set EnergyMin 0.5
152 set EnergySignificanceMin 1.0
153
154 set SmearTowerCenter true
155
156 set pi [expr {acos(-1)}]
157
158 # lists of the edges of each tower in eta and phi
159 # each list starts with the lower edge of the first tower
160 # the list ends with the higher edged of the last tower
161
162 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
163
164 set PhiBins {}
165 for {set i -180} {$i <= 180} {incr i} {
166 add PhiBins [expr {$i * $pi/180.0}]
167 }
168
169 # 0.02 unit in eta up to eta = 1.5 (barrel)
170 for {set i -85} {$i <= 86} {incr i} {
171 set eta [expr {$i * 0.0174}]
172 add EtaPhiBins $eta $PhiBins
173 }
174
175 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
176
177 set PhiBins {}
178 for {set i -180} {$i <= 180} {incr i} {
179 add PhiBins [expr {$i * $pi/180.0}]
180 }
181
182 # 0.02 unit in eta up to eta = 3
183 for {set i 1} {$i <= 84} {incr i} {
184 set eta [expr { -2.958 + $i * 0.0174}]
185 add EtaPhiBins $eta $PhiBins
186 }
187
188 for {set i 1} {$i <= 84} {incr i} {
189 set eta [expr { 1.4964 + $i * 0.0174}]
190 add EtaPhiBins $eta $PhiBins
191 }
192
193 # take present CMS granularity for HF
194
195 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
196 set PhiBins {}
197 for {set i -18} {$i <= 18} {incr i} {
198 add PhiBins [expr {$i * $pi/18.0}]
199 }
200
201 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} {
202 add EtaPhiBins $eta $PhiBins
203 }
204
205 add EnergyFraction {0} {0.0}
206 # energy fractions for e, gamma and pi0
207 add EnergyFraction {11} {1.0}
208 add EnergyFraction {22} {1.0}
209 add EnergyFraction {111} {1.0}
210 # energy fractions for muon, neutrinos and neutralinos
211 add EnergyFraction {12} {0.0}
212 add EnergyFraction {13} {0.0}
213 add EnergyFraction {14} {0.0}
214 add EnergyFraction {16} {0.0}
215 add EnergyFraction {1000022} {0.0}
216 add EnergyFraction {1000023} {0.0}
217 add EnergyFraction {1000025} {0.0}
218 add EnergyFraction {1000035} {0.0}
219 add EnergyFraction {1000045} {0.0}
220 # energy fractions for K0short and Lambda
221 add EnergyFraction {310} {0.3}
222 add EnergyFraction {3122} {0.3}
223
224 # set ResolutionFormula {resolution formula as a function of eta and energy}
225
226 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
227
228 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
229 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
230 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
231 (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) +
232 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
233
234}
235
236
237#############
238# HCAL
239#############
240
241module SimpleCalorimeter HCal {
242 set ParticleInputArray ParticlePropagator/stableParticles
243 set TrackInputArray ECal/eflowTracks
244
245 set TowerOutputArray hcalTowers
246 set EFlowTrackOutputArray eflowTracks
247 set EFlowTowerOutputArray eflowNeutralHadrons
248
249 set IsEcal false
250
251 set EnergyMin 1.0
252 set EnergySignificanceMin 1.0
253
254 set SmearTowerCenter true
255
256 set pi [expr {acos(-1)}]
257
258 # lists of the edges of each tower in eta and phi
259 # each list starts with the lower edge of the first tower
260 # the list ends with the higher edged of the last tower
261
262 # 5 degrees towers
263 set PhiBins {}
264 for {set i -36} {$i <= 36} {incr i} {
265 add PhiBins [expr {$i * $pi/36.0}]
266 }
267 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} {
268 add EtaPhiBins $eta $PhiBins
269 }
270
271 # 10 degrees towers
272 set PhiBins {}
273 for {set i -18} {$i <= 18} {incr i} {
274 add PhiBins [expr {$i * $pi/18.0}]
275 }
276 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} {
277 add EtaPhiBins $eta $PhiBins
278 }
279
280 # 20 degrees towers
281 set PhiBins {}
282 for {set i -9} {$i <= 9} {incr i} {
283 add PhiBins [expr {$i * $pi/9.0}]
284 }
285 foreach eta {-5 -4.7 -4.525 4.7 5} {
286 add EtaPhiBins $eta $PhiBins
287 }
288
289 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
290 add EnergyFraction {0} {1.0}
291 # energy fractions for e, gamma and pi0
292 add EnergyFraction {11} {0.0}
293 add EnergyFraction {22} {0.0}
294 add EnergyFraction {111} {0.0}
295 # energy fractions for muon, neutrinos and neutralinos
296 add EnergyFraction {12} {0.0}
297 add EnergyFraction {13} {0.0}
298 add EnergyFraction {14} {0.0}
299 add EnergyFraction {16} {0.0}
300 add EnergyFraction {1000022} {0.0}
301 add EnergyFraction {1000023} {0.0}
302 add EnergyFraction {1000025} {0.0}
303 add EnergyFraction {1000035} {0.0}
304 add EnergyFraction {1000045} {0.0}
305 # energy fractions for K0short and Lambda
306 add EnergyFraction {310} {0.7}
307 add EnergyFraction {3122} {0.7}
308
309 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
310 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
311 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
312
313}
314
315
316###################################################
317# Tower Merger (in case not using e-flow algorithm)
318###################################################
319
320module Merger Calorimeter {
321 add InputArray ECal/ecalTowers
322 add InputArray HCal/hcalTowers
323 set OutputArray towers
324}
325
326
327####################
328# Energy flow merger
329####################
330module Merger EFlowMerger {
331 add InputArray HCal/eflowTracks
332 add InputArray ECal/eflowPhotons
333 add InputArray HCal/eflowNeutralHadrons
334 set OutputArray eflow
335}
336
337
338######################
339# EFlowFilter
340######################
341module PdgCodeFilter EFlowFilter {
342 set InputArray EFlowMerger/eflow
343 set OutputArray eflow
344 add PdgCode {11}
345 add PdgCode {-11}
346 add PdgCode {13}
347 add PdgCode {-13}
348}
349
350
351#################
352# Electron filter
353#################
354module PdgCodeFilter ElectronFilter {
355 set InputArray HCal/eflowTracks
356 set OutputArray electrons
357 set Invert true
358 add PdgCode {11}
359 add PdgCode {-11}
360}
361
362
363######################
364# ChargedHadronFilter
365######################
366
367module PdgCodeFilter ChargedHadronFilter {
368 set InputArray HCal/eflowTracks
369 set OutputArray chargedHadrons
370
371 add PdgCode {11}
372 add PdgCode {-11}
373 add PdgCode {13}
374 add PdgCode {-13}
375}
376
377
378###################
379# Photon efficiency
380###################
381module Efficiency PhotonEfficiency {
382 set InputArray ECal/eflowPhotons
383 set OutputArray photons
384 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
385 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
386 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
387 (abs(eta) > 2.5) * (0.00)}
388}
389
390
391##############################
392# Electron tracking efficiency
393##############################
394
395module Efficiency ElectronEfficiency {
396 set InputArray ElectronFilter/electrons
397 set OutputArray electrons
398
399 set EfficiencyFormula {
400 (pt <= 5.0) * (0.00) +
401 (abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.336) +
402 (abs(eta) <= 0.8) * (pt > 10.0 && pt <= 15.0) * (0.412) +
403 (abs(eta) <= 0.8) * (pt > 15.0 && pt <= 20.0) * (0.465) +
404 (abs(eta) <= 0.8) * (pt > 20.0 && pt <= 25.0) * (0.496) +
405 (abs(eta) <= 0.8) * (pt > 25.0 ) * (0.503) +
406 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 5.0 && pt <= 10.0) * (0.344) +
407 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 10.0 && pt <= 15.0) * (0.402) +
408 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 15.0 && pt <= 20.0) * (0.448) +
409 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 20.0 && pt <= 25.0) * (0.476) +
410 (abs(eta) > 0.8 && abs(eta) <= 1.442) * (pt > 25.0) * (0.482) +
411 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 5.0 && pt <= 10.0) * (0.233) +
412 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 10.0 && pt <= 15.0) * (0.229) +
413 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 15.0 && pt <= 20.0) * (0.250) +
414 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 20.0 && pt <= 25.0) * (0.261) +
415 (abs(eta) > 1.442 && abs(eta) <= 1.556) * (pt > 25.0) * (0.255) +
416 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.309) +
417 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 15.0) * (0.359) +
418 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 15.0 && pt <= 20.0) * (0.394) +
419 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 25.0) * (0.408) +
420 (abs(eta) > 1.556 && abs(eta) <= 2.0) * (pt > 25.0 ) * (0.418) +
421 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.243) +
422 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 15.0) * (0.287) +
423 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 15.0 && pt <= 20.0) * (0.327) +
424 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 25.0) * (0.341) +
425 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 25.0 ) * (0.352) +
426 (abs(eta) > 2.5) * (0.00)}
427}
428
429##########################
430# Muon tracking efficiency
431##########################
432
433module Efficiency MuonEfficiency {
434 set InputArray MuonMomentumSmearing/muons
435 set OutputArray muons
436
437 set EfficiencyFormula {
438 (pt <= 3.5) * (0.00) +
439 (abs(eta) <= 0.9) * (pt > 3.5 && pt <= 10.0) * (0.647) +
440 (abs(eta) <= 0.9) * (pt > 10.0 && pt <= 15.0) * (0.718) +
441 (abs(eta) <= 0.9) * (pt > 15.0 && pt <= 20.0) * (0.739) +
442 (abs(eta) <= 0.9) * (pt > 20.0 && pt <= 25.0) * (0.760) +
443 (abs(eta) <= 0.9) * (pt > 25.0 ) * (0.763) +
444 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 3.5 && pt <= 10.0) * (0.627) +
445 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 10.0 && pt <= 15.0) * (0.662) +
446 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 15.0 && pt <= 20.0) * (0.694) +
447 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 20.0 && pt <= 25.0) * (0.725) +
448 (abs(eta) > 0.9 && abs(eta) <= 1.2) * (pt > 25.0 ) * (0.733) +
449 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 3.5 && pt <= 10.0) * (0.610) +
450 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 10.0 && pt <= 15.0) * (0.660) +
451 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 15.0 && pt <= 20.0) * (0.678) +
452 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 20.0 && pt <= 25.0) * (0.685) +
453 (abs(eta) > 1.2 && abs(eta) <= 2.1) * (pt > 25.0 ) * (0.723) +
454 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 3.5 && pt <= 10.0) * (0.566) +
455 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 15.0) * (0.629) +
456 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 15.0 && pt <= 20.0) * (0.655) +
457 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 25.0) * (0.670) +
458 (abs(eta) > 2.1 && abs(eta) <= 2.4) * (pt > 25.0 ) * (0.696) +
459 (abs(eta) > 2.4) * (0.00)}
460}
461
462
463
464###################
465# Missing ET merger
466###################
467module Merger MissingET {
468 add InputArray EFlowMerger/eflow
469 set MomentumOutputArray momentum
470}
471
472
473#####################
474# Neutrino Filter
475#####################
476module PdgCodeFilter NeutrinoFilter {
477
478 set InputArray Delphes/stableParticles
479 set OutputArray filteredParticles
480 set PTMin 0.0
481 add PdgCode {12}
482 add PdgCode {14}
483 add PdgCode {16}
484 add PdgCode {-12}
485 add PdgCode {-14}
486 add PdgCode {-16}
487 add PdgCode {1000022}
488}
489
490
491#####################
492# MC truth jet finder
493#####################
494
495module FastJetFinder GenJetFinder {
496 set InputArray NeutrinoFilter/filteredParticles
497 set OutputArray jets
498 set JetAlgorithm 6
499 set ParameterR 0.4
500 set JetPTMin 10.0
501}
502
503#########################
504# Gen Missing ET merger
505########################
506module Merger GenMissingET {
507 add InputArray NeutrinoFilter/filteredParticles
508 set MomentumOutputArray momentum
509}
510
511
512#################
513# Reco jet finder
514#################
515module FastJetFinder FastJetFinder {
516 set InputArray EFlowMerger/eflow
517 set OutputArray jets
518 set JetAlgorithm 6
519 set ParameterR 0.4
520 set JetPTMin 10.0
521}
522
523
524##################
525# Jet Energy Scale
526##################
527
528module EnergyScale JetEnergyScale {
529 set InputArray FastJetFinder/jets
530 set OutputArray jets
531 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
532}
533
534########################
535# Jet Flavor Association
536########################
537
538module JetFlavorAssociation JetFlavorAssociation {
539
540 set PartonInputArray Delphes/partons
541 set ParticleInputArray Delphes/allParticles
542 set ParticleLHEFInputArray Delphes/allParticlesLHEF
543 set JetInputArray JetEnergyScale/jets
544
545 set DeltaR 0.4
546 set PartonPTMin 1.0
547 set PartonEtaMax 2.5
548
549}
550
551
552###########
553# b-tagging
554###########
555
556module BTagging BTagging {
557 set JetInputArray JetEnergyScale/jets
558
559 set BitNumber 0
560
561 add EfficiencyFormula {0} { (pt >= 20.0 && pt < 150.0 ) * (0.2450 - 0.005400*pt + 6.92*10^(-5)*pt^2 - 3.89*10^(-7)*pt^3 + 1.021*10^(-9)*pt^4 - 1.007*10^(-12)*pt^5) +
562 (pt >= 150.0 && pt < 1000.0) * (0.0558 + 0.000428*pt - 1.00*10^(-7)*pt^2)}
563
564 add EfficiencyFormula {4} { (pt >= 20.0 && pt < 300.0 ) * (0.398 - 0.000182*pt + 2.53*10^(-6)*pt^2 - 6.796*10^(-9)*pt^3 + 8.66*10^(-12)*pt^4 - 4.42*10^(-15)*pt^5) +
565 (pt >= 300.0 && pt < 1000.0) * (0.350 + 0.000374*pt - 1.81*10^(-7)*pt^2)}
566
567 add EfficiencyFormula {5} {
568 (pt >= 20.0 && pt < 160.0 ) * (0.4344 + 0.02069*pt - 0.0004429*pt^2 + 5.137*10^(-6)*pt^3 - 3.406*10^(-8)*pt^4 + 1.285*10^(-10)*pt^5 - 2.559*10^(-13)*pt^6 + 2.084*10^(-16)*pt^7) +
569 (pt >=160.0 && pt < 300.0 ) * (0.7140 + 0.002617*pt - 1.656*10^(-5)*pt^2 + 4.767*10^(-8)*pt^3 - 6.431*10^(-11)*pt^4 + 3.287*10^(-14)*pt^5) +
570 (pt >=300.0 && pt < 1000.0) * (0.8720 - 6.885*10^(-5)*pt + 4.34*10^(-8)*pt^2) }
571}
572
573module TauTagging TauTagging {
574 set ParticleInputArray Delphes/allParticles
575 set PartonInputArray Delphes/partons
576 set JetInputArray JetEnergyScale/jets
577
578 set DeltaR 0.4
579 set TauPTMin 1.0
580 set TauEtaMax 2.5
581
582 add EfficiencyFormula {0} {0.01}
583 add EfficiencyFormula {15} {0.6}
584}
585
586
587
588
589##################
590# ROOT tree writer
591##################
592module TreeWriter TreeWriter {
593 add Branch Delphes/allParticles Particle GenParticle
594
595 add Branch TrackMerger/tracks Track Track
596 add Branch Calorimeter/towers Tower Tower
597 add Branch HCal/eflowTracks EFlowTrack Track
598 add Branch ECal/eflowPhotons EFlowPhoton Tower
599 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
600
601 add Branch GenJetFinder/jets GenJet Jet
602 add Branch GenMissingET/momentum GenMissingET MissingET
603
604 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
605 add Branch MuonEfficiency/muons MuonMA5 Muon
606 add Branch PhotonEfficiency/photons PhotonMA5 Photon
607 add Branch JetEnergyScale/jets JetMA5 Jet
608
609 add Branch MissingET/momentum MissingET MissingET
610}