Fork me on GitHub

source: git/cards/CMS_PhaseII_0PU.tcl@ bd8c884

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

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