Fork me on GitHub

source: git/cards/CMS_PhaseII_0PU.tcl@ 61c0925

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

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