Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ 1bf8604

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 1bf8604 was 90ce654, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

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