Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ 61f2714

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 61f2714 was 61f2714, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

correct path pile-up file
  • Property mode set to 100644
File size: 75.0 KB
Line 
1#
2# Phase II - Pile-Up
3#
4# Main authors: Michele Selvaggi (UCL)
5#
6# Released on:
7#
8# Version: v01
9#
10#
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16
17 PileUpMerger
18
19 ParticlePropagator
20
21 ChargedHadronTrackingEfficiency
22 ElectronTrackingEfficiency
23 MuonTrackingEfficiency
24
25 ChargedHadronMomentumSmearing
26 ElectronMomentumSmearing
27 MuonMomentumSmearing
28
29 TrackMerger
30
31 ECal
32 HCal
33
34 ElectronFilter
35 TrackPileUpSubtractor
36
37 TowerMerger
38 NeutralEFlowMerger
39 EFlowMergerAllTracks
40 EFlowMerger
41
42 RunPUPPI
43
44 PhotonEfficiency
45 PhotonIsolation
46
47 ElectronEfficiency
48 ElectronIsolation
49
50 MuonEfficiency
51 MuonIsolation
52
53 NeutrinoFilter
54
55 MissingET
56 GenMissingET
57 GenPileUpMissingET
58
59 GenJetFinder
60 FastJetFinder
61
62 ScalarHT
63
64 JetEnergyScale
65
66 JetFlavorAssociation
67
68 BTaggingLoose
69 BTaggingMedium
70 BTaggingTight
71
72 TauTagging
73
74 GenParticleFilter
75
76 TreeWriter
77}
78
79
80
81###############
82# PileUp Merger
83###############
84
85module PileUpMerger PileUpMerger {
86 set InputArray Delphes/stableParticles
87
88 set ParticleOutputArray stableParticles
89 set VertexOutputArray vertices
90
91 # pre-generated minbias input file
92 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
93
94 # average expected pile up
95 set MeanPileUp 200
96
97 # maximum spread in the beam direction in m
98 set ZVertexSpread 0.15
99
100 # maximum spread in time in s
101 set TVertexSpread 1.5E-09
102
103 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
104
105 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
106
107
108 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
109 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
110 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
111 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
112
113}
114
115
116#################################
117# Propagate particles in cylinder
118#################################
119
120module ParticlePropagator ParticlePropagator {
121 set InputArray PileUpMerger/stableParticles
122
123 set OutputArray stableParticles
124 set ChargedHadronOutputArray chargedHadrons
125 set ElectronOutputArray electrons
126 set MuonOutputArray muons
127
128 # radius of the magnetic field coverage, in m
129 set Radius 1.29
130 # half-length of the magnetic field coverage, in m
131 set HalfLength 3.0
132
133 # magnetic field
134 set Bz 3.8
135}
136
137
138####################################
139# Charged hadron tracking efficiency
140####################################
141
142module Efficiency ChargedHadronTrackingEfficiency {
143 ## particles after propagation
144 set InputArray ParticlePropagator/chargedHadrons
145 set OutputArray chargedHadrons
146 # tracking efficiency formula for charged hadrons
147 set EfficiencyFormula {
148 (pt <= 0.2) * (0.00) + \
149 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
150 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
151 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
152 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
153 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
154 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
155 (abs(eta) > 4.0) * (0.00)
156 }
157}
158
159
160#####################################
161# Electron tracking efficiency - ID
162####################################
163
164module Efficiency ElectronTrackingEfficiency {
165 set InputArray ParticlePropagator/electrons
166 set OutputArray electrons
167 # tracking efficiency formula for electrons
168 set EfficiencyFormula {
169 (pt <= 0.2) * (0.00) + \
170 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
171 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
172 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
173 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
174 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
175 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
176 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
177 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
178 (abs(eta) > 4.0) * (0.00)
179
180 }
181}
182
183##########################
184# Muon tracking efficiency
185##########################
186
187module Efficiency MuonTrackingEfficiency {
188 set InputArray ParticlePropagator/muons
189 set OutputArray muons
190 # tracking efficiency formula for muons
191 set EfficiencyFormula {
192 (pt <= 0.2) * (0.00) + \
193 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
194 (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
195 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
197 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
198 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
199 (abs(eta) > 4.0) * (0.00)
200
201 }
202}
203
204########################################
205# Momentum resolution for charged tracks
206########################################
207
208module MomentumSmearing ChargedHadronMomentumSmearing {
209 ## hadrons after having applied the tracking efficiency
210 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
211 set OutputArray chargedHadrons
212 # resolution formula for charged hadrons
213 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
214 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
215 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
216 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
217 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
218 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
219 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
220 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
221 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
222 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
223 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
224 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
225 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
226 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
227 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
228 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
229 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
230 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
231 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
232 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
233 }
234}
235
236########################################
237# Momentum resolution for electrons
238########################################
239
240module MomentumSmearing ElectronMomentumSmearing {
241 ## hadrons after having applied the tracking efficiency
242 set InputArray ElectronTrackingEfficiency/electrons
243 set OutputArray electrons
244 # resolution formula for electrons
245 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
246 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
247 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
248 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
249 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
250 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
251 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
252 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
253 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
254 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
255 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
256 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
257 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
258 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
259 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
260 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
261 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
262 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
263 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
264 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
265 }
266}
267
268
269###############################
270# Momentum resolution for muons
271###############################
272
273module MomentumSmearing MuonMomentumSmearing {
274 set InputArray MuonTrackingEfficiency/muons
275 set OutputArray muons
276 # resolution formula for muons
277 set ResolutionFormula {
278 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.012) + \
279 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
280 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.012) + \
281 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.030) + \
282 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.011) + \
283 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.011) + \
284 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.022) + \
285 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.030) + \
286 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
287 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
288 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
289 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
290 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
291 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
292 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
293 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
294 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
295 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
296 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
297 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
298
299 }
300}
301
302
303##############
304# Track merger
305##############
306
307module Merger TrackMerger {
308# add InputArray InputArray
309 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
310 add InputArray ElectronMomentumSmearing/electrons
311 add InputArray MuonMomentumSmearing/muons
312 set OutputArray tracks
313}
314
315#############
316# ECAL
317#############
318
319module SimpleCalorimeter ECal {
320 set ParticleInputArray ParticlePropagator/stableParticles
321 set TrackInputArray TrackMerger/tracks
322
323 set TowerOutputArray ecalTowers
324 set EFlowTrackOutputArray eflowTracks
325 set EFlowTowerOutputArray eflowPhotons
326
327 set IsEcal true
328
329 set EnergyMin 0.5
330 set EnergySignificanceMin 1.0
331
332 set SmearTowerCenter true
333
334 set pi [expr {acos(-1)}]
335
336 # lists of the edges of each tower in eta and phi
337 # each list starts with the lower edge of the first tower
338 # the list ends with the higher edged of the last tower
339
340 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
341
342 set PhiBins {}
343 for {set i -157} {$i <= 157} {incr i} {
344 add PhiBins [expr {$i * $pi/157.0}]
345 }
346
347 # 0.02 unit in eta up to eta = 1.5 (barrel)
348 for {set i -74} {$i <= 75} {incr i} {
349 set eta [expr {$i * 0.02}]
350 add EtaPhiBins $eta $PhiBins
351 }
352
353 # assume 0.05 x 0.05 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL)
354
355 set PhiBins {}
356 for {set i -65} {$i <= 65} {incr i} {
357 add PhiBins [expr {$i * $pi/65.0}]
358 }
359
360 # 0.05 unit in eta up to eta = 3
361 for {set i 1} {$i <= 30} {incr i} {
362 set eta [expr { -3.00 + $i * 0.05}]
363 add EtaPhiBins $eta $PhiBins
364 }
365
366 for {set i 1} {$i <= 30} {incr i} {
367 set eta [expr { 1.5 + $i * 0.05}]
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 # take present CMS granularity for HF
372
373 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
374 set PhiBins {}
375 for {set i -18} {$i <= 18} {incr i} {
376 add PhiBins [expr {$i * $pi/18.0}]
377 }
378
379 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -3 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
380 add EtaPhiBins $eta $PhiBins
381 }
382
383
384 add EnergyFraction {0} {0.0}
385 # energy fractions for e, gamma and pi0
386 add EnergyFraction {11} {1.0}
387 add EnergyFraction {22} {1.0}
388 add EnergyFraction {111} {1.0}
389 # energy fractions for muon, neutrinos and neutralinos
390 add EnergyFraction {12} {0.0}
391 add EnergyFraction {13} {0.0}
392 add EnergyFraction {14} {0.0}
393 add EnergyFraction {16} {0.0}
394 add EnergyFraction {1000022} {0.0}
395 add EnergyFraction {1000023} {0.0}
396 add EnergyFraction {1000025} {0.0}
397 add EnergyFraction {1000035} {0.0}
398 add EnergyFraction {1000045} {0.0}
399 # energy fractions for K0short and Lambda
400 add EnergyFraction {310} {0.3}
401 add EnergyFraction {3122} {0.3}
402
403 # set ResolutionFormula {resolution formula as a function of eta and energy}
404
405 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
406 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
407
408 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.003^2 + energy*0.028^2 + 0.12^2) + \
409 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
410 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
411 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
412 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
413
414}
415
416#############
417# HCAL
418#############
419
420module SimpleCalorimeter HCal {
421 set ParticleInputArray ParticlePropagator/stableParticles
422 set TrackInputArray ECal/eflowTracks
423
424 set TowerOutputArray hcalTowers
425 set EFlowTrackOutputArray eflowTracks
426 set EFlowTowerOutputArray eflowNeutralHadrons
427
428 set IsEcal false
429
430 set EnergyMin 1.0
431 set EnergySignificanceMin 1.0
432
433 set SmearTowerCenter true
434
435 set pi [expr {acos(-1)}]
436
437 # lists of the edges of each tower in eta and phi
438 # each list starts with the lower edge of the first tower
439 # the list ends with the higher edged of the last tower
440
441 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
442
443 set PhiBins {}
444 for {set i -36} {$i <= 36} {incr i} {
445 add PhiBins [expr {$i * $pi/36.0}]
446 }
447 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.65} {
448 add EtaPhiBins $eta $PhiBins
449 }
450
451 # assume 0.05 x 0.05 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL)
452
453 set PhiBins {}
454 for {set i -65} {$i <= 65} {incr i} {
455 add PhiBins [expr {$i * $pi/65.0}]
456 }
457
458 # 0.05 unit in eta up to eta = 3
459 for {set i 1} {$i <= 27} {incr i} {
460 set eta [expr { -3.00 + $i * 0.05}]
461 add EtaPhiBins $eta $PhiBins
462 }
463
464 for {set i 4} {$i <= 30} {incr i} {
465 set eta [expr { 1.5 + $i * 0.05}]
466 add EtaPhiBins $eta $PhiBins
467 }
468
469 # take present CMS granularity for HF
470
471 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
472 set PhiBins {}
473 for {set i -18} {$i <= 18} {incr i} {
474 add PhiBins [expr {$i * $pi/18.0}]
475 }
476
477 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -3 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
478 add EtaPhiBins $eta $PhiBins
479 }
480
481 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
482 add EnergyFraction {0} {1.0}
483 # energy fractions for e, gamma and pi0
484 add EnergyFraction {11} {0.0}
485 add EnergyFraction {22} {0.0}
486 add EnergyFraction {111} {0.0}
487 # energy fractions for muon, neutrinos and neutralinos
488 add EnergyFraction {12} {0.0}
489 add EnergyFraction {13} {0.0}
490 add EnergyFraction {14} {0.0}
491 add EnergyFraction {16} {0.0}
492 add EnergyFraction {1000022} {0.0}
493 add EnergyFraction {1000023} {0.0}
494 add EnergyFraction {1000025} {0.0}
495 add EnergyFraction {1000035} {0.0}
496 add EnergyFraction {1000045} {0.0}
497 # energy fractions for K0short and Lambda
498 add EnergyFraction {310} {0.7}
499 add EnergyFraction {3122} {0.7}
500
501# set ResolutionFormula {resolution formula as a function of eta and energy}
502 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.52^2 + 1.59^2) + \
503 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.05^2 + energy*0.706^2) + \
504 (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
505
506}
507
508#################
509# Electron filter
510#################
511
512module PdgCodeFilter ElectronFilter {
513 set InputArray HCal/eflowTracks
514 set OutputArray electrons
515 set Invert true
516 add PdgCode {11}
517 add PdgCode {-11}
518}
519
520##########################
521# Track pile-up subtractor
522##########################
523
524module TrackPileUpSubtractor TrackPileUpSubtractor {
525# add InputArray InputArray OutputArray
526 add InputArray HCal/eflowTracks eflowTracks
527 add InputArray ElectronFilter/electrons electrons
528 add InputArray MuonMomentumSmearing/muons muons
529
530 set VertexInputArray PileUpMerger/vertices
531 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
532 # Z vertex resolution in m
533 set ZVertexResolution 0.0001
534}
535
536
537###################################################
538# Tower Merger (in case not using e-flow algorithm)
539###################################################
540
541module Merger TowerMerger {
542# add InputArray InputArray
543 add InputArray ECal/ecalTowers
544 add InputArray HCal/hcalTowers
545 set OutputArray towers
546}
547
548
549####################
550# Neutral eflow erger
551####################
552
553module Merger NeutralEFlowMerger {
554# add InputArray InputArray
555 add InputArray ECal/eflowPhotons
556 add InputArray HCal/eflowNeutralHadrons
557 set OutputArray eflowTowers
558}
559
560
561####################
562# Energy flow merger
563####################
564
565module Merger EFlowMerger {
566# add InputArray InputArray
567 add InputArray HCal/eflowTracks
568 add InputArray ECal/eflowPhotons
569 add InputArray HCal/eflowNeutralHadrons
570 set OutputArray eflow
571}
572
573##################################
574# Energy flow merger (all tracks)
575##################################
576
577module Merger EFlowMergerAllTracks {
578# add InputArray InputArray
579 add InputArray TrackMerger/tracks
580 add InputArray ECal/eflowPhotons
581 add InputArray HCal/eflowNeutralHadrons
582 set OutputArray eflow
583}
584
585#########################################
586### Run the puppi code (to be tuned) ###
587#########################################
588
589module RunPUPPI RunPUPPI {
590 ## input information
591 set TrackInputArray TrackMerger/tracks
592 set NeutralInputArray NeutralEFlowMerger/eflowTowers
593 set PVInputArray PileUpMerger/vertices
594 set MinPuppiWeight 0.05
595 set UseExp false
596
597 ## define puppi algorithm parameters (more than one for the same eta region is possible)
598 add EtaMinBin 0. 2.5 2.5 3.0 3.0
599 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0
600 add PtMinBin 0. 0.5 0.5 0.5 0.5
601 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25
602 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5
603 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0
604 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5
605 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02
606 add ApplyCHS true true true true true
607 add UseCharged true false false false false
608 add ApplyLowPUCorr true true true true true
609 add MetricId 5 5 1 5 1
610
611 ## output name
612 set OutputArray PuppiParticles
613 set OutputArrayTracks puppiTracks
614 set OutputArrayNeutrals puppiNeutrals
615}
616
617
618
619###################
620# Missing ET merger
621###################
622
623module Merger MissingET {
624# add InputArray InputArray
625# add InputArray RunPUPPI/PuppiParticles
626 add InputArray EFlowMerger/eflow
627 set MomentumOutputArray momentum
628}
629
630###################
631# Ger PileUp Missing ET
632###################
633
634module Merger GenPileUpMissingET {
635# add InputArray InputArray
636# add InputArray RunPUPPI/PuppiParticles
637 add InputArray ParticlePropagator/stableParticles
638 set MomentumOutputArray momentum
639}
640
641##################
642# Scalar HT merger
643##################
644
645module Merger ScalarHT {
646# add InputArray InputArray
647 add InputArray RunPUPPI/PuppiParticles
648 set EnergyOutputArray energy
649}
650
651#################
652# Neutrino Filter
653#################
654
655module PdgCodeFilter NeutrinoFilter {
656
657 set InputArray Delphes/stableParticles
658 set OutputArray filteredParticles
659
660 set PTMin 0.0
661
662 add PdgCode {12}
663 add PdgCode {14}
664 add PdgCode {16}
665 add PdgCode {-12}
666 add PdgCode {-14}
667 add PdgCode {-16}
668
669}
670
671
672#####################
673# MC truth jet finder
674#####################
675
676module FastJetFinder GenJetFinder {
677 set InputArray NeutrinoFilter/filteredParticles
678
679 set OutputArray jets
680
681 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
682 set JetAlgorithm 6
683 set ParameterR 0.4
684
685 set JetPTMin 15.0
686}
687
688#########################
689# Gen Missing ET merger
690########################
691
692module Merger GenMissingET {
693
694# add InputArray InputArray
695 add InputArray NeutrinoFilter/filteredParticles
696 set MomentumOutputArray momentum
697}
698
699
700
701############
702# Jet finder
703############
704
705module FastJetFinder FastJetFinder {
706# set InputArray TowerMerger/towers
707 set InputArray RunPUPPI/PuppiParticles
708
709 set OutputArray jets
710
711 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
712 set JetAlgorithm 6
713 set ParameterR 0.4
714
715 set JetPTMin 15.0
716}
717
718##################
719# Jet Energy Scale
720##################
721
722module EnergyScale JetEnergyScale {
723 set InputArray FastJetFinder/jets
724 set OutputArray jets
725
726 # scale formula for jets
727 set ScaleFormula {1.00}
728}
729
730
731#####################
732# Photon efficiency #
733#####################
734
735module Efficiency PhotonEfficiency {
736
737 ## input particles
738 set InputArray ECal/eflowPhotons
739 ## output particles
740 set OutputArray photons
741 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
742 # efficiency formula for photons
743 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
744 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
745 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
746 (abs(eta) > 4.0) * (0.00)}
747
748}
749
750####################
751# Photon isolation #
752####################
753
754module Isolation PhotonIsolation {
755
756 # particle for which calculate the isolation
757 set CandidateInputArray PhotonEfficiency/photons
758
759 # isolation collection
760 set IsolationInputArray RunPUPPI/PuppiParticles
761
762 # output array
763 set OutputArray photons
764
765 # isolation cone
766 set DeltaRMax 0.3
767
768 # minimum pT
769 set PTMin 1.0
770
771 # iso ratio to cut
772 set PTRatioMax 9999.
773}
774
775
776#######################
777# Electron efficiency #
778#######################
779
780module Efficiency ElectronEfficiency {
781
782 set InputArray ElectronFilter/electrons
783 set OutputArray electrons
784
785 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
786 # efficiency formula for electrons
787 set EfficiencyFormula {
788 (pt <= 4.0) * (0.00) + \
789 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
790 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
791 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
792 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
793 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
794 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
795 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
796 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
797 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
798 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
799 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
800 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
801 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
802 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
803 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
804 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
805 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
806 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
807 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
808 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
809 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
810 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
811 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
812 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
813 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
814 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
815 (abs(eta) > 4.0) * (0.00)
816
817 }
818}
819
820######################
821# Electron isolation #
822######################
823
824module Isolation ElectronIsolation {
825
826 set CandidateInputArray ElectronEfficiency/electrons
827
828 # isolation collection
829 set IsolationInputArray RunPUPPI/PuppiParticles
830
831 set OutputArray electrons
832
833 set DeltaRMax 0.3
834 set PTMin 1.0
835 set PTRatioMax 9999.
836
837}
838
839
840
841###################
842# Muon efficiency #
843###################
844
845module Efficiency MuonEfficiency {
846
847 set InputArray MuonMomentumSmearing/muons
848
849 set OutputArray muons
850 # set EfficiencyFormula {efficiency as a function of eta and pt}
851 # efficiency formula for muons
852 set EfficiencyFormula { (pt <= 2.0) * (0.00) + \
853 (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \
854 (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \
855 (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \
856 (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \
857 (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \
858 (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \
859 (abs(eta) > 4.00) * (0.00)
860 }
861
862}
863
864
865
866##################
867# Muon isolation #
868##################
869
870module Isolation MuonIsolation {
871 set CandidateInputArray MuonEfficiency/muons
872
873 # isolation collection
874 set IsolationInputArray RunPUPPI/PuppiParticles
875
876 set OutputArray muons
877
878 set DeltaRMax 0.3
879 set PTMin 1.0
880 set PTRatioMax 9999.
881
882}
883
884
885
886########################
887# Jet Flavor Association
888########################
889
890module JetFlavorAssociation JetFlavorAssociation {
891
892 set PartonInputArray Delphes/partons
893 set ParticleInputArray Delphes/allParticles
894 set ParticleLHEFInputArray Delphes/allParticlesLHEF
895 set JetInputArray JetEnergyScale/jets
896
897 set DeltaR 0.5
898 set PartonPTMin 10.0
899 set PartonEtaMax 4.0
900
901}
902
903
904#############
905# b-tagging #
906#############
907module BTagging BTaggingLoose {
908
909 set JetInputArray JetEnergyScale/jets
910
911 set BitNumber 0
912
913 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
914 (pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
915 (abs(eta) > 3.4) * (0.000)
916 }
917
918
919 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
920 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
921 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
922 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
923 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
924 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
925 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
926 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
927 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
928 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
929 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
930 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
931 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
932 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
933 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
934 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
935 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
936 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
937 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
938 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
939 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
940 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
941 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
942 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
943 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
944 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
945
946 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
947 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
948 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
949 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
950 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
951 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
952 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
953 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
954 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
955 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
956 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
957 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
958 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
959 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
960 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
961 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
962 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
963 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
964 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
965 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
966 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
967 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
968
969 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
970 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
971 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
972 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
973 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
974 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
975 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
976 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
977 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
978 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
979 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
980 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
981 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
982 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
983 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
984 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
985 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
986 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
987
988 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
989 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
990 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
991 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
992 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
993 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
994 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
995 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
996 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
997 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
998 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
999 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
1000 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
1001 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
1002 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
1003 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
1004 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
1005 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1006 (abs(eta) > 3.4) * (0.000)
1007 }
1008
1009 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1010 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
1011 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
1012 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1013 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
1014 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
1015 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
1016 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
1017 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
1018 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
1019 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
1020 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
1021 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
1022 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
1023 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
1024 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
1025 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
1026 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
1027 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
1028 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
1029 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
1030 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
1031 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
1032 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
1033 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
1034 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1035
1036 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1037 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
1038 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
1039 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
1040 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
1041 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
1042 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
1043 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
1044 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
1045 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
1046 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
1047 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
1048 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
1049 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
1050 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
1051 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
1052 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
1053 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
1054 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
1055 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
1056 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
1057 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1058
1059 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1060 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
1061 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
1062 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
1063 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
1064 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
1065 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
1066 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
1067 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
1068 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
1069 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
1070 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
1071 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
1072 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
1073 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
1074 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
1075 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
1076 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1077
1078 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1079 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
1080 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
1081 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
1082 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
1083 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
1084 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
1085 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
1086 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
1087 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
1088 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
1089 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
1090 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
1091 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
1092 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
1093 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
1094 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
1095 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1096 (abs(eta) > 3.4) * (0.000)
1097 }
1098}
1099
1100
1101#############
1102# b-tagging #
1103#############
1104module BTagging BTaggingMedium {
1105
1106 set JetInputArray JetEnergyScale/jets
1107
1108 set BitNumber 1
1109
1110 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1111 (pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
1112 (abs(eta) > 3.4) * (0.000)
1113 }
1114
1115 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1116 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
1117 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
1118 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
1119 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
1120 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
1121 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
1122 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
1123 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
1124 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
1125 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
1126 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
1127 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
1128 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
1129 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
1130 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
1131 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
1132 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
1133 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
1134 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
1135 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
1136 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
1137 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
1138 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
1139 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
1140 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1141
1142 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1143 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
1144 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
1145 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
1146 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
1147 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
1148 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
1149 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
1150 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
1151 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
1152 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
1153 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
1154 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
1155 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
1156 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
1157 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
1158 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
1159 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
1160 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
1161 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
1162 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
1163 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1164
1165 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1166 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
1167 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
1168 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
1169 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
1170 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
1171 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
1172 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
1173 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
1174 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
1175 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
1176 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
1177 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
1178 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
1179 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
1180 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
1181 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
1182 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1183
1184 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1185 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
1186 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
1187 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
1188 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
1189 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
1190 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
1191 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
1192 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
1193 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
1194 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
1195 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
1196 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
1197 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
1198 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
1199 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
1200 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
1201 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1202 (abs(eta) > 3.4) * (0.000)
1203 }
1204
1205
1206 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1207 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
1208 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
1209 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
1210 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
1211 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
1212 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
1213 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
1214 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
1215 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
1216 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
1217 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
1218 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
1219 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
1220 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
1221 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
1222 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
1223 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
1224 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
1225 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
1226 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
1227 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
1228 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
1229 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
1230 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
1231 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1232
1233 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1234 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
1235 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
1236 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
1237 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
1238 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
1239 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
1240 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
1241 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
1242 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
1243 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
1244 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
1245 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
1246 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
1247 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
1248 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
1249 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
1250 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
1251 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
1252 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
1253 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
1254 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1255
1256 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1257 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
1258 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
1259 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
1260 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
1261 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
1262 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
1263 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
1264 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
1265 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
1266 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
1267 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
1268 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
1269 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
1270 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
1271 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
1272 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
1273 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1274
1275 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1276 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
1277 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
1278 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
1279 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
1280 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
1281 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
1282 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
1283 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
1284 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
1285 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
1286 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
1287 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
1288 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
1289 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
1290 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
1291 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
1292 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1293 (abs(eta) > 3.4) * (0.000)
1294 }
1295}
1296
1297
1298
1299#############
1300# b-tagging #
1301#############
1302module BTagging BTaggingTight {
1303
1304 set JetInputArray JetEnergyScale/jets
1305
1306 set BitNumber 2
1307
1308 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1309 (pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
1310 (abs(eta) > 3.4) * (0.000)
1311 }
1312
1313
1314 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1315 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
1316 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
1317 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
1318 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
1319 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
1320 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
1321 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
1322 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
1323 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
1324 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
1325 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
1326 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
1327 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
1328 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
1329 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
1330 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
1331 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
1332 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
1333 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
1334 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
1335 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
1336 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
1337 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
1338 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
1339 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1340
1341 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1342 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
1343 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
1344 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
1345 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
1346 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
1347 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
1348 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
1349 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
1350 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
1351 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
1352 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
1353 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
1354 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
1355 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
1356 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
1357 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
1358 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
1359 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
1360 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
1361 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
1362 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1363
1364 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1365 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
1366 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
1367 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
1368 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
1369 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
1370 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
1371 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
1372 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
1373 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
1374 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
1375 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
1376 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
1377 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
1378 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
1379 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
1380 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
1381 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1382
1383 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1384 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
1385 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
1386 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
1387 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
1388 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
1389 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
1390 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
1391 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
1392 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
1393 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
1394 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
1395 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
1396 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
1397 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
1398 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
1399 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
1400 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1401 (abs(eta) > 3.4) * (0.000)
1402 }
1403
1404
1405
1406 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1407 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
1408 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
1409 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
1410 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
1411 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
1412 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
1413 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
1414 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
1415 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
1416 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
1417 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
1418 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
1419 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
1420 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
1421 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
1422 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
1423 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
1424 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
1425 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
1426 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
1427 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
1428 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
1429 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
1430 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
1431 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1432
1433 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1434 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
1435 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
1436 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
1437 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
1438 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
1439 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
1440 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
1441 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
1442 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
1443 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
1444 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
1445 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
1446 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
1447 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
1448 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
1449 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
1450 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
1451 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
1452 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
1453 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
1454 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1455
1456 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1457 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
1458 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
1459 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
1460 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
1461 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
1462 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
1463 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
1464 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
1465 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
1466 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
1467 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
1468 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
1469 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
1470 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
1471 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
1472 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
1473 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1474
1475 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1476 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
1477 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
1478 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
1479 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
1480 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
1481 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
1482 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
1483 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
1484 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
1485 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
1486 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
1487 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
1488 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
1489 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
1490 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
1491 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
1492 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1493 (abs(eta) > 3.4) * (0.000)
1494 }
1495}
1496
1497
1498
1499#############
1500# tau-tagging
1501#############
1502
1503
1504module TauTagging TauTagging {
1505 set ParticleInputArray Delphes/allParticles
1506 set PartonInputArray Delphes/partons
1507 set JetInputArray JetEnergyScale/jets
1508
1509 set DeltaR 0.5
1510
1511 set TauPTMin 20.0
1512
1513 set TauEtaMax 2.3
1514
1515 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1516
1517 add EfficiencyFormula {0} { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097*pt-2.19642e-5*pt^2+1.49393e-7*pt^3-4.58972e-10*pt^4+5.27983e-13*pt^5 )) * (pt<250) + 0.0032*(pt>250)) + \
1518 (abs(eta) > 2.3) * (0.000)
1519 }
1520 add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
1521 (abs(eta) > 2.3) * (0.000)
1522 }
1523}
1524
1525
1526###############################################################################################################
1527# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1528###############################################################################################################
1529
1530module StatusPidFilter GenParticleFilter {
1531
1532 set InputArray Delphes/allParticles
1533 set OutputArray filteredParticles
1534 set PTMin 5.0
1535
1536}
1537
1538
1539##################
1540# ROOT tree writer
1541##################
1542
1543module TreeWriter TreeWriter {
1544# add Branch InputArray BranchName BranchClass
1545 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1546 add Branch PileUpMerger/vertices Vertex Vertex
1547
1548 add Branch GenJetFinder/jets GenJet Jet
1549 add Branch GenMissingET/momentum GenMissingET MissingET
1550
1551# add Branch HCal/eflowTracks EFlowTrack Track
1552# add Branch ECal/eflowPhotons EFlowPhoton Tower
1553# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1554
1555 add Branch PhotonIsolation/photons Photon Photon
1556 add Branch ElectronIsolation/electrons Electron Electron
1557 add Branch MuonIsolation/muons Muon Muon
1558
1559 add Branch JetEnergyScale/jets Jet Jet
1560
1561 add Branch MissingET/momentum MissingET MissingET
1562 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1563 add Branch ScalarHT/energy ScalarHT ScalarHT
1564}
1565
Note: See TracBrowser for help on using the repository browser.