Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ a0236d5

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

updated muons resolutions in PhaseII cards
  • Property mode set to 100644
File size: 79.8 KB
Line 
1#
2# Phase II - Pile-Up
3#
4# Main authors: Michele Selvaggi (UCL)
5#
6# Released on:
7#
8# Version: v01
9#
10#
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16
17 PileUpMerger
18
19 ParticlePropagator
20
21 ChargedHadronTrackingEfficiency
22 ElectronTrackingEfficiency
23 MuonTrackingEfficiency
24
25 ChargedHadronMomentumSmearing
26 ElectronMomentumSmearing
27 MuonMomentumSmearing
28
29 TrackMerger
30
31 ECal
32 HCal
33
34 ElectronFilter
35 TrackPileUpSubtractor
36
37 TowerMerger
38 NeutralEFlowMerger
39 EFlowMergerAllTracks
40 EFlowMerger
41
42 RunPUPPI
43
44 PhotonEfficiency
45 PhotonIsolation
46
47 ElectronEfficiency
48 ElectronIsolation
49
50 MuonEfficiency
51 MuonIsolation
52
53 NeutrinoFilter
54
55 MissingET
56 GenMissingET
57 GenPileUpMissingET
58
59 GenJetFinder
60 FastJetFinder
61
62 ScalarHT
63
64 JetEnergyScale
65
66 JetFlavorAssociation
67
68 BTaggingLoose
69 BTaggingMedium
70 BTaggingTight
71
72 TauTagging
73
74 GenParticleFilter
75
76 TreeWriter
77}
78
79
80
81###############
82# PileUp Merger
83###############
84
85module PileUpMerger PileUpMerger {
86 set InputArray Delphes/stableParticles
87
88 set ParticleOutputArray stableParticles
89 set VertexOutputArray vertices
90
91 # pre-generated minbias input file
92 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
93
94 # average expected pile up
95 set MeanPileUp 200
96
97 # maximum spread in the beam direction in m
98 set ZVertexSpread 0.15
99
100 # maximum spread in time in s
101 set TVertexSpread 1.5E-09
102
103 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
104
105 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
106
107
108 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
109 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
110 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
111 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
112
113}
114
115
116#################################
117# Propagate particles in cylinder
118#################################
119
120module ParticlePropagator ParticlePropagator {
121 set InputArray PileUpMerger/stableParticles
122
123 set OutputArray stableParticles
124 set ChargedHadronOutputArray chargedHadrons
125 set ElectronOutputArray electrons
126 set MuonOutputArray muons
127
128 # radius of the magnetic field coverage, in m
129 set Radius 1.29
130 # half-length of the magnetic field coverage, in m
131 set HalfLength 3.0
132
133 # magnetic field
134 set Bz 3.8
135}
136
137
138####################################
139# Charged hadron tracking efficiency
140####################################
141
142module Efficiency ChargedHadronTrackingEfficiency {
143 ## particles after propagation
144 set InputArray ParticlePropagator/chargedHadrons
145 set OutputArray chargedHadrons
146 # tracking efficiency formula for charged hadrons
147 set EfficiencyFormula {
148 (pt <= 0.2) * (0.00) + \
149 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
150 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
151 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
152 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
153 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
154 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
155 (abs(eta) > 4.0) * (0.00)
156 }
157}
158
159
160#####################################
161# Electron tracking efficiency - ID
162####################################
163
164module Efficiency ElectronTrackingEfficiency {
165 set InputArray ParticlePropagator/electrons
166 set OutputArray electrons
167 # tracking efficiency formula for electrons
168 set EfficiencyFormula {
169 (pt <= 0.2) * (0.00) + \
170 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
171 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
172 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
173 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
174 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
175 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
176 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
177 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
178 (abs(eta) > 4.0) * (0.00)
179
180 }
181}
182
183##########################
184# Muon tracking efficiency
185##########################
186
187module Efficiency MuonTrackingEfficiency {
188 set InputArray ParticlePropagator/muons
189 set OutputArray muons
190 # tracking efficiency formula for muons
191 set EfficiencyFormula {
192 (pt <= 0.2) * (0.00) + \
193 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
194 (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
195 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
197 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
198 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
199 (abs(eta) > 4.0) * (0.00)
200
201 }
202}
203
204########################################
205# Momentum resolution for charged tracks
206########################################
207
208module MomentumSmearing ChargedHadronMomentumSmearing {
209 ## hadrons after having applied the tracking efficiency
210 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
211 set OutputArray chargedHadrons
212 # resolution formula for charged hadrons
213 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
214 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
215 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
216 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
217 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
218 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
219 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
220 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
221 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
222 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
223 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
224 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
225 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
226 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
227 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
228 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
229 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
230 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
231 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
232 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
233 }
234}
235
236########################################
237# Momentum resolution for electrons
238########################################
239
240module MomentumSmearing ElectronMomentumSmearing {
241 ## hadrons after having applied the tracking efficiency
242 set InputArray ElectronTrackingEfficiency/electrons
243 set OutputArray electrons
244 # resolution formula for electrons
245 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
246 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
247 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
248 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
249 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
250 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
251 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
252 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
253 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
254 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
255 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
256 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
257 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
258 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
259 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
260 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
261 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
262 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
263 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
264 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
265 }
266}
267
268
269###############################
270# Momentum resolution for muons
271###############################
272
273module MomentumSmearing MuonMomentumSmearing {
274 set InputArray MuonTrackingEfficiency/muons
275 set OutputArray muons
276 # 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 efficiency #
808#####################
809
810module Efficiency PhotonEfficiency {
811
812 ## input particles
813 set InputArray ECal/eflowPhotons
814 ## output particles
815 set OutputArray photons
816 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
817 # efficiency formula for photons
818 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
819 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
820 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
821 (abs(eta) > 4.0) * (0.00)}
822
823}
824
825####################
826# Photon isolation #
827####################
828
829module Isolation PhotonIsolation {
830
831 # particle for which calculate the isolation
832 set CandidateInputArray PhotonEfficiency/photons
833
834 # isolation collection
835 set IsolationInputArray RunPUPPI/PuppiParticles
836
837 # output array
838 set OutputArray photons
839
840 # isolation cone
841 set DeltaRMax 0.3
842
843 # minimum pT
844 set PTMin 1.0
845
846 # iso ratio to cut
847 set PTRatioMax 9999.
848}
849
850
851#######################
852# Electron efficiency #
853#######################
854
855module Efficiency ElectronEfficiency {
856
857 set InputArray ElectronFilter/electrons
858 set OutputArray electrons
859
860 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
861 # efficiency formula for electrons
862 set EfficiencyFormula {
863 (pt <= 4.0) * (0.00) + \
864 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
865 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
866 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
867 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
868 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
869 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
870 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
871 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
872 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
873 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
874 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
875 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
876 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
877 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
878 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
879 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
880 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
881 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
882 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
883 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
884 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
885 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
886 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
887 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
888 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
889 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
890 (abs(eta) > 4.0) * (0.00)
891
892 }
893}
894
895######################
896# Electron isolation #
897######################
898
899module Isolation ElectronIsolation {
900
901 set CandidateInputArray ElectronEfficiency/electrons
902
903 # isolation collection
904 set IsolationInputArray RunPUPPI/PuppiParticles
905
906 set OutputArray electrons
907
908 set DeltaRMax 0.3
909 set PTMin 1.0
910 set PTRatioMax 9999.
911
912}
913
914
915
916###################
917# Muon efficiency #
918###################
919
920module Efficiency MuonEfficiency {
921
922 set InputArray MuonMomentumSmearing/muons
923
924 set OutputArray muons
925 # set EfficiencyFormula {efficiency as a function of eta and pt}
926 # efficiency formula for muons
927 set EfficiencyFormula { (pt <= 2.0) * (0.00) + \
928 (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \
929 (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \
930 (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \
931 (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \
932 (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \
933 (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \
934 (abs(eta) > 4.00) * (0.00)
935 }
936
937}
938
939
940
941##################
942# Muon isolation #
943##################
944
945module Isolation MuonIsolation {
946 set CandidateInputArray MuonEfficiency/muons
947
948 # isolation collection
949 set IsolationInputArray RunPUPPI/PuppiParticles
950
951 set OutputArray muons
952
953 set DeltaRMax 0.3
954 set PTMin 1.0
955 set PTRatioMax 9999.
956
957}
958
959
960
961########################
962# Jet Flavor Association
963########################
964
965module JetFlavorAssociation JetFlavorAssociation {
966
967 set PartonInputArray Delphes/partons
968 set ParticleInputArray Delphes/allParticles
969 set ParticleLHEFInputArray Delphes/allParticlesLHEF
970 set JetInputArray JetEnergyScale/jets
971
972 set DeltaR 0.5
973 set PartonPTMin 10.0
974 set PartonEtaMax 4.0
975
976}
977
978
979#############
980# b-tagging #
981#############
982module BTagging BTaggingLoose {
983
984 set JetInputArray JetEnergyScale/jets
985
986 set BitNumber 0
987
988 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
989 (pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
990 (abs(eta) > 3.4) * (0.000)
991 }
992
993
994 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
995 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
996 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
997 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
998 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
999 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
1000 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
1001 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
1002 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
1003 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
1004 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
1005 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
1006 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
1007 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
1008 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
1009 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
1010 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
1011 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
1012 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
1013 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
1014 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
1015 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
1016 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
1017 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
1018 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
1019 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1020
1021 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1022 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
1023 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
1024 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
1025 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
1026 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
1027 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
1028 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
1029 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
1030 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
1031 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
1032 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
1033 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
1034 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
1035 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
1036 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
1037 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
1038 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
1039 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
1040 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
1041 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
1042 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1043
1044 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1045 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
1046 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
1047 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
1048 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
1049 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
1050 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
1051 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
1052 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
1053 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
1054 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
1055 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
1056 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
1057 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
1058 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
1059 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
1060 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
1061 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1062
1063 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1064 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
1065 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
1066 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
1067 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
1068 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
1069 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
1070 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
1071 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
1072 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
1073 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
1074 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
1075 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
1076 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
1077 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
1078 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
1079 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
1080 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1081 (abs(eta) > 3.4) * (0.000)
1082 }
1083
1084 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1085 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
1086 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
1087 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1088 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
1089 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
1090 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
1091 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
1092 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
1093 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
1094 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
1095 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
1096 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
1097 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
1098 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
1099 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
1100 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
1101 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
1102 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
1103 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
1104 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
1105 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
1106 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
1107 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
1108 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
1109 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1110
1111 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1112 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
1113 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
1114 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
1115 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
1116 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
1117 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
1118 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
1119 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
1120 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
1121 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
1122 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
1123 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
1124 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
1125 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
1126 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
1127 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
1128 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
1129 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
1130 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
1131 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
1132 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1133
1134 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1135 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
1136 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
1137 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
1138 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
1139 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
1140 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
1141 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
1142 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
1143 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
1144 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
1145 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
1146 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
1147 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
1148 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
1149 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
1150 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
1151 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1152
1153 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1154 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
1155 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
1156 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
1157 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
1158 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
1159 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
1160 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
1161 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
1162 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
1163 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
1164 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
1165 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
1166 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
1167 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
1168 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
1169 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
1170 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1171 (abs(eta) > 3.4) * (0.000)
1172 }
1173}
1174
1175
1176#############
1177# b-tagging #
1178#############
1179module BTagging BTaggingMedium {
1180
1181 set JetInputArray JetEnergyScale/jets
1182
1183 set BitNumber 1
1184
1185 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1186 (pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
1187 (abs(eta) > 3.4) * (0.000)
1188 }
1189
1190 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1191 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
1192 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
1193 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
1194 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
1195 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
1196 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
1197 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
1198 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
1199 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
1200 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
1201 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
1202 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
1203 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
1204 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
1205 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
1206 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
1207 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
1208 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
1209 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
1210 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
1211 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
1212 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
1213 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
1214 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
1215 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1216
1217 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1218 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
1219 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
1220 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
1221 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
1222 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
1223 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
1224 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
1225 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
1226 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
1227 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
1228 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
1229 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
1230 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
1231 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
1232 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
1233 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
1234 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
1235 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
1236 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
1237 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
1238 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1239
1240 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1241 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
1242 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
1243 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
1244 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
1245 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
1246 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
1247 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
1248 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
1249 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
1250 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
1251 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
1252 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
1253 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
1254 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
1255 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
1256 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
1257 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1258
1259 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1260 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
1261 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
1262 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
1263 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
1264 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
1265 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
1266 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
1267 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
1268 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
1269 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
1270 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
1271 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
1272 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
1273 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
1274 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
1275 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
1276 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1277 (abs(eta) > 3.4) * (0.000)
1278 }
1279
1280
1281 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1282 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
1283 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
1284 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
1285 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
1286 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
1287 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
1288 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
1289 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
1290 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
1291 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
1292 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
1293 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
1294 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
1295 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
1296 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
1297 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
1298 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
1299 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
1300 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
1301 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
1302 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
1303 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
1304 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
1305 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
1306 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1307
1308 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1309 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
1310 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
1311 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
1312 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
1313 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
1314 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
1315 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
1316 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
1317 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
1318 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
1319 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
1320 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
1321 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
1322 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
1323 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
1324 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
1325 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
1326 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
1327 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
1328 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
1329 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1330
1331 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1332 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
1333 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
1334 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
1335 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
1336 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
1337 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
1338 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
1339 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
1340 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
1341 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
1342 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
1343 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
1344 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
1345 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
1346 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
1347 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
1348 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1349
1350 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1351 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
1352 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
1353 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
1354 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
1355 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
1356 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
1357 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
1358 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
1359 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
1360 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
1361 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
1362 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
1363 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
1364 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
1365 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
1366 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
1367 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1368 (abs(eta) > 3.4) * (0.000)
1369 }
1370}
1371
1372
1373
1374#############
1375# b-tagging #
1376#############
1377module BTagging BTaggingTight {
1378
1379 set JetInputArray JetEnergyScale/jets
1380
1381 set BitNumber 2
1382
1383 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1384 (pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
1385 (abs(eta) > 3.4) * (0.000)
1386 }
1387
1388
1389 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1390 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
1391 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
1392 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
1393 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
1394 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
1395 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
1396 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
1397 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
1398 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
1399 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
1400 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
1401 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
1402 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
1403 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
1404 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
1405 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
1406 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
1407 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
1408 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
1409 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
1410 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
1411 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
1412 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
1413 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
1414 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1415
1416 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1417 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
1418 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
1419 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
1420 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
1421 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
1422 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
1423 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
1424 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
1425 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
1426 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
1427 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
1428 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
1429 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
1430 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
1431 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
1432 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
1433 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
1434 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
1435 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
1436 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
1437 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1438
1439 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1440 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
1441 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
1442 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
1443 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
1444 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
1445 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
1446 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
1447 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
1448 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
1449 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
1450 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
1451 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
1452 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
1453 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
1454 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
1455 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
1456 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1457
1458 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1459 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
1460 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
1461 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
1462 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
1463 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
1464 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
1465 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
1466 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
1467 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
1468 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
1469 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
1470 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
1471 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
1472 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
1473 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
1474 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
1475 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1476 (abs(eta) > 3.4) * (0.000)
1477 }
1478
1479
1480
1481 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1482 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
1483 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
1484 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
1485 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
1486 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
1487 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
1488 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
1489 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
1490 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
1491 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
1492 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
1493 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
1494 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
1495 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
1496 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
1497 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
1498 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
1499 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
1500 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
1501 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
1502 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
1503 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
1504 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
1505 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
1506 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1507
1508 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1509 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
1510 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
1511 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
1512 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
1513 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
1514 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
1515 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
1516 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
1517 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
1518 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
1519 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
1520 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
1521 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
1522 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
1523 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
1524 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
1525 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
1526 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
1527 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
1528 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
1529 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1530
1531 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1532 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
1533 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
1534 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
1535 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
1536 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
1537 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
1538 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
1539 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
1540 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
1541 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
1542 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
1543 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
1544 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
1545 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
1546 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
1547 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
1548 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1549
1550 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1551 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
1552 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
1553 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
1554 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
1555 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
1556 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
1557 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
1558 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
1559 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
1560 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
1561 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
1562 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
1563 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
1564 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
1565 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
1566 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
1567 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1568 (abs(eta) > 3.4) * (0.000)
1569 }
1570}
1571
1572
1573
1574#############
1575# tau-tagging
1576#############
1577
1578
1579module TauTagging TauTagging {
1580 set ParticleInputArray Delphes/allParticles
1581 set PartonInputArray Delphes/partons
1582 set JetInputArray JetEnergyScale/jets
1583
1584 set DeltaR 0.5
1585
1586 set TauPTMin 20.0
1587
1588 set TauEtaMax 2.3
1589
1590 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1591
1592 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)) + \
1593 (abs(eta) > 2.3) * (0.000)
1594 }
1595 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) ) + \
1596 (abs(eta) > 2.3) * (0.000)
1597 }
1598}
1599
1600
1601###############################################################################################################
1602# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1603###############################################################################################################
1604
1605module StatusPidFilter GenParticleFilter {
1606
1607 set InputArray Delphes/allParticles
1608 set OutputArray filteredParticles
1609 set PTMin 5.0
1610
1611}
1612
1613
1614##################
1615# ROOT tree writer
1616##################
1617
1618module TreeWriter TreeWriter {
1619# add Branch InputArray BranchName BranchClass
1620 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1621 add Branch PileUpMerger/vertices Vertex Vertex
1622
1623 add Branch GenJetFinder/jets GenJet Jet
1624 add Branch GenMissingET/momentum GenMissingET MissingET
1625
1626# add Branch HCal/eflowTracks EFlowTrack Track
1627# add Branch ECal/eflowPhotons EFlowPhoton Tower
1628# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1629
1630 add Branch PhotonIsolation/photons Photon Photon
1631 add Branch ElectronIsolation/electrons Electron Electron
1632 add Branch MuonIsolation/muons Muon Muon
1633
1634 add Branch JetEnergyScale/jets Jet Jet
1635
1636 add Branch MissingET/momentum MissingET MissingET
1637 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1638 add Branch ScalarHT/energy ScalarHT ScalarHT
1639}
1640
Note: See TracBrowser for help on using the repository browser.