Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ 576428a

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

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