Fork me on GitHub

source: git/cards/delphes_card_MuonColliderDet.tcl@ b8a6aa3

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

fixed forward muon collection
  • Property mode set to 100644
File size: 98.9 KB
RevLine 
[a02a49e]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.
[4b3df0d]13# Hybrid between FCC-hh and CLIC performance.
[a02a49e]14#
15#
16#######################################
17
18#######################################
19# Order of execution of various modules
20#######################################
21
22set ExecutionPath {
23 ParticlePropagator
24 TrackMergerProp
25
26 DenseProp
27 DenseMergeTracks
28 DenseTrackFilter
29
30 ChargedHadronTrackingEfficiency
31 ElectronTrackingEfficiency
32 MuonTrackingEfficiency
[e01b141]33 ForwardMuonEfficiency
[a02a49e]34
35 ChargedHadronMomentumSmearing
36 ElectronMomentumSmearing
37 MuonMomentumSmearing
[e01b141]38 ForwardMuonMomentumSmearing
[a02a49e]39
40 TrackMerger
41
42 ECal
43 HCal
44
45 Calorimeter
46 EFlowMerger
47 EFlowFilter
48
49 PhotonEfficiency
50 PhotonIsolation
51
52 ElectronFilter
53 ElectronEfficiency
54 ElectronIsolation
55
56 ChargedHadronFilter
57
58 MuonEfficiency
59 MuonIsolation
60
61 UniqueObjectFinder
62
63 NeutrinoFilter
64 GenJetFinder
65
66
67 FastJetFinderKt
68 FastJetFinderVLC_R02_N2
69 FastJetFinderVLC_R02_N3
70 FastJetFinderVLC_R02_N4
71 FastJetFinderVLC_R02_N5
72 FastJetFinderVLC_R02_N6
73 FastJetFinderVLC_R05_N2
74 FastJetFinderVLC_R05_N3
75 FastJetFinderVLC_R05_N4
76 FastJetFinderVLC_R05_N5
77 FastJetFinderVLC_R05_N6
78 FastJetFinderVLC_R07_N2
79 FastJetFinderVLC_R07_N3
80 FastJetFinderVLC_R07_N4
81 FastJetFinderVLC_R07_N5
82 FastJetFinderVLC_R07_N6
83 FastJetFinderVLC_R10_N2
84 FastJetFinderVLC_R10_N3
85 FastJetFinderVLC_R10_N4
86 FastJetFinderVLC_R10_N5
87 FastJetFinderVLC_R10_N6
88 FastJetFinderVLC_R12_N2
89 FastJetFinderVLC_R12_N3
90 FastJetFinderVLC_R12_N4
91 FastJetFinderVLC_R12_N5
92 FastJetFinderVLC_R12_N6
93 FastJetFinderVLC_R15_N2
94 FastJetFinderVLC_R15_N3
95 FastJetFinderVLC_R15_N4
96 FastJetFinderVLC_R15_N5
97 FastJetFinderVLC_R15_N6
98
99 FastJetFinderVLC_R05_inclusive
100 FastJetFinderVLC_R02_inclusive
101 FastJetFinderVLC_R07_inclusive
102 FastJetFinderVLC_R10_inclusive
103 FastJetFinderVLC_R12_inclusive
104 FastJetFinderVLC_R15_inclusive
105
106 MissingET
107 GenMissingET
108
109 JetMomentumSmearing_VLCR02N2
110 JetMomentumSmearing_VLCR02N3
111 JetMomentumSmearing_VLCR02N4
112 JetMomentumSmearing_VLCR02N5
113 JetMomentumSmearing_VLCR02N6
114 JetMomentumSmearing_VLCR02_inclusive
115 JetMomentumSmearing_VLCR05N2
116 JetMomentumSmearing_VLCR05N3
117 JetMomentumSmearing_VLCR05N4
118 JetMomentumSmearing_VLCR05N5
119 JetMomentumSmearing_VLCR05N6
120 JetMomentumSmearing_VLCR05_inclusive
121 JetMomentumSmearing_VLCR07N2
122 JetMomentumSmearing_VLCR07N3
123 JetMomentumSmearing_VLCR07N4
124 JetMomentumSmearing_VLCR07N5
125 JetMomentumSmearing_VLCR07N6
126 JetMomentumSmearing_VLCR07_inclusive
127 JetMomentumSmearing_VLCR10N2
128 JetMomentumSmearing_VLCR10N3
129 JetMomentumSmearing_VLCR10N4
130 JetMomentumSmearing_VLCR10N5
131 JetMomentumSmearing_VLCR10N6
132 JetMomentumSmearing_VLCR10_inclusive
133 JetMomentumSmearing_VLCR12N2
134 JetMomentumSmearing_VLCR12N3
135 JetMomentumSmearing_VLCR12N4
136 JetMomentumSmearing_VLCR12N5
137 JetMomentumSmearing_VLCR12N6
138 JetMomentumSmearing_VLCR12_inclusive
139 JetMomentumSmearing_VLCR15N2
140 JetMomentumSmearing_VLCR15N3
141 JetMomentumSmearing_VLCR15N4
142 JetMomentumSmearing_VLCR15N5
143 JetMomentumSmearing_VLCR15N6
144 JetMomentumSmearing_VLCR15_inclusive
145
146
147 JetFlavorAssociation_R02N2
148 JetFlavorAssociation_R02N3
149 JetFlavorAssociation_R02N4
150 JetFlavorAssociation_R02N5
151 JetFlavorAssociation_R02N6
152 JetFlavorAssociation_R05N2
153 JetFlavorAssociation_R05N3
154 JetFlavorAssociation_R05N4
155 JetFlavorAssociation_R05N5
156 JetFlavorAssociation_R05N6
157
158 JetFlavorAssociation_R07N2
159 JetFlavorAssociation_R07N3
160 JetFlavorAssociation_R07N4
161 JetFlavorAssociation_R07N5
162 JetFlavorAssociation_R07N6
163
164 JetFlavorAssociation_R10N2
165 JetFlavorAssociation_R10N3
166 JetFlavorAssociation_R10N4
167 JetFlavorAssociation_R10N5
168 JetFlavorAssociation_R10N6
169
170 JetFlavorAssociation_R12N2
171 JetFlavorAssociation_R12N3
172 JetFlavorAssociation_R12N4
173 JetFlavorAssociation_R12N5
174 JetFlavorAssociation_R12N6
175
176 JetFlavorAssociation_R15N2
177 JetFlavorAssociation_R15N3
178 JetFlavorAssociation_R15N4
179 JetFlavorAssociation_R15N5
180 JetFlavorAssociation_R15N6
181
182 JetFlavorAssociation_R02_inclusive
183 JetFlavorAssociation_R05_inclusive
184 JetFlavorAssociation_R07_inclusive
185 JetFlavorAssociation_R10_inclusive
186 JetFlavorAssociation_R12_inclusive
187 JetFlavorAssociation_R15_inclusive
188
189
190 BTagging_WP50_R02N2
191 BTagging_WP70_R02N2
192 BTagging_WP90_R02N2
193 BTagging_WP50_R02N3
194 BTagging_WP70_R02N3
195 BTagging_WP90_R02N3
196 BTagging_WP50_R02N4
197 BTagging_WP70_R02N4
198 BTagging_WP90_R02N4
199 BTagging_WP50_R02N5
200 BTagging_WP70_R02N5
201 BTagging_WP90_R02N5
202 BTagging_WP50_R02N6
203 BTagging_WP70_R02N6
204 BTagging_WP90_R02N6
205 BTagging_WP50_R05N2
206 BTagging_WP70_R05N2
207 BTagging_WP90_R05N2
208 BTagging_WP50_R05N3
209 BTagging_WP70_R05N3
210 BTagging_WP90_R05N3
211 BTagging_WP50_R05N4
212 BTagging_WP70_R05N4
213 BTagging_WP90_R05N4
214 BTagging_WP50_R05N5
215 BTagging_WP70_R05N5
216 BTagging_WP90_R05N5
217 BTagging_WP50_R05N6
218 BTagging_WP70_R05N6
219 BTagging_WP90_R05N6
220 BTagging_WP50_R07N2
221 BTagging_WP70_R07N2
222 BTagging_WP90_R07N2
223 BTagging_WP50_R07N3
224 BTagging_WP70_R07N3
225 BTagging_WP90_R07N3
226 BTagging_WP50_R07N4
227 BTagging_WP70_R07N4
228 BTagging_WP90_R07N4
229 BTagging_WP50_R07N5
230 BTagging_WP70_R07N5
231 BTagging_WP90_R07N5
232 BTagging_WP50_R07N6
233 BTagging_WP70_R07N6
234 BTagging_WP90_R07N6
235 BTagging_WP50_R10N2
236 BTagging_WP70_R10N2
237 BTagging_WP90_R10N2
238 BTagging_WP50_R10N3
239 BTagging_WP70_R10N3
240 BTagging_WP90_R10N3
241 BTagging_WP50_R10N4
242 BTagging_WP70_R10N4
243 BTagging_WP90_R10N4
244 BTagging_WP50_R10N5
245 BTagging_WP70_R10N5
246 BTagging_WP90_R10N5
247 BTagging_WP50_R10N6
248 BTagging_WP70_R10N6
249 BTagging_WP90_R10N6
250 BTagging_WP50_R12N2
251 BTagging_WP70_R12N2
252 BTagging_WP90_R12N2
253 BTagging_WP50_R12N3
254 BTagging_WP70_R12N3
255 BTagging_WP90_R12N3
256 BTagging_WP50_R12N4
257 BTagging_WP70_R12N4
258 BTagging_WP90_R12N4
259 BTagging_WP50_R12N5
260 BTagging_WP70_R12N5
261 BTagging_WP90_R12N5
262 BTagging_WP50_R12N6
263 BTagging_WP70_R12N6
264 BTagging_WP90_R12N6
265 BTagging_WP50_R15N2
266 BTagging_WP70_R15N2
267 BTagging_WP90_R15N2
268 BTagging_WP50_R15N3
269 BTagging_WP70_R15N3
270 BTagging_WP90_R15N3
271 BTagging_WP50_R15N4
272 BTagging_WP70_R15N4
273 BTagging_WP90_R15N4
274 BTagging_WP50_R15N5
275 BTagging_WP70_R15N5
276 BTagging_WP90_R15N5
277 BTagging_WP50_R15N6
278 BTagging_WP70_R15N6
279 BTagging_WP90_R15N6
280 BTagging_WP50_R02_inclusive
281 BTagging_WP70_R02_inclusive
282 BTagging_WP90_R02_inclusive
283 BTagging_WP50_R05_inclusive
284 BTagging_WP70_R05_inclusive
285 BTagging_WP90_R05_inclusive
286 BTagging_WP50_R07_inclusive
287 BTagging_WP70_R07_inclusive
288 BTagging_WP90_R07_inclusive
289 BTagging_WP50_R10_inclusive
290 BTagging_WP70_R10_inclusive
291 BTagging_WP90_R10_inclusive
292 BTagging_WP50_R12_inclusive
293 BTagging_WP70_R12_inclusive
294 BTagging_WP90_R12_inclusive
295 BTagging_WP50_R15_inclusive
296 BTagging_WP70_R15_inclusive
297 BTagging_WP90_R15_inclusive
298
299
300 TauTagging_R02N2
301 TauTagging_R02N3
302 TauTagging_R02N4
303 TauTagging_R02N5
304 TauTagging_R02N6
305 TauTagging_R05N2
306 TauTagging_R05N3
307 TauTagging_R05N4
308 TauTagging_R05N5
309 TauTagging_R05N6
310 TauTagging_R07N2
311 TauTagging_R07N3
312 TauTagging_R07N4
313 TauTagging_R07N5
314 TauTagging_R07N6
315 TauTagging_R10N2
316 TauTagging_R10N3
317 TauTagging_R10N4
318 TauTagging_R10N5
319 TauTagging_R10N6
320 TauTagging_R12N2
321 TauTagging_R12N3
322 TauTagging_R12N4
323 TauTagging_R12N5
324 TauTagging_R12N6
325 TauTagging_R15N2
326 TauTagging_R15N3
327 TauTagging_R15N4
328 TauTagging_R15N5
329 TauTagging_R15N6
330 TauTagging_R02_inclusive
331 TauTagging_R05_inclusive
332 TauTagging_R07_inclusive
333 TauTagging_R10_inclusive
334 TauTagging_R12_inclusive
335 TauTagging_R15_inclusive
336
337 ScalarHT
338
339 TreeWriter
340}
341
342#################################
343# Propagate particles in cylinder
344#################################
345
346module ParticlePropagator ParticlePropagator {
347 set InputArray Delphes/stableParticles
348
349 set OutputArray stableParticles
350 set ChargedHadronOutputArray chargedHadrons
351 set ElectronOutputArray electrons
352 set MuonOutputArray muons
353
354 # radius of the magnetic field coverage in the calorimeter, in m
355 set Radius 1.5
356 # half-length of the magnetic field coverage in the calorimeter, in m
357 set HalfLength 2.31
358
359 # magnetic field, in T
360 set Bz 4.0
361}
362
363
364##############
365# Track merger
366##############
367
368module Merger TrackMergerProp {
369# add InputArray InputArray
370 add InputArray ParticlePropagator/chargedHadrons
371 add InputArray ParticlePropagator/electrons
372 add InputArray ParticlePropagator/muons
373 set OutputArray tracks
374}
375
376
377
378####################################
379# Track propagation to pseudo-pixel
380####################################
381
382module ParticlePropagator DenseProp {
383
384 set InputArray TrackMergerProp/tracks
385
386 # radius of the magnetic field coverage, in m
387 set Radius 0.45
388 set RadiusMax 1.5
389 # half-length of the magnetic field coverage, in m
390 set HalfLength 0.8
391 set HalfLengthMax 2.31
392
393 # magnetic field
394 set Bz 4.0
395}
396
397#####################
398# Dense Track merger
399#####################
400
401module Merger DenseMergeTracks {
402# add InputArray InputArray
403 add InputArray DenseProp/chargedHadrons
404 add InputArray DenseProp/electrons
405 add InputArray DenseProp/muons
406 set OutputArray tracks
407}
408
409
410######################
411# Dense Track Filter
412######################
413
414module DenseTrackFilter DenseTrackFilter {
415
416 set TrackInputArray DenseMergeTracks/tracks
417
418 set TrackOutputArray tracks
419 set ChargedHadronOutputArray chargedHadrons
420 set ElectronOutputArray electrons
421 set MuonOutputArray muons
422
[e69eea2]423 set EtaPhiRes 0.003
[a02a49e]424 set EtaMax 2.5
425
426 set pi [expr {acos(-1)}]
427
428 set nbins_phi [expr {$pi/$EtaPhiRes} ]
429 set nbins_phi [expr {int($nbins_phi)} ]
430
431 set PhiBins {}
432 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
433 add PhiBins [expr {$i * $pi/$nbins_phi}]
434 }
435
436 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
437 set nbins_eta [expr {int($nbins_eta)} ]
438
439 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
440 set eta [expr {$i * $EtaPhiRes}]
441 add EtaPhiBins $eta $PhiBins
442 }
443}
444
445
446
447
448####################################
449# Charged hadron tracking efficiency
450####################################
451
452module Efficiency ChargedHadronTrackingEfficiency {
453 set InputArray DenseTrackFilter/chargedHadrons
454 set OutputArray chargedHadrons
455 # tracking efficiency formula for charged hadrons
456
457 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
458 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
459 (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
460 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
461 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
462 (abs(eta) > 2.5 ) * (0.00)
463 }
464}
465
466##############################
467# Electron tracking efficiency
468##############################
469
470module Efficiency ElectronTrackingEfficiency {
471 set InputArray DenseTrackFilter/electrons
472 set OutputArray electrons
473
474 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
475 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
476 (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
477 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
478 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
479 (abs(eta) > 2.5 ) * (0.00)
480 }
481}
482
483##########################
484# Muon tracking efficiency
485##########################
486
487module Efficiency MuonTrackingEfficiency {
488 set InputArray DenseTrackFilter/muons
489 set OutputArray muons
490
491 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
492 (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.95) +
493 (abs(eta) <= 2.0) * (pt > 1) * (0.99) +
494 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.90) +
495 (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.95) +
496 (abs(eta) > 2.5 ) * (0.00)
497 }
498}
499
[e01b141]500##########################
501# Forward Muon efficiency
502##########################
503
504## hypothetical forward muon spectrometer
505module Efficiency ForwardMuonEfficiency {
506 set InputArray ParticlePropagator/muons
507 set OutputArray muons
508
509 set EfficiencyFormula { (pt <= 0.5) * (0.00) +
[55f821a]510 (abs(eta) < 2.5 ) * (0.00) +
511 (abs(eta) > 2.5 && abs(eta) < 8.0) * (pt > 0.5 && pt <= 1) * (0.90) +
512 (abs(eta) > 2.5 && abs(eta) < 8.0) * (pt > 1.0) * (0.95) +
513 (abs(eta) > 8.0 ) * (0.00)
[e01b141]514 }
515}
516
517
[a02a49e]518########################################
519# Momentum resolution for charged tracks
520########################################
521
522module MomentumSmearing ChargedHadronMomentumSmearing {
523 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
524 set OutputArray chargedHadrons
525
526
[e01b141]527 # Resolution given in dpT/pT (from FCC-hh)
528 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
529 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
530 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
531 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
532 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
533 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
534 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
535 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
536 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
537 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
538 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
539 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
540 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
541 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
542 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
543 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
544 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
545 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
546 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
547 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
548 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
549 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
550 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
551 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
552 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
553 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
554 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
555 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
556 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
557 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
558 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
559 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
560 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
561 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
562 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
563 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
564 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
565 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
566 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
567 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
568 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
569 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
570 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
571 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
572 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
573 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
574 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
575 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
576 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
577 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
578 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
579 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
580 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
581 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
582 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
583 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
584 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
585 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
586 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
587 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
588 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
589 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
590 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
591 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
592 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
593 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
594 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
595 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
596 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
597 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
598 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
599 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
600 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
601 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
602 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
603 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
604 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
605 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
606 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
607 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
608 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
609 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
610 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
611 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
612 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
613 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
614 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
615 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
616 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
617 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
618 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
619 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
620 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
621 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
622 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
623 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
624 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
625 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
626 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
627 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
628 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
629 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
630 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
631 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
632 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
633 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
634 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
635 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
636 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
637 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
638 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
639 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
640 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
641 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
642 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
643 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
644 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
645 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
646 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
647 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
648 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
649 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
650 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
651 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
652 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
653 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
654 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
655 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
656 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
657 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
658 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
659 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
660 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
661 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
662 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
663 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
664 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
665 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
666 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
667 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
668 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
669 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
670 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
671 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
672 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
673 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
674 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
675 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
676 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
677 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
678 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
679 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
680 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
681 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
682 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
683 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
684 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
685 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
686 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
687 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
688 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
689 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
690 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
691 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
692 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
693 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
694 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
695 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
696 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
697 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
698 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
699 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
700 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
701 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
702 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
703 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
704 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
705 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
706 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
707 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
708 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
709 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
710 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
711 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
712 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
713 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
714 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
715 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
716 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
717 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
718 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
719 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
720 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
721 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
722 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
723 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
724 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
725 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
726 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
727 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
[a02a49e]728 }
729}
730
731###################################
732# Momentum resolution for electrons
733###################################
734
735module MomentumSmearing ElectronMomentumSmearing {
736 set InputArray ElectronTrackingEfficiency/electrons
737 set OutputArray electrons
738
[e01b141]739 # Resolution given in dpT/pT (from FCC-hh)
740 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
741 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
742 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
743 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
744 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
745 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
746 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
747 (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
748 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
749 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
750 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
751 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
752 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
753 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
754 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
755 (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
756 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
757 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
758 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
759 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
760 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
761 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
762 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
763 (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
764 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
765 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
766 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
767 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
768 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
769 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
770 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
771 (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
772 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
773 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
774 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
775 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
776 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
777 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
778 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
779 (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
780 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
781 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
782 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
783 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
784 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
785 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
786 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
787 (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
788 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
789 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
790 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
791 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
792 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
793 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
794 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
795 (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
796 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
797 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
798 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
799 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
800 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
801 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
802 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
803 (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
804 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
805 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
806 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
807 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
808 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
809 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
810 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
811 (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
812 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
813 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
814 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
815 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
816 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
817 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
818 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
819 (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
820 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
821 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
822 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
823 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
824 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
825 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
826 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
827 (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
828 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
829 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
830 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
831 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
832 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
833 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
834 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
835 (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
836 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
837 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
838 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
839 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
840 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
841 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
842 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
843 (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
844 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
845 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
846 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
847 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
848 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
849 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
850 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
851 (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
852 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
853 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
854 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
855 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
856 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
857 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
858 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
859 (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
860 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
861 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
862 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
863 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
864 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
865 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
866 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
867 (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
868 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
869 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
870 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
871 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
872 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
873 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
874 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
875 (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
876 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
877 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
878 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
879 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
880 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
881 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
882 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
883 (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
884 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
885 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
886 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
887 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
888 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
889 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
890 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
891 (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
892 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
893 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
894 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
895 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
896 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
897 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
898 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
899 (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
900 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
901 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
902 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
903 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
904 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
905 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
906 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
907 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
908 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
909 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
910 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
911 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
912 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
913 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
914 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
915 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
916 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
917 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
918 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
919 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
920 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
921 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
922 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
923 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
924 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
925 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
926 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
927 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
928 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
929 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
930 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
931 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
932 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
933 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
934 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
935 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
936 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
937 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
938 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
939 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
[a02a49e]940 }
941}
942
943###############################
944# Momentum resolution for muons
945###############################
946
947module MomentumSmearing MuonMomentumSmearing {
948 set InputArray MuonTrackingEfficiency/muons
949 set OutputArray muons
950
951 # Resolution given in dpT/pT.
952
953 set ResolutionFormula {
954
[e01b141]955 ( abs(eta) < 1.35 ) *
956
957 (sqrt(0.0000364164 + (
[4b3df0d]958 3*9.06262e-8 *pt^2* cosh(
[e01b141]959 eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + (
960 504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/(
961 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
962 sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/(
963 0.00516429/sin(2*atan(exp(-abs(eta))))^2 + (
[4b3df0d]964 96868.8 *(1/400000000 + 5*(0.117945 * 1/cosh(eta)^2)/(
[e01b141]965 pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
966 sin(2*atan(exp(-abs(eta))))^2))
967
968 ) +
969
970 ( abs(eta) > 1.35 && abs(eta) < 2.00) *
971
972 ( 1.53846*(
973
974 (
975 (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
976 (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
977 (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
978 (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
979 (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
980 (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
981 (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
982 (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
983 (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
984 (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
985 (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
986 (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
987 (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
988 (energy >= 10000.0000) * (0.062437*energy/10000.000000)
989 )
990
991
992 - 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) +
993
994 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)))
995
996 - 2.07692*
997 (
998
999 (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
1000 (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
1001 (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
1002 (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
1003 (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
1004 (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
1005 (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
1006 (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
1007 (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
1008 (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
1009 (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
1010 (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
1011 (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
1012 (energy >= 10000.0000) * (0.062437*energy/10000.000000)
1013
1014 )
1015
1016
1017 ) +
1018
1019 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01062416) + \
1020 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010624 + (energy-1.000000)* -0.005532) + \
1021 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005092 + (energy-2.000000)* -0.000326) + \
1022 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004115 + (energy-5.000000)* 0.000006) + \
1023 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 20.0000) * (0.004143 + (energy-10.000000)* 0.000007) + \
1024 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 20.0000 && energy < 50.0000) * (0.004209 + (energy-20.000000)* 0.000007) + \
1025 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 50.0000 && energy < 100.0000) * (0.004413 + (energy-50.000000)* 0.000005) + \
1026 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 200.0000) * (0.004681 + (energy-100.000000)* 0.000004) + \
1027 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 200.0000 && energy < 500.0000) * (0.005076 + (energy-200.000000)* 0.000004) + \
1028 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006270 + (energy-500.000000)* 0.000005) + \
1029 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.008960 + (energy-1000.000000)* 0.000007) + \
1030 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.015510 + (energy-2000.000000)* 0.000007) + \
1031 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.036867 + (energy-5000.000000)* 0.000007) + \
1032 (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073168*energy/10000.000000) + \
1033 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01007098) + \
1034 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010071 + (energy-1.000000)* -0.004627) + \
1035 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000322) + \
1036 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004478 + (energy-5.000000)* 0.000005) + \
1037 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 20.0000) * (0.004501 + (energy-10.000000)* 0.000006) + \
1038 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 20.0000 && energy < 50.0000) * (0.004558 + (energy-20.000000)* 0.000005) + \
1039 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 50.0000 && energy < 100.0000) * (0.004701 + (energy-50.000000)* 0.000004) + \
1040 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 200.0000) * (0.004888 + (energy-100.000000)* 0.000003) + \
1041 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 200.0000 && energy < 500.0000) * (0.005213 + (energy-200.000000)* 0.000004) + \
1042 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006429 + (energy-500.000000)* 0.000006) + \
1043 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.009343 + (energy-1000.000000)* 0.000007) + \
1044 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.016410 + (energy-2000.000000)* 0.000008) + \
1045 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.039265 + (energy-5000.000000)* 0.000008) + \
1046 (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078014*energy/10000.000000) + \
1047 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01095892) + \
1048 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010959 + (energy-1.000000)* -0.005458) + \
1049 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005501 + (energy-2.000000)* -0.000281) + \
1050 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004660 + (energy-5.000000)* 0.000005) + \
1051 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 20.0000) * (0.004686 + (energy-10.000000)* 0.000007) + \
1052 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 20.0000 && energy < 50.0000) * (0.004757 + (energy-20.000000)* 0.000006) + \
1053 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 50.0000 && energy < 100.0000) * (0.004937 + (energy-50.000000)* 0.000004) + \
1054 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 200.0000) * (0.005143 + (energy-100.000000)* 0.000004) + \
1055 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 200.0000 && energy < 500.0000) * (0.005505 + (energy-200.000000)* 0.000005) + \
1056 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006975 + (energy-500.000000)* 0.000007) + \
1057 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.010462 + (energy-1000.000000)* 0.000008) + \
1058 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.018731 + (energy-2000.000000)* 0.000009) + \
1059 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.045159 + (energy-5000.000000)* 0.000009) + \
1060 (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.089830*energy/10000.000000) + \
1061 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01279214) + \
1062 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012792 + (energy-1.000000)* -0.005763) + \
1063 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.007029 + (energy-2.000000)* -0.000513) + \
1064 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005489 + (energy-5.000000)* 0.000003) + \
1065 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 20.0000) * (0.005503 + (energy-10.000000)* 0.000003) + \
1066 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 20.0000 && energy < 50.0000) * (0.005537 + (energy-20.000000)* 0.000003) + \
1067 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 50.0000 && energy < 100.0000) * (0.005636 + (energy-50.000000)* 0.000003) + \
1068 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 200.0000) * (0.005773 + (energy-100.000000)* 0.000003) + \
1069 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 200.0000 && energy < 500.0000) * (0.006060 + (energy-200.000000)* 0.000005) + \
1070 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 500.0000 && energy < 1000.0000) * (0.007489 + (energy-500.000000)* 0.000007) + \
1071 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.011117 + (energy-1000.000000)* 0.000009) + \
1072 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.019824 + (energy-2000.000000)* 0.000009) + \
1073 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.047732 + (energy-5000.000000)* 0.000009) + \
1074 (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.094931*energy/10000.000000) + \
1075 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01502671) + \
1076 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015027 + (energy-1.000000)* -0.007177) + \
1077 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007850 + (energy-2.000000)* -0.000651) + \
1078 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005898 + (energy-5.000000)* 0.000003) + \
1079 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 20.0000) * (0.005913 + (energy-10.000000)* 0.000003) + \
1080 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 20.0000 && energy < 50.0000) * (0.005947 + (energy-20.000000)* 0.000003) + \
1081 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 50.0000 && energy < 100.0000) * (0.006039 + (energy-50.000000)* 0.000003) + \
1082 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 200.0000) * (0.006170 + (energy-100.000000)* 0.000003) + \
1083 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 200.0000 && energy < 500.0000) * (0.006485 + (energy-200.000000)* 0.000006) + \
1084 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 500.0000 && energy < 1000.0000) * (0.008140 + (energy-500.000000)* 0.000008) + \
1085 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.012304 + (energy-1000.000000)* 0.000010) + \
1086 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.022168 + (energy-2000.000000)* 0.000010) + \
1087 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.053585 + (energy-5000.000000)* 0.000011) + \
1088 (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.106635*energy/10000.000000)
[a02a49e]1089 }
1090}
1091
[e01b141]1092###############################
1093# Momentum resolution for muons
1094###############################
1095
1096module MomentumSmearing ForwardMuonMomentumSmearing {
1097 set InputArray ForwardMuonEfficiency/muons
1098 set OutputArray muons
1099
[55f821a]1100 # Resolution given in dpT/pT assume constant 10% regardless of energy
[e01b141]1101 set ResolutionFormula {
[55f821a]1102 (abs(eta) >= 2.5 && abs(eta) < 8.0) * (0.10)
[e01b141]1103 }
1104}
1105
1106
[a02a49e]1107##############
1108# Track merger
1109##############
1110
1111module Merger TrackMerger {
1112 # add InputArray InputArray
1113 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
1114 add InputArray ElectronMomentumSmearing/electrons
1115 add InputArray MuonMomentumSmearing/muons
1116 set OutputArray tracks
1117}
1118
1119#############
1120# ECAL
1121#############
1122
1123module SimpleCalorimeter ECal {
1124 set ParticleInputArray ParticlePropagator/stableParticles
1125 set TrackInputArray TrackMerger/tracks
1126
1127 set TowerOutputArray ecalTowers
1128 set EFlowTrackOutputArray eflowTracks
1129 set EFlowTowerOutputArray eflowPhotons
1130
1131 set IsEcal true
1132
1133 set EnergyMin 0.5
1134 set EnergySignificanceMin 1.0
1135
1136 set SmearTowerCenter true
1137
1138 set pi [expr {acos(-1)}]
1139
1140 # lists of the edges of each tower in eta and phi
1141 # each list starts with the lower edge of the first tower
1142 # the list ends with the higher edged of the last tower
1143
1144 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
1145 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
1146
1147 #barrel:
1148 #dphi = 0.2 degree towers up to eta <=1.2
1149 set PhiBins {}
1150 for {set i -900} {$i <= 900} {incr i} {
1151 add PhiBins [expr {$i * $pi/900.0 }]
1152 }
1153 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
1154 for {set i -400} {$i <=400} {incr i} {
1155 set eta [expr {$i * 0.003}]
1156 add EtaPhiBins $eta $PhiBins
1157 }
1158
1159 #endcaps:
1160 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
1161 set PhiBins {}
1162 for {set i -225} {$i <= 225} {incr i} {
1163 add PhiBins [expr {$i * $pi/225.}]
1164 }
1165 #deta=0.02 units for 1.2 < |eta| <=2.5
1166 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
1167 for {set i 1} {$i <=66} {incr i} {
1168 set eta [expr {-2.52 + $i * 0.02}]
1169 add EtaPhiBins $eta $PhiBins
1170 }
1171 #same for 1.2 to 2.5
1172 for {set i 1} {$i <=66} {incr i} {
1173 set eta [expr {1.18 + $i*0.02}]
1174 add EtaPhiBins $eta $PhiBins
1175 }
1176
1177
1178 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
1179
1180 add EnergyFraction {0} {0.0}
1181 # energy fractions for e, gamma and pi0
1182 add EnergyFraction {11} {1.0}
1183 add EnergyFraction {22} {1.0}
1184 add EnergyFraction {111} {1.0}
1185 # energy fractions for muon, neutrinos and neutralinos
1186 add EnergyFraction {12} {0.0}
1187 add EnergyFraction {13} {0.0}
1188 add EnergyFraction {14} {0.0}
1189 add EnergyFraction {16} {0.0}
1190 add EnergyFraction {1000022} {0.0}
1191 add EnergyFraction {1000023} {0.0}
1192 add EnergyFraction {1000025} {0.0}
1193 add EnergyFraction {1000035} {0.0}
1194 add EnergyFraction {1000045} {0.0}
1195 # energy fractions for K0short and Lambda
1196 add EnergyFraction {310} {0.3}
1197 add EnergyFraction {3122} {0.3}
1198
1199 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
1200 set ResolutionFormula {
1201 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
1202 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
1203 (abs(eta) <= 2.5 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
1204}
1205
1206#############
1207# HCAL
1208#############
1209
1210module SimpleCalorimeter HCal {
1211 set ParticleInputArray ParticlePropagator/stableParticles
1212 set TrackInputArray ECal/eflowTracks
1213
1214 set TowerOutputArray hcalTowers
1215 set EFlowTrackOutputArray eflowTracks
1216 set EFlowTowerOutputArray eflowNeutralHadrons
1217
1218 set IsEcal false
1219
1220 set EnergyMin 1.0
1221 set EnergySignificanceMin 1.0
1222
1223 set SmearTowerCenter true
1224
1225 set pi [expr {acos(-1)}]
1226
1227 # lists of the edges of each tower in eta and phi
1228 # each list starts with the lower edge of the first tower
1229 # the list ends with the higher edged of the last tower
1230
1231
1232 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
1233 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
1234 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=2.5
1235 #HCAL cell sizes always 30x30 mm^2
1236
1237 #barrel and ring:
1238 #dphi = 1 degree up to |eta| <=0.9
1239 set PhiBins {}
1240 for {set i -180} {$i <=180} {incr i} {
1241 add PhiBins [expr {$i * $pi/180.0}]
1242 }
1243 #deta= 0.02 towers up to |eta| <=0.9
1244 for {set i -45} {$i <=45} {incr i} {
1245 set eta [expr {$i * 0.02}]
1246 add EtaPhiBins $eta $PhiBins
1247 }
1248
1249 #endcaps:
1250 # dphi = 6 degree
1251 set PhiBins {}
1252 for {set i -30} {$i <=30} {incr i} {
1253 add PhiBins [expr {$i * $pi/30.0}]
1254 }
1255 # deta =0.1 for 0.9 < |eta| <=2.5
1256 #for -2.5 to -0.9, 21 segments
1257 for {set i 1} {$i <=17} {incr i} {
1258 set eta [expr {-2.5 + $i * 0.1}]
1259 add EtaPhiBins $eta $PhiBins
1260 }
1261 #same for 0.9 to 2.5
1262 for {set i 1} {$i <=17} {incr i} {
1263 set eta [expr {0.8 + $i * 0.1 }]
1264 add EtaPhiBins $eta $PhiBins
1265 }
1266
1267 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
1268 add EnergyFraction {0} {1.0}
1269 # energy fractions for e, gamma and pi0
1270 add EnergyFraction {11} {0.0}
1271 add EnergyFraction {22} {0.0}
1272 add EnergyFraction {111} {0.0}
1273 # energy fractions for muon, neutrinos and neutralinos
1274 add EnergyFraction {12} {0.0}
1275 add EnergyFraction {13} {0.0}
1276 add EnergyFraction {14} {0.0}
1277 add EnergyFraction {16} {0.0}
1278 add EnergyFraction {1000022} {0.0}
1279 add EnergyFraction {1000023} {0.0}
1280 add EnergyFraction {1000025} {0.0}
1281 add EnergyFraction {1000035} {0.0}
1282 add EnergyFraction {1000045} {0.0}
1283 # energy fractions for K0short and Lambda
1284 add EnergyFraction {310} {0.7}
1285 add EnergyFraction {3122} {0.7}
1286
1287 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
1288 #CLICdet internal studies
1289 set ResolutionFormula {
1290 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
1291 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
1292 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
1293 (abs(eta)<=2.5 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
1294 }
1295
1296}
1297
1298#################
1299# Electron filter
1300#################
1301
1302module PdgCodeFilter ElectronFilter {
1303 set InputArray HCal/eflowTracks
1304 set OutputArray electrons
1305 set Invert true
1306 add PdgCode {11}
1307 add PdgCode {-11}
1308}
1309
1310######################
1311# ChargedHadronFilter
1312######################
1313
1314module PdgCodeFilter ChargedHadronFilter {
1315 set InputArray HCal/eflowTracks
1316 set OutputArray chargedHadrons
1317
1318 add PdgCode {11}
1319 add PdgCode {-11}
1320 add PdgCode {13}
1321 add PdgCode {-13}
1322}
1323
1324
1325
1326###################################################
1327# Tower Merger (in case not using e-flow algorithm)
1328###################################################
1329
1330module Merger Calorimeter {
1331 # add InputArray InputArray
1332 add InputArray ECal/ecalTowers
1333 add InputArray HCal/hcalTowers
1334 add InputArray MuonMomentumSmearing/muons
1335 set OutputArray towers
1336}
1337
1338
1339####################
1340# Energy flow merger
1341####################
1342
1343module Merger EFlowMerger {
1344 # add InputArray InputArray
1345 add InputArray HCal/eflowTracks
1346 add InputArray ECal/eflowPhotons
1347 add InputArray HCal/eflowNeutralHadrons
1348 set OutputArray eflow
1349}
1350
1351######################
1352# EFlowFilter
1353######################
1354
1355module PdgCodeFilter EFlowFilter {
1356 set InputArray EFlowMerger/eflow
1357 set OutputArray eflow
1358
1359 add PdgCode {11}
1360 add PdgCode {-11}
1361 add PdgCode {13}
1362 add PdgCode {-13}
1363}
1364
1365
1366###################
1367# Photon efficiency
1368###################
1369
1370module Efficiency PhotonEfficiency {
1371 set InputArray ECal/eflowPhotons
1372 set OutputArray photons
1373
1374 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1375
1376 # efficiency formula for photons
1377 set EfficiencyFormula {
1378 (energy < 2.0 ) * (0.000) +
1379 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
1380 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=2.5) * (0.9) }
1381
1382}
1383
1384
1385##################
1386# Photon isolation
1387##################
1388
1389module Isolation PhotonIsolation {
1390 set CandidateInputArray PhotonEfficiency/photons
1391 set IsolationInputArray EFlowMerger/eflow
1392
1393 set OutputArray photons
1394
1395 set DeltaRMax 0.1
1396
1397 set PTMin 0.5
1398
1399 set PTRatioMax 0.2
1400}
1401
1402
1403#####################
1404# Electron efficiency
1405#####################
1406
1407module Efficiency ElectronEfficiency {
1408 set InputArray ElectronFilter/electrons
1409 set OutputArray electrons
1410
1411 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1412
1413 set EfficiencyFormula {
1414 (energy < 3.0 ) * ( 0.00 ) +
1415 (abs(eta) > 2.50) * ( 0.00 ) +
1416 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.58 ) +
1417 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
1418 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
1419 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
1420 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
1421 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
1422 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.6 ) +
1423 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
1424 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
1425 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
1426 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
1427 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
1428 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.6 ) +
1429 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
1430 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
1431 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
1432 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
1433 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
1434 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.64 ) +
1435 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
1436 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
1437 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
1438 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
1439 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
1440 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.64 ) +
1441 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
1442 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
1443 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
1444 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
1445 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
1446 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.67 ) +
1447 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
1448 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
1449 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
1450 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
1451 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
1452 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1453 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
1454 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
1455 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
1456 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
1457 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
1458 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.7 ) +
1459 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
1460 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
1461 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
1462 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
1463 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
1464 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1465 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
1466 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
1467 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
1468 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
1469 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
1470 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.68 ) +
1471 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
1472 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
1473 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
1474 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
1475 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
1476 ( energy >=400 ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
1477 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
1478 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
1479 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
1480 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
1481 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
1482 }
1483}
1484
1485####################
1486# Electron isolation
1487####################
1488
1489module Isolation ElectronIsolation {
1490 set CandidateInputArray ElectronEfficiency/electrons
1491 set IsolationInputArray EFlowMerger/eflow
1492
1493 set OutputArray electrons
1494
1495 set DeltaRMax 0.1
1496
1497 set PTMin 0.5
1498
1499 set PTRatioMax 0.2
1500}
1501
1502#################
1503# Muon efficiency
1504#################
1505
1506module Efficiency MuonEfficiency {
1507 set InputArray MuonMomentumSmearing/muons
1508 set OutputArray muons
1509
1510 # set EfficiencyFormula {efficiency as a function of eta and pt}
1511
1512 set EfficiencyFormula {
1513 (energy < 2.5 ) * (0.00) +
1514 (energy>=2.5 ) * (0.999)
1515 }
1516}
1517
1518################
1519# Muon isolation
1520################
1521
1522module Isolation MuonIsolation {
1523 set CandidateInputArray MuonEfficiency/muons
1524 set IsolationInputArray EFlowMerger/eflow
1525
1526 set OutputArray muons
1527
1528 set DeltaRMax 0.1
1529
1530 set PTMin 0.5
1531
1532 set PTRatioMax 0.2
1533}
1534
1535
1536
1537###################
1538# Missing ET merger
1539###################
1540
1541module Merger MissingET {
1542 # add InputArray InputArray
1543 add InputArray EFlowMerger/eflow
1544 set MomentumOutputArray momentum
1545}
1546
1547
1548##################
1549# Scalar HT merger
1550##################
1551
1552module Merger ScalarHT {
1553 # add InputArray InputArray
1554 add InputArray EFlowMerger/eflow
1555 set EnergyOutputArray energy
1556}
1557######################
1558# EFlowFilter (UniqueObjectFinder)
1559######################
1560module UniqueObjectFinder UniqueObjectFinder {
1561 add InputArray PhotonIsolation/photons photons
1562 add InputArray ElectronIsolation/electrons electrons
1563 add InputArray MuonIsolation/muons muons
1564 add InputArray EFlowMerger/eflow eflow
1565}
1566
1567#################
1568# Neutrino Filter
1569#################
1570
1571module PdgCodeFilter NeutrinoFilter {
1572
1573 set InputArray Delphes/stableParticles
1574 set OutputArray filteredParticles
1575
1576 set PTMin 0.0
1577
1578 add PdgCode {12}
1579 add PdgCode {14}
1580 add PdgCode {16}
1581 add PdgCode {-12}
1582 add PdgCode {-14}
1583 add PdgCode {-16}
1584
1585}
1586
1587
1588#####################
1589# MC truth jet finder
1590#####################
1591
1592module FastJetFinder GenJetFinder {
1593 set InputArray NeutrinoFilter/filteredParticles
1594
1595 set OutputArray jets
1596
1597 # 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
1598 set JetAlgorithm 9
1599 set ParameterR 0.5
1600
1601 set JetPTMin 20.0
1602}
1603
1604#########################
1605# Gen Missing ET merger
1606########################
1607
1608module Merger GenMissingET {
1609 # add InputArray InputArray
1610 add InputArray NeutrinoFilter/filteredParticles
1611 set MomentumOutputArray momentum
1612}
1613
1614
1615
1616############
1617# Jet finder
1618############
1619
1620module FastJetFinder FastJetFinderKt {
1621 # set InputArray Calorimeter/towers
1622 set InputArray EFlowMerger/eflow
1623
1624 set OutputArray KTjets
1625
1626 # 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
1627 set JetAlgorithm 4
1628 set ParameterR 0.5
1629
1630 set JetPTMin 20.0
1631}
1632
1633################
1634# Jet finder VLC
1635################
1636
1637
1638source MuonCollider/MuonColliderDet_JetReco.tcl
1639
1640
1641#########################################
1642# Jet Momentum Smearing to mimick overlay
1643#########################################
1644
1645
1646source MuonCollider/MuonColliderDet_JetSmearing.tcl
1647
1648
1649
1650########################
1651# Jet Flavor Association
1652########################
1653
1654source MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl
1655
1656###########
1657# b-tagging
1658###########
1659# based on CLICdp-Note-2014-002
1660
1661source MuonCollider/MuonColliderDet_BTagging.tcl
1662
1663
1664#############
1665# tau-tagging
1666#############
1667
1668source MuonCollider/MuonColliderDet_TauTagging.tcl
1669
1670##################
1671# ROOT tree writer
1672##################
1673
1674module TreeWriter TreeWriter {
1675 # add Branch InputArray BranchName BranchClass
1676 add Branch Delphes/allParticles Particle GenParticle
1677
1678 add Branch GenJetFinder/jets GenJet Jet
1679
1680 add Branch FastJetFinderKt/KTjets KTjet Jet
1681
1682 add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 VLCjetR05N2 Jet
1683 add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 VLCjetR05N3 Jet
1684 add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 VLCjetR05N4 Jet
1685 add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 VLCjetR05N5 Jet
1686 add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 VLCjetR05N6 Jet
1687
1688 add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 VLCjetR07N2 Jet
1689 add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 VLCjetR07N3 Jet
1690 add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 VLCjetR07N4 Jet
1691 add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 VLCjetR07N5 Jet
1692 add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 VLCjetR07N6 Jet
1693
1694 add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 VLCjetR10N2 Jet
1695 add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 VLCjetR10N3 Jet
1696 add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 VLCjetR10N4 Jet
1697 add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 VLCjetR10N5 Jet
1698 add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 VLCjetR10N6 Jet
1699
1700 add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 VLCjetR12N2 Jet
1701 add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 VLCjetR12N3 Jet
1702 add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 VLCjetR12N4 Jet
1703 add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 VLCjetR12N5 Jet
1704 add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 VLCjetR12N6 Jet
1705
1706 add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 VLCjetR15N2 Jet
1707 add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 VLCjetR15N3 Jet
1708 add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 VLCjetR15N4 Jet
1709 add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 VLCjetR15N5 Jet
1710 add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 VLCjetR15N6 Jet
1711
[e01b141]1712 add Branch JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive VLCjetR02_inclusive Jet
[a02a49e]1713 add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive VLCjetR05_inclusive Jet
1714 add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive VLCjetR07_inclusive Jet
1715 add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive VLCjetR10_inclusive Jet
1716 add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive VLCjetR12_inclusive Jet
1717 add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive VLCjetR15_inclusive Jet
1718
1719
1720 ####
1721
1722 add Branch GenMissingET/momentum GenMissingET MissingET
1723
1724 add Branch TrackMerger/tracks Track Track
1725 add Branch Calorimeter/towers Tower Tower
1726
1727 add Branch HCal/eflowTracks EFlowTrack Track
1728 add Branch ECal/eflowPhotons EFlowPhoton Tower
1729 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1730
1731 add Branch UniqueObjectFinder/photons Photon Photon
1732 add Branch UniqueObjectFinder/electrons Electron Electron
1733 add Branch UniqueObjectFinder/muons Muon Muon
[e01b141]1734 add Branch ForwardMuonMomentumSmearing/muons ForwardMuon Muon
[a02a49e]1735
1736 add Branch MissingET/momentum MissingET MissingET
1737 add Branch ScalarHT/energy ScalarHT ScalarHT
1738}
Note: See TracBrowser for help on using the repository browser.