Fork me on GitHub

Ignore:
Timestamp:
Jan 15, 2018, 11:36:43 AM (7 years ago)
Author:
Ulrike Schnoor <schnooru@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
1f4c2c3, cd60feb
Parents:
5a69ba1d
Message:

added tracking efficiency numbers and momentum resolutions for e, mu, pions

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CLICdet.tcl

    r5a69ba1d r3037765  
    8383
    8484    JetFlavorAssociation
    85 
    8685    BTagging
    8786   
    … …  
    123122    set InputArray ParticlePropagator/chargedHadrons
    124123    set OutputArray chargedHadrons
    125 
    126     # FIXME currently uses tracking efficiency from muon simulation
     124    # Current full simulation with CLICdet provides for pions:
     125
    127126    set EfficiencyFormula {
    128         (abs(eta) > 2.66)                                               * (0.000) +
    129         (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1)               * (0.997) +
    130         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
    131         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174)             * (1.000) +
    132         (abs(eta) < 2.25)                    * (pt > 0.1)               * (1.000)                     }
     127(abs(eta) > 2.54) * (0.000) +
     128(energy >= 100) * (abs(eta) < 2.54)  * (1.000) +
     129(energy < 100 && energy >= 10) * (abs(eta) <=2.54 && abs(eta) > 2.34)  * (0.994) +
     130(energy < 100 && energy >= 10) * (abs(eta) <= 2.34) * (1.000) +
     131(energy < 10  && energy >= 1) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) +
     132(energy < 10  && energy >= 1) * (abs(eta) <= 0.55 ) * (1.000)
     133    }
    133134}
    134135
    … …  
    141142    set OutputArray electrons
    142143
    143     # FIXME currently uses tracking efficiency from muon simulation
     144
     145    # Current full simulation with CLICdet provides for electrons:
    144146    set EfficiencyFormula {
    145         (abs(eta) > 2.66)                                               * (0.000) +
    146         (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1)               * (0.997) +
    147         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
    148         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174)             * (1.000) +
    149         (abs(eta) < 2.25)                    * (pt > 0.1)               * (1.000)
     147                (pt <= 1)                                                               * (0.000) +
     148                (abs(eta) > 2.54)                                                       * (0.000) +
     149                (energy > 100)                 * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
     150                (energy > 100)                 * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
     151                (energy > 100)                 * (abs(eta) <= 2.34  )                   * (1.000) +
     152                (energy <= 100 && energy > 10) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
     153                (energy <= 100 && energy > 10) * (abs(eta) <= 2.17)                     * (1.000) +
     154                (energy <= 10 && energy > 1)   * (abs(eta) > 2.34 )                     * (0.000) +
     155                (energy <= 10 && energy > 1)   * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
     156                (energy <= 10 && energy > 1)   * (abs(eta) <= 0.76)                     * (0.999)
    150157     }
    151158}
    … …  
    161168    # Current full simulation with CLICdet provides for muons:
    162169    set EfficiencyFormula {
    163         (abs(eta) > 2.66)                                               * (0.000) +
    164         (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1)               * (0.997) +
    165         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
    166         (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174)             * (1.000) +
    167         (abs(eta) < 2.25)                    * (pt > 0.1)               * (1.000)
    168      }
     170        (pt < 1) * (0.000)+
     171        (abs(eta) > 2.54) * (0.0000) +
     172        (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 100)                * (0.994) +
     173        (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 10 && energy < 100) * (0.997) +
     174        (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 1  && energy < 10)  * (0.996) +
     175        (abs(eta) <= 2.44 )                    * (energy >= 10)                 * (1.000) +
     176        (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy >= 1 && energy < 10)   * (0.999) +
     177        (abs(eta) <= 2.25 )                    * (energy >= 1)                  * (1.000)
     178    }
    169179}
    170180
    … …  
    177187    set OutputArray chargedHadrons
    178188
    179     # Using eta mid-points between evaluated resolutions from full simulation
     189
    180190    # Resolution given in dpT/pT.
    181     # FIXME: currently uses momentum resolution from single muon simulation!
    182191    set ResolutionFormula {
    183                            (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
    184         (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
    185         (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
    186         (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
    187         (abs(eta) > 2.75)                     * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
     192        (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2    ) +
     193        (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2   ) +
     194        (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2    ) +
     195        (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2   ) +
     196        (abs(eta) < 0.18)                      * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2   )
     197       
    188198    }
    189199}
    … …  
    197207    set OutputArray electrons
    198208
    199     # Using eta mid-points between evaluated resolutions from full simulation
    200209    # Resolution given in dpT/pT.
    201     # FIXME: currently uses momentum resolution from single muon simulation!
    202210    set ResolutionFormula {
    203                            (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
    204         (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
    205         (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
    206         (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
    207         (abs(eta) > 2.75)                     * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
     211        (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt( 8.62283e-05^2 * pt^2  + 0.0177556^2   ) +
     212        (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt( 1.0915e-05 ^2 * pt^2  + 0.00663766^2  ) +
     213        (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt( 1.15518e-05^2 * pt^2  + 0.00398644^2  ) +
     214        (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2  + 0.00317807^2  ) +
     215        (abs(eta) < 0.18)                      * sqrt( 1.40722e-05^2 * pt^2  + 0.00292138^2  )
     216
    208217    }
    209218}
    … …  
    217226    set OutputArray muons
    218227
    219     # Using eta mid-points between evaluated resolutions from full simulation
    220228    # Resolution given in dpT/pT.
    221229    set ResolutionFormula {
    222                            (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
    223         (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
    224         (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
    225         (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
    226         (abs(eta) > 2.75)                     * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
     230
     231        (abs(eta) < 2.66 && abs(eta) >= 2.03 ) * sqrt(4.57439e-05^2 * pt^2*   + 0.0149328^2        ) +
     232        (abs(eta) < 2.03 && abs(eta) >= 1.01 ) * sqrt(9.81626e-06^2 * pt^2*   + 0.00379895^2  ) +
     233        (abs(eta) < 1.01 && abs(eta) >= 0.55 ) * sqrt(1.1959e-05^2 * pt^2*   +  0.00242417^2 ) +
     234        (abs(eta) < 0.55 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2  + 0.00219291^2  ) +
     235        (abs(eta) < 0.18)                      * sqrt(1.29686e-05^2 * pt^2  + 0.0020392^2      )
     236
    227237    }
    228238}
    … …  
    10581068    # gluon's PDG code has the lowest priority
    10591069
    1060     # based on arXiv:1211.4462
     1070    # based on arXiv:1211.4462 (CMS)
    10611071   
    10621072    # default efficiency formula (misidentification rate)
Note: See TracChangeset for help on using the changeset viewer.