Fork me on GitHub

Changeset 679a87e in git


Ignore:
Timestamp:
Oct 9, 2015, 12:16:13 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
ad71e62
Parents:
221abf3
Message:

put sqrt in mom. reso. formulae

Location:
cards
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • cards/delphes_card_FCC.tcl

    r221abf3 r679a87e  
    145145
    146146  # resolution formula for charged hadrons
    147   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    148                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5) +
    149                              (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*5.e-5)}
     147  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.e-5^2) +
     148                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.e-5^2) +
     149                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * sqrt(0.03^2 + pt^2*5.e-5^2)}
    150150
    151151
     
    163163
    164164  # resolution formula for electrons
    165   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    166                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5) +
    167                              (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*5.e-5)}
     165 set ResolutionFormula {     (abs(eta) <= 1.5)                   * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.e-5^2) +
     166                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.e-5^2) +
     167                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * sqrt(0.03^2 + pt^2*5.e-5^2)}
    168168
    169169}
     
    180180
    181181  # resolution formula for muons
    182   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.0e-6) +
    183                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.0e-5) +
    184                              (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*1.5e-5)}
     182  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * sqrt(0.01^2 + pt^2*5.0e-6^2) +
     183                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * sqrt(0.02^2 + pt^2*1.0e-5^2) +
     184                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.5e-5^2)}
    185185
    186186
  • cards/delphes_card_ILD.tcl

    r221abf3 r679a87e  
    124124
    125125  # resolution formula for charged hadrons
    126   set ResolutionFormula {    (abs(eta) <= 1.0)                   * (0.001 + pt*1.e-5) +
    127                              (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
     126  set ResolutionFormula {    (abs(eta) <= 1.0)                   * sqrt(0.001^2 + pt^2*1.e-5^2) +
     127                             (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
    128128
    129129
     
    140140  # set ResolutionFormula {resolution formula as a function of eta and energy}
    141141
    142   # resolution formula for electrons (we keep nominal tracking reso since ECAL can't beat this, even at high E)
    143   set ResolutionFormula {    (abs(eta) <= 1.0)                   * (0.001 + pt*1.e-5) +
    144                              (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
     142   # resolution formula for charged hadrons
     143  set ResolutionFormula {    (abs(eta) <= 1.0)                   * sqrt(0.001^2 + pt^2*1.e-5^2) +
     144                             (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
    145145}
    146146
     
    155155  # set ResolutionFormula {resolution formula as a function of eta and pt}
    156156
    157   # resolution formula for muons
    158   set ResolutionFormula {    (abs(eta) <= 1.0)                   * (0.001 + pt*1.e-5) +
    159                              (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
    160 
     157   # resolution formula for charged hadrons
     158  set ResolutionFormula {    (abs(eta) <= 1.0)                   * sqrt(0.001^2 + pt^2*1.e-5^2) +
     159                             (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
    161160
    162161}
Note: See TracChangeset for help on using the changeset viewer.