Fork me on GitHub

Changeset 934d037 in git for cards/delphes_card_ATLAS.tcl


Ignore:
Timestamp:
Oct 6, 2015, 4:26:29 PM (9 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
c4b2aae
Parents:
b3c0535
Message:

replace ElectronEnergySmearing with ElectronMomentumSmearing and add ElectronFilter

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_ATLAS.tcl

    rb3c0535 r934d037  
    1111
    1212  ChargedHadronMomentumSmearing
    13   ElectronEnergySmearing
     13  ElectronMomentumSmearing
    1414  MuonMomentumSmearing
    1515
    … …  
    2121  PhotonIsolation
    2222
     23  ElectronFilter
    2324  ElectronEfficiency
    2425  ElectronIsolation
    … …  
    138139
    139140  # resolution formula for charged hadrons
    140   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
    141                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
    142                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
    143                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    145                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
    146                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    147                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    148 }
    149 
    150 #################################
    151 # Energy resolution for electrons
    152 #################################
    153 
    154 module EnergySmearing ElectronEnergySmearing {
     141  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     142                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     143                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     144}
     145
     146###################################
     147# Momentum resolution for electrons
     148###################################
     149
     150module MomentumSmearing ElectronMomentumSmearing {
    155151  set InputArray ElectronTrackingEfficiency/electrons
    156152  set OutputArray electrons
    … …  
    158154  # set ResolutionFormula {resolution formula as a function of eta and energy}
    159155
    160   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) +
    161                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    162                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    163                          (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    164 
     156  # resolution formula for electrons
     157  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     158                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     159                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    165160}
    166161
    … …  
    176171
    177172  # resolution formula for muons
    178   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    179                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) +
    180                                            (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) +
    181                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) +
    182                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) +
    183                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) +
    184                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) +
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
     173  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*2.0e-4^2) +
     174                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*3.0e-4^2) +
     175                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*6.0e-4^2)}
    186176}
    187177
    … …  
    193183# add InputArray InputArray
    194184  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    195   add InputArray ElectronEnergySmearing/electrons
     185  add InputArray ElectronMomentumSmearing/electrons
    196186  add InputArray MuonMomentumSmearing/muons
    197187  set OutputArray tracks
    … …  
    326316}
    327317
     318#################
     319# Electron filter
     320#################
     321
     322module PdgCodeFilter ElectronFilter {
     323  set InputArray Calorimeter/eflowTracks
     324  set OutputArray electrons
     325  set Invert true
     326  add PdgCode {11}
     327  add PdgCode {-11}
     328}
     329
    328330#####################
    329331# Electron efficiency
    … …  
    331333
    332334module Efficiency ElectronEfficiency {
    333   set InputArray ElectronEnergySmearing/electrons
     335  set InputArray ElectronFilter/electrons
    334336  set OutputArray electrons
    335337
Note: See TracChangeset for help on using the changeset viewer.