Fork me on GitHub

source: git/cards/delphes_card_MuonColliderDet.tcl@ e01b141

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

added forward muon collection, jet substructure to inclusive jets
  • Property mode set to 100644
File size: 134.5 KB
Line 
1#######################################
2# Muon Collider Detector TARGET model
3#
4# Michele Selvaggi michele.selvaggi@cern.ch
5# Ulrike Schnoor ulrike.schnoor@cern.ch
6#
7#
8# !!! DISCLAIMER !!!
9#
10# The parameterisation of the Muon Collider
11# has to be intended as a target performance.
12# This has not been validated by full simulation.
13# Very similar to CLIC performance, added DenseTrackFilter
14#
15# Comments:
16# - fix: replace momentum resolution by FCC-hh
17# - fix: angle param in DenseTrackFilter
18# - fix: add FWD muon collection
19# - what to do with Eta acceptance?
20# - for now everything reduced at 2.5
21# - add jet substructure to valencia jets ?
22# - added R02 jets
23# - added electron misId for taus, and make flat eff
24#
25#######################################
26
27#######################################
28# Order of execution of various modules
29#######################################
30
31set ExecutionPath {
32 ParticlePropagator
33 TrackMergerProp
34
35 DenseProp
36 DenseMergeTracks
37 DenseTrackFilter
38
39 ChargedHadronTrackingEfficiency
40 ElectronTrackingEfficiency
41 MuonTrackingEfficiency
42 ForwardMuonEfficiency
43
44 ChargedHadronMomentumSmearing
45 ElectronMomentumSmearing
46 MuonMomentumSmearing
47 ForwardMuonMomentumSmearing
48
49 TrackMerger
50
51 ECal
52 HCal
53
54 Calorimeter
55 EFlowMerger
56 EFlowFilter
57
58 PhotonEfficiency
59 PhotonIsolation
60
61 ElectronFilter
62 ElectronEfficiency
63 ElectronIsolation
64
65 ChargedHadronFilter
66
67 MuonEfficiency
68 MuonIsolation
69
70 UniqueObjectFinder
71
72 NeutrinoFilter
73 GenJetFinder
74
75
76 FastJetFinderKt
77 FastJetFinderVLC_R02_N2
78 FastJetFinderVLC_R02_N3
79 FastJetFinderVLC_R02_N4
80 FastJetFinderVLC_R02_N5
81 FastJetFinderVLC_R02_N6
82 FastJetFinderVLC_R05_N2
83 FastJetFinderVLC_R05_N3
84 FastJetFinderVLC_R05_N4
85 FastJetFinderVLC_R05_N5
86 FastJetFinderVLC_R05_N6
87 FastJetFinderVLC_R07_N2
88 FastJetFinderVLC_R07_N3
89 FastJetFinderVLC_R07_N4
90 FastJetFinderVLC_R07_N5
91 FastJetFinderVLC_R07_N6
92 FastJetFinderVLC_R10_N2
93 FastJetFinderVLC_R10_N3
94 FastJetFinderVLC_R10_N4
95 FastJetFinderVLC_R10_N5
96 FastJetFinderVLC_R10_N6
97 FastJetFinderVLC_R12_N2
98 FastJetFinderVLC_R12_N3
99 FastJetFinderVLC_R12_N4
100 FastJetFinderVLC_R12_N5
101 FastJetFinderVLC_R12_N6
102 FastJetFinderVLC_R15_N2
103 FastJetFinderVLC_R15_N3
104 FastJetFinderVLC_R15_N4
105 FastJetFinderVLC_R15_N5
106 FastJetFinderVLC_R15_N6
107
108 FastJetFinderVLC_R05_inclusive
109 FastJetFinderVLC_R02_inclusive
110 FastJetFinderVLC_R07_inclusive
111 FastJetFinderVLC_R10_inclusive
112 FastJetFinderVLC_R12_inclusive
113 FastJetFinderVLC_R15_inclusive
114
115 MissingET
116 GenMissingET
117
118 JetMomentumSmearing_VLCR02N2
119 JetMomentumSmearing_VLCR02N3
120 JetMomentumSmearing_VLCR02N4
121 JetMomentumSmearing_VLCR02N5
122 JetMomentumSmearing_VLCR02N6
123 JetMomentumSmearing_VLCR02_inclusive
124 JetMomentumSmearing_VLCR05N2
125 JetMomentumSmearing_VLCR05N3
126 JetMomentumSmearing_VLCR05N4
127 JetMomentumSmearing_VLCR05N5
128 JetMomentumSmearing_VLCR05N6
129 JetMomentumSmearing_VLCR05_inclusive
130 JetMomentumSmearing_VLCR07N2
131 JetMomentumSmearing_VLCR07N3
132 JetMomentumSmearing_VLCR07N4
133 JetMomentumSmearing_VLCR07N5
134 JetMomentumSmearing_VLCR07N6
135 JetMomentumSmearing_VLCR07_inclusive
136 JetMomentumSmearing_VLCR10N2
137 JetMomentumSmearing_VLCR10N3
138 JetMomentumSmearing_VLCR10N4
139 JetMomentumSmearing_VLCR10N5
140 JetMomentumSmearing_VLCR10N6
141 JetMomentumSmearing_VLCR10_inclusive
142 JetMomentumSmearing_VLCR12N2
143 JetMomentumSmearing_VLCR12N3
144 JetMomentumSmearing_VLCR12N4
145 JetMomentumSmearing_VLCR12N5
146 JetMomentumSmearing_VLCR12N6
147 JetMomentumSmearing_VLCR12_inclusive
148 JetMomentumSmearing_VLCR15N2
149 JetMomentumSmearing_VLCR15N3
150 JetMomentumSmearing_VLCR15N4
151 JetMomentumSmearing_VLCR15N5
152 JetMomentumSmearing_VLCR15N6
153 JetMomentumSmearing_VLCR15_inclusive
154
155
156 JetFlavorAssociation_R02N2
157 JetFlavorAssociation_R02N3
158 JetFlavorAssociation_R02N4
159 JetFlavorAssociation_R02N5
160 JetFlavorAssociation_R02N6
161 JetFlavorAssociation_R05N2
162 JetFlavorAssociation_R05N3
163 JetFlavorAssociation_R05N4
164 JetFlavorAssociation_R05N5
165 JetFlavorAssociation_R05N6
166
167 JetFlavorAssociation_R07N2
168 JetFlavorAssociation_R07N3
169 JetFlavorAssociation_R07N4
170 JetFlavorAssociation_R07N5
171 JetFlavorAssociation_R07N6
172
173 JetFlavorAssociation_R10N2
174 JetFlavorAssociation_R10N3
175 JetFlavorAssociation_R10N4
176 JetFlavorAssociation_R10N5
177 JetFlavorAssociation_R10N6
178
179 JetFlavorAssociation_R12N2
180 JetFlavorAssociation_R12N3
181 JetFlavorAssociation_R12N4
182 JetFlavorAssociation_R12N5
183 JetFlavorAssociation_R12N6
184
185 JetFlavorAssociation_R15N2
186 JetFlavorAssociation_R15N3
187 JetFlavorAssociation_R15N4
188 JetFlavorAssociation_R15N5
189 JetFlavorAssociation_R15N6
190
191 JetFlavorAssociation_R02_inclusive
192 JetFlavorAssociation_R05_inclusive
193 JetFlavorAssociation_R07_inclusive
194 JetFlavorAssociation_R10_inclusive
195 JetFlavorAssociation_R12_inclusive
196 JetFlavorAssociation_R15_inclusive
197
198
199 BTagging_WP50_R02N2
200 BTagging_WP70_R02N2
201 BTagging_WP90_R02N2
202 BTagging_WP50_R02N3
203 BTagging_WP70_R02N3
204 BTagging_WP90_R02N3
205 BTagging_WP50_R02N4
206 BTagging_WP70_R02N4
207 BTagging_WP90_R02N4
208 BTagging_WP50_R02N5
209 BTagging_WP70_R02N5
210 BTagging_WP90_R02N5
211 BTagging_WP50_R02N6
212 BTagging_WP70_R02N6
213 BTagging_WP90_R02N6
214 BTagging_WP50_R05N2
215 BTagging_WP70_R05N2
216 BTagging_WP90_R05N2
217 BTagging_WP50_R05N3
218 BTagging_WP70_R05N3
219 BTagging_WP90_R05N3
220 BTagging_WP50_R05N4
221 BTagging_WP70_R05N4
222 BTagging_WP90_R05N4
223 BTagging_WP50_R05N5
224 BTagging_WP70_R05N5
225 BTagging_WP90_R05N5
226 BTagging_WP50_R05N6
227 BTagging_WP70_R05N6
228 BTagging_WP90_R05N6
229 BTagging_WP50_R07N2
230 BTagging_WP70_R07N2
231 BTagging_WP90_R07N2
232 BTagging_WP50_R07N3
233 BTagging_WP70_R07N3
234 BTagging_WP90_R07N3
235 BTagging_WP50_R07N4
236 BTagging_WP70_R07N4
237 BTagging_WP90_R07N4
238 BTagging_WP50_R07N5
239 BTagging_WP70_R07N5
240 BTagging_WP90_R07N5
241 BTagging_WP50_R07N6
242 BTagging_WP70_R07N6
243 BTagging_WP90_R07N6
244 BTagging_WP50_R10N2
245 BTagging_WP70_R10N2
246 BTagging_WP90_R10N2
247 BTagging_WP50_R10N3
248 BTagging_WP70_R10N3
249 BTagging_WP90_R10N3
250 BTagging_WP50_R10N4
251 BTagging_WP70_R10N4
252 BTagging_WP90_R10N4
253 BTagging_WP50_R10N5
254 BTagging_WP70_R10N5
255 BTagging_WP90_R10N5
256 BTagging_WP50_R10N6
257 BTagging_WP70_R10N6
258 BTagging_WP90_R10N6
259 BTagging_WP50_R12N2
260 BTagging_WP70_R12N2
261 BTagging_WP90_R12N2
262 BTagging_WP50_R12N3
263 BTagging_WP70_R12N3
264 BTagging_WP90_R12N3
265 BTagging_WP50_R12N4
266 BTagging_WP70_R12N4
267 BTagging_WP90_R12N4
268 BTagging_WP50_R12N5
269 BTagging_WP70_R12N5
270 BTagging_WP90_R12N5
271 BTagging_WP50_R12N6
272 BTagging_WP70_R12N6
273 BTagging_WP90_R12N6
274 BTagging_WP50_R15N2
275 BTagging_WP70_R15N2
276 BTagging_WP90_R15N2
277 BTagging_WP50_R15N3
278 BTagging_WP70_R15N3
279 BTagging_WP90_R15N3
280 BTagging_WP50_R15N4
281 BTagging_WP70_R15N4
282 BTagging_WP90_R15N4
283 BTagging_WP50_R15N5
284 BTagging_WP70_R15N5
285 BTagging_WP90_R15N5
286 BTagging_WP50_R15N6
287 BTagging_WP70_R15N6
288 BTagging_WP90_R15N6
289 BTagging_WP50_R02_inclusive
290 BTagging_WP70_R02_inclusive
291 BTagging_WP90_R02_inclusive
292 BTagging_WP50_R05_inclusive
293 BTagging_WP70_R05_inclusive
294 BTagging_WP90_R05_inclusive
295 BTagging_WP50_R07_inclusive
296 BTagging_WP70_R07_inclusive
297 BTagging_WP90_R07_inclusive
298 BTagging_WP50_R10_inclusive
299 BTagging_WP70_R10_inclusive
300 BTagging_WP90_R10_inclusive
301 BTagging_WP50_R12_inclusive
302 BTagging_WP70_R12_inclusive
303 BTagging_WP90_R12_inclusive
304 BTagging_WP50_R15_inclusive
305 BTagging_WP70_R15_inclusive
306 BTagging_WP90_R15_inclusive
307
308
309 TauTagging_R02N2
310 TauTagging_R02N3
311 TauTagging_R02N4
312 TauTagging_R02N5
313 TauTagging_R02N6
314 TauTagging_R05N2
315 TauTagging_R05N3
316 TauTagging_R05N4
317 TauTagging_R05N5
318 TauTagging_R05N6
319 TauTagging_R07N2
320 TauTagging_R07N3
321 TauTagging_R07N4
322 TauTagging_R07N5
323 TauTagging_R07N6
324 TauTagging_R10N2
325 TauTagging_R10N3
326 TauTagging_R10N4
327 TauTagging_R10N5
328 TauTagging_R10N6
329 TauTagging_R12N2
330 TauTagging_R12N3
331 TauTagging_R12N4
332 TauTagging_R12N5
333 TauTagging_R12N6
334 TauTagging_R15N2
335 TauTagging_R15N3
336 TauTagging_R15N4
337 TauTagging_R15N5
338 TauTagging_R15N6
339 TauTagging_R02_inclusive
340 TauTagging_R05_inclusive
341 TauTagging_R07_inclusive
342 TauTagging_R10_inclusive
343 TauTagging_R12_inclusive
344 TauTagging_R15_inclusive
345
346 ScalarHT
347
348 TreeWriter
349}
350
351#################################
352# Propagate particles in cylinder
353#################################
354
355module ParticlePropagator ParticlePropagator {
356 set InputArray Delphes/stableParticles
357
358 set OutputArray stableParticles
359 set ChargedHadronOutputArray chargedHadrons
360 set ElectronOutputArray electrons
361 set MuonOutputArray muons
362
363 # radius of the magnetic field coverage in the calorimeter, in m
364 set Radius 1.5
365 # half-length of the magnetic field coverage in the calorimeter, in m
366 set HalfLength 2.31
367
368 # magnetic field, in T
369 set Bz 4.0
370}
371
372
373##############
374# Track merger
375##############
376
377module Merger TrackMergerProp {
378# add InputArray InputArray
379 add InputArray ParticlePropagator/chargedHadrons
380 add InputArray ParticlePropagator/electrons
381 add InputArray ParticlePropagator/muons
382 set OutputArray tracks
383}
384
385
386
387####################################
388# Track propagation to pseudo-pixel
389####################################
390
391module ParticlePropagator DenseProp {
392
393 set InputArray TrackMergerProp/tracks
394
395 # radius of the magnetic field coverage, in m
396 set Radius 0.45
397 set RadiusMax 1.5
398 # half-length of the magnetic field coverage, in m
399 set HalfLength 0.8
400 set HalfLengthMax 2.31
401
402 # magnetic field
403 set Bz 4.0
404}
405
406#####################
407# Dense Track merger
408#####################
409
410module Merger DenseMergeTracks {
411# add InputArray InputArray
412 add InputArray DenseProp/chargedHadrons
413 add InputArray DenseProp/electrons
414 add InputArray DenseProp/muons
415 set OutputArray tracks
416}
417
418
419######################
420# Dense Track Filter
421######################
422
423module DenseTrackFilter DenseTrackFilter {
424
425 set TrackInputArray DenseMergeTracks/tracks
426
427 set TrackOutputArray tracks
428 set ChargedHadronOutputArray chargedHadrons
429 set ElectronOutputArray electrons
430 set MuonOutputArray muons
431
432 set EtaPhiRes 0.01
433 set EtaMax 2.5
434
435 set pi [expr {acos(-1)}]
436
437 set nbins_phi [expr {$pi/$EtaPhiRes} ]
438 set nbins_phi [expr {int($nbins_phi)} ]
439
440 set PhiBins {}
441 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
442 add PhiBins [expr {$i * $pi/$nbins_phi}]
443 }
444
445 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
446 set nbins_eta [expr {int($nbins_eta)} ]
447
448 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
449 set eta [expr {$i * $EtaPhiRes}]
450 add EtaPhiBins $eta $PhiBins
451 }
452}
453
454
455
456
457####################################
458# Charged hadron tracking efficiency
459####################################
460
461module Efficiency ChargedHadronTrackingEfficiency {
462 set InputArray DenseTrackFilter/chargedHadrons
463 set OutputArray chargedHadrons
464 # tracking efficiency formula for charged hadrons
465
466 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
467 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
468 (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
469 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
470 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
471 (abs(eta) > 2.5 ) * (0.00)
472 }
473}
474
475##############################
476# Electron tracking efficiency
477##############################
478
479module Efficiency ElectronTrackingEfficiency {
480 set InputArray DenseTrackFilter/electrons
481 set OutputArray electrons
482
483 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
484 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
485 (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
486 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
487 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
488 (abs(eta) > 2.5 ) * (0.00)
489 }
490}
491
492##########################
493# Muon tracking efficiency
494##########################
495
496module Efficiency MuonTrackingEfficiency {
497 set InputArray DenseTrackFilter/muons
498 set OutputArray muons
499
500 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
501 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.95) +
502 (abs(eta) <= 2.0) * (pt > 1) * (0.99) +
503 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.90) +
504 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.95) +
505 (abs(eta) > 2.5 ) * (0.00)
506 }
507}
508
509##########################
510# Forward Muon efficiency
511##########################
512
513## hypothetical forward muon spectrometer
514module Efficiency ForwardMuonEfficiency {
515 set InputArray ParticlePropagator/muons
516 set OutputArray muons
517
518 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
519 (abs(eta) > 2.5 && abs(eta) < 6.0) * (pt > 0.5 && pt <= 1) * (0.90) +
520 (abs(eta) > 2.5 && abs(eta) < 6.0) * (pt > 1.0) * (0.95) +
521 (abs(eta) > 6.0 ) * (0.00)
522 }
523}
524
525
526########################################
527# Momentum resolution for charged tracks
528########################################
529
530module MomentumSmearing ChargedHadronMomentumSmearing {
531 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
532 set OutputArray chargedHadrons
533
534
535 # Resolution given in dpT/pT (from FCC-hh)
536 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
537 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
538 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
539 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
540 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
541 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
542 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
543 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
544 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
545 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
546 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
547 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
548 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
549 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
550 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
551 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
552 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
553 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
554 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
555 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
556 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
557 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
558 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
559 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
560 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
561 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
562 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
563 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
564 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
565 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
566 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
567 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
568 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
569 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
570 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
571 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
572 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
573 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
574 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
575 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
576 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
577 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
578 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
579 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
580 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
581 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
582 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
583 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
584 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
585 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
586 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
587 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
588 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
589 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
590 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
591 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
592 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
593 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
594 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
595 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
596 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
597 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
598 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
599 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
600 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
601 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
602 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
603 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
604 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
605 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
606 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
607 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
608 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
609 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
610 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
611 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
612 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
613 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
614 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
615 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
616 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
617 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
618 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
619 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
620 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
621 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
622 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
623 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
624 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
625 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
626 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
627 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
628 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
629 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
630 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
631 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
632 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
633 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
634 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
635 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
636 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
637 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
638 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
639 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
640 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
641 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
642 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
643 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
644 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
645 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
646 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
647 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
648 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
649 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
650 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
651 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
652 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
653 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
654 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
655 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
656 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
657 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
658 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
659 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
660 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
661 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
662 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
663 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
664 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
665 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
666 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
667 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
668 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
669 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
670 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
671 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
672 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
673 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
674 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
675 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
676 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
677 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
678 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
679 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
680 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
681 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
682 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
683 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
684 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
685 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
686 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
687 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
688 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
689 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
690 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
691 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
692 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
693 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
694 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
695 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
696 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
697 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
698 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
699 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
700 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
701 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
702 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
703 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
704 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
705 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
706 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
707 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
708 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
709 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
710 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
711 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
712 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
713 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
714 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
715 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
716 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
717 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
718 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
719 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
720 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
721 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
722 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
723 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
724 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
725 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
726 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
727 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
728 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
729 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
730 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
731 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
732 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
733 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
734 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
735 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
736 }
737}
738
739###################################
740# Momentum resolution for electrons
741###################################
742
743module MomentumSmearing ElectronMomentumSmearing {
744 set InputArray ElectronTrackingEfficiency/electrons
745 set OutputArray electrons
746
747 # Resolution given in dpT/pT (from FCC-hh)
748 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
749 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
750 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
751 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
752 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
753 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
754 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
755 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
756 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
757 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
758 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
759 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
760 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
761 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
762 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
763 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
764 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
765 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
766 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
767 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
768 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
769 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
770 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
771 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
772 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
773 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
774 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
775 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
776 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
777 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
778 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
779 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
780 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
781 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
782 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
783 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
784 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
785 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
786 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
787 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
788 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
789 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
790 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
791 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
792 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
793 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
794 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
795 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
796 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
797 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
798 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
799 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
800 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
801 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
802 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
803 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
804 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
805 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
806 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
807 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
808 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
809 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
810 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
811 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
812 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
813 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
814 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
815 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
816 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
817 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
818 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
819 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
820 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
821 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
822 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
823 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
824 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
825 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
826 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
827 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
828 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
829 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
830 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
831 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
832 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
833 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
834 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
835 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
836 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
837 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
838 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
839 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
840 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
841 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
842 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
843 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
844 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
845 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
846 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
847 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
848 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
849 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
850 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
851 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
852 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
853 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
854 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
855 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
856 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
857 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
858 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
859 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
860 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
861 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
862 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
863 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
864 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
865 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
866 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
867 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
868 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
869 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
870 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
871 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
872 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
873 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
874 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
875 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
876 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
877 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
878 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
879 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
880 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
881 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
882 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
883 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
884 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
885 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
886 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
887 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
888 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
889 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
890 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
891 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
892 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
893 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
894 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
895 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
896 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
897 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
898 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
899 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
900 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
901 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
902 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
903 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
904 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
905 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
906 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
907 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
908 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
909 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
910 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
911 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
912 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
913 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
914 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
915 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
916 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
917 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
918 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
919 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
920 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
921 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
922 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
923 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
924 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
925 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
926 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
927 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
928 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
929 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
930 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
931 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
932 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
933 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
934 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
935 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
936 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
937 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
938 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
939 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
940 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
941 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
942 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
943 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
944 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
945 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
946 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
947 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
948 }
949}
950
951###############################
952# Momentum resolution for muons
953###############################
954
955module MomentumSmearing MuonMomentumSmearing {
956 set InputArray MuonTrackingEfficiency/muons
957 set OutputArray muons
958
959 # Resolution given in dpT/pT.
960
961 set ResolutionFormula {
962
963 ( abs(eta) < 1.35 ) *
964
965 (sqrt(0.0000364164 + (
966 9.06262e-8 *pt^2* cosh(
967 eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + (
968 504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/(
969 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
970 sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/(
971 0.00516429/sin(2*atan(exp(-abs(eta))))^2 + (
972 96868.8 *(1/400000000 + (0.117945 * 1/cosh(eta)^2)/(
973 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
974 sin(2*atan(exp(-abs(eta))))^2))
975
976 ) +
977
978 ( abs(eta) > 1.35 && abs(eta) < 2.00) *
979
980 ( 1.53846*(
981
982 (
983 (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
984 (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
985 (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
986 (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
987 (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
988 (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
989 (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
990 (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
991 (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
992 (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
993 (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
994 (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
995 (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
996 (energy >= 10000.0000) * (0.062437*energy/10000.000000)
997 )
998
999
1000 - sqrt(0.0000364164 + (9.06262e-8 *(1.19507e-6 + 2137.54 *(1/400000000 + 0.155982/pt^2)) *pt^2)/(0.0218797 + 410407. *(1/400000000 + 0.155982/pt^2))))*abs(eta) +
1001
1002 3.07692*sqrt(0.0000364164 + (9.06262e-8 *(1.19507e-6 + 2137.54 *(1/400000000 + 0.155982/pt^2)) *pt^2)/(0.0218797 + 410407. *(1/400000000 + 0.155982/pt^2)))
1003
1004 - 2.07692*
1005 (
1006
1007 (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
1008 (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
1009 (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
1010 (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
1011 (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
1012 (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
1013 (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
1014 (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
1015 (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
1016 (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
1017 (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
1018 (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
1019 (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
1020 (energy >= 10000.0000) * (0.062437*energy/10000.000000)
1021
1022 )
1023
1024
1025 ) +
1026
1027 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01062416) + \
1028 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010624 + (energy-1.000000)* -0.005532) + \
1029 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005092 + (energy-2.000000)* -0.000326) + \
1030 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004115 + (energy-5.000000)* 0.000006) + \
1031 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 20.0000) * (0.004143 + (energy-10.000000)* 0.000007) + \
1032 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 20.0000 && energy < 50.0000) * (0.004209 + (energy-20.000000)* 0.000007) + \
1033 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 50.0000 && energy < 100.0000) * (0.004413 + (energy-50.000000)* 0.000005) + \
1034 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 200.0000) * (0.004681 + (energy-100.000000)* 0.000004) + \
1035 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 200.0000 && energy < 500.0000) * (0.005076 + (energy-200.000000)* 0.000004) + \
1036 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006270 + (energy-500.000000)* 0.000005) + \
1037 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.008960 + (energy-1000.000000)* 0.000007) + \
1038 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.015510 + (energy-2000.000000)* 0.000007) + \
1039 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.036867 + (energy-5000.000000)* 0.000007) + \
1040 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073168*energy/10000.000000) + \
1041 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01007098) + \
1042 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010071 + (energy-1.000000)* -0.004627) + \
1043 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000322) + \
1044 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004478 + (energy-5.000000)* 0.000005) + \
1045 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 20.0000) * (0.004501 + (energy-10.000000)* 0.000006) + \
1046 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 20.0000 && energy < 50.0000) * (0.004558 + (energy-20.000000)* 0.000005) + \
1047 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 50.0000 && energy < 100.0000) * (0.004701 + (energy-50.000000)* 0.000004) + \
1048 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 200.0000) * (0.004888 + (energy-100.000000)* 0.000003) + \
1049 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 200.0000 && energy < 500.0000) * (0.005213 + (energy-200.000000)* 0.000004) + \
1050 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006429 + (energy-500.000000)* 0.000006) + \
1051 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.009343 + (energy-1000.000000)* 0.000007) + \
1052 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.016410 + (energy-2000.000000)* 0.000008) + \
1053 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.039265 + (energy-5000.000000)* 0.000008) + \
1054 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078014*energy/10000.000000) + \
1055 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01095892) + \
1056 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010959 + (energy-1.000000)* -0.005458) + \
1057 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005501 + (energy-2.000000)* -0.000281) + \
1058 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004660 + (energy-5.000000)* 0.000005) + \
1059 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 20.0000) * (0.004686 + (energy-10.000000)* 0.000007) + \
1060 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 20.0000 && energy < 50.0000) * (0.004757 + (energy-20.000000)* 0.000006) + \
1061 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 50.0000 && energy < 100.0000) * (0.004937 + (energy-50.000000)* 0.000004) + \
1062 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 200.0000) * (0.005143 + (energy-100.000000)* 0.000004) + \
1063 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 200.0000 && energy < 500.0000) * (0.005505 + (energy-200.000000)* 0.000005) + \
1064 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006975 + (energy-500.000000)* 0.000007) + \
1065 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.010462 + (energy-1000.000000)* 0.000008) + \
1066 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.018731 + (energy-2000.000000)* 0.000009) + \
1067 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.045159 + (energy-5000.000000)* 0.000009) + \
1068 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.089830*energy/10000.000000) + \
1069 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01279214) + \
1070 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012792 + (energy-1.000000)* -0.005763) + \
1071 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.007029 + (energy-2.000000)* -0.000513) + \
1072 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005489 + (energy-5.000000)* 0.000003) + \
1073 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 20.0000) * (0.005503 + (energy-10.000000)* 0.000003) + \
1074 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 20.0000 && energy < 50.0000) * (0.005537 + (energy-20.000000)* 0.000003) + \
1075 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 50.0000 && energy < 100.0000) * (0.005636 + (energy-50.000000)* 0.000003) + \
1076 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 200.0000) * (0.005773 + (energy-100.000000)* 0.000003) + \
1077 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 200.0000 && energy < 500.0000) * (0.006060 + (energy-200.000000)* 0.000005) + \
1078 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 500.0000 && energy < 1000.0000) * (0.007489 + (energy-500.000000)* 0.000007) + \
1079 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.011117 + (energy-1000.000000)* 0.000009) + \
1080 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.019824 + (energy-2000.000000)* 0.000009) + \
1081 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.047732 + (energy-5000.000000)* 0.000009) + \
1082 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.094931*energy/10000.000000) + \
1083 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01502671) + \
1084 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015027 + (energy-1.000000)* -0.007177) + \
1085 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007850 + (energy-2.000000)* -0.000651) + \
1086 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005898 + (energy-5.000000)* 0.000003) + \
1087 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 20.0000) * (0.005913 + (energy-10.000000)* 0.000003) + \
1088 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 20.0000 && energy < 50.0000) * (0.005947 + (energy-20.000000)* 0.000003) + \
1089 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 50.0000 && energy < 100.0000) * (0.006039 + (energy-50.000000)* 0.000003) + \
1090 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 200.0000) * (0.006170 + (energy-100.000000)* 0.000003) + \
1091 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 200.0000 && energy < 500.0000) * (0.006485 + (energy-200.000000)* 0.000006) + \
1092 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 500.0000 && energy < 1000.0000) * (0.008140 + (energy-500.000000)* 0.000008) + \
1093 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.012304 + (energy-1000.000000)* 0.000010) + \
1094 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.022168 + (energy-2000.000000)* 0.000010) + \
1095 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.053585 + (energy-5000.000000)* 0.000011) + \
1096 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.106635*energy/10000.000000)
1097 }
1098}
1099
1100###############################
1101# Momentum resolution for muons
1102###############################
1103
1104module MomentumSmearing ForwardMuonMomentumSmearing {
1105 set InputArray ForwardMuonEfficiency/muons
1106 set OutputArray muons
1107
1108 # Resolution given in dpT/pT (assume forward tracking spectrometer from FCC-hh)
1109
1110 set ResolutionFormula {
1111 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.01778926) + \
1112 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.017789 + (energy-1.000000)* -0.008811) + \
1113 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.008978 + (energy-2.000000)* -0.000844) + \
1114 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.006446 + (energy-5.000000)* -0.000195) + \
1115 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.005473 + (energy-10.000000)* 0.000003) + \
1116 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005782 + (energy-100.000000)* 0.000005) + \
1117 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010655 + (energy-1000.000000)* 0.000009) + \
1118 (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 10000.0000) * (0.087190*energy/10000.000000) + \
1119 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.02061617) + \
1120 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.020616 + (energy-1.000000)* -0.010732) + \
1121 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.009884 + (energy-2.000000)* -0.001349) + \
1122 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.005836 + (energy-5.000000)* -0.000484) + \
1123 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003417 + (energy-10.000000)* 0.000003) + \
1124 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003655 + (energy-100.000000)* 0.000002) + \
1125 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005164 + (energy-1000.000000)* 0.000003) + \
1126 (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 10000.0000) * (0.029564*energy/10000.000000) + \
1127 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.02328940) + \
1128 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.023289 + (energy-1.000000)* -0.012327) + \
1129 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.010962 + (energy-2.000000)* -0.001777) + \
1130 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.005632 + (energy-5.000000)* -0.000436) + \
1131 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003452 + (energy-10.000000)* 0.000001) + \
1132 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003539 + (energy-100.000000)* 0.000001) + \
1133 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004876 + (energy-1000.000000)* 0.000002) + \
1134 (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 10000.0000) * (0.026660*energy/10000.000000) + \
1135 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.02303055) + \
1136 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.023031 + (energy-1.000000)* -0.011212) + \
1137 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.011819 + (energy-2.000000)* -0.002216) + \
1138 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.005171 + (energy-5.000000)* -0.000352) + \
1139 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003412 + (energy-10.000000)* 0.000001) + \
1140 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003516 + (energy-100.000000)* 0.000001) + \
1141 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004565 + (energy-1000.000000)* 0.000002) + \
1142 (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 10000.0000) * (0.025158*energy/10000.000000) + \
1143 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.02611889) + \
1144 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.026119 + (energy-1.000000)* -0.012835) + \
1145 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.013284 + (energy-2.000000)* -0.002589) + \
1146 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.005516 + (energy-5.000000)* -0.000388) + \
1147 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003574 + (energy-10.000000)* 0.000001) + \
1148 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003690 + (energy-100.000000)* 0.000001) + \
1149 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004493 + (energy-1000.000000)* 0.000002) + \
1150 (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 10000.0000) * (0.022856*energy/10000.000000) + \
1151 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.02991154) + \
1152 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.029912 + (energy-1.000000)* -0.015226) + \
1153 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.014686 + (energy-2.000000)* -0.002867) + \
1154 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.006086 + (energy-5.000000)* -0.000457) + \
1155 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003803 + (energy-10.000000)* 0.000001) + \
1156 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003932 + (energy-100.000000)* 0.000001) + \
1157 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004619 + (energy-1000.000000)* 0.000002) + \
1158 (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 10000.0000) * (0.020643*energy/10000.000000) + \
1159 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.03556765) + \
1160 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.035568 + (energy-1.000000)* -0.019079) + \
1161 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.016488 + (energy-2.000000)* -0.003242) + \
1162 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.006763 + (energy-5.000000)* -0.000526) + \
1163 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004133 + (energy-10.000000)* 0.000002) + \
1164 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004270 + (energy-100.000000)* 0.000001) + \
1165 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004991 + (energy-1000.000000)* 0.000002) + \
1166 (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 10000.0000) * (0.021331*energy/10000.000000) + \
1167 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.03803690) + \
1168 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.038037 + (energy-1.000000)* -0.020041) + \
1169 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.017996 + (energy-2.000000)* -0.003460) + \
1170 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.007617 + (energy-5.000000)* -0.000724) + \
1171 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003999 + (energy-10.000000)* 0.000002) + \
1172 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004159 + (energy-100.000000)* 0.000001) + \
1173 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005300 + (energy-1000.000000)* 0.000002) + \
1174 (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 10000.0000) * (0.023674*energy/10000.000000) + \
1175 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.03760924) + \
1176 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.037609 + (energy-1.000000)* -0.018328) + \
1177 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.019281 + (energy-2.000000)* -0.003682) + \
1178 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.008234 + (energy-5.000000)* -0.000772) + \
1179 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004372 + (energy-10.000000)* 0.000002) + \
1180 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004545 + (energy-100.000000)* 0.000001) + \
1181 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005786 + (energy-1000.000000)* 0.000002) + \
1182 (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 10000.0000) * (0.026924*energy/10000.000000) + \
1183 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.03763696) + \
1184 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.037637 + (energy-1.000000)* -0.016962) + \
1185 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.020675 + (energy-2.000000)* -0.003893) + \
1186 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.008995 + (energy-5.000000)* -0.000848) + \
1187 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004755 + (energy-10.000000)* 0.000002) + \
1188 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004949 + (energy-100.000000)* 0.000002) + \
1189 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.006336 + (energy-1000.000000)* 0.000003) + \
1190 (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 10000.0000) * (0.030130*energy/10000.000000) + \
1191 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.04368466) + \
1192 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.043685 + (energy-1.000000)* -0.020430) + \
1193 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.023254 + (energy-2.000000)* -0.004385) + \
1194 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.010100 + (energy-5.000000)* -0.001016) + \
1195 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.005020 + (energy-10.000000)* 0.000002) + \
1196 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005206 + (energy-100.000000)* 0.000002) + \
1197 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.006679 + (energy-1000.000000)* 0.000003) + \
1198 (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 10000.0000) * (0.033540*energy/10000.000000) + \
1199 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.05055352) + \
1200 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.050554 + (energy-1.000000)* -0.024554) + \
1201 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.025999 + (energy-2.000000)* -0.004964) + \
1202 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.011106 + (energy-5.000000)* -0.001134) + \
1203 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.005436 + (energy-10.000000)* 0.000002) + \
1204 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005582 + (energy-100.000000)* 0.000002) + \
1205 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.007047 + (energy-1000.000000)* 0.000003) + \
1206 (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 10000.0000) * (0.036533*energy/10000.000000) + \
1207 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.05790950) + \
1208 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.057910 + (energy-1.000000)* -0.028845) + \
1209 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.029065 + (energy-2.000000)* -0.005581) + \
1210 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.012322 + (energy-5.000000)* -0.001294) + \
1211 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.005853 + (energy-10.000000)* 0.000002) + \
1212 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006005 + (energy-100.000000)* 0.000002) + \
1213 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.007594 + (energy-1000.000000)* 0.000004) + \
1214 (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 10000.0000) * (0.040788*energy/10000.000000) + \
1215 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.06787795) + \
1216 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.067878 + (energy-1.000000)* -0.035271) + \
1217 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.032607 + (energy-2.000000)* -0.006346) + \
1218 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.013570 + (energy-5.000000)* -0.001453) + \
1219 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.006304 + (energy-10.000000)* 0.000002) + \
1220 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006489 + (energy-100.000000)* 0.000002) + \
1221 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008302 + (energy-1000.000000)* 0.000004) + \
1222 (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 10000.0000) * (0.045331*energy/10000.000000) + \
1223 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.08769401) + \
1224 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.087694 + (energy-1.000000)* -0.050535) + \
1225 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.037159 + (energy-2.000000)* -0.007293) + \
1226 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.015278 + (energy-5.000000)* -0.001678) + \
1227 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.006886 + (energy-10.000000)* 0.000003) + \
1228 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.007116 + (energy-100.000000)* 0.000002) + \
1229 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009229 + (energy-1000.000000)* 0.000005) + \
1230 (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 10000.0000) * (0.051460*energy/10000.000000) + \
1231 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.12378717) + \
1232 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.123787 + (energy-1.000000)* -0.081534) + \
1233 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.042253 + (energy-2.000000)* -0.008504) + \
1234 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.016742 + (energy-5.000000)* -0.001830) + \
1235 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.007595 + (energy-10.000000)* 0.000003) + \
1236 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.007850 + (energy-100.000000)* 0.000003) + \
1237 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010277 + (energy-1000.000000)* 0.000005) + \
1238 (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 10000.0000) * (0.059714*energy/10000.000000) + \
1239 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.14293546) + \
1240 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.142935 + (energy-1.000000)* -0.095632) + \
1241 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.047303 + (energy-2.000000)* -0.009550) + \
1242 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.018654 + (energy-5.000000)* -0.002035) + \
1243 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.008480 + (energy-10.000000)* 0.000003) + \
1244 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.008767 + (energy-100.000000)* 0.000003) + \
1245 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011459 + (energy-1000.000000)* 0.000006) + \
1246 (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 10000.0000) * (0.066014*energy/10000.000000) + \
1247 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.16175578) + \
1248 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.161756 + (energy-1.000000)* -0.108957) + \
1249 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.052799 + (energy-2.000000)* -0.010753) + \
1250 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.020541 + (energy-5.000000)* -0.002235) + \
1251 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.009366 + (energy-10.000000)* 0.000004) + \
1252 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.009686 + (energy-100.000000)* 0.000003) + \
1253 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012672 + (energy-1000.000000)* 0.000007) + \
1254 (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 10000.0000) * (0.072840*energy/10000.000000) + \
1255 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.17768991) + \
1256 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.177690 + (energy-1.000000)* -0.119322) + \
1257 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.058368 + (energy-2.000000)* -0.011801) + \
1258 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.022966 + (energy-5.000000)* -0.002532) + \
1259 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.010308 + (energy-10.000000)* 0.000004) + \
1260 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.010654 + (energy-100.000000)* 0.000004) + \
1261 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013982 + (energy-1000.000000)* 0.000007) + \
1262 (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 10000.0000) * (0.081143*energy/10000.000000) + \
1263 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.19156432) + \
1264 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.191564 + (energy-1.000000)* -0.125357) + \
1265 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.066207 + (energy-2.000000)* -0.013680) + \
1266 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.025168 + (energy-5.000000)* -0.002755) + \
1267 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.011394 + (energy-10.000000)* 0.000004) + \
1268 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.011796 + (energy-100.000000)* 0.000004) + \
1269 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015576 + (energy-1000.000000)* 0.000009) + \
1270 (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 10000.0000) * (0.092193*energy/10000.000000) + \
1271 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1272 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.075652) + \
1273 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.075652 + (energy-2.000000)* -0.015621) + \
1274 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.028790 + (energy-5.000000)* -0.003236) + \
1275 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.012611 + (energy-10.000000)* 0.000005) + \
1276 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.013037 + (energy-100.000000)* 0.000005) + \
1277 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017269 + (energy-1000.000000)* 0.000010) + \
1278 (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 10000.0000) * (0.103252*energy/10000.000000) + \
1279 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1280 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.088768) + \
1281 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.088768 + (energy-2.000000)* -0.018880) + \
1282 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.032126 + (energy-5.000000)* -0.003627) + \
1283 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.013990 + (energy-10.000000)* 0.000005) + \
1284 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.014468 + (energy-100.000000)* 0.000005) + \
1285 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019295 + (energy-1000.000000)* 0.000011) + \
1286 (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 10000.0000) * (0.117161*energy/10000.000000) + \
1287 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1288 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.099752) + \
1289 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.099752 + (energy-2.000000)* -0.021799) + \
1290 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.034356 + (energy-5.000000)* -0.003793) + \
1291 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.015389 + (energy-10.000000)* 0.000006) + \
1292 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.015906 + (energy-100.000000)* 0.000006) + \
1293 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.021315 + (energy-1000.000000)* 0.000012) + \
1294 (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 10000.0000) * (0.133463*energy/10000.000000) + \
1295 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1296 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.129119) + \
1297 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.129119 + (energy-2.000000)* -0.029895) + \
1298 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.039433 + (energy-5.000000)* -0.004396) + \
1299 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.017452 + (energy-10.000000)* 0.000006) + \
1300 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.018028 + (energy-100.000000)* 0.000007) + \
1301 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.024225 + (energy-1000.000000)* 0.000014) + \
1302 (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 10000.0000) * (0.151204*energy/10000.000000) + \
1303 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1304 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.148587) + \
1305 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.148587 + (energy-2.000000)* -0.035133) + \
1306 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.043189 + (energy-5.000000)* -0.004829) + \
1307 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.019045 + (energy-10.000000)* 0.000006) + \
1308 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.019625 + (energy-100.000000)* 0.000007) + \
1309 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.026335 + (energy-1000.000000)* 0.000015) + \
1310 (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 10000.0000) * (0.165614*energy/10000.000000) + \
1311 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1312 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.230194) + \
1313 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.230194 + (energy-2.000000)* -0.059949) + \
1314 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.050347 + (energy-5.000000)* -0.005793) + \
1315 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.021380 + (energy-10.000000)* 0.000007) + \
1316 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.022044 + (energy-100.000000)* 0.000009) + \
1317 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.029961 + (energy-1000.000000)* 0.000018) + \
1318 (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 10000.0000) * (0.188507*energy/10000.000000) + \
1319 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1320 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.256226) + \
1321 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.256226 + (energy-2.000000)* -0.067581) + \
1322 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.053485 + (energy-5.000000)* -0.006077) + \
1323 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.023102 + (energy-10.000000)* 0.000008) + \
1324 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.023822 + (energy-100.000000)* 0.000010) + \
1325 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.032699 + (energy-1000.000000)* 0.000020) + \
1326 (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 10000.0000) * (0.212432*energy/10000.000000) + \
1327 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1328 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.428061) + \
1329 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.428061 + (energy-2.000000)* -0.122293) + \
1330 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.061183 + (energy-5.000000)* -0.007029) + \
1331 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.026038 + (energy-10.000000)* 0.000009) + \
1332 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.026834 + (energy-100.000000)* 0.000011) + \
1333 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.037168 + (energy-1000.000000)* 0.000023) + \
1334 (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 10000.0000) * (0.240865*energy/10000.000000) + \
1335 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1336 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.302517) + \
1337 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.302517 + (energy-2.000000)* -0.078082) + \
1338 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.068271 + (energy-5.000000)* -0.007915) + \
1339 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.028698 + (energy-10.000000)* 0.000010) + \
1340 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.029589 + (energy-100.000000)* 0.000014) + \
1341 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.042319 + (energy-1000.000000)* 0.000027) + \
1342 (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 10000.0000) * (0.284210*energy/10000.000000) + \
1343 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1344 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1345 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.025424) + \
1346 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.076273 + (energy-5.000000)* -0.008946) + \
1347 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.031543 + (energy-10.000000)* 0.000010) + \
1348 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.032442 + (energy-100.000000)* 0.000017) + \
1349 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.047712 + (energy-1000.000000)* 0.000031) + \
1350 (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 10000.0000) * (0.326399*energy/10000.000000) + \
1351 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1352 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1353 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.029080) + \
1354 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.087241 + (energy-5.000000)* -0.010397) + \
1355 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.035256 + (energy-10.000000)* 0.000011) + \
1356 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.036252 + (energy-100.000000)* 0.000019) + \
1357 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.053456 + (energy-1000.000000)* 0.000034) + \
1358 (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 10000.0000) * (0.357698*energy/10000.000000) + \
1359 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1360 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1361 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.033344) + \
1362 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.100031 + (energy-5.000000)* -0.012242) + \
1363 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.038819 + (energy-10.000000)* 0.000015) + \
1364 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.040162 + (energy-100.000000)* 0.000023) + \
1365 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.061269 + (energy-1000.000000)* 0.000040) + \
1366 (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 10000.0000) * (0.420216*energy/10000.000000) + \
1367 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1368 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1369 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.039257) + \
1370 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.117770 + (energy-5.000000)* -0.014993) + \
1371 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.042808 + (energy-10.000000)* 0.000025) + \
1372 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.045061 + (energy-100.000000)* 0.000029) + \
1373 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.071538 + (energy-1000.000000)* 0.000049) + \
1374 (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 10000.0000) * (0.510223*energy/10000.000000) + \
1375 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1376 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1377 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.041967) + \
1378 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.125902 + (energy-5.000000)* -0.016273) + \
1379 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.044536 + (energy-10.000000)* 0.000023) + \
1380 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.046571 + (energy-100.000000)* 0.000043) + \
1381 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.085014 + (energy-1000.000000)* 0.000066) + \
1382 (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 10000.0000) * (0.680639*energy/10000.000000) + \
1383 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
1384 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
1385 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.000000) + \
1386 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.000000 + (energy-5.000000)* 0.013098) + \
1387 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.065492 + (energy-10.000000)* 0.000107) + \
1388 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.075097 + (energy-100.000000)* 0.000243) + \
1389 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.294198 + (energy-1000.000000)* 0.000280) + \
1390 (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 10000.0000) * (2.814894*energy/10000.000000)
1391 }
1392}
1393
1394
1395##############
1396# Track merger
1397##############
1398
1399module Merger TrackMerger {
1400 # add InputArray InputArray
1401 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
1402 add InputArray ElectronMomentumSmearing/electrons
1403 add InputArray MuonMomentumSmearing/muons
1404 set OutputArray tracks
1405}
1406
1407#############
1408# ECAL
1409#############
1410
1411module SimpleCalorimeter ECal {
1412 set ParticleInputArray ParticlePropagator/stableParticles
1413 set TrackInputArray TrackMerger/tracks
1414
1415 set TowerOutputArray ecalTowers
1416 set EFlowTrackOutputArray eflowTracks
1417 set EFlowTowerOutputArray eflowPhotons
1418
1419 set IsEcal true
1420
1421 set EnergyMin 0.5
1422 set EnergySignificanceMin 1.0
1423
1424 set SmearTowerCenter true
1425
1426 set pi [expr {acos(-1)}]
1427
1428 # lists of the edges of each tower in eta and phi
1429 # each list starts with the lower edge of the first tower
1430 # the list ends with the higher edged of the last tower
1431
1432 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
1433 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
1434
1435 #barrel:
1436 #dphi = 0.2 degree towers up to eta <=1.2
1437 set PhiBins {}
1438 for {set i -900} {$i <= 900} {incr i} {
1439 add PhiBins [expr {$i * $pi/900.0 }]
1440 }
1441 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
1442 for {set i -400} {$i <=400} {incr i} {
1443 set eta [expr {$i * 0.003}]
1444 add EtaPhiBins $eta $PhiBins
1445 }
1446
1447 #endcaps:
1448 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
1449 set PhiBins {}
1450 for {set i -225} {$i <= 225} {incr i} {
1451 add PhiBins [expr {$i * $pi/225.}]
1452 }
1453 #deta=0.02 units for 1.2 < |eta| <=2.5
1454 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
1455 for {set i 1} {$i <=66} {incr i} {
1456 set eta [expr {-2.52 + $i * 0.02}]
1457 add EtaPhiBins $eta $PhiBins
1458 }
1459 #same for 1.2 to 2.5
1460 for {set i 1} {$i <=66} {incr i} {
1461 set eta [expr {1.18 + $i*0.02}]
1462 add EtaPhiBins $eta $PhiBins
1463 }
1464
1465
1466 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
1467
1468 add EnergyFraction {0} {0.0}
1469 # energy fractions for e, gamma and pi0
1470 add EnergyFraction {11} {1.0}
1471 add EnergyFraction {22} {1.0}
1472 add EnergyFraction {111} {1.0}
1473 # energy fractions for muon, neutrinos and neutralinos
1474 add EnergyFraction {12} {0.0}
1475 add EnergyFraction {13} {0.0}
1476 add EnergyFraction {14} {0.0}
1477 add EnergyFraction {16} {0.0}
1478 add EnergyFraction {1000022} {0.0}
1479 add EnergyFraction {1000023} {0.0}
1480 add EnergyFraction {1000025} {0.0}
1481 add EnergyFraction {1000035} {0.0}
1482 add EnergyFraction {1000045} {0.0}
1483 # energy fractions for K0short and Lambda
1484 add EnergyFraction {310} {0.3}
1485 add EnergyFraction {3122} {0.3}
1486
1487 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
1488 set ResolutionFormula {
1489 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
1490 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
1491 (abs(eta) <= 2.5 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
1492}
1493
1494#############
1495# HCAL
1496#############
1497
1498module SimpleCalorimeter HCal {
1499 set ParticleInputArray ParticlePropagator/stableParticles
1500 set TrackInputArray ECal/eflowTracks
1501
1502 set TowerOutputArray hcalTowers
1503 set EFlowTrackOutputArray eflowTracks
1504 set EFlowTowerOutputArray eflowNeutralHadrons
1505
1506 set IsEcal false
1507
1508 set EnergyMin 1.0
1509 set EnergySignificanceMin 1.0
1510
1511 set SmearTowerCenter true
1512
1513 set pi [expr {acos(-1)}]
1514
1515 # lists of the edges of each tower in eta and phi
1516 # each list starts with the lower edge of the first tower
1517 # the list ends with the higher edged of the last tower
1518
1519
1520 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
1521 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
1522 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=2.5
1523 #HCAL cell sizes always 30x30 mm^2
1524
1525 #barrel and ring:
1526 #dphi = 1 degree up to |eta| <=0.9
1527 set PhiBins {}
1528 for {set i -180} {$i <=180} {incr i} {
1529 add PhiBins [expr {$i * $pi/180.0}]
1530 }
1531 #deta= 0.02 towers up to |eta| <=0.9
1532 for {set i -45} {$i <=45} {incr i} {
1533 set eta [expr {$i * 0.02}]
1534 add EtaPhiBins $eta $PhiBins
1535 }
1536
1537 #endcaps:
1538 # dphi = 6 degree
1539 set PhiBins {}
1540 for {set i -30} {$i <=30} {incr i} {
1541 add PhiBins [expr {$i * $pi/30.0}]
1542 }
1543 # deta =0.1 for 0.9 < |eta| <=2.5
1544 #for -2.5 to -0.9, 21 segments
1545 for {set i 1} {$i <=17} {incr i} {
1546 set eta [expr {-2.5 + $i * 0.1}]
1547 add EtaPhiBins $eta $PhiBins
1548 }
1549 #same for 0.9 to 2.5
1550 for {set i 1} {$i <=17} {incr i} {
1551 set eta [expr {0.8 + $i * 0.1 }]
1552 add EtaPhiBins $eta $PhiBins
1553 }
1554
1555 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
1556 add EnergyFraction {0} {1.0}
1557 # energy fractions for e, gamma and pi0
1558 add EnergyFraction {11} {0.0}
1559 add EnergyFraction {22} {0.0}
1560 add EnergyFraction {111} {0.0}
1561 # energy fractions for muon, neutrinos and neutralinos
1562 add EnergyFraction {12} {0.0}
1563 add EnergyFraction {13} {0.0}
1564 add EnergyFraction {14} {0.0}
1565 add EnergyFraction {16} {0.0}
1566 add EnergyFraction {1000022} {0.0}
1567 add EnergyFraction {1000023} {0.0}
1568 add EnergyFraction {1000025} {0.0}
1569 add EnergyFraction {1000035} {0.0}
1570 add EnergyFraction {1000045} {0.0}
1571 # energy fractions for K0short and Lambda
1572 add EnergyFraction {310} {0.7}
1573 add EnergyFraction {3122} {0.7}
1574
1575 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
1576 #CLICdet internal studies
1577 set ResolutionFormula {
1578 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
1579 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
1580 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
1581 (abs(eta)<=2.5 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
1582 }
1583
1584}
1585
1586#################
1587# Electron filter
1588#################
1589
1590module PdgCodeFilter ElectronFilter {
1591 set InputArray HCal/eflowTracks
1592 set OutputArray electrons
1593 set Invert true
1594 add PdgCode {11}
1595 add PdgCode {-11}
1596}
1597
1598######################
1599# ChargedHadronFilter
1600######################
1601
1602module PdgCodeFilter ChargedHadronFilter {
1603 set InputArray HCal/eflowTracks
1604 set OutputArray chargedHadrons
1605
1606 add PdgCode {11}
1607 add PdgCode {-11}
1608 add PdgCode {13}
1609 add PdgCode {-13}
1610}
1611
1612
1613
1614###################################################
1615# Tower Merger (in case not using e-flow algorithm)
1616###################################################
1617
1618module Merger Calorimeter {
1619 # add InputArray InputArray
1620 add InputArray ECal/ecalTowers
1621 add InputArray HCal/hcalTowers
1622 add InputArray MuonMomentumSmearing/muons
1623 set OutputArray towers
1624}
1625
1626
1627####################
1628# Energy flow merger
1629####################
1630
1631module Merger EFlowMerger {
1632 # add InputArray InputArray
1633 add InputArray HCal/eflowTracks
1634 add InputArray ECal/eflowPhotons
1635 add InputArray HCal/eflowNeutralHadrons
1636 set OutputArray eflow
1637}
1638
1639######################
1640# EFlowFilter
1641######################
1642
1643module PdgCodeFilter EFlowFilter {
1644 set InputArray EFlowMerger/eflow
1645 set OutputArray eflow
1646
1647 add PdgCode {11}
1648 add PdgCode {-11}
1649 add PdgCode {13}
1650 add PdgCode {-13}
1651}
1652
1653
1654###################
1655# Photon efficiency
1656###################
1657
1658module Efficiency PhotonEfficiency {
1659 set InputArray ECal/eflowPhotons
1660 set OutputArray photons
1661
1662 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1663
1664 # efficiency formula for photons
1665 set EfficiencyFormula {
1666 (energy < 2.0 ) * (0.000) +
1667 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
1668 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=2.5) * (0.9) }
1669
1670}
1671
1672
1673##################
1674# Photon isolation
1675##################
1676
1677module Isolation PhotonIsolation {
1678 set CandidateInputArray PhotonEfficiency/photons
1679 set IsolationInputArray EFlowMerger/eflow
1680
1681 set OutputArray photons
1682
1683 set DeltaRMax 0.1
1684
1685 set PTMin 0.5
1686
1687 set PTRatioMax 0.2
1688}
1689
1690
1691#####################
1692# Electron efficiency
1693#####################
1694
1695module Efficiency ElectronEfficiency {
1696 set InputArray ElectronFilter/electrons
1697 set OutputArray electrons
1698
1699 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1700
1701 set EfficiencyFormula {
1702 (energy < 3.0 ) * ( 0.00 ) +
1703 (abs(eta) > 2.50) * ( 0.00 ) +
1704 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.58 ) +
1705 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
1706 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
1707 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
1708 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
1709 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
1710 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.6 ) +
1711 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
1712 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
1713 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
1714 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
1715 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
1716 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.6 ) +
1717 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
1718 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
1719 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
1720 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
1721 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
1722 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.64 ) +
1723 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
1724 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
1725 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
1726 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
1727 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
1728 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.64 ) +
1729 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
1730 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
1731 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
1732 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
1733 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
1734 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.67 ) +
1735 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
1736 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
1737 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
1738 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
1739 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
1740 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1741 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
1742 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
1743 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
1744 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
1745 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
1746 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.7 ) +
1747 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
1748 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
1749 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
1750 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
1751 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
1752 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1753 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
1754 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
1755 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
1756 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
1757 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
1758 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.68 ) +
1759 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
1760 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
1761 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
1762 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
1763 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
1764 ( energy >=400 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1765 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
1766 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
1767 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
1768 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
1769 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
1770 }
1771}
1772
1773####################
1774# Electron isolation
1775####################
1776
1777module Isolation ElectronIsolation {
1778 set CandidateInputArray ElectronEfficiency/electrons
1779 set IsolationInputArray EFlowMerger/eflow
1780
1781 set OutputArray electrons
1782
1783 set DeltaRMax 0.1
1784
1785 set PTMin 0.5
1786
1787 set PTRatioMax 0.2
1788}
1789
1790#################
1791# Muon efficiency
1792#################
1793
1794module Efficiency MuonEfficiency {
1795 set InputArray MuonMomentumSmearing/muons
1796 set OutputArray muons
1797
1798 # set EfficiencyFormula {efficiency as a function of eta and pt}
1799
1800 set EfficiencyFormula {
1801 (energy < 2.5 ) * (0.00) +
1802 (energy>=2.5 ) * (0.999)
1803 }
1804}
1805
1806################
1807# Muon isolation
1808################
1809
1810module Isolation MuonIsolation {
1811 set CandidateInputArray MuonEfficiency/muons
1812 set IsolationInputArray EFlowMerger/eflow
1813
1814 set OutputArray muons
1815
1816 set DeltaRMax 0.1
1817
1818 set PTMin 0.5
1819
1820 set PTRatioMax 0.2
1821}
1822
1823
1824
1825###################
1826# Missing ET merger
1827###################
1828
1829module Merger MissingET {
1830 # add InputArray InputArray
1831 add InputArray EFlowMerger/eflow
1832 set MomentumOutputArray momentum
1833}
1834
1835
1836##################
1837# Scalar HT merger
1838##################
1839
1840module Merger ScalarHT {
1841 # add InputArray InputArray
1842 add InputArray EFlowMerger/eflow
1843 set EnergyOutputArray energy
1844}
1845######################
1846# EFlowFilter (UniqueObjectFinder)
1847######################
1848module UniqueObjectFinder UniqueObjectFinder {
1849 add InputArray PhotonIsolation/photons photons
1850 add InputArray ElectronIsolation/electrons electrons
1851 add InputArray MuonIsolation/muons muons
1852 add InputArray EFlowMerger/eflow eflow
1853}
1854
1855#################
1856# Neutrino Filter
1857#################
1858
1859module PdgCodeFilter NeutrinoFilter {
1860
1861 set InputArray Delphes/stableParticles
1862 set OutputArray filteredParticles
1863
1864 set PTMin 0.0
1865
1866 add PdgCode {12}
1867 add PdgCode {14}
1868 add PdgCode {16}
1869 add PdgCode {-12}
1870 add PdgCode {-14}
1871 add PdgCode {-16}
1872
1873}
1874
1875
1876#####################
1877# MC truth jet finder
1878#####################
1879
1880module FastJetFinder GenJetFinder {
1881 set InputArray NeutrinoFilter/filteredParticles
1882
1883 set OutputArray jets
1884
1885 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1886 set JetAlgorithm 9
1887 set ParameterR 0.5
1888
1889 set JetPTMin 20.0
1890}
1891
1892#########################
1893# Gen Missing ET merger
1894########################
1895
1896module Merger GenMissingET {
1897 # add InputArray InputArray
1898 add InputArray NeutrinoFilter/filteredParticles
1899 set MomentumOutputArray momentum
1900}
1901
1902
1903
1904############
1905# Jet finder
1906############
1907
1908module FastJetFinder FastJetFinderKt {
1909 # set InputArray Calorimeter/towers
1910 set InputArray EFlowMerger/eflow
1911
1912 set OutputArray KTjets
1913
1914 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1915 set JetAlgorithm 4
1916 set ParameterR 0.5
1917
1918 set JetPTMin 20.0
1919}
1920
1921################
1922# Jet finder VLC
1923################
1924
1925
1926source MuonCollider/MuonColliderDet_JetReco.tcl
1927
1928
1929#########################################
1930# Jet Momentum Smearing to mimick overlay
1931#########################################
1932
1933
1934source MuonCollider/MuonColliderDet_JetSmearing.tcl
1935
1936
1937
1938########################
1939# Jet Flavor Association
1940########################
1941
1942source MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl
1943
1944###########
1945# b-tagging
1946###########
1947# based on CLICdp-Note-2014-002
1948
1949source MuonCollider/MuonColliderDet_BTagging.tcl
1950
1951
1952#############
1953# tau-tagging
1954#############
1955
1956source MuonCollider/MuonColliderDet_TauTagging.tcl
1957
1958##################
1959# ROOT tree writer
1960##################
1961
1962module TreeWriter TreeWriter {
1963 # add Branch InputArray BranchName BranchClass
1964 add Branch Delphes/allParticles Particle GenParticle
1965
1966 add Branch GenJetFinder/jets GenJet Jet
1967
1968 add Branch FastJetFinderKt/KTjets KTjet Jet
1969
1970 add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 VLCjetR05N2 Jet
1971 add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 VLCjetR05N3 Jet
1972 add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 VLCjetR05N4 Jet
1973 add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 VLCjetR05N5 Jet
1974 add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 VLCjetR05N6 Jet
1975
1976 add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 VLCjetR07N2 Jet
1977 add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 VLCjetR07N3 Jet
1978 add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 VLCjetR07N4 Jet
1979 add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 VLCjetR07N5 Jet
1980 add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 VLCjetR07N6 Jet
1981
1982 add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 VLCjetR10N2 Jet
1983 add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 VLCjetR10N3 Jet
1984 add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 VLCjetR10N4 Jet
1985 add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 VLCjetR10N5 Jet
1986 add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 VLCjetR10N6 Jet
1987
1988 add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 VLCjetR12N2 Jet
1989 add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 VLCjetR12N3 Jet
1990 add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 VLCjetR12N4 Jet
1991 add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 VLCjetR12N5 Jet
1992 add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 VLCjetR12N6 Jet
1993
1994 add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 VLCjetR15N2 Jet
1995 add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 VLCjetR15N3 Jet
1996 add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 VLCjetR15N4 Jet
1997 add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 VLCjetR15N5 Jet
1998 add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 VLCjetR15N6 Jet
1999
2000 add Branch JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive VLCjetR02_inclusive Jet
2001 add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive VLCjetR05_inclusive Jet
2002 add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive VLCjetR07_inclusive Jet
2003 add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive VLCjetR10_inclusive Jet
2004 add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive VLCjetR12_inclusive Jet
2005 add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive VLCjetR15_inclusive Jet
2006
2007
2008 ####
2009
2010 add Branch GenMissingET/momentum GenMissingET MissingET
2011
2012 add Branch TrackMerger/tracks Track Track
2013 add Branch Calorimeter/towers Tower Tower
2014
2015 add Branch HCal/eflowTracks EFlowTrack Track
2016 add Branch ECal/eflowPhotons EFlowPhoton Tower
2017 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
2018
2019 add Branch UniqueObjectFinder/photons Photon Photon
2020 add Branch UniqueObjectFinder/electrons Electron Electron
2021 add Branch UniqueObjectFinder/muons Muon Muon
2022 add Branch ForwardMuonMomentumSmearing/muons ForwardMuon Muon
2023
2024 add Branch MissingET/momentum MissingET MissingET
2025 add Branch ScalarHT/energy ScalarHT ScalarHT
2026}
Note: See TracBrowser for help on using the repository browser.