Fork me on GitHub

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_PileUp.tcl

    rb3c0535 r934d037  
    1313
    1414  ChargedHadronMomentumSmearing
    15   ElectronEnergySmearing
     15  ElectronMomentumSmearing
    1616  MuonMomentumSmearing
    1717
    1818  TrackMerger
    1919  Calorimeter
     20  ElectronFilter
    2021  TrackPileUpSubtractor
    2122  NeutralTowerMerger
     
    180181
    181182  # resolution formula for charged hadrons
    182   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
    183                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
    184                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
    185                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
    188                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    189                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    190 }
    191 
    192 #################################
    193 # Energy resolution for electrons
    194 #################################
    195 
    196 module EnergySmearing ElectronEnergySmearing {
     183  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     184                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     185                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     186}
     187
     188###################################
     189# Momentum resolution for electrons
     190###################################
     191
     192module MomentumSmearing ElectronMomentumSmearing {
    197193  set InputArray ElectronTrackingEfficiency/electrons
    198194  set OutputArray electrons
     
    200196  # set ResolutionFormula {resolution formula as a function of eta and energy}
    201197
    202   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) +
    203                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    204                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    205                          (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    206 
     198  # resolution formula for electrons
     199  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     200                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     201                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    207202}
    208203
     
    218213
    219214  # resolution formula for muons
    220   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    221                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) +
    222                                            (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) +
    223                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) +
    224                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) +
    225                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) +
    226                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) +
    227                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
     215  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*2.0e-4^2) +
     216                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*3.0e-4^2) +
     217                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
    228218}
    229219
     
    235225# add InputArray InputArray
    236226  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    237   add InputArray ElectronEnergySmearing/electrons
     227  add InputArray ElectronMomentumSmearing/electrons
    238228  add InputArray MuonMomentumSmearing/muons
    239229  set OutputArray tracks
     
    322312}
    323313
     314#################
     315# Electron filter
     316#################
     317
     318module PdgCodeFilter ElectronFilter {
     319  set InputArray Calorimeter/eflowTracks
     320  set OutputArray electrons
     321  set Invert true
     322  add PdgCode {11}
     323  add PdgCode {-11}
     324}
     325
    324326##########################
    325327# Track pile-up subtractor
     
    329331# add InputArray InputArray OutputArray
    330332  add InputArray Calorimeter/eflowTracks eflowTracks
    331   add InputArray ElectronEnergySmearing/electrons electrons
     333  add InputArray ElectronFilter/electrons electrons
    332334  add InputArray MuonMomentumSmearing/muons muons
    333335
Note: See TracChangeset for help on using the changeset viewer.