Fork me on GitHub

Changes in / [fa7729a:4f52352] in git


Ignore:
Location:
cards
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • cards/MuonCollider/MuonColliderDet_BTagging.tcl

    rfa7729a r4f52352  
    1212        source MuonCollider/MuonColliderDet_BTag_50.tcl
    1313}
    14 
    15 module BTagging BTagging_WP50_R02N2 {
     14module BTagging BTagging_WP70_R02N2 {
    1615        set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
     16        set BitNumber 1
     17        source MuonCollider/MuonColliderDet_BTag_70.tcl
     18}
     19module BTagging BTagging_WP90_R02N2 {
     20        set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
     21        set BitNumber 2
     22        source MuonCollider/MuonColliderDet_BTag_90.tcl
     23}
     24module BTagging BTagging_WP50_R02N3 {
     25        set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
     26        set BitNumber 0
     27        source MuonCollider/MuonColliderDet_BTag_50.tcl
     28}
     29module BTagging BTagging_WP70_R02N3 {
     30        set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
     31        set BitNumber 1
     32        source MuonCollider/MuonColliderDet_BTag_70.tcl
     33}
     34module BTagging BTagging_WP90_R02N3 {
     35        set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
     36        set BitNumber 2
     37        source MuonCollider/MuonColliderDet_BTag_90.tcl
     38}
     39module BTagging BTagging_WP50_R02N4 {
     40        set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
     41        set BitNumber 0
     42        source MuonCollider/MuonColliderDet_BTag_50.tcl
     43}
     44module BTagging BTagging_WP70_R02N4 {
     45        set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
     46        set BitNumber 1
     47        source MuonCollider/MuonColliderDet_BTag_70.tcl
     48}
     49module BTagging BTagging_WP90_R02N4 {
     50        set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
     51        set BitNumber 2
     52        source MuonCollider/MuonColliderDet_BTag_90.tcl
     53}
     54module BTagging BTagging_WP50_R02N5 {
     55        set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
     56        set BitNumber 0
     57        source MuonCollider/MuonColliderDet_BTag_50.tcl
     58}
     59module BTagging BTagging_WP70_R02N5 {
     60        set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
     61        set BitNumber 1
     62        source MuonCollider/MuonColliderDet_BTag_70.tcl
     63}
     64module BTagging BTagging_WP90_R02N5 {
     65        set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
     66        set BitNumber 2
     67        source MuonCollider/MuonColliderDet_BTag_90.tcl
     68}
     69module BTagging BTagging_WP50_R02N6 {
     70        set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
     71        set BitNumber 0
     72        source MuonCollider/MuonColliderDet_BTag_50.tcl
     73}
     74module BTagging BTagging_WP70_R02N6 {
     75        set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
     76        set BitNumber 1
     77        source MuonCollider/MuonColliderDet_BTag_70.tcl
     78}
     79module BTagging BTagging_WP90_R02N6 {
     80        set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
     81        set BitNumber 2
     82        source MuonCollider/MuonColliderDet_BTag_90.tcl
     83}
     84
     85module BTagging BTagging_WP50_R05N2 {
     86        set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
    1787        set BitNumber 0
    1888
     
    2696        source MuonCollider/MuonColliderDet_BTag_50.tcl
    2797}
    28 module BTagging BTagging_WP70_R02N2 {
    29         set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
    30         set BitNumber 1
    31         source MuonCollider/MuonColliderDet_BTag_70.tcl
    32 }
    33 module BTagging BTagging_WP90_R02N2 {
    34         set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
    35         set BitNumber 2
    36         source MuonCollider/MuonColliderDet_BTag_90.tcl
    37 }
    38 module BTagging BTagging_WP50_R02N3 {
    39         set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
    40         set BitNumber 0
    41         source MuonCollider/MuonColliderDet_BTag_50.tcl
    42 }
    43 module BTagging BTagging_WP70_R02N3 {
    44         set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
    45         set BitNumber 1
    46         source MuonCollider/MuonColliderDet_BTag_70.tcl
    47 }
    48 module BTagging BTagging_WP90_R02N3 {
    49         set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
    50         set BitNumber 2
    51         source MuonCollider/MuonColliderDet_BTag_90.tcl
    52 }
    53 module BTagging BTagging_WP50_R02N4 {
    54         set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
    55         set BitNumber 0
    56         source MuonCollider/MuonColliderDet_BTag_50.tcl
    57 }
    58 module BTagging BTagging_WP70_R02N4 {
    59         set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
    60         set BitNumber 1
    61         source MuonCollider/MuonColliderDet_BTag_70.tcl
    62 }
    63 module BTagging BTagging_WP90_R02N4 {
    64         set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
    65         set BitNumber 2
    66         source MuonCollider/MuonColliderDet_BTag_90.tcl
    67 }
    68 module BTagging BTagging_WP50_R02N5 {
    69         set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
    70         set BitNumber 0
    71         source MuonCollider/MuonColliderDet_BTag_50.tcl
    72 }
    73 module BTagging BTagging_WP70_R02N5 {
    74         set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
    75         set BitNumber 1
    76         source MuonCollider/MuonColliderDet_BTag_70.tcl
    77 }
    78 module BTagging BTagging_WP90_R02N5 {
    79         set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
    80         set BitNumber 2
    81         source MuonCollider/MuonColliderDet_BTag_90.tcl
    82 }
    83 module BTagging BTagging_WP50_R02N6 {
    84         set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
    85         set BitNumber 0
    86         source MuonCollider/MuonColliderDet_BTag_50.tcl
    87 }
    88 module BTagging BTagging_WP70_R02N6 {
    89         set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
    90         set BitNumber 1
    91         source MuonCollider/MuonColliderDet_BTag_70.tcl
    92 }
    93 module BTagging BTagging_WP90_R02N6 {
    94         set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
    95         set BitNumber 2
    96         source MuonCollider/MuonColliderDet_BTag_90.tcl
    97 }
    98 module BTagging BTagging_WP50_R07N2 {
    99         set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
    100         set BitNumber 0
    101         source MuonCollider/MuonColliderDet_BTag_50.tcl
    102 }
    103 
    104 module BTagging BTagging_WP50_R05N2 {
    105         set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
    106         set BitNumber 0
    107 
    108         # 50% efficiency working point
    109     # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    110     # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
    111     # gluon's PDG code has the lowest priority
    112 
    113         # based on CLICdp-Note-2014-002
    114 
    115         source MuonCollider/MuonColliderDet_BTag_50.tcl
    116 }
    11798module BTagging BTagging_WP70_R05N2 {
    11899        set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
     
    185166        source MuonCollider/MuonColliderDet_BTag_90.tcl
    186167}
     168
    187169module BTagging BTagging_WP50_R07N2 {
    188170        set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
  • cards/MuonCollider/MuonColliderDet_JetReco.tcl

    rfa7729a r4f52352  
    390390    set Beta 1.0
    391391    set Gamma 1.0
    392 
    393     set ComputeNsubjettiness 1
    394     set AxisMode 4
    395 
    396     set ComputeSoftDrop 1
    397     set BetaSoftDrop 0.0
    398     set SymmetryCutSoftDrop 0.1
    399     set R0SoftDrop 0.2
    400 
    401     set JetPTMin 20.0
    402 
    403 
    404 
    405 
    406 
     392    set JetPTMin 20.0
    407393}
    408394#R05
     
    415401    set Beta 1.0
    416402    set Gamma 1.0
    417     set ComputeNsubjettiness 1
    418     set AxisMode 4
    419 
    420     set ComputeSoftDrop 1
    421     set BetaSoftDrop 0.0
    422     set SymmetryCutSoftDrop 0.1
    423     set R0SoftDrop 0.5
    424403    set JetPTMin 20.0
    425404}
     
    433412    set Beta 1.0
    434413    set Gamma 1.0
    435     set ComputeNsubjettiness 1
    436     set AxisMode 4
    437 
    438     set ComputeSoftDrop 1
    439     set BetaSoftDrop 0.0
    440     set SymmetryCutSoftDrop 0.1
    441     set R0SoftDrop 0.7
    442414    set JetPTMin 20.0
    443415}
     
    451423    set Beta 1.0
    452424    set Gamma 1.0
    453     set ComputeNsubjettiness 1
    454     set AxisMode 4
    455 
    456     set ComputeSoftDrop 1
    457     set BetaSoftDrop 0.0
    458     set SymmetryCutSoftDrop 0.1
    459     set R0SoftDrop 1.0
    460425    set JetPTMin 20.0
    461426}
     
    469434    set Beta 1.0
    470435    set Gamma 1.0
    471     set ComputeNsubjettiness 1
    472     set AxisMode 4
    473 
    474     set ComputeSoftDrop 1
    475     set BetaSoftDrop 0.0
    476     set SymmetryCutSoftDrop 0.1
    477     set R0SoftDrop 1.2
    478436    set JetPTMin 20.0
    479437}
     
    487445    set Beta 1.0
    488446    set Gamma 1.0
    489     set ComputeNsubjettiness 1
    490     set AxisMode 4
    491 
    492     set ComputeSoftDrop 1
    493     set BetaSoftDrop 0.0
    494     set SymmetryCutSoftDrop 0.1
    495     set R0SoftDrop 1.5
    496 
    497     set JetPTMin 20.0
    498 }
     447    set JetPTMin 20.0
     448}
  • cards/delphes_card_MuonColliderDet.tcl

    rfa7729a r4f52352  
    1414#
    1515# Comments:
    16 #   - fix: replace momentum resolution by FCC-hh
    1716#   - fix: angle param in DenseTrackFilter
    1817#   - fix: add FWD muon collection
     
    4039    ElectronTrackingEfficiency
    4140    MuonTrackingEfficiency
    42     ForwardMuonEfficiency
    4341
    4442    ChargedHadronMomentumSmearing
    4543    ElectronMomentumSmearing
    4644    MuonMomentumSmearing
    47     ForwardMuonMomentumSmearing
    4845
    4946    TrackMerger
     
    507504}
    508505
    509 ##########################
    510 # Forward Muon efficiency
    511 ##########################
    512 
    513 ## hypothetical forward muon spectrometer
    514 module Efficiency ForwardMuonEfficiency {
    515     set InputArray ParticlePropagator/muons
    516     set OutputArray muons
    517 
    518     set EfficiencyFormula { (pt <= 0.5) * (0.00) +
    519     (abs(eta) > 2.5 && abs(eta) < 6.0) * (pt > 0.5 && pt <= 1) * (0.90) +
    520     (abs(eta) > 2.5 && abs(eta) < 6.0) * (pt > 1.0) * (0.95) +
    521     (abs(eta) > 6.0 ) * (0.00)
    522    }
    523 }
    524 
    525 
    526506########################################
    527507# Momentum resolution for charged tracks
     
    533513
    534514
    535     # Resolution given in dpT/pT (from FCC-hh)
    536     set ResolutionFormula {    (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
    537        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
    538        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
    539        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
    540        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
    541        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
    542        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
    543        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
    544        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
    545        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
    546        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
    547        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
    548        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
    549        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
    550        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
    551        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
    552        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
    553        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
    554        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
    555        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
    556        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
    557        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
    558        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
    559        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
    560        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
    561        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
    562        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
    563        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
    564        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
    565        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
    566        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
    567        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
    568        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
    569        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
    570        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
    571        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
    572        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
    573        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
    574        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
    575        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
    576        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
    577        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
    578        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
    579        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
    580        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
    581        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
    582        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
    583        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
    584        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
    585        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
    586        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
    587        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
    588        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
    589        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
    590        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
    591        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
    592        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
    593        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
    594        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
    595        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
    596        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
    597        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
    598        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
    599        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
    600        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
    601        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
    602        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
    603        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
    604        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
    605        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
    606        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
    607        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
    608        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
    609        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
    610        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
    611        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
    612        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
    613        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
    614        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
    615        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
    616        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
    617        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
    618        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
    619        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
    620        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
    621        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
    622        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
    623        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
    624        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
    625        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
    626        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
    627        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
    628        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
    629        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
    630        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
    631        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
    632        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
    633        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
    634        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
    635        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
    636        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
    637        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
    638        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
    639        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
    640        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
    641        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
    642        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
    643        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
    644        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
    645        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
    646        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
    647        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
    648        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
    649        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
    650        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
    651        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
    652        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
    653        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
    654        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
    655        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
    656        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
    657        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
    658        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
    659        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
    660        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
    661        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
    662        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
    663        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
    664        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
    665        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
    666        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
    667        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
    668        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
    669        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
    670        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
    671        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
    672        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
    673        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
    674        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
    675        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
    676        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
    677        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
    678        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
    679        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
    680        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
    681        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
    682        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
    683        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
    684        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
    685        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
    686        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
    687        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
    688        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
    689        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
    690        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
    691        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
    692        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
    693        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
    694        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
    695        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
    696        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
    697        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
    698        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
    699        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
    700        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
    701        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
    702        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
    703        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
    704        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
    705        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
    706        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
    707        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
    708        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
    709        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
    710        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
    711        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
    712        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
    713        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
    714        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
    715        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
    716        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
    717        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
    718        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
    719        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
    720        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
    721        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
    722        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
    723        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
    724        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
    725        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
    726        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
    727        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
    728        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
    729        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
    730        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
    731        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
    732        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
    733        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
    734        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
    735        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
     515    # Resolution given in dpT/pT.
     516    # CLICdet internal studies
     517    set ResolutionFormula {
     518        (abs(eta) < 2.50 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2    ) +
     519        (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2    ) +
     520        (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2   ) +
     521        (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2    ) +
     522        (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2   ) +
     523        (abs(eta) < 0.18)                      * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2   )
     524
    736525    }
    737526}
     
    745534    set OutputArray electrons
    746535
    747     # Resolution given in dpT/pT (from FCC-hh)
    748     set ResolutionFormula {    (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00315864) + \
    749        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.003159 + (energy-1.000000)* 0.000007) + \
    750        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.003166 + (energy-2.000000)* 0.000011) + \
    751        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000012) + \
    752        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003259 + (energy-10.000000)* 0.000010) + \
    753        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004173 + (energy-100.000000)* 0.000019) + \
    754        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020916 + (energy-1000.000000)* 0.000021) + \
    755        (abs(eta) >= 0.0000 && abs(eta) < 0.1000) * (energy >= 10000.0000) * (0.205876*energy/10000.000000) + \
    756        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00316278) + \
    757        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.003163 + (energy-1.000000)* 0.000006) + \
    758        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.003169 + (energy-2.000000)* 0.000010) + \
    759        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.003198 + (energy-5.000000)* 0.000011) + \
    760        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.003255 + (energy-10.000000)* 0.000010) + \
    761        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004165 + (energy-100.000000)* 0.000019) + \
    762        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020917 + (energy-1000.000000)* 0.000021) + \
    763        (abs(eta) >= 0.1000 && abs(eta) < 0.2000) * (energy >= 10000.0000) * (0.205952*energy/10000.000000) + \
    764        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00320482) + \
    765        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.003205 + (energy-1.000000)* 0.000006) + \
    766        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.003211 + (energy-2.000000)* 0.000009) + \
    767        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.003238 + (energy-5.000000)* 0.000011) + \
    768        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003294 + (energy-10.000000)* 0.000010) + \
    769        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004176 + (energy-100.000000)* 0.000018) + \
    770        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.020586 + (energy-1000.000000)* 0.000020) + \
    771        (abs(eta) >= 0.2000 && abs(eta) < 0.3000) * (energy >= 10000.0000) * (0.202528*energy/10000.000000) + \
    772        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00325680) + \
    773        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.003257 + (energy-1.000000)* 0.000001) + \
    774        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.003257 + (energy-2.000000)* 0.000009) + \
    775        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.003286 + (energy-5.000000)* 0.000011) + \
    776        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.003342 + (energy-10.000000)* 0.000009) + \
    777        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004195 + (energy-100.000000)* 0.000017) + \
    778        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019873 + (energy-1000.000000)* 0.000019) + \
    779        (abs(eta) >= 0.3000 && abs(eta) < 0.4000) * (energy >= 10000.0000) * (0.195142*energy/10000.000000) + \
    780        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00354020) + \
    781        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.003540 + (energy-1.000000)* -0.000201) + \
    782        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.003340 + (energy-2.000000)* 0.000009) + \
    783        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.003366 + (energy-5.000000)* 0.000011) + \
    784        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.003422 + (energy-10.000000)* 0.000009) + \
    785        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004265 + (energy-100.000000)* 0.000017) + \
    786        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019240 + (energy-1000.000000)* 0.000019) + \
    787        (abs(eta) >= 0.4000 && abs(eta) < 0.5000) * (energy >= 10000.0000) * (0.188429*energy/10000.000000) + \
    788        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00362672) + \
    789        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.003627 + (energy-1.000000)* -0.000223) + \
    790        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.003403 + (energy-2.000000)* 0.000008) + \
    791        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.003428 + (energy-5.000000)* 0.000010) + \
    792        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.003479 + (energy-10.000000)* 0.000009) + \
    793        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004256 + (energy-100.000000)* 0.000016) + \
    794        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.018471 + (energy-1000.000000)* 0.000018) + \
    795        (abs(eta) >= 0.5000 && abs(eta) < 0.6000) * (energy >= 10000.0000) * (0.180531*energy/10000.000000) + \
    796        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00386864) + \
    797        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.003869 + (energy-1.000000)* -0.000354) + \
    798        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.003515 + (energy-2.000000)* 0.000007) + \
    799        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.003536 + (energy-5.000000)* 0.000009) + \
    800        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003583 + (energy-10.000000)* 0.000009) + \
    801        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004349 + (energy-100.000000)* 0.000015) + \
    802        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017442 + (energy-1000.000000)* 0.000017) + \
    803        (abs(eta) >= 0.6000 && abs(eta) < 0.7000) * (energy >= 10000.0000) * (0.169559*energy/10000.000000) + \
    804        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00418203) + \
    805        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.004182 + (energy-1.000000)* -0.000556) + \
    806        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.003626 + (energy-2.000000)* 0.000007) + \
    807        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.003645 + (energy-5.000000)* 0.000008) + \
    808        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003687 + (energy-10.000000)* 0.000008) + \
    809        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004399 + (energy-100.000000)* 0.000013) + \
    810        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016509 + (energy-1000.000000)* 0.000016) + \
    811        (abs(eta) >= 0.7000 && abs(eta) < 0.8000) * (energy >= 10000.0000) * (0.159676*energy/10000.000000) + \
    812        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00436103) + \
    813        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.004361 + (energy-1.000000)* -0.000597) + \
    814        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.003764 + (energy-2.000000)* 0.000006) + \
    815        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.003781 + (energy-5.000000)* 0.000008) + \
    816        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003821 + (energy-10.000000)* 0.000007) + \
    817        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004495 + (energy-100.000000)* 0.000012) + \
    818        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015532 + (energy-1000.000000)* 0.000015) + \
    819        (abs(eta) >= 0.8000 && abs(eta) < 0.9000) * (energy >= 10000.0000) * (0.149090*energy/10000.000000) + \
    820        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00488279) + \
    821        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.004883 + (energy-1.000000)* -0.000969) + \
    822        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.003914 + (energy-2.000000)* 0.000006) + \
    823        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.003930 + (energy-5.000000)* 0.000007) + \
    824        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003967 + (energy-10.000000)* 0.000007) + \
    825        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004593 + (energy-100.000000)* 0.000011) + \
    826        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.014592 + (energy-1000.000000)* 0.000014) + \
    827        (abs(eta) >= 0.9000 && abs(eta) < 1.0000) * (energy >= 10000.0000) * (0.138764*energy/10000.000000) + \
    828        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00513716) + \
    829        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.005137 + (energy-1.000000)* -0.001026) + \
    830        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.004111 + (energy-2.000000)* 0.000005) + \
    831        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004125 + (energy-5.000000)* 0.000007) + \
    832        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004159 + (energy-10.000000)* 0.000007) + \
    833        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004752 + (energy-100.000000)* 0.000010) + \
    834        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013718 + (energy-1000.000000)* 0.000013) + \
    835        (abs(eta) >= 1.0000 && abs(eta) < 1.1000) * (energy >= 10000.0000) * (0.128750*energy/10000.000000) + \
    836        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00572019) + \
    837        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.005720 + (energy-1.000000)* -0.001362) + \
    838        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.004359 + (energy-2.000000)* 0.000005) + \
    839        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004372 + (energy-5.000000)* 0.000006) + \
    840        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004405 + (energy-10.000000)* 0.000007) + \
    841        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005024 + (energy-100.000000)* 0.000010) + \
    842        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013902 + (energy-1000.000000)* 0.000013) + \
    843        (abs(eta) >= 1.1000 && abs(eta) < 1.2000) * (energy >= 10000.0000) * (0.129437*energy/10000.000000) + \
    844        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00613558) + \
    845        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.006136 + (energy-1.000000)* -0.001331) + \
    846        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.004805 + (energy-2.000000)* 0.000005) + \
    847        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004818 + (energy-5.000000)* 0.000007) + \
    848        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004851 + (energy-10.000000)* 0.000007) + \
    849        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005512 + (energy-100.000000)* 0.000012) + \
    850        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.016144 + (energy-1000.000000)* 0.000015) + \
    851        (abs(eta) >= 1.2000 && abs(eta) < 1.3000) * (energy >= 10000.0000) * (0.151739*energy/10000.000000) + \
    852        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00655464) + \
    853        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.006555 + (energy-1.000000)* -0.001843) + \
    854        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.004711 + (energy-2.000000)* 0.000004) + \
    855        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.004724 + (energy-5.000000)* 0.000005) + \
    856        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004748 + (energy-10.000000)* 0.000006) + \
    857        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005326 + (energy-100.000000)* 0.000009) + \
    858        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013686 + (energy-1000.000000)* 0.000012) + \
    859        (abs(eta) >= 1.3000 && abs(eta) < 1.4000) * (energy >= 10000.0000) * (0.125361*energy/10000.000000) + \
    860        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00714442) + \
    861        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.007144 + (energy-1.000000)* -0.002674) + \
    862        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.004470 + (energy-2.000000)* -0.000070) + \
    863        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.004259 + (energy-5.000000)* 0.000006) + \
    864        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004287 + (energy-10.000000)* 0.000009) + \
    865        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005101 + (energy-100.000000)* 0.000008) + \
    866        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012303 + (energy-1000.000000)* 0.000011) + \
    867        (abs(eta) >= 1.4000 && abs(eta) < 1.5000) * (energy >= 10000.0000) * (0.110091*energy/10000.000000) + \
    868        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00680449) + \
    869        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.006804 + (energy-1.000000)* -0.002108) + \
    870        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.004696 + (energy-2.000000)* -0.000116) + \
    871        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.004348 + (energy-5.000000)* 0.000009) + \
    872        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.004391 + (energy-10.000000)* 0.000009) + \
    873        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005232 + (energy-100.000000)* 0.000007) + \
    874        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011392 + (energy-1000.000000)* 0.000010) + \
    875        (abs(eta) >= 1.5000 && abs(eta) < 1.6000) * (energy >= 10000.0000) * (0.097988*energy/10000.000000) + \
    876        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00763793) + \
    877        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.007638 + (energy-1.000000)* -0.003061) + \
    878        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.004577 + (energy-2.000000)* -0.000168) + \
    879        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.004074 + (energy-5.000000)* 0.000009) + \
    880        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.004121 + (energy-10.000000)* 0.000009) + \
    881        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004910 + (energy-100.000000)* 0.000006) + \
    882        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010635 + (energy-1000.000000)* 0.000009) + \
    883        (abs(eta) >= 1.6000 && abs(eta) < 1.7000) * (energy >= 10000.0000) * (0.090704*energy/10000.000000) + \
    884        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00913948) + \
    885        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.009139 + (energy-1.000000)* -0.004405) + \
    886        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.004735 + (energy-2.000000)* -0.000231) + \
    887        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.004043 + (energy-5.000000)* 0.000005) + \
    888        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.004066 + (energy-10.000000)* 0.000007) + \
    889        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004739 + (energy-100.000000)* 0.000006) + \
    890        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009891 + (energy-1000.000000)* 0.000008) + \
    891        (abs(eta) >= 1.7000 && abs(eta) < 1.8000) * (energy >= 10000.0000) * (0.083216*energy/10000.000000) + \
    892        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00956747) + \
    893        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.009567 + (energy-1.000000)* -0.004497) + \
    894        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.005070 + (energy-2.000000)* -0.000331) + \
    895        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.004078 + (energy-5.000000)* 0.000004) + \
    896        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.004100 + (energy-10.000000)* 0.000006) + \
    897        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004682 + (energy-100.000000)* 0.000005) + \
    898        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009157 + (energy-1000.000000)* 0.000007) + \
    899        (abs(eta) >= 1.8000 && abs(eta) < 1.9000) * (energy >= 10000.0000) * (0.074702*energy/10000.000000) + \
    900        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00964334) + \
    901        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.009643 + (energy-1.000000)* -0.003950) + \
    902        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.005694 + (energy-2.000000)* -0.000479) + \
    903        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.004256 + (energy-5.000000)* 0.000005) + \
    904        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.004281 + (energy-10.000000)* 0.000006) + \
    905        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004862 + (energy-100.000000)* 0.000004) + \
    906        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008501 + (energy-1000.000000)* 0.000006) + \
    907        (abs(eta) >= 1.9000 && abs(eta) < 2.0000) * (energy >= 10000.0000) * (0.062525*energy/10000.000000) + \
    908        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01045039) + \
    909        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010450 + (energy-1.000000)* -0.005379) + \
    910        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005072 + (energy-2.000000)* -0.000321) + \
    911        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004109 + (energy-5.000000)* 0.000006) + \
    912        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.004137 + (energy-10.000000)* 0.000006) + \
    913        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004666 + (energy-100.000000)* 0.000005) + \
    914        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008951 + (energy-1000.000000)* 0.000007) + \
    915        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073400*energy/10000.000000) + \
    916        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01046694) + \
    917        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010467 + (energy-1.000000)* -0.005023) + \
    918        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000330) + \
    919        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004455 + (energy-5.000000)* 0.000005) + \
    920        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004479 + (energy-10.000000)* 0.000004) + \
    921        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004883 + (energy-100.000000)* 0.000005) + \
    922        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009382 + (energy-1000.000000)* 0.000008) + \
    923        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078852*energy/10000.000000) + \
    924        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01090933) + \
    925        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010909 + (energy-1.000000)* -0.005299) + \
    926        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005610 + (energy-2.000000)* -0.000302) + \
    927        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004704 + (energy-5.000000)* 0.000005) + \
    928        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.004730 + (energy-10.000000)* 0.000005) + \
    929        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005146 + (energy-100.000000)* 0.000006) + \
    930        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010323 + (energy-1000.000000)* 0.000009) + \
    931        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.088469*energy/10000.000000) + \
    932        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01271833) + \
    933        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012718 + (energy-1.000000)* -0.005764) + \
    934        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.006954 + (energy-2.000000)* -0.000492) + \
    935        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005479 + (energy-5.000000)* 0.000003) + \
    936        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.005495 + (energy-10.000000)* 0.000003) + \
    937        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005774 + (energy-100.000000)* 0.000006) + \
    938        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011265 + (energy-1000.000000)* 0.000009) + \
    939        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.096592*energy/10000.000000) + \
    940        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01515272) + \
    941        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015153 + (energy-1.000000)* -0.007272) + \
    942        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007881 + (energy-2.000000)* -0.000660) + \
    943        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005900 + (energy-5.000000)* 0.000003) + \
    944        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.005914 + (energy-10.000000)* 0.000003) + \
    945        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006174 + (energy-100.000000)* 0.000007) + \
    946        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012486 + (energy-1000.000000)* 0.000011) + \
    947        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.108659*energy/10000.000000)
     536    # Resolution given in dpT/pT.
     537    # CLICdet internal studies
     538    set ResolutionFormula {
     539        (abs(eta) < 2.50 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2  + 0.0177556^2   ) +
     540        (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.62283e-05^2 * pt^2  + 0.0177556^2   ) +
     541        (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.0915e-05 ^2 * pt^2  + 0.00663766^2  ) +
     542        (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.15518e-05^2 * pt^2  + 0.00398644^2  ) +
     543        (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2  + 0.00317807^2  ) +
     544        (abs(eta) < 0.18)                      * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2  )
     545
    948546    }
    949547}
     
    959557    # Resolution given in dpT/pT.
    960558
     559    # CLICdet internal studies
    961560    set ResolutionFormula {
    962561
    963     ( abs(eta) < 1.35 ) *
    964 
    965     (sqrt(0.0000364164 + (
    966        9.06262e-8 *pt^2* cosh(
    967          eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + (
    968           504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/(
    969              pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
    970           sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/(
    971        0.00516429/sin(2*atan(exp(-abs(eta))))^2 + (
    972         96868.8 *(1/400000000 + (0.117945 * 1/cosh(eta)^2)/(
    973            pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
    974         sin(2*atan(exp(-abs(eta))))^2))
    975 
    976     ) +
    977 
    978     ( abs(eta) > 1.35 && abs(eta) < 2.00) *
    979 
    980     ( 1.53846*(
    981 
    982      (
    983        (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
    984     (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
    985     (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
    986     (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
    987     (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
    988     (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
    989     (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
    990     (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
    991     (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
    992     (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
    993     (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
    994     (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
    995     (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
    996     (energy >= 10000.0000) * (0.062437*energy/10000.000000)
    997      )
    998 
    999 
    1000      - sqrt(0.0000364164 + (9.06262e-8 *(1.19507e-6 + 2137.54 *(1/400000000 + 0.155982/pt^2)) *pt^2)/(0.0218797 + 410407. *(1/400000000 + 0.155982/pt^2))))*abs(eta) +
    1001 
    1002     3.07692*sqrt(0.0000364164 + (9.06262e-8 *(1.19507e-6 + 2137.54 *(1/400000000 + 0.155982/pt^2)) *pt^2)/(0.0218797 + 410407. *(1/400000000 + 0.155982/pt^2)))
    1003 
    1004      - 2.07692*
    1005      (
    1006 
    1007       (energy >= 0.0000 && energy < 1.0000) * (0.00953577) +
    1008     (energy >= 1.0000 && energy < 2.0000) * (0.009536 + (energy-1.000000)* -0.003793) +
    1009     (energy >= 2.0000 && energy < 5.0000) * (0.005742 + (energy-2.000000)* -0.000489) +
    1010     (energy >= 5.0000 && energy < 10.0000) * (0.004277 + (energy-5.000000)* 0.000005) +
    1011     (energy >= 10.0000 && energy < 20.0000) * (0.004302 + (energy-10.000000)* 0.000007) +
    1012     (energy >= 20.0000 && energy < 50.0000) * (0.004368 + (energy-20.000000)* 0.000007) +
    1013     (energy >= 50.0000 && energy < 100.0000) * (0.004581 + (energy-50.000000)* 0.000006) +
    1014     (energy >= 100.0000 && energy < 200.0000) * (0.004875 + (energy-100.000000)* 0.000005) +
    1015     (energy >= 200.0000 && energy < 500.0000) * (0.005344 + (energy-200.000000)* 0.000004) +
    1016     (energy >= 500.0000 && energy < 1000.0000) * (0.006395 + (energy-500.000000)* 0.000004) +
    1017     (energy >= 1000.0000 && energy < 2000.0000) * (0.008441 + (energy-1000.000000)* 0.000005) +
    1018     (energy >= 2000.0000 && energy < 5000.0000) * (0.013700 + (energy-2000.000000)* 0.000006) +
    1019     (energy >= 5000.0000 && energy < 10000.0000) * (0.031615 + (energy-5000.000000)* 0.000006) +
    1020     (energy >= 10000.0000) * (0.062437*energy/10000.000000)
    1021 
    1022      )
    1023 
    1024 
    1025     ) +
    1026 
    1027        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.01062416) + \
    1028        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.010624 + (energy-1.000000)* -0.005532) + \
    1029        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.005092 + (energy-2.000000)* -0.000326) + \
    1030        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.004115 + (energy-5.000000)* 0.000006) + \
    1031        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10.0000 && energy < 20.0000) * (0.004143 + (energy-10.000000)* 0.000007) + \
    1032        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 20.0000 && energy < 50.0000) * (0.004209 + (energy-20.000000)* 0.000007) + \
    1033        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 50.0000 && energy < 100.0000) * (0.004413 + (energy-50.000000)* 0.000005) + \
    1034        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 100.0000 && energy < 200.0000) * (0.004681 + (energy-100.000000)* 0.000004) + \
    1035        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 200.0000 && energy < 500.0000) * (0.005076 + (energy-200.000000)* 0.000004) + \
    1036        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006270 + (energy-500.000000)* 0.000005) + \
    1037        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.008960 + (energy-1000.000000)* 0.000007) + \
    1038        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.015510 + (energy-2000.000000)* 0.000007) + \
    1039        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.036867 + (energy-5000.000000)* 0.000007) + \
    1040        (abs(eta) >= 2.0000 && abs(eta) < 2.1000) * (energy >= 10000.0000) * (0.073168*energy/10000.000000) + \
    1041        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.01007098) + \
    1042        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.010071 + (energy-1.000000)* -0.004627) + \
    1043        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.005444 + (energy-2.000000)* -0.000322) + \
    1044        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.004478 + (energy-5.000000)* 0.000005) + \
    1045        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10.0000 && energy < 20.0000) * (0.004501 + (energy-10.000000)* 0.000006) + \
    1046        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 20.0000 && energy < 50.0000) * (0.004558 + (energy-20.000000)* 0.000005) + \
    1047        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 50.0000 && energy < 100.0000) * (0.004701 + (energy-50.000000)* 0.000004) + \
    1048        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 100.0000 && energy < 200.0000) * (0.004888 + (energy-100.000000)* 0.000003) + \
    1049        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 200.0000 && energy < 500.0000) * (0.005213 + (energy-200.000000)* 0.000004) + \
    1050        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006429 + (energy-500.000000)* 0.000006) + \
    1051        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.009343 + (energy-1000.000000)* 0.000007) + \
    1052        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.016410 + (energy-2000.000000)* 0.000008) + \
    1053        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.039265 + (energy-5000.000000)* 0.000008) + \
    1054        (abs(eta) >= 2.1000 && abs(eta) < 2.2000) * (energy >= 10000.0000) * (0.078014*energy/10000.000000) + \
    1055        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.01095892) + \
    1056        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.010959 + (energy-1.000000)* -0.005458) + \
    1057        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.005501 + (energy-2.000000)* -0.000281) + \
    1058        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.004660 + (energy-5.000000)* 0.000005) + \
    1059        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10.0000 && energy < 20.0000) * (0.004686 + (energy-10.000000)* 0.000007) + \
    1060        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 20.0000 && energy < 50.0000) * (0.004757 + (energy-20.000000)* 0.000006) + \
    1061        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 50.0000 && energy < 100.0000) * (0.004937 + (energy-50.000000)* 0.000004) + \
    1062        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 100.0000 && energy < 200.0000) * (0.005143 + (energy-100.000000)* 0.000004) + \
    1063        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 200.0000 && energy < 500.0000) * (0.005505 + (energy-200.000000)* 0.000005) + \
    1064        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 500.0000 && energy < 1000.0000) * (0.006975 + (energy-500.000000)* 0.000007) + \
    1065        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.010462 + (energy-1000.000000)* 0.000008) + \
    1066        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.018731 + (energy-2000.000000)* 0.000009) + \
    1067        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.045159 + (energy-5000.000000)* 0.000009) + \
    1068        (abs(eta) >= 2.2000 && abs(eta) < 2.3000) * (energy >= 10000.0000) * (0.089830*energy/10000.000000) + \
    1069        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.01279214) + \
    1070        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.012792 + (energy-1.000000)* -0.005763) + \
    1071        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.007029 + (energy-2.000000)* -0.000513) + \
    1072        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.005489 + (energy-5.000000)* 0.000003) + \
    1073        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10.0000 && energy < 20.0000) * (0.005503 + (energy-10.000000)* 0.000003) + \
    1074        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 20.0000 && energy < 50.0000) * (0.005537 + (energy-20.000000)* 0.000003) + \
    1075        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 50.0000 && energy < 100.0000) * (0.005636 + (energy-50.000000)* 0.000003) + \
    1076        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 100.0000 && energy < 200.0000) * (0.005773 + (energy-100.000000)* 0.000003) + \
    1077        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 200.0000 && energy < 500.0000) * (0.006060 + (energy-200.000000)* 0.000005) + \
    1078        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 500.0000 && energy < 1000.0000) * (0.007489 + (energy-500.000000)* 0.000007) + \
    1079        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.011117 + (energy-1000.000000)* 0.000009) + \
    1080        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.019824 + (energy-2000.000000)* 0.000009) + \
    1081        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.047732 + (energy-5000.000000)* 0.000009) + \
    1082        (abs(eta) >= 2.3000 && abs(eta) < 2.4000) * (energy >= 10000.0000) * (0.094931*energy/10000.000000) + \
    1083        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.01502671) + \
    1084        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.015027 + (energy-1.000000)* -0.007177) + \
    1085        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.007850 + (energy-2.000000)* -0.000651) + \
    1086        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.005898 + (energy-5.000000)* 0.000003) + \
    1087        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10.0000 && energy < 20.0000) * (0.005913 + (energy-10.000000)* 0.000003) + \
    1088        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 20.0000 && energy < 50.0000) * (0.005947 + (energy-20.000000)* 0.000003) + \
    1089        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 50.0000 && energy < 100.0000) * (0.006039 + (energy-50.000000)* 0.000003) + \
    1090        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 100.0000 && energy < 200.0000) * (0.006170 + (energy-100.000000)* 0.000003) + \
    1091        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 200.0000 && energy < 500.0000) * (0.006485 + (energy-200.000000)* 0.000006) + \
    1092        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 500.0000 && energy < 1000.0000) * (0.008140 + (energy-500.000000)* 0.000008) + \
    1093        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 1000.0000 && energy < 2000.0000) * (0.012304 + (energy-1000.000000)* 0.000010) + \
    1094        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 2000.0000 && energy < 5000.0000) * (0.022168 + (energy-2000.000000)* 0.000010) + \
    1095        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 5000.0000 && energy < 10000.0000) * (0.053585 + (energy-5000.000000)* 0.000011) + \
    1096        (abs(eta) >= 2.4000 && abs(eta) < 2.5000) * (energy >= 10000.0000) * (0.106635*energy/10000.000000)
    1097     }
    1098 }
    1099 
    1100 ###############################
    1101 # Momentum resolution for muons
    1102 ###############################
    1103 
    1104 module MomentumSmearing ForwardMuonMomentumSmearing {
    1105     set InputArray ForwardMuonEfficiency/muons
    1106     set OutputArray muons
    1107 
    1108     # Resolution given in dpT/pT (assume forward tracking spectrometer from FCC-hh)
    1109 
    1110     set ResolutionFormula {
    1111     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.01778926) + \
    1112     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.017789 + (energy-1.000000)* -0.008811) + \
    1113     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.008978 + (energy-2.000000)* -0.000844) + \
    1114     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.006446 + (energy-5.000000)* -0.000195) + \
    1115     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.005473 + (energy-10.000000)* 0.000003) + \
    1116     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005782 + (energy-100.000000)* 0.000005) + \
    1117     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010655 + (energy-1000.000000)* 0.000009) + \
    1118     (abs(eta) >= 2.5000 && abs(eta) < 2.6000) * (energy >= 10000.0000) * (0.087190*energy/10000.000000) + \
    1119     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.02061617) + \
    1120     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.020616 + (energy-1.000000)* -0.010732) + \
    1121     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.009884 + (energy-2.000000)* -0.001349) + \
    1122     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.005836 + (energy-5.000000)* -0.000484) + \
    1123     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.003417 + (energy-10.000000)* 0.000003) + \
    1124     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003655 + (energy-100.000000)* 0.000002) + \
    1125     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005164 + (energy-1000.000000)* 0.000003) + \
    1126     (abs(eta) >= 2.6000 && abs(eta) < 2.7000) * (energy >= 10000.0000) * (0.029564*energy/10000.000000) + \
    1127     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.02328940) + \
    1128     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.023289 + (energy-1.000000)* -0.012327) + \
    1129     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.010962 + (energy-2.000000)* -0.001777) + \
    1130     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.005632 + (energy-5.000000)* -0.000436) + \
    1131     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.003452 + (energy-10.000000)* 0.000001) + \
    1132     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003539 + (energy-100.000000)* 0.000001) + \
    1133     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004876 + (energy-1000.000000)* 0.000002) + \
    1134     (abs(eta) >= 2.7000 && abs(eta) < 2.8000) * (energy >= 10000.0000) * (0.026660*energy/10000.000000) + \
    1135     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.02303055) + \
    1136     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.023031 + (energy-1.000000)* -0.011212) + \
    1137     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.011819 + (energy-2.000000)* -0.002216) + \
    1138     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.005171 + (energy-5.000000)* -0.000352) + \
    1139     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.003412 + (energy-10.000000)* 0.000001) + \
    1140     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003516 + (energy-100.000000)* 0.000001) + \
    1141     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004565 + (energy-1000.000000)* 0.000002) + \
    1142     (abs(eta) >= 2.8000 && abs(eta) < 2.9000) * (energy >= 10000.0000) * (0.025158*energy/10000.000000) + \
    1143     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.02611889) + \
    1144     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.026119 + (energy-1.000000)* -0.012835) + \
    1145     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.013284 + (energy-2.000000)* -0.002589) + \
    1146     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.005516 + (energy-5.000000)* -0.000388) + \
    1147     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.003574 + (energy-10.000000)* 0.000001) + \
    1148     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003690 + (energy-100.000000)* 0.000001) + \
    1149     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004493 + (energy-1000.000000)* 0.000002) + \
    1150     (abs(eta) >= 2.9000 && abs(eta) < 3.0000) * (energy >= 10000.0000) * (0.022856*energy/10000.000000) + \
    1151     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.02991154) + \
    1152     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.029912 + (energy-1.000000)* -0.015226) + \
    1153     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.014686 + (energy-2.000000)* -0.002867) + \
    1154     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.006086 + (energy-5.000000)* -0.000457) + \
    1155     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.003803 + (energy-10.000000)* 0.000001) + \
    1156     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.003932 + (energy-100.000000)* 0.000001) + \
    1157     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004619 + (energy-1000.000000)* 0.000002) + \
    1158     (abs(eta) >= 3.0000 && abs(eta) < 3.1000) * (energy >= 10000.0000) * (0.020643*energy/10000.000000) + \
    1159     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.03556765) + \
    1160     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.035568 + (energy-1.000000)* -0.019079) + \
    1161     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.016488 + (energy-2.000000)* -0.003242) + \
    1162     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.006763 + (energy-5.000000)* -0.000526) + \
    1163     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.004133 + (energy-10.000000)* 0.000002) + \
    1164     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004270 + (energy-100.000000)* 0.000001) + \
    1165     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.004991 + (energy-1000.000000)* 0.000002) + \
    1166     (abs(eta) >= 3.1000 && abs(eta) < 3.2000) * (energy >= 10000.0000) * (0.021331*energy/10000.000000) + \
    1167     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.03803690) + \
    1168     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.038037 + (energy-1.000000)* -0.020041) + \
    1169     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.017996 + (energy-2.000000)* -0.003460) + \
    1170     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.007617 + (energy-5.000000)* -0.000724) + \
    1171     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.003999 + (energy-10.000000)* 0.000002) + \
    1172     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004159 + (energy-100.000000)* 0.000001) + \
    1173     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005300 + (energy-1000.000000)* 0.000002) + \
    1174     (abs(eta) >= 3.2000 && abs(eta) < 3.3000) * (energy >= 10000.0000) * (0.023674*energy/10000.000000) + \
    1175     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.03760924) + \
    1176     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.037609 + (energy-1.000000)* -0.018328) + \
    1177     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.019281 + (energy-2.000000)* -0.003682) + \
    1178     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.008234 + (energy-5.000000)* -0.000772) + \
    1179     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.004372 + (energy-10.000000)* 0.000002) + \
    1180     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004545 + (energy-100.000000)* 0.000001) + \
    1181     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.005786 + (energy-1000.000000)* 0.000002) + \
    1182     (abs(eta) >= 3.3000 && abs(eta) < 3.4000) * (energy >= 10000.0000) * (0.026924*energy/10000.000000) + \
    1183     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.03763696) + \
    1184     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.037637 + (energy-1.000000)* -0.016962) + \
    1185     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.020675 + (energy-2.000000)* -0.003893) + \
    1186     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.008995 + (energy-5.000000)* -0.000848) + \
    1187     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.004755 + (energy-10.000000)* 0.000002) + \
    1188     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.004949 + (energy-100.000000)* 0.000002) + \
    1189     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.006336 + (energy-1000.000000)* 0.000003) + \
    1190     (abs(eta) >= 3.4000 && abs(eta) < 3.5000) * (energy >= 10000.0000) * (0.030130*energy/10000.000000) + \
    1191     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.04368466) + \
    1192     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.043685 + (energy-1.000000)* -0.020430) + \
    1193     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.023254 + (energy-2.000000)* -0.004385) + \
    1194     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.010100 + (energy-5.000000)* -0.001016) + \
    1195     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.005020 + (energy-10.000000)* 0.000002) + \
    1196     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005206 + (energy-100.000000)* 0.000002) + \
    1197     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.006679 + (energy-1000.000000)* 0.000003) + \
    1198     (abs(eta) >= 3.5000 && abs(eta) < 3.6000) * (energy >= 10000.0000) * (0.033540*energy/10000.000000) + \
    1199     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.05055352) + \
    1200     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.050554 + (energy-1.000000)* -0.024554) + \
    1201     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.025999 + (energy-2.000000)* -0.004964) + \
    1202     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.011106 + (energy-5.000000)* -0.001134) + \
    1203     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.005436 + (energy-10.000000)* 0.000002) + \
    1204     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.005582 + (energy-100.000000)* 0.000002) + \
    1205     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.007047 + (energy-1000.000000)* 0.000003) + \
    1206     (abs(eta) >= 3.6000 && abs(eta) < 3.7000) * (energy >= 10000.0000) * (0.036533*energy/10000.000000) + \
    1207     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.05790950) + \
    1208     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.057910 + (energy-1.000000)* -0.028845) + \
    1209     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.029065 + (energy-2.000000)* -0.005581) + \
    1210     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.012322 + (energy-5.000000)* -0.001294) + \
    1211     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.005853 + (energy-10.000000)* 0.000002) + \
    1212     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006005 + (energy-100.000000)* 0.000002) + \
    1213     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.007594 + (energy-1000.000000)* 0.000004) + \
    1214     (abs(eta) >= 3.7000 && abs(eta) < 3.8000) * (energy >= 10000.0000) * (0.040788*energy/10000.000000) + \
    1215     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.06787795) + \
    1216     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.067878 + (energy-1.000000)* -0.035271) + \
    1217     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.032607 + (energy-2.000000)* -0.006346) + \
    1218     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.013570 + (energy-5.000000)* -0.001453) + \
    1219     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.006304 + (energy-10.000000)* 0.000002) + \
    1220     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.006489 + (energy-100.000000)* 0.000002) + \
    1221     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.008302 + (energy-1000.000000)* 0.000004) + \
    1222     (abs(eta) >= 3.8000 && abs(eta) < 3.9000) * (energy >= 10000.0000) * (0.045331*energy/10000.000000) + \
    1223     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.08769401) + \
    1224     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.087694 + (energy-1.000000)* -0.050535) + \
    1225     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.037159 + (energy-2.000000)* -0.007293) + \
    1226     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.015278 + (energy-5.000000)* -0.001678) + \
    1227     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.006886 + (energy-10.000000)* 0.000003) + \
    1228     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.007116 + (energy-100.000000)* 0.000002) + \
    1229     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.009229 + (energy-1000.000000)* 0.000005) + \
    1230     (abs(eta) >= 3.9000 && abs(eta) < 4.0000) * (energy >= 10000.0000) * (0.051460*energy/10000.000000) + \
    1231     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.12378717) + \
    1232     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.123787 + (energy-1.000000)* -0.081534) + \
    1233     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.042253 + (energy-2.000000)* -0.008504) + \
    1234     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.016742 + (energy-5.000000)* -0.001830) + \
    1235     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.007595 + (energy-10.000000)* 0.000003) + \
    1236     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.007850 + (energy-100.000000)* 0.000003) + \
    1237     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.010277 + (energy-1000.000000)* 0.000005) + \
    1238     (abs(eta) >= 4.0000 && abs(eta) < 4.1000) * (energy >= 10000.0000) * (0.059714*energy/10000.000000) + \
    1239     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.14293546) + \
    1240     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.142935 + (energy-1.000000)* -0.095632) + \
    1241     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.047303 + (energy-2.000000)* -0.009550) + \
    1242     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.018654 + (energy-5.000000)* -0.002035) + \
    1243     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.008480 + (energy-10.000000)* 0.000003) + \
    1244     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.008767 + (energy-100.000000)* 0.000003) + \
    1245     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.011459 + (energy-1000.000000)* 0.000006) + \
    1246     (abs(eta) >= 4.1000 && abs(eta) < 4.2000) * (energy >= 10000.0000) * (0.066014*energy/10000.000000) + \
    1247     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.16175578) + \
    1248     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.161756 + (energy-1.000000)* -0.108957) + \
    1249     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.052799 + (energy-2.000000)* -0.010753) + \
    1250     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.020541 + (energy-5.000000)* -0.002235) + \
    1251     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.009366 + (energy-10.000000)* 0.000004) + \
    1252     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.009686 + (energy-100.000000)* 0.000003) + \
    1253     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.012672 + (energy-1000.000000)* 0.000007) + \
    1254     (abs(eta) >= 4.2000 && abs(eta) < 4.3000) * (energy >= 10000.0000) * (0.072840*energy/10000.000000) + \
    1255     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.17768991) + \
    1256     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.177690 + (energy-1.000000)* -0.119322) + \
    1257     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.058368 + (energy-2.000000)* -0.011801) + \
    1258     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.022966 + (energy-5.000000)* -0.002532) + \
    1259     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.010308 + (energy-10.000000)* 0.000004) + \
    1260     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.010654 + (energy-100.000000)* 0.000004) + \
    1261     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.013982 + (energy-1000.000000)* 0.000007) + \
    1262     (abs(eta) >= 4.3000 && abs(eta) < 4.4000) * (energy >= 10000.0000) * (0.081143*energy/10000.000000) + \
    1263     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.19156432) + \
    1264     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.191564 + (energy-1.000000)* -0.125357) + \
    1265     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.066207 + (energy-2.000000)* -0.013680) + \
    1266     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.025168 + (energy-5.000000)* -0.002755) + \
    1267     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.011394 + (energy-10.000000)* 0.000004) + \
    1268     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.011796 + (energy-100.000000)* 0.000004) + \
    1269     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.015576 + (energy-1000.000000)* 0.000009) + \
    1270     (abs(eta) >= 4.4000 && abs(eta) < 4.5000) * (energy >= 10000.0000) * (0.092193*energy/10000.000000) + \
    1271     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1272     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.075652) + \
    1273     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.075652 + (energy-2.000000)* -0.015621) + \
    1274     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.028790 + (energy-5.000000)* -0.003236) + \
    1275     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.012611 + (energy-10.000000)* 0.000005) + \
    1276     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.013037 + (energy-100.000000)* 0.000005) + \
    1277     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.017269 + (energy-1000.000000)* 0.000010) + \
    1278     (abs(eta) >= 4.5000 && abs(eta) < 4.6000) * (energy >= 10000.0000) * (0.103252*energy/10000.000000) + \
    1279     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1280     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.088768) + \
    1281     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.088768 + (energy-2.000000)* -0.018880) + \
    1282     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.032126 + (energy-5.000000)* -0.003627) + \
    1283     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.013990 + (energy-10.000000)* 0.000005) + \
    1284     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.014468 + (energy-100.000000)* 0.000005) + \
    1285     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.019295 + (energy-1000.000000)* 0.000011) + \
    1286     (abs(eta) >= 4.6000 && abs(eta) < 4.7000) * (energy >= 10000.0000) * (0.117161*energy/10000.000000) + \
    1287     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1288     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.099752) + \
    1289     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.099752 + (energy-2.000000)* -0.021799) + \
    1290     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.034356 + (energy-5.000000)* -0.003793) + \
    1291     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.015389 + (energy-10.000000)* 0.000006) + \
    1292     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.015906 + (energy-100.000000)* 0.000006) + \
    1293     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.021315 + (energy-1000.000000)* 0.000012) + \
    1294     (abs(eta) >= 4.7000 && abs(eta) < 4.8000) * (energy >= 10000.0000) * (0.133463*energy/10000.000000) + \
    1295     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1296     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.129119) + \
    1297     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.129119 + (energy-2.000000)* -0.029895) + \
    1298     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.039433 + (energy-5.000000)* -0.004396) + \
    1299     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.017452 + (energy-10.000000)* 0.000006) + \
    1300     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.018028 + (energy-100.000000)* 0.000007) + \
    1301     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.024225 + (energy-1000.000000)* 0.000014) + \
    1302     (abs(eta) >= 4.8000 && abs(eta) < 4.9000) * (energy >= 10000.0000) * (0.151204*energy/10000.000000) + \
    1303     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1304     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.148587) + \
    1305     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.148587 + (energy-2.000000)* -0.035133) + \
    1306     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.043189 + (energy-5.000000)* -0.004829) + \
    1307     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.019045 + (energy-10.000000)* 0.000006) + \
    1308     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.019625 + (energy-100.000000)* 0.000007) + \
    1309     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.026335 + (energy-1000.000000)* 0.000015) + \
    1310     (abs(eta) >= 4.9000 && abs(eta) < 5.0000) * (energy >= 10000.0000) * (0.165614*energy/10000.000000) + \
    1311     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1312     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.230194) + \
    1313     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 2.0000 && energy < 5.0000) * (0.230194 + (energy-2.000000)* -0.059949) + \
    1314     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 5.0000 && energy < 10.0000) * (0.050347 + (energy-5.000000)* -0.005793) + \
    1315     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 10.0000 && energy < 100.0000) * (0.021380 + (energy-10.000000)* 0.000007) + \
    1316     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 100.0000 && energy < 1000.0000) * (0.022044 + (energy-100.000000)* 0.000009) + \
    1317     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.029961 + (energy-1000.000000)* 0.000018) + \
    1318     (abs(eta) >= 5.0000 && abs(eta) < 5.1000) * (energy >= 10000.0000) * (0.188507*energy/10000.000000) + \
    1319     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1320     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.256226) + \
    1321     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 2.0000 && energy < 5.0000) * (0.256226 + (energy-2.000000)* -0.067581) + \
    1322     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 5.0000 && energy < 10.0000) * (0.053485 + (energy-5.000000)* -0.006077) + \
    1323     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 10.0000 && energy < 100.0000) * (0.023102 + (energy-10.000000)* 0.000008) + \
    1324     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 100.0000 && energy < 1000.0000) * (0.023822 + (energy-100.000000)* 0.000010) + \
    1325     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.032699 + (energy-1000.000000)* 0.000020) + \
    1326     (abs(eta) >= 5.1000 && abs(eta) < 5.2000) * (energy >= 10000.0000) * (0.212432*energy/10000.000000) + \
    1327     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1328     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.428061) + \
    1329     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 2.0000 && energy < 5.0000) * (0.428061 + (energy-2.000000)* -0.122293) + \
    1330     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 5.0000 && energy < 10.0000) * (0.061183 + (energy-5.000000)* -0.007029) + \
    1331     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 10.0000 && energy < 100.0000) * (0.026038 + (energy-10.000000)* 0.000009) + \
    1332     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 100.0000 && energy < 1000.0000) * (0.026834 + (energy-100.000000)* 0.000011) + \
    1333     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.037168 + (energy-1000.000000)* 0.000023) + \
    1334     (abs(eta) >= 5.2000 && abs(eta) < 5.3000) * (energy >= 10000.0000) * (0.240865*energy/10000.000000) + \
    1335     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1336     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.302517) + \
    1337     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 2.0000 && energy < 5.0000) * (0.302517 + (energy-2.000000)* -0.078082) + \
    1338     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 5.0000 && energy < 10.0000) * (0.068271 + (energy-5.000000)* -0.007915) + \
    1339     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 10.0000 && energy < 100.0000) * (0.028698 + (energy-10.000000)* 0.000010) + \
    1340     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 100.0000 && energy < 1000.0000) * (0.029589 + (energy-100.000000)* 0.000014) + \
    1341     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.042319 + (energy-1000.000000)* 0.000027) + \
    1342     (abs(eta) >= 5.3000 && abs(eta) < 5.4000) * (energy >= 10000.0000) * (0.284210*energy/10000.000000) + \
    1343     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1344     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1345     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.025424) + \
    1346     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 5.0000 && energy < 10.0000) * (0.076273 + (energy-5.000000)* -0.008946) + \
    1347     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 10.0000 && energy < 100.0000) * (0.031543 + (energy-10.000000)* 0.000010) + \
    1348     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 100.0000 && energy < 1000.0000) * (0.032442 + (energy-100.000000)* 0.000017) + \
    1349     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.047712 + (energy-1000.000000)* 0.000031) + \
    1350     (abs(eta) >= 5.4000 && abs(eta) < 5.5000) * (energy >= 10000.0000) * (0.326399*energy/10000.000000) + \
    1351     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1352     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1353     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.029080) + \
    1354     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 5.0000 && energy < 10.0000) * (0.087241 + (energy-5.000000)* -0.010397) + \
    1355     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 10.0000 && energy < 100.0000) * (0.035256 + (energy-10.000000)* 0.000011) + \
    1356     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 100.0000 && energy < 1000.0000) * (0.036252 + (energy-100.000000)* 0.000019) + \
    1357     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.053456 + (energy-1000.000000)* 0.000034) + \
    1358     (abs(eta) >= 5.5000 && abs(eta) < 5.6000) * (energy >= 10000.0000) * (0.357698*energy/10000.000000) + \
    1359     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1360     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1361     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.033344) + \
    1362     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 5.0000 && energy < 10.0000) * (0.100031 + (energy-5.000000)* -0.012242) + \
    1363     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 10.0000 && energy < 100.0000) * (0.038819 + (energy-10.000000)* 0.000015) + \
    1364     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 100.0000 && energy < 1000.0000) * (0.040162 + (energy-100.000000)* 0.000023) + \
    1365     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.061269 + (energy-1000.000000)* 0.000040) + \
    1366     (abs(eta) >= 5.6000 && abs(eta) < 5.7000) * (energy >= 10000.0000) * (0.420216*energy/10000.000000) + \
    1367     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1368     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1369     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.039257) + \
    1370     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 5.0000 && energy < 10.0000) * (0.117770 + (energy-5.000000)* -0.014993) + \
    1371     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 10.0000 && energy < 100.0000) * (0.042808 + (energy-10.000000)* 0.000025) + \
    1372     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 100.0000 && energy < 1000.0000) * (0.045061 + (energy-100.000000)* 0.000029) + \
    1373     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.071538 + (energy-1000.000000)* 0.000049) + \
    1374     (abs(eta) >= 5.7000 && abs(eta) < 5.8000) * (energy >= 10000.0000) * (0.510223*energy/10000.000000) + \
    1375     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1376     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1377     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.041967) + \
    1378     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 5.0000 && energy < 10.0000) * (0.125902 + (energy-5.000000)* -0.016273) + \
    1379     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 10.0000 && energy < 100.0000) * (0.044536 + (energy-10.000000)* 0.000023) + \
    1380     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 100.0000 && energy < 1000.0000) * (0.046571 + (energy-100.000000)* 0.000043) + \
    1381     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.085014 + (energy-1000.000000)* 0.000066) + \
    1382     (abs(eta) >= 5.8000 && abs(eta) < 5.9000) * (energy >= 10000.0000) * (0.680639*energy/10000.000000) + \
    1383     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 0.0000 && energy < 1.0000) * (0.00000000) + \
    1384     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 1.0000 && energy < 2.0000) * (0.000000 + (energy-1.000000)* 0.000000) + \
    1385     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 2.0000 && energy < 5.0000) * (0.000000 + (energy-2.000000)* 0.000000) + \
    1386     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 5.0000 && energy < 10.0000) * (0.000000 + (energy-5.000000)* 0.013098) + \
    1387     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 10.0000 && energy < 100.0000) * (0.065492 + (energy-10.000000)* 0.000107) + \
    1388     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 100.0000 && energy < 1000.0000) * (0.075097 + (energy-100.000000)* 0.000243) + \
    1389     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 1000.0000 && energy < 10000.0000) * (0.294198 + (energy-1000.000000)* 0.000280) + \
    1390     (abs(eta) >= 5.9000 && abs(eta) < 6.0000) * (energy >= 10000.0000) * (2.814894*energy/10000.000000)
    1391  }
    1392 }
    1393 
     562        (abs(eta) < 2.50 && abs(eta) >= 1.74 )   * 2 * sqrt(4.57439e-05^2 * pt^2   + 0.0149328^2           ) +
     563        (abs(eta) < 1.74 && abs(eta) >= 1.32 )   * 1.5 * sqrt(4.57439e-05^2 * pt^2   + 0.0149328^2         ) +
     564        (abs(eta) < 1.32 && abs(eta) >= 0.76 )   * sqrt(9.81626e-06^2 * pt^2   + 0.00379895^2  ) +
     565        (abs(eta) < 0.76 && abs(eta) >= 0.36 )   * sqrt(1.1959e-05^2 * pt^2   +  0.00242417^2 ) +
     566        (abs(eta) < 0.36 && abs(eta) >= 0.18 )   * sqrt(1.20149e-05^2 * pt^2  + 0.00219291^2  ) +
     567        (abs(eta) < 0.18)                        * sqrt(1.29686e-05^2 * pt^2  + 0.0020392^2      )
     568
     569    }
     570}
    1394571
    1395572##############
     
    19981175    add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 VLCjetR15N6 Jet
    19991176
    2000     add Branch JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive VLCjetR02_inclusive Jet
    20011177    add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive VLCjetR05_inclusive Jet
    20021178    add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive VLCjetR07_inclusive Jet
     
    20201196    add Branch UniqueObjectFinder/electrons Electron Electron
    20211197    add Branch UniqueObjectFinder/muons Muon Muon
    2022     add Branch ForwardMuonMomentumSmearing/muons ForwardMuon Muon
    20231198
    20241199    add Branch MissingET/momentum MissingET MissingET
Note: See TracChangeset for help on using the changeset viewer.