Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ 667a02a

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

Merge branch 'master' of github.com:delphes/delphes
  • 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
765 add InputArray ECal/ecalTowers
766 add InputArray HCal/hcalTowers
767 add InputArray MuonMomentumSmearing/muons
768 set OutputArray towers
769
770}
771
772####################
773# Energy flow merger
774####################
775
776module Merger EFlowMerger {
777# add InputArray InputArray
778 add InputArray HCal/eflowTracks
779 add InputArray ECal/eflowPhotons
780 add InputArray HCal/eflowNeutralHadrons
781 set OutputArray eflow
782}
783
784######################
785# EFlowFilter
786######################
787
788module PdgCodeFilter EFlowFilter {
789 set InputArray EFlowMerger/eflow
790 set OutputArray eflow
791
792 add PdgCode {11}
793 add PdgCode {-11}
794 add PdgCode {13}
795 add PdgCode {-13}
796}
797
798
799###################
800# Missing ET merger
801###################
802
803module Merger MissingET {
804# add InputArray InputArray
805 add InputArray EFlowMerger/eflow
806 set MomentumOutputArray momentum
807}
808
809
810
811##################
812# Scalar HT merger
813##################
814
815module Merger ScalarHT {
816# add InputArray InputArray
817 add InputArray EFlowMerger/eflow
818 set EnergyOutputArray energy
819}
820
821#################
822# Neutrino Filter
823#################
824
825module PdgCodeFilter NeutrinoFilter {
826
827 set InputArray Delphes/stableParticles
828 set OutputArray filteredParticles
829
830 set PTMin 0.0
831
832 add PdgCode {12}
833 add PdgCode {14}
834 add PdgCode {16}
835 add PdgCode {-12}
836 add PdgCode {-14}
837 add PdgCode {-16}
838
839}
840
841
842#####################
843# MC truth jet finder
844#####################
845
846module FastJetFinder GenJetFinder {
847 set InputArray NeutrinoFilter/filteredParticles
848
849 set OutputArray jets
850
851 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
852 set JetAlgorithm 6
853 set ParameterR 0.4
854
855 set JetPTMin 15.0
856}
857
858#########################
859# Gen Missing ET merger
860########################
861
862module Merger GenMissingET {
863
864# add InputArray InputArray
865 add InputArray NeutrinoFilter/filteredParticles
866 set MomentumOutputArray momentum
867}
868
869
870
871#####################
872# MC truth jet finder
873#####################
874
875# TBC: is jet radius fine?
876
877module FastJetFinder GenJetFinder02 {
878 set InputArray NeutrinoFilter/filteredParticles
879
880 set OutputArray jets
881
882 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
883 set JetAlgorithm 6
884 set ParameterR 0.2
885
886 set ComputeNsubjettiness 1
887 set Beta 1.0
888 set AxisMode 4
889
890 set ComputeSoftDrop 1
891 set BetaSoftDrop 0.0
892 set SymmetryCutSoftDrop 0.1
893 set R0SoftDrop 0.2
894
895 set JetPTMin 25.0
896}
897
898
899#####################
900# MC truth jet finder
901#####################
902
903# TBC: is jet radius fine?
904
905module FastJetFinder GenJetFinder04 {
906 set InputArray NeutrinoFilter/filteredParticles
907
908 set OutputArray jets
909
910 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
911 set JetAlgorithm 6
912 set ParameterR 0.4
913
914 set ComputeNsubjettiness 1
915 set Beta 1.0
916 set AxisMode 4
917
918 set ComputeSoftDrop 1
919 set BetaSoftDrop 0.0
920 set SymmetryCutSoftDrop 0.1
921 set R0SoftDrop 0.4
922
923 set JetPTMin 25.0
924}
925#####################
926# MC truth jet finder
927#####################
928
929# TBC: is jet radius fine?
930
931module FastJetFinder GenJetFinder08 {
932 set InputArray NeutrinoFilter/filteredParticles
933
934 set OutputArray jets
935
936 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
937 set JetAlgorithm 6
938 set ParameterR 0.8
939
940 set ComputeNsubjettiness 1
941 set Beta 1.0
942 set AxisMode 4
943
944 set ComputeSoftDrop 1
945 set BetaSoftDrop 0.0
946 set SymmetryCutSoftDrop 0.1
947 set R0SoftDrop 0.8
948
949 set JetPTMin 25.0
950}
951
952#####################
953# MC truth jet finder
954#####################
955
956# TBC: is jet radius fine?
957
958module FastJetFinder GenJetFinder15 {
959 set InputArray NeutrinoFilter/filteredParticles
960
961 set OutputArray jets
962
963 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
964 set JetAlgorithm 6
965 set ParameterR 1.5
966
967 set ComputeNsubjettiness 1
968 set Beta 1.0
969 set AxisMode 4
970
971 set ComputeSoftDrop 1
972 set BetaSoftDrop 0.0
973 set SymmetryCutSoftDrop 0.1
974 set R0SoftDrop 1.5
975
976 set JetPTMin 25.0
977}
978
979
980##################
981# Fast Jet finder
982##################
983
984module FastJetFinder FastJetFinder02 {
985 set InputArray EFlowMerger/eflow
986
987 set OutputArray jets
988
989 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
990 set JetAlgorithm 6
991 set ParameterR 0.2
992
993 set ComputeNsubjettiness 1
994 set Beta 1.0
995 set AxisMode 4
996
997 set ComputeSoftDrop 1
998 set BetaSoftDrop 0.0
999 set SymmetryCutSoftDrop 0.1
1000 set R0SoftDrop 0.2
1001
1002 set JetPTMin 25.0
1003}
1004
1005##################
1006# Fast Jet finder
1007##################
1008
1009module FastJetFinder FastJetFinder04 {
1010 set InputArray EFlowMerger/eflow
1011
1012 set OutputArray jets
1013
1014 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1015 set JetAlgorithm 6
1016 set ParameterR 0.4
1017
1018 set ComputeNsubjettiness 1
1019 set Beta 1.0
1020 set AxisMode 4
1021
1022 set ComputeSoftDrop 1
1023 set BetaSoftDrop 0.0
1024 set SymmetryCutSoftDrop 0.1
1025 set R0SoftDrop 0.4
1026
1027 set JetPTMin 25.0
1028}
1029
1030
1031##################
1032# Fast Jet finder
1033##################
1034
1035module FastJetFinder FastJetFinder08 {
1036 set InputArray EFlowMerger/eflow
1037
1038 set OutputArray jets
1039
1040 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1041 set JetAlgorithm 6
1042 set ParameterR 0.8
1043
1044 set ComputeNsubjettiness 1
1045 set Beta 1.0
1046 set AxisMode 4
1047
1048 set ComputeSoftDrop 1
1049 set BetaSoftDrop 0.0
1050 set SymmetryCutSoftDrop 0.1
1051 set R0SoftDrop 0.8
1052
1053 set JetPTMin 25.0
1054}
1055
1056##################
1057# Fast Jet finder
1058##################
1059
1060module FastJetFinder FastJetFinder15 {
1061 set InputArray EFlowMerger/eflow
1062
1063 set OutputArray jets
1064
1065 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1066 set JetAlgorithm 6
1067 set ParameterR 1.5
1068
1069 set ComputeNsubjettiness 1
1070 set Beta 1.0
1071 set AxisMode 4
1072
1073 set ComputeSoftDrop 1
1074 set BetaSoftDrop 0.0
1075 set SymmetryCutSoftDrop 0.1
1076 set R0SoftDrop 1.5
1077
1078 set JetPTMin 25.0
1079}
1080
1081
1082##################
1083# Fast Jet finder
1084##################
1085
1086module FastJetFinder CaloJetFinder02 {
1087 set InputArray Calorimeter/towers
1088
1089 set OutputArray jets
1090
1091 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1092 set JetAlgorithm 6
1093 set ParameterR 0.2
1094
1095 set ComputeNsubjettiness 1
1096 set Beta 1.0
1097 set AxisMode 4
1098
1099 set ComputeSoftDrop 1
1100 set BetaSoftDrop 0.0
1101 set SymmetryCutSoftDrop 0.1
1102 set R0SoftDrop 0.2
1103
1104 set JetPTMin 25.0
1105}
1106
1107##################
1108# Fast Jet finder
1109##################
1110
1111module FastJetFinder CaloJetFinder04 {
1112 set InputArray Calorimeter/towers
1113
1114 set OutputArray jets
1115
1116 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1117 set JetAlgorithm 6
1118 set ParameterR 0.4
1119
1120 set ComputeNsubjettiness 1
1121 set Beta 1.0
1122 set AxisMode 4
1123
1124 set ComputeSoftDrop 1
1125 set BetaSoftDrop 0.0
1126 set SymmetryCutSoftDrop 0.1
1127 set R0SoftDrop 0.4
1128
1129 set JetPTMin 25.0
1130}
1131
1132
1133##################
1134# Fast Jet finder
1135##################
1136
1137module FastJetFinder CaloJetFinder08 {
1138 set InputArray Calorimeter/towers
1139
1140 set OutputArray jets
1141
1142 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1143 set JetAlgorithm 6
1144 set ParameterR 0.8
1145
1146 set ComputeNsubjettiness 1
1147 set Beta 1.0
1148 set AxisMode 4
1149
1150 set ComputeSoftDrop 1
1151 set BetaSoftDrop 0.0
1152 set SymmetryCutSoftDrop 0.1
1153 set R0SoftDrop 0.8
1154
1155 set JetPTMin 25.0
1156}
1157
1158##################
1159# Fast Jet finder
1160##################
1161
1162module FastJetFinder CaloJetFinder15 {
1163 set InputArray Calorimeter/towers
1164
1165 set OutputArray jets
1166
1167 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1168 set JetAlgorithm 6
1169 set ParameterR 1.5
1170
1171 set ComputeNsubjettiness 1
1172 set Beta 1.0
1173 set AxisMode 4
1174
1175 set ComputeSoftDrop 1
1176 set BetaSoftDrop 0.0
1177 set SymmetryCutSoftDrop 0.1
1178 set R0SoftDrop 1.5
1179
1180 set JetPTMin 25.0
1181}
1182
1183
1184##################
1185# Fast Jet finder
1186##################
1187
1188module FastJetFinder TrackJetFinder02 {
1189 set InputArray TrackMerger/tracks
1190
1191 set OutputArray jets
1192
1193 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1194 set JetAlgorithm 6
1195 set ParameterR 0.2
1196
1197 set ComputeNsubjettiness 1
1198 set Beta 1.0
1199 set AxisMode 4
1200
1201 set ComputeSoftDrop 1
1202 set BetaSoftDrop 0.0
1203 set SymmetryCutSoftDrop 0.1
1204 set R0SoftDrop 0.2
1205
1206 set JetPTMin 25.0
1207}
1208
1209##################
1210# Fast Jet finder
1211##################
1212
1213module FastJetFinder TrackJetFinder04 {
1214 set InputArray TrackMerger/tracks
1215
1216 set OutputArray jets
1217
1218 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1219 set JetAlgorithm 6
1220 set ParameterR 0.4
1221
1222 set ComputeNsubjettiness 1
1223 set Beta 1.0
1224 set AxisMode 4
1225
1226 set ComputeSoftDrop 1
1227 set BetaSoftDrop 0.0
1228 set SymmetryCutSoftDrop 0.1
1229 set R0SoftDrop 0.4
1230
1231 set JetPTMin 25.0
1232}
1233
1234
1235##################
1236# Fast Jet finder
1237##################
1238
1239module FastJetFinder TrackJetFinder08 {
1240 set InputArray TrackMerger/tracks
1241
1242 set OutputArray jets
1243
1244 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1245 set JetAlgorithm 6
1246 set ParameterR 0.8
1247
1248 set ComputeNsubjettiness 1
1249 set Beta 1.0
1250 set AxisMode 4
1251
1252 set ComputeSoftDrop 1
1253 set BetaSoftDrop 0.0
1254 set SymmetryCutSoftDrop 0.1
1255 set R0SoftDrop 0.8
1256
1257 set JetPTMin 25.0
1258}
1259
1260##################
1261# Fast Jet finder
1262##################
1263
1264module FastJetFinder TrackJetFinder15 {
1265 set InputArray TrackMerger/tracks
1266
1267 set OutputArray jets
1268
1269 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1270 set JetAlgorithm 6
1271 set ParameterR 1.5
1272
1273 set ComputeNsubjettiness 1
1274 set Beta 1.0
1275 set AxisMode 4
1276
1277 set ComputeSoftDrop 1
1278 set BetaSoftDrop 0.0
1279 set SymmetryCutSoftDrop 0.1
1280 set R0SoftDrop 1.5
1281
1282 set JetPTMin 25.0
1283}
1284
1285
1286##################
1287# Jet Energy Scale
1288##################
1289
1290module EnergyScale JetEnergyScale {
1291 set InputArray CaloJetFinder04/jets
1292 set OutputArray jets
1293
1294 # scale formula for jets
1295 set ScaleFormula {1.00}
1296}
1297
1298
1299########################
1300# Jet Flavor Association
1301########################
1302
1303module JetFlavorAssociation JetFlavorAssociation {
1304
1305 set PartonInputArray Delphes/partons
1306 set ParticleInputArray Delphes/allParticles
1307 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1308 set JetInputArray JetEnergyScale/jets
1309
1310 set DeltaR 0.5
1311 set PartonPTMin 5.0
1312 set PartonEtaMax 4.0
1313
1314}
1315
1316##################
1317# Photon isolation
1318##################
1319
1320# TBC: check values for iso cuts
1321
1322module Isolation PhotonIsolation {
1323 set CandidateInputArray ECal/eflowPhotons
1324 set IsolationInputArray EFlowFilter/eflow
1325
1326 set OutputArray photons
1327
1328 set DeltaRMax 0.3
1329 set PTMin 0.0
1330 set PTRatioMax 0.1
1331}
1332
1333
1334
1335###################
1336# Photon efficiency
1337###################
1338
1339module Efficiency PhotonEfficiency {
1340 set InputArray PhotonIsolation/photons
1341 set OutputArray photons
1342
1343 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1344
1345 set EfficiencyFormula {
1346 (pt <= 1.0) * (0.00) + \
1347 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
1348 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
1349 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
1350
1351 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
1352 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
1353 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
1354
1355 (abs(eta) > 4.0) * (0.00)}
1356
1357}
1358
1359
1360
1361####################
1362# Electron isolation
1363####################
1364
1365# TBC: check values for iso cuts
1366
1367module Isolation ElectronIsolation {
1368 set CandidateInputArray ElectronFilter/electrons
1369 set IsolationInputArray EFlowFilter/eflow
1370
1371 set OutputArray electrons
1372
1373 set DeltaRMax 0.3
1374
1375 set PTMin 0.0
1376
1377 set PTRatioMax 0.1
1378}
1379
1380
1381#######################
1382# Electron efficiency #
1383#######################
1384
1385module Efficiency ElectronEfficiency {
1386
1387 set InputArray ElectronIsolation/electrons
1388 set OutputArray electrons
1389
1390 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1391 # efficiency formula for electrons (taken from CMS PhaseII card)
1392 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
1393 set EfficiencyFormula {
1394 (pt <= 2.0) * (0.00) +
1395 (abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1396 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1397 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1398 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1399 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1400 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1401 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1402 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1403 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1404 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1405 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1406 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1407 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1408 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1409 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1410 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1411 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1412 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1413 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1414 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1415 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1416 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1417 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1418 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1419 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1420 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1421 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1422 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1423 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1424 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1425 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1426 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1427 (abs(eta) > 4.0) * (0.00)
1428
1429 }
1430
1431}
1432
1433
1434
1435################
1436# Muon isolation
1437################
1438
1439# TBC: check values for iso cuts
1440
1441module Isolation MuonIsolation {
1442 set CandidateInputArray MuonMomentumSmearing/muons
1443 set IsolationInputArray EFlowFilter/eflow
1444
1445 set OutputArray muons
1446
1447 set DeltaRMax 0.3
1448
1449 set PTMin 0.0
1450
1451 set PTRatioMax 0.2
1452}
1453
1454####################
1455# Muon Efficiency #
1456####################
1457
1458module Efficiency MuonEfficiency {
1459 set InputArray MuonIsolation/muons
1460 set OutputArray muons
1461 # tracking + TightID efficiency formula for muons
1462 set EfficiencyFormula {
1463 (pt <= 2.0) * (0.00) +
1464 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1465 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1466 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1467 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1468 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1469 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1470 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1471 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1472 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1473 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1474 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1475 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
1476 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1477 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
1478 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
1479 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1480 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
1481 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
1482 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1483 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
1484 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1485 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1486 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1487 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1488 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
1489 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1490 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
1491 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
1492 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
1493 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1494 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1495 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1496 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1497 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
1498 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
1499 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1500 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1501 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1502 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1503 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
1504 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1505 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1506 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1507 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1508 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1509 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1510 (abs(eta) > 4.0) * (0.00)
1511 }
1512}
1513
1514
1515
1516###########
1517# b-tagging
1518###########
1519
1520module BTagging BTagging {
1521 set JetInputArray JetEnergyScale/jets
1522
1523 set BitNumber 0
1524
1525 add EfficiencyFormula {0} {
1526
1527 (pt <= 10.0) * (0.00) +
1528 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
1529 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
1530 (abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
1531 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
1532 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
1533 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1534 (abs(eta) > 4.0) * (0.00)}
1535
1536 add EfficiencyFormula {4} {
1537
1538 (pt <= 10.0) * (0.00) +
1539 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
1540 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
1541 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1542 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
1543 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
1544 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1545 (abs(eta) > 4.0) * (0.00)}
1546
1547 add EfficiencyFormula {5} {
1548
1549 (pt <= 10.0) * (0.00) +
1550 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
1551 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
1552 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
1553 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
1554 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
1555 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
1556 (abs(eta) >= 4.0) * (0.00)}
1557
1558}
1559
1560
1561#############
1562# tau-tagging
1563#############
1564
1565
1566module TauTagging TauTagging {
1567 set ParticleInputArray Delphes/allParticles
1568 set PartonInputArray Delphes/partons
1569 set JetInputArray JetEnergyScale/jets
1570
1571 set DeltaR 0.5
1572
1573 set TauPTMin 1.0
1574
1575 set TauEtaMax 4.0
1576
1577 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1578 add EfficiencyFormula {0} {
1579
1580 (pt <= 10.0) * (0.00) +
1581 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
1582 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
1583 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1584 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
1585 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
1586 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1587 (abs(eta) > 4.0) * (0.00)}
1588
1589 add EfficiencyFormula {11} {
1590
1591 (pt <= 10.0) * (0.00) +
1592 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
1593 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
1594 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1595 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
1596 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
1597 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1598 (abs(eta) > 4.0) * (0.00)}
1599
1600 add EfficiencyFormula {15} {
1601
1602 (pt <= 10.0) * (0.00) +
1603 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
1604 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
1605 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1606 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
1607 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
1608 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1609 (abs(eta) > 4.0) * (0.00)}
1610
1611}
1612
1613#####################################################
1614# Find uniquely identified photons/electrons/tau/jets
1615#####################################################
1616
1617module UniqueObjectFinder UniqueObjectFinder {
1618# earlier arrays take precedence over later ones
1619# add InputArray InputArray OutputArray
1620 add InputArray PhotonEfficiency/photons photons
1621 add InputArray ElectronEfficiency/electrons electrons
1622 add InputArray MuonEfficiency/muons muons
1623 add InputArray JetEnergyScale/jets jets
1624
1625}
1626
1627
1628##################
1629# ROOT tree writer
1630##################
1631
1632module TreeWriter TreeWriter {
1633# add Branch InputArray BranchName BranchClass
1634 add Branch Delphes/allParticles Particle GenParticle
1635
1636 add Branch GenJetFinder/jets GenJet Jet
1637 add Branch GenMissingET/momentum GenMissingET MissingET
1638
1639 add Branch TrackMerger/tracks Track Track
1640 add Branch Calorimeter/towers Tower Tower
1641
1642 add Branch HCal/eflowTracks EFlowTrack Track
1643 add Branch ECal/eflowPhotons EFlowPhoton Tower
1644 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1645
1646 add Branch UniqueObjectFinder/photons Photon Photon
1647 add Branch UniqueObjectFinder/electrons Electron Electron
1648 add Branch UniqueObjectFinder/muons Muon Muon
1649 add Branch UniqueObjectFinder/jets Jet Jet
1650
1651 add Branch GenJetFinder02/jets GenJet02 Jet
1652 add Branch GenJetFinder04/jets GenJet04 Jet
1653 add Branch GenJetFinder08/jets GenJet08 Jet
1654 add Branch GenJetFinder15/jets GenJet15 Jet
1655
1656 add Branch FastJetFinder02/jets ParticleFlowJet02 Jet
1657 add Branch FastJetFinder04/jets ParticleFlowJet04 Jet
1658 add Branch FastJetFinder08/jets ParticleFlowJet08 Jet
1659 add Branch FastJetFinder15/jets ParticleFlowJet15 Jet
1660
1661 add Branch CaloJetFinder02/jets CaloJet02 Jet
1662 add Branch CaloJetFinder04/jets CaloJet04 Jet
1663 add Branch CaloJetFinder08/jets CaloJet08 Jet
1664 add Branch CaloJetFinder15/jets CaloJet15 Jet
1665
1666 add Branch TrackJetFinder02/jets TrackJet02 Jet
1667 add Branch TrackJetFinder04/jets TrackJet04 Jet
1668 add Branch TrackJetFinder08/jets TrackJet08 Jet
1669 add Branch TrackJetFinder15/jets TrackJet15 Jet
1670
1671 add Branch MissingET/momentum MissingET MissingET
1672 add Branch ScalarHT/energy ScalarHT ScalarHT
1673}
Note: See TracBrowser for help on using the repository browser.