Fork me on GitHub

source: git/cards/CMS_PhaseII_0PU.tcl@ 730bc64

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

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