Fork me on GitHub

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