Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ a7bd051

Last change on this file since a7bd051 was e57c062, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

fixed large R jet collections
  • Property mode set to 100644
File size: 66.6 KB
Line 
1#
2# Beta card for HL-LHC and HE-LHC studies
3#
4# Main authors: Michele Selvaggi (CERN)
5#
6# Released on: Dec. 6th, 2017
7#
8#
9#######################################
10# Order of execution of various modules
11#######################################
12
13set ExecutionPath {
14
15 ParticlePropagator
16 TrackMergerProp
17
18 DenseProp
19 DenseMergeTracks
20 DenseTrackFilter
21
22 ChargedHadronTrackingEfficiency
23 ElectronTrackingEfficiency
24 MuonTrackingEfficiency
25
26 ChargedHadronMomentumSmearing
27 ElectronMomentumSmearing
28 MuonMomentumSmearing
29
30 TrackMerger
31
32 ECal
33 HCal
34
35 Calorimeter
36 EFlowMerger
37 EFlowFilter
38
39 PhotonIsolation
40 PhotonEfficiency
41
42 ElectronFilter
43 ElectronIsolation
44 ElectronEfficiency
45
46 ChargedHadronFilter
47
48 MuonIsolation
49 MuonEfficiency
50
51 MissingET
52
53 NeutrinoFilter
54 GenJetFinder
55 GenMissingET
56
57 GenJetFinder02
58 GenJetFinder04
59 GenJetFinder08
60 GenJetFinder15
61
62 FastJetFinder02
63 FastJetFinder04
64 FastJetFinder08
65 FastJetFinder15
66
67 CaloJetFinder02
68 CaloJetFinder04
69 CaloJetFinder08
70 CaloJetFinder15
71
72 TrackJetFinder02
73 TrackJetFinder04
74 TrackJetFinder08
75 TrackJetFinder15
76
77 JetEnergyScale
78
79 JetFlavorAssociation
80
81 BTagging
82 TauTagging
83
84 UniqueObjectFinder
85
86 ScalarHT
87
88 TreeWriter
89}
90
91
92#####################################
93# Track propagation to calorimeters
94#####################################
95
96module ParticlePropagator ParticlePropagator {
97 set InputArray Delphes/stableParticles
98
99 set OutputArray stableParticles
100 set ChargedHadronOutputArray chargedHadrons
101 set ElectronOutputArray electrons
102 set MuonOutputArray muons
103
104 # radius of the magnetic field coverage, in m
105 set Radius 1.2
106 # half-length of the magnetic field coverage, in m
107 set HalfLength 3.25
108
109 # magnetic field
110 set Bz 3.0
111}
112
113##############
114# Track merger
115##############
116
117module Merger TrackMergerProp {
118# add InputArray InputArray
119 add InputArray ParticlePropagator/chargedHadrons
120 add InputArray ParticlePropagator/electrons
121 add InputArray ParticlePropagator/muons
122 set OutputArray tracks
123}
124
125####################################
126# Track propagation to pseudo-pixel
127####################################
128
129module ParticlePropagator DenseProp {
130
131 set InputArray TrackMergerProp/tracks
132
133 # radius of the magnetic field coverage, in m
134 set Radius 0.3
135 set RadiusMax 1.2
136 # half-length of the magnetic field coverage, in m
137 set HalfLength 0.7
138 set HalfLengthMax 3.25
139
140 # magnetic field
141 set Bz 3.0
142}
143
144#####################
145# Dense Track merger
146#####################
147
148module Merger DenseMergeTracks {
149# add InputArray InputArray
150 add InputArray DenseProp/chargedHadrons
151 add InputArray DenseProp/electrons
152 add InputArray DenseProp/muons
153 set OutputArray tracks
154}
155
156
157######################
158# Dense Track Filter
159######################
160
161module DenseTrackFilter DenseTrackFilter {
162
163 set TrackInputArray DenseMergeTracks/tracks
164
165 set TrackOutputArray tracks
166 set ChargedHadronOutputArray chargedHadrons
167 set ElectronOutputArray electrons
168 set MuonOutputArray muons
169
170 set EtaPhiRes 0.003
171 set EtaMax 4.0
172
173 set pi [expr {acos(-1)}]
174
175 set nbins_phi [expr {$pi/$EtaPhiRes} ]
176 set nbins_phi [expr {int($nbins_phi)} ]
177
178 set PhiBins {}
179 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
180 add PhiBins [expr {$i * $pi/$nbins_phi}]
181 }
182
183 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
184 set nbins_eta [expr {int($nbins_eta)} ]
185
186 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
187 set eta [expr {$i * $EtaPhiRes}]
188 add EtaPhiBins $eta $PhiBins
189 }
190}
191
192
193
194####################################
195# Charged hadron tracking efficiency
196####################################
197
198module Efficiency ChargedHadronTrackingEfficiency {
199 set InputArray DenseTrackFilter/chargedHadrons
200 set OutputArray chargedHadrons
201
202 # TBC (which eta_max ? which pT min?)
203
204 # tracking efficiency formula for charged hadrons
205
206 set EfficiencyFormula {
207 (pt <= 0.5) * (0.00) + \
208 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
209 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
210 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
211 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
212 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
213 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
214 (abs(eta) > 4.0) * (0.00)
215 }
216
217}
218
219##############################
220# Electron tracking efficiency
221##############################
222
223module Efficiency ElectronTrackingEfficiency {
224 set InputArray DenseTrackFilter/electrons
225 set OutputArray electrons
226
227# TBC (which eta_max ?)
228# putting same as charged hadrons for now...
229
230 set EfficiencyFormula {
231 (pt <= 0.5) * (0.00) + \
232 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
233 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
234 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
235 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90)/(-0.033*log10(pt) + 1.033) + \
236 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
237 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85)/(-0.033*log10(pt) + 1.033) + \
238 (abs(eta) > 4.0) * (0.00)
239 }
240
241}
242##########################
243# Muon tracking efficiency
244##########################
245
246module Efficiency MuonTrackingEfficiency {
247 set InputArray DenseTrackFilter/muons
248 set OutputArray muons
249
250 set EfficiencyFormula {
251 (pt <= 0.5) * (0.00) + \
252 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 1.00) + \
253 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
254 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.5 && pt <= 1.0) * (pt*1.00) + \
255 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
256 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.95) + \
257 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
258 (abs(eta) > 4.0) * (0.00)
259
260 }
261
262}
263
264########################################
265# Momentum resolution for charged tracks
266########################################
267
268module MomentumSmearing ChargedHadronMomentumSmearing {
269 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
270 set OutputArray chargedHadrons
271
272 # --- CMS resolution Phase II ---
273
274 set ResolutionFormula { 2*((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
275 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
276 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
277 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
278 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
279 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
280 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
281 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
282 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
283 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
284 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
285 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
286 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
287 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
288 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
289 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
290 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
291 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
292 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
293 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
294 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
295 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
296 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
297 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
298 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
299 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
300 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
301 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
302 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
303 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
304 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
305 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
306 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
307 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
308 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
309 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
310 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
311 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
312 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
313 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
314 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
315 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
316 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
317 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
318 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
319 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
320 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
321 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
322 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
323 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
324 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
325 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
326 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
327 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
328 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
329 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
330 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
331 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
332 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
333 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
334 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
335 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
336 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
337 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
338 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
339 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
340 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
341 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
342 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
343 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
344 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
345 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
346 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
347 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
348 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
349 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
350 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
351 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
352 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
353 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
354 }
355
356}
357
358
359###################################
360# Momentum resolution for electrons
361###################################
362
363module MomentumSmearing ElectronMomentumSmearing {
364 set InputArray ElectronTrackingEfficiency/electrons
365 set OutputArray electrons
366
367 set ResolutionFormula { 2*((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
368 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
369 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
370 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
371 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
372 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
373 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
374 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
375 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
376 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
377 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
378 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
379 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
380 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
381 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
382 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
383 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
384 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
385 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
386 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
387 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
388 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
389 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
390 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
391 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
392 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
393 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
394 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
395 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
396 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
397 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
398 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
399 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
400 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
401 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
402 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
403 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
404 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
405 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
406 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
407 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
408 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
409 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
410 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
411 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
412 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
413 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
414 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
415 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
416 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
417 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
418 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
419 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
420 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
421 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
422 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
423 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
424 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
425 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
426 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
427 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
428 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
429 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
430 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
431 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
432 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
433 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
434 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
435 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
436 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
437 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
438 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
439 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
440 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
441 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
442 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
443 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
444 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
445 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
446 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
447 }
448}
449
450
451###############################
452# Momentum resolution for muons
453###############################
454
455module MomentumSmearing MuonMomentumSmearing {
456 set InputArray MuonTrackingEfficiency/muons
457 set OutputArray muons
458
459 # Equal to Tk resolution up to pT = 100 GeV, then tk slope/5.
460
461 set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
462 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
463 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
464 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244/5.*pt/100.000000 + 4./5.*0.009244) + \
465 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
466 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
467 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
468 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172/5.*pt/100.000000 + 4./5.*0.009172) + \
469 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
470 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
471 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
472 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077/5.*pt/100.000000 + 4./5.*0.009077) + \
473 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
474 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
475 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
476 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177/5.*pt/100.000000 + 4./5.*0.009177) + \
477 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
478 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
479 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
480 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429/5.*pt/100.000000 + 4./5.*0.010429) + \
481 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
482 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
483 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
484 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867/5.*pt/100.000000 + 4./5.*0.013867) + \
485 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
486 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
487 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
488 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800/5.*pt/100.000000 + 4./5.*0.016800) + \
489 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
490 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
491 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
492 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161/5.*pt/100.000000 + 4./5.*0.018161) + \
493 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
494 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
495 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
496 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039/5.*pt/100.000000 + 4./5.*0.018039) + \
497 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
498 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
499 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
500 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098/5.*pt/100.000000 + 4./5.*0.020098) + \
501 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
502 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
503 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
504 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289/5.*pt/100.000000 + 4./5.*0.020289) + \
505 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
506 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
507 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
508 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524/5.*pt/100.000000 + 4./5.*0.021524) + \
509 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
510 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
511 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
512 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477/5.*pt/100.000000 + 4./5.*0.024477) + \
513 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
514 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
515 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
516 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392/5.*pt/100.000000 + 4./5.*0.030392) + \
517 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
518 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
519 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
520 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445/5.*pt/100.000000 + 4./5.*0.027445) + \
521 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
522 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
523 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
524 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053/5.*pt/100.000000 + 4./5.*0.025053) + \
525 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
526 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
527 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
528 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710/5.*pt/100.000000 + 4./5.*0.030710) + \
529 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
530 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
531 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
532 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629/5.*pt/100.000000 + 4./5.*0.040629) + \
533 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
534 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
535 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
536 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345/5.*pt/100.000000 + 4./5.*0.052345) + \
537 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
538 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
539 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
540 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617/5.*pt/100.000000 + 4./5.*0.069617))
541 }
542
543}
544
545
546
547##############
548# Track merger
549##############
550
551module Merger TrackMerger {
552# add InputArray InputArray
553 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
554 add InputArray ElectronMomentumSmearing/electrons
555 add InputArray MuonMomentumSmearing/muons
556 set OutputArray tracks
557}
558
559
560#############
561# ECAL
562#############
563
564module SimpleCalorimeter ECal {
565 set ParticleInputArray ParticlePropagator/stableParticles
566 set TrackInputArray TrackMerger/tracks
567
568 set TowerOutputArray ecalTowers
569 set EFlowTrackOutputArray eflowTracks
570 set EFlowTowerOutputArray eflowPhotons
571
572 set IsEcal true
573
574 set EnergyMin 0.5
575 set EnergySignificanceMin 1.0
576
577 set SmearTowerCenter true
578
579 set pi [expr {acos(-1)}]
580
581 # lists of the edges of each tower in eta and phi
582 # each list starts with the lower edge of the first tower
583 # the list ends with the higher edged of the last tower
584
585 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
586
587 set PhiBins {}
588 for {set i -180} {$i <= 180} {incr i} {
589 add PhiBins [expr {$i * $pi/180.0}]
590 }
591
592 # 0.02 unit in eta up to eta = 1.5 (barrel)
593 for {set i -85} {$i <= 86} {incr i} {
594 set eta [expr {$i * 0.0174}]
595 add EtaPhiBins $eta $PhiBins
596 }
597
598 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
599
600 set PhiBins {}
601 for {set i -180} {$i <= 180} {incr i} {
602 add PhiBins [expr {$i * $pi/180.0}]
603 }
604
605 # 0.02 unit in eta up to eta = 3
606 for {set i 1} {$i <= 84} {incr i} {
607 set eta [expr { -2.958 + $i * 0.0174}]
608 add EtaPhiBins $eta $PhiBins
609 }
610
611 for {set i 1} {$i <= 84} {incr i} {
612 set eta [expr { 1.4964 + $i * 0.0174}]
613 add EtaPhiBins $eta $PhiBins
614 }
615
616 # take present CMS granularity for HF
617
618 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
619 set PhiBins {}
620 for {set i -18} {$i <= 18} {incr i} {
621 add PhiBins [expr {$i * $pi/18.0}]
622 }
623
624 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} {
625 add EtaPhiBins $eta $PhiBins
626 }
627
628
629 add EnergyFraction {0} {0.0}
630 # energy fractions for e, gamma and pi0
631 add EnergyFraction {11} {1.0}
632 add EnergyFraction {22} {1.0}
633 add EnergyFraction {111} {1.0}
634 # energy fractions for muon, neutrinos and neutralinos
635 add EnergyFraction {12} {0.0}
636 add EnergyFraction {13} {0.0}
637 add EnergyFraction {14} {0.0}
638 add EnergyFraction {16} {0.0}
639 add EnergyFraction {1000022} {0.0}
640 add EnergyFraction {1000023} {0.0}
641 add EnergyFraction {1000025} {0.0}
642 add EnergyFraction {1000035} {0.0}
643 add EnergyFraction {1000045} {0.0}
644 # energy fractions for K0short and Lambda
645 add EnergyFraction {310} {0.3}
646 add EnergyFraction {3122} {0.3}
647
648
649 # ECAL is CMS
650 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
651 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
652
653 set ResolutionFormula { 2*((abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
654 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
655 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
656 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2)) + \
657 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
658
659
660}
661
662#############
663# HCAL
664#############
665module SimpleCalorimeter HCal {
666 set ParticleInputArray ParticlePropagator/stableParticles
667 set TrackInputArray ECal/eflowTracks
668
669 set TowerOutputArray hcalTowers
670 set EFlowTrackOutputArray eflowTracks
671 set EFlowTowerOutputArray eflowNeutralHadrons
672
673 set IsEcal false
674
675 set EnergyMin 1.0
676 set EnergySignificanceMin 2.0
677
678 set SmearTowerCenter true
679
680 set pi [expr {acos(-1)}]
681
682 # lists of the edges of each tower in eta and phi
683 # each list starts with the lower edge of the first tower
684 # the list ends with the higher edged of the last tower
685
686 # 10 degrees towers
687 set PhiBins {}
688 for {set i -18} {$i <= 18} {incr i} {
689 add PhiBins [expr {$i * $pi/18.0}]
690 }
691 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
692 add EtaPhiBins $eta $PhiBins
693 }
694
695 # 20 degrees towers
696 set PhiBins {}
697 for {set i -9} {$i <= 9} {incr i} {
698 add PhiBins [expr {$i * $pi/9.0}]
699 }
700 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
701 add EtaPhiBins $eta $PhiBins
702 }
703
704 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
705 add EnergyFraction {0} {1.0}
706 # energy fractions for e, gamma and pi0
707 add EnergyFraction {11} {0.0}
708 add EnergyFraction {22} {0.0}
709 add EnergyFraction {111} {0.0}
710 # energy fractions for muon, neutrinos and neutralinos
711 add EnergyFraction {12} {0.0}
712 add EnergyFraction {13} {0.0}
713 add EnergyFraction {14} {0.0}
714 add EnergyFraction {16} {0.0}
715 add EnergyFraction {1000022} {0.0}
716 add EnergyFraction {1000023} {0.0}
717 add EnergyFraction {1000025} {0.0}
718 add EnergyFraction {1000035} {0.0}
719 add EnergyFraction {1000045} {0.0}
720 # energy fractions for K0short and Lambda
721 add EnergyFraction {310} {0.7}
722 add EnergyFraction {3122} {0.7}
723
724 # HCAL is ALTAS
725 # http://arxiv.org/pdf/hep-ex/0004009v1
726 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
727 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
728 set ResolutionFormula { 2*((abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
729 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2)) +
730 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
731}
732
733
734#################
735# Electron filter
736#################
737
738module PdgCodeFilter ElectronFilter {
739 set InputArray HCal/eflowTracks
740 set OutputArray electrons
741 set Invert true
742 add PdgCode {11}
743 add PdgCode {-11}
744}
745
746
747######################
748# ChargedHadronFilter
749######################
750
751module PdgCodeFilter ChargedHadronFilter {
752 set InputArray HCal/eflowTracks
753 set OutputArray chargedHadrons
754
755 add PdgCode {11}
756 add PdgCode {-11}
757 add PdgCode {13}
758 add PdgCode {-13}
759}
760
761
762###################################################
763# Tower Merger (in case not using e-flow algorithm)
764###################################################
765
766module Merger Calorimeter {
767# add InputArray InputArray
768
769 add InputArray ECal/ecalTowers
770 add InputArray HCal/hcalTowers
771 add InputArray MuonMomentumSmearing/muons
772 set OutputArray towers
773
774}
775
776####################
777# Energy flow merger
778####################
779
780module Merger EFlowMerger {
781# add InputArray InputArray
782 add InputArray HCal/eflowTracks
783 add InputArray ECal/eflowPhotons
784 add InputArray HCal/eflowNeutralHadrons
785 set OutputArray eflow
786}
787
788######################
789# EFlowFilter
790######################
791
792module PdgCodeFilter EFlowFilter {
793 set InputArray EFlowMerger/eflow
794 set OutputArray eflow
795
796 add PdgCode {11}
797 add PdgCode {-11}
798 add PdgCode {13}
799 add PdgCode {-13}
800}
801
802
803###################
804# Missing ET merger
805###################
806
807module Merger MissingET {
808# add InputArray InputArray
809 add InputArray EFlowMerger/eflow
810 set MomentumOutputArray momentum
811}
812
813
814
815##################
816# Scalar HT merger
817##################
818
819module Merger ScalarHT {
820# add InputArray InputArray
821 add InputArray EFlowMerger/eflow
822 set EnergyOutputArray energy
823}
824
825#################
826# Neutrino Filter
827#################
828
829module PdgCodeFilter NeutrinoFilter {
830
831 set InputArray Delphes/stableParticles
832 set OutputArray filteredParticles
833
834 set PTMin 0.0
835
836 add PdgCode {12}
837 add PdgCode {14}
838 add PdgCode {16}
839 add PdgCode {-12}
840 add PdgCode {-14}
841 add PdgCode {-16}
842
843}
844
845
846#####################
847# MC truth jet finder
848#####################
849
850module FastJetFinder GenJetFinder {
851 set InputArray NeutrinoFilter/filteredParticles
852
853 set OutputArray jets
854
855 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
856 set JetAlgorithm 6
857 set ParameterR 0.4
858
859 set JetPTMin 15.0
860}
861
862#########################
863# Gen Missing ET merger
864########################
865
866module Merger GenMissingET {
867
868# add InputArray InputArray
869 add InputArray NeutrinoFilter/filteredParticles
870 set MomentumOutputArray momentum
871}
872
873
874
875#####################
876# MC truth jet finder
877#####################
878
879# TBC: is jet radius fine?
880
881module FastJetFinder GenJetFinder02 {
882 set InputArray NeutrinoFilter/filteredParticles
883
884 set OutputArray jets
885
886 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
887 set JetAlgorithm 6
888 set ParameterR 0.2
889
890 set ComputeNsubjettiness 1
891 set Beta 1.0
892 set AxisMode 4
893
894 set ComputeSoftDrop 1
895 set BetaSoftDrop 0.0
896 set SymmetryCutSoftDrop 0.1
897 set R0SoftDrop 0.2
898
899 set JetPTMin 25.0
900}
901
902
903#####################
904# MC truth jet finder
905#####################
906
907# TBC: is jet radius fine?
908
909module FastJetFinder GenJetFinder04 {
910 set InputArray NeutrinoFilter/filteredParticles
911
912 set OutputArray jets
913
914 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
915 set JetAlgorithm 6
916 set ParameterR 0.4
917
918 set ComputeNsubjettiness 1
919 set Beta 1.0
920 set AxisMode 4
921
922 set ComputeSoftDrop 1
923 set BetaSoftDrop 0.0
924 set SymmetryCutSoftDrop 0.1
925 set R0SoftDrop 0.4
926
927 set JetPTMin 25.0
928}
929#####################
930# MC truth jet finder
931#####################
932
933# TBC: is jet radius fine?
934
935module FastJetFinder GenJetFinder08 {
936 set InputArray NeutrinoFilter/filteredParticles
937
938 set OutputArray jets
939
940 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
941 set JetAlgorithm 6
942 set ParameterR 0.8
943
944 set ComputeNsubjettiness 1
945 set Beta 1.0
946 set AxisMode 4
947
948 set ComputeSoftDrop 1
949 set BetaSoftDrop 0.0
950 set SymmetryCutSoftDrop 0.1
951 set R0SoftDrop 0.8
952
953 set JetPTMin 25.0
954}
955
956#####################
957# MC truth jet finder
958#####################
959
960# TBC: is jet radius fine?
961
962module FastJetFinder GenJetFinder15 {
963 set InputArray NeutrinoFilter/filteredParticles
964
965 set OutputArray jets
966
967 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
968 set JetAlgorithm 6
969 set ParameterR 1.5
970
971 set ComputeNsubjettiness 1
972 set Beta 1.0
973 set AxisMode 4
974
975 set ComputeSoftDrop 1
976 set BetaSoftDrop 0.0
977 set SymmetryCutSoftDrop 0.1
978 set R0SoftDrop 1.5
979
980 set JetPTMin 25.0
981}
982
983
984##################
985# Fast Jet finder
986##################
987
988module FastJetFinder FastJetFinder02 {
989 set InputArray EFlowMerger/eflow
990
991 set OutputArray jets
992
993 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
994 set JetAlgorithm 6
995 set ParameterR 0.2
996
997 set ComputeNsubjettiness 1
998 set Beta 1.0
999 set AxisMode 4
1000
1001 set ComputeSoftDrop 1
1002 set BetaSoftDrop 0.0
1003 set SymmetryCutSoftDrop 0.1
1004 set R0SoftDrop 0.2
1005
1006 set JetPTMin 25.0
1007}
1008
1009##################
1010# Fast Jet finder
1011##################
1012
1013module FastJetFinder FastJetFinder04 {
1014 set InputArray EFlowMerger/eflow
1015
1016 set OutputArray jets
1017
1018 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1019 set JetAlgorithm 6
1020 set ParameterR 0.4
1021
1022 set ComputeNsubjettiness 1
1023 set Beta 1.0
1024 set AxisMode 4
1025
1026 set ComputeSoftDrop 1
1027 set BetaSoftDrop 0.0
1028 set SymmetryCutSoftDrop 0.1
1029 set R0SoftDrop 0.4
1030
1031 set JetPTMin 25.0
1032}
1033
1034
1035##################
1036# Fast Jet finder
1037##################
1038
1039module FastJetFinder FastJetFinder08 {
1040 set InputArray EFlowMerger/eflow
1041
1042 set OutputArray jets
1043
1044 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1045 set JetAlgorithm 6
1046 set ParameterR 0.8
1047
1048 set ComputeNsubjettiness 1
1049 set Beta 1.0
1050 set AxisMode 4
1051
1052 set ComputeSoftDrop 1
1053 set BetaSoftDrop 0.0
1054 set SymmetryCutSoftDrop 0.1
1055 set R0SoftDrop 0.8
1056
1057 set JetPTMin 25.0
1058}
1059
1060##################
1061# Fast Jet finder
1062##################
1063
1064module FastJetFinder FastJetFinder15 {
1065 set InputArray EFlowMerger/eflow
1066
1067 set OutputArray jets
1068
1069 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1070 set JetAlgorithm 6
1071 set ParameterR 1.5
1072
1073 set ComputeNsubjettiness 1
1074 set Beta 1.0
1075 set AxisMode 4
1076
1077 set ComputeSoftDrop 1
1078 set BetaSoftDrop 0.0
1079 set SymmetryCutSoftDrop 0.1
1080 set R0SoftDrop 1.5
1081
1082 set JetPTMin 25.0
1083}
1084
1085
1086##################
1087# Fast Jet finder
1088##################
1089
1090module FastJetFinder CaloJetFinder02 {
1091 set InputArray Calorimeter/towers
1092
1093 set OutputArray jets
1094
1095 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1096 set JetAlgorithm 6
1097 set ParameterR 0.2
1098
1099 set ComputeNsubjettiness 1
1100 set Beta 1.0
1101 set AxisMode 4
1102
1103 set ComputeSoftDrop 1
1104 set BetaSoftDrop 0.0
1105 set SymmetryCutSoftDrop 0.1
1106 set R0SoftDrop 0.2
1107
1108 set JetPTMin 25.0
1109}
1110
1111##################
1112# Fast Jet finder
1113##################
1114
1115module FastJetFinder CaloJetFinder04 {
1116 set InputArray Calorimeter/towers
1117
1118 set OutputArray jets
1119
1120 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1121 set JetAlgorithm 6
1122 set ParameterR 0.4
1123
1124 set ComputeNsubjettiness 1
1125 set Beta 1.0
1126 set AxisMode 4
1127
1128 set ComputeSoftDrop 1
1129 set BetaSoftDrop 0.0
1130 set SymmetryCutSoftDrop 0.1
1131 set R0SoftDrop 0.4
1132
1133 set JetPTMin 25.0
1134}
1135
1136
1137##################
1138# Fast Jet finder
1139##################
1140
1141module FastJetFinder CaloJetFinder08 {
1142 set InputArray Calorimeter/towers
1143
1144 set OutputArray jets
1145
1146 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1147 set JetAlgorithm 6
1148 set ParameterR 0.8
1149
1150 set ComputeNsubjettiness 1
1151 set Beta 1.0
1152 set AxisMode 4
1153
1154 set ComputeSoftDrop 1
1155 set BetaSoftDrop 0.0
1156 set SymmetryCutSoftDrop 0.1
1157 set R0SoftDrop 0.8
1158
1159 set JetPTMin 25.0
1160}
1161
1162##################
1163# Fast Jet finder
1164##################
1165
1166module FastJetFinder CaloJetFinder15 {
1167 set InputArray Calorimeter/towers
1168
1169 set OutputArray jets
1170
1171 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1172 set JetAlgorithm 6
1173 set ParameterR 1.5
1174
1175 set ComputeNsubjettiness 1
1176 set Beta 1.0
1177 set AxisMode 4
1178
1179 set ComputeSoftDrop 1
1180 set BetaSoftDrop 0.0
1181 set SymmetryCutSoftDrop 0.1
1182 set R0SoftDrop 1.5
1183
1184 set JetPTMin 25.0
1185}
1186
1187
1188##################
1189# Fast Jet finder
1190##################
1191
1192module FastJetFinder TrackJetFinder02 {
1193 set InputArray TrackMerger/tracks
1194
1195 set OutputArray jets
1196
1197 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1198 set JetAlgorithm 6
1199 set ParameterR 0.2
1200
1201 set ComputeNsubjettiness 1
1202 set Beta 1.0
1203 set AxisMode 4
1204
1205 set ComputeSoftDrop 1
1206 set BetaSoftDrop 0.0
1207 set SymmetryCutSoftDrop 0.1
1208 set R0SoftDrop 0.2
1209
1210 set JetPTMin 25.0
1211}
1212
1213##################
1214# Fast Jet finder
1215##################
1216
1217module FastJetFinder TrackJetFinder04 {
1218 set InputArray TrackMerger/tracks
1219
1220 set OutputArray jets
1221
1222 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1223 set JetAlgorithm 6
1224 set ParameterR 0.4
1225
1226 set ComputeNsubjettiness 1
1227 set Beta 1.0
1228 set AxisMode 4
1229
1230 set ComputeSoftDrop 1
1231 set BetaSoftDrop 0.0
1232 set SymmetryCutSoftDrop 0.1
1233 set R0SoftDrop 0.4
1234
1235 set JetPTMin 25.0
1236}
1237
1238
1239##################
1240# Fast Jet finder
1241##################
1242
1243module FastJetFinder TrackJetFinder08 {
1244 set InputArray TrackMerger/tracks
1245
1246 set OutputArray jets
1247
1248 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1249 set JetAlgorithm 6
1250 set ParameterR 0.8
1251
1252 set ComputeNsubjettiness 1
1253 set Beta 1.0
1254 set AxisMode 4
1255
1256 set ComputeSoftDrop 1
1257 set BetaSoftDrop 0.0
1258 set SymmetryCutSoftDrop 0.1
1259 set R0SoftDrop 0.8
1260
1261 set JetPTMin 25.0
1262}
1263
1264##################
1265# Fast Jet finder
1266##################
1267
1268module FastJetFinder TrackJetFinder15 {
1269 set InputArray TrackMerger/tracks
1270
1271 set OutputArray jets
1272
1273 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1274 set JetAlgorithm 6
1275 set ParameterR 1.5
1276
1277 set ComputeNsubjettiness 1
1278 set Beta 1.0
1279 set AxisMode 4
1280
1281 set ComputeSoftDrop 1
1282 set BetaSoftDrop 0.0
1283 set SymmetryCutSoftDrop 0.1
1284 set R0SoftDrop 1.5
1285
1286 set JetPTMin 25.0
1287}
1288
1289
1290##################
1291# Jet Energy Scale
1292##################
1293
1294module EnergyScale JetEnergyScale {
1295 set InputArray CaloJetFinder04/jets
1296 set OutputArray jets
1297
1298 # scale formula for jets
1299 set ScaleFormula {1.00}
1300}
1301
1302
1303########################
1304# Jet Flavor Association
1305########################
1306
1307module JetFlavorAssociation JetFlavorAssociation {
1308
1309 set PartonInputArray Delphes/partons
1310 set ParticleInputArray Delphes/allParticles
1311 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1312 set JetInputArray JetEnergyScale/jets
1313
1314 set DeltaR 0.5
1315 set PartonPTMin 5.0
1316 set PartonEtaMax 4.0
1317
1318}
1319
1320##################
1321# Photon isolation
1322##################
1323
1324# TBC: check values for iso cuts
1325
1326module Isolation PhotonIsolation {
1327 set CandidateInputArray ECal/eflowPhotons
1328 set IsolationInputArray EFlowFilter/eflow
1329
1330 set OutputArray photons
1331
1332 set DeltaRMax 0.3
1333 set PTMin 0.0
1334 set PTRatioMax 0.1
1335}
1336
1337
1338
1339###################
1340# Photon efficiency
1341###################
1342
1343module Efficiency PhotonEfficiency {
1344 set InputArray PhotonIsolation/photons
1345 set OutputArray photons
1346
1347 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1348
1349 set EfficiencyFormula {
1350 (pt <= 1.0) * (0.00) + \
1351 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
1352 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
1353 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
1354
1355 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
1356 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
1357 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
1358
1359 (abs(eta) > 4.0) * (0.00)}
1360
1361}
1362
1363
1364
1365####################
1366# Electron isolation
1367####################
1368
1369# TBC: check values for iso cuts
1370
1371module Isolation ElectronIsolation {
1372 set CandidateInputArray ElectronFilter/electrons
1373 set IsolationInputArray EFlowFilter/eflow
1374
1375 set OutputArray electrons
1376
1377 set DeltaRMax 0.3
1378
1379 set PTMin 0.0
1380
1381 set PTRatioMax 0.1
1382}
1383
1384
1385#######################
1386# Electron efficiency #
1387#######################
1388
1389module Efficiency ElectronEfficiency {
1390
1391 set InputArray ElectronIsolation/electrons
1392 set OutputArray electrons
1393
1394 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1395 # efficiency formula for electrons (taken from CMS PhaseII card)
1396 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
1397 set EfficiencyFormula {
1398 (pt <= 2.0) * (0.00) +
1399 (abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1400 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1401 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1402 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1403 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1404 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1405 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1406 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1407 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1408 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1409 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1410 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1411 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1412 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1413 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1414 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1415 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1416 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1417 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1418 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1419 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1420 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1421 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1422 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1423 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1424 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1425 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1426 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1427 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1428 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1429 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1430 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1431 (abs(eta) > 4.0) * (0.00)
1432
1433 }
1434
1435}
1436
1437
1438
1439################
1440# Muon isolation
1441################
1442
1443# TBC: check values for iso cuts
1444
1445module Isolation MuonIsolation {
1446 set CandidateInputArray MuonMomentumSmearing/muons
1447 set IsolationInputArray EFlowFilter/eflow
1448
1449 set OutputArray muons
1450
1451 set DeltaRMax 0.3
1452
1453 set PTMin 0.0
1454
1455 set PTRatioMax 0.2
1456}
1457
1458####################
1459# Muon Efficiency #
1460####################
1461
1462module Efficiency MuonEfficiency {
1463 set InputArray MuonIsolation/muons
1464 set OutputArray muons
1465 # tracking + TightID efficiency formula for muons
1466 set EfficiencyFormula {
1467 (pt <= 2.0) * (0.00) +
1468 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1469 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1470 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1471 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1472 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1473 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1474 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1475 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1476 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1477 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1478 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1479 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
1480 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1481 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
1482 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
1483 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1484 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
1485 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
1486 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1487 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
1488 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1489 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1490 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1491 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1492 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
1493 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1494 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
1495 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
1496 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
1497 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1498 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1499 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1500 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1501 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
1502 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
1503 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1504 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1505 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1506 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1507 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
1508 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1509 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1510 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1511 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1512 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1513 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1514 (abs(eta) > 4.0) * (0.00)
1515 }
1516}
1517
1518
1519
1520###########
1521# b-tagging
1522###########
1523
1524module BTagging BTagging {
1525 set JetInputArray JetEnergyScale/jets
1526
1527 set BitNumber 0
1528
1529 add EfficiencyFormula {0} {
1530
1531 (pt <= 10.0) * (0.00) +
1532 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
1533 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
1534 (abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
1535 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
1536 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
1537 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1538 (abs(eta) > 4.0) * (0.00)}
1539
1540 add EfficiencyFormula {4} {
1541
1542 (pt <= 10.0) * (0.00) +
1543 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
1544 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
1545 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1546 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
1547 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
1548 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1549 (abs(eta) > 4.0) * (0.00)}
1550
1551 add EfficiencyFormula {5} {
1552
1553 (pt <= 10.0) * (0.00) +
1554 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
1555 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
1556 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
1557 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
1558 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
1559 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
1560 (abs(eta) >= 4.0) * (0.00)}
1561
1562}
1563
1564
1565#############
1566# tau-tagging
1567#############
1568
1569
1570module TauTagging TauTagging {
1571 set ParticleInputArray Delphes/allParticles
1572 set PartonInputArray Delphes/partons
1573 set JetInputArray JetEnergyScale/jets
1574
1575 set DeltaR 0.5
1576
1577 set TauPTMin 1.0
1578
1579 set TauEtaMax 4.0
1580
1581 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1582 add EfficiencyFormula {0} {
1583
1584 (pt <= 10.0) * (0.00) +
1585 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
1586 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
1587 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1588 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
1589 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
1590 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1591 (abs(eta) > 4.0) * (0.00)}
1592
1593 add EfficiencyFormula {11} {
1594
1595 (pt <= 10.0) * (0.00) +
1596 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
1597 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
1598 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1599 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
1600 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
1601 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1602 (abs(eta) > 4.0) * (0.00)}
1603
1604 add EfficiencyFormula {15} {
1605
1606 (pt <= 10.0) * (0.00) +
1607 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
1608 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
1609 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1610 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
1611 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
1612 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1613 (abs(eta) > 4.0) * (0.00)}
1614
1615}
1616
1617#####################################################
1618# Find uniquely identified photons/electrons/tau/jets
1619#####################################################
1620
1621module UniqueObjectFinder UniqueObjectFinder {
1622# earlier arrays take precedence over later ones
1623# add InputArray InputArray OutputArray
1624 add InputArray PhotonEfficiency/photons photons
1625 add InputArray ElectronEfficiency/electrons electrons
1626 add InputArray MuonEfficiency/muons muons
1627 add InputArray JetEnergyScale/jets jets
1628
1629}
1630
1631
1632##################
1633# ROOT tree writer
1634##################
1635
1636module TreeWriter TreeWriter {
1637# add Branch InputArray BranchName BranchClass
1638 add Branch Delphes/allParticles Particle GenParticle
1639
1640 add Branch GenJetFinder/jets GenJet Jet
1641 add Branch GenMissingET/momentum GenMissingET MissingET
1642
1643 add Branch TrackMerger/tracks Track Track
1644 add Branch Calorimeter/towers Tower Tower
1645
1646 add Branch HCal/eflowTracks EFlowTrack Track
1647 add Branch ECal/eflowPhotons EFlowPhoton Tower
1648 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1649
1650 add Branch UniqueObjectFinder/photons Photon Photon
1651 add Branch UniqueObjectFinder/electrons Electron Electron
1652 add Branch UniqueObjectFinder/muons Muon Muon
1653 add Branch UniqueObjectFinder/jets Jet Jet
1654
1655 add Branch GenJetFinder02/jets GenJet02 Jet
1656 add Branch GenJetFinder04/jets GenJet04 Jet
1657 add Branch GenJetFinder08/jets GenJet08 Jet
1658 add Branch GenJetFinder15/jets GenJet15 Jet
1659
1660 add Branch FastJetFinder02/jets ParticleFlowJet02 Jet
1661 add Branch FastJetFinder04/jets ParticleFlowJet04 Jet
1662 add Branch FastJetFinder08/jets ParticleFlowJet08 Jet
1663 add Branch FastJetFinder15/jets ParticleFlowJet15 Jet
1664
1665 add Branch CaloJetFinder02/jets CaloJet02 Jet
1666 add Branch CaloJetFinder04/jets CaloJet04 Jet
1667 add Branch CaloJetFinder08/jets CaloJet08 Jet
1668 add Branch CaloJetFinder15/jets CaloJet15 Jet
1669
1670 add Branch TrackJetFinder02/jets TrackJet02 Jet
1671 add Branch TrackJetFinder04/jets TrackJet04 Jet
1672 add Branch TrackJetFinder08/jets TrackJet08 Jet
1673 add Branch TrackJetFinder15/jets TrackJet15 Jet
1674
1675 add Branch MissingET/momentum MissingET MissingET
1676 add Branch ScalarHT/energy ScalarHT ScalarHT
1677}
Note: See TracBrowser for help on using the repository browser.