Fork me on GitHub

Changeset 934d037 in git


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

Location:
cards
Files:
8 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

  • 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

  • cards/delphes_card_CMS.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  # based on arXiv:1405.6569
     142  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     143                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     145}
     146
     147###################################
     148# Momentum resolution for electrons
     149###################################
     150
     151module MomentumSmearing ElectronMomentumSmearing {
    155152  set InputArray ElectronTrackingEfficiency/electrons
    156153  set OutputArray electrons
    … …  
    159156
    160157  # resolution formula for electrons
    161   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
    162                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    163                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    164                          (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    165 
     158  # based on arXiv:1405.6569
     159  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     160                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     161                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    166162}
    167163
    … …  
    177173
    178174  # resolution formula for muons
    179   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
    180                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
    181                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
    182                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    183                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
    184                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
    185                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
    186                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
    188                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
    189                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    190                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
     175  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
     176                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
     177                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
    191178}
    192179
    … …  
    198185# add InputArray InputArray
    199186  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    200   add InputArray ElectronEnergySmearing/electrons
     187  add InputArray ElectronMomentumSmearing/electrons
    201188  add InputArray MuonMomentumSmearing/muons
    202189  set OutputArray tracks
    … …  
    334321}
    335322
     323#################
     324# Electron filter
     325#################
     326
     327module PdgCodeFilter ElectronFilter {
     328  set InputArray Calorimeter/eflowTracks
     329  set OutputArray electrons
     330  set Invert true
     331  add PdgCode {11}
     332  add PdgCode {-11}
     333}
     334
    336335#####################
    337336# Electron efficiency
    … …  
    339338
    340339module Efficiency ElectronEfficiency {
    341   set InputArray ElectronEnergySmearing/electrons
     340  set InputArray ElectronFilter/electrons
    342341  set OutputArray electrons
    343342
    … …  
    500499
    501500module JetFlavorAssociation JetFlavorAssociation {
    502  
     501
    503502  set PartonInputArray Delphes/partons
    504503  set ParticleInputArray Delphes/allParticles
    505504  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    506505  set JetInputArray JetEnergyScale/jets
    507  
     506
    508507  set DeltaR 0.5
    509508  set PartonPTMin 1.0

  • cards/delphes_card_CMS_NoFastJet.tcl

    rb3c0535 r934d037  
    1111
    1212  ChargedHadronMomentumSmearing
    13   ElectronEnergySmearing
     13  ElectronMomentumSmearing
    1414  MuonMomentumSmearing
    1515
    … …  
    110110
    111111  # resolution formula for charged hadrons
    112   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
    113                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
    114                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
    115                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
    116                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    117                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
    118                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    119                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    120 }
    121 
    122 #################################
    123 # Energy resolution for electrons
    124 #################################
    125 
    126 module EnergySmearing ElectronEnergySmearing {
     112  # based on arXiv:1405.6569
     113  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     114                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     115                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     116}
     117
     118###################################
     119# Momentum resolution for electrons
     120###################################
     121
     122module MomentumSmearing ElectronMomentumSmearing {
    127123  set InputArray ElectronTrackingEfficiency/electrons
    128124  set OutputArray electrons
    … …  
    131127
    132128  # resolution formula for electrons
    133   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
    134                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    135                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    136                          (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    137 
     129  # based on arXiv:1405.6569
     130  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     131                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     132                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    138133}
    139134
    … …  
    149144
    150145  # resolution formula for muons
    151   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
    152                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
    153                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
    154                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    155                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
    156                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
    157                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
    158                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    159                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
    160                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
    161                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    162                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
     146  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
     147                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
     148                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
    163149}
    164150
    … …  
    170156# add InputArray InputArray
    171157  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    172   add InputArray ElectronEnergySmearing/electrons
     158  add InputArray ElectronMomentumSmearing/electrons
    173159  add InputArray MuonMomentumSmearing/muons
    174160  set OutputArray tracks

  • cards/delphes_card_CMS_PileUp.tcl

    rb3c0535 r934d037  
    1313
    1414  ChargedHadronMomentumSmearing
    15   ElectronEnergySmearing
     15  ElectronMomentumSmearing
    1616  MuonMomentumSmearing
    1717
    1818  TrackMerger
    1919  Calorimeter
     20  ElectronFilter
    2021  TrackPileUpSubtractor
    2122  NeutralTowerMerger
    … …  
    181182
    182183  # resolution formula for charged hadrons
    183   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
    184                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
    185                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
    186                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
    188                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
    189                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    190                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    191 }
    192 
    193 #################################
    194 # Energy resolution for electrons
    195 #################################
    196 
    197 module EnergySmearing ElectronEnergySmearing {
     184  # based on arXiv:1405.6569
     185  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     186                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     188}
     189
     190###################################
     191# Momentum resolution for electrons
     192###################################
     193
     194module MomentumSmearing ElectronMomentumSmearing {
    198195  set InputArray ElectronTrackingEfficiency/electrons
    199196  set OutputArray electrons
    … …  
    202199
    203200  # resolution formula for electrons
    204   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
    205                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    206                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    207                          (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    208 
     201  # based on arXiv:1405.6569
     202  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     203                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     204                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    209205}
    210206
    … …  
    220216
    221217  # resolution formula for muons
    222   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
    223                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
    224                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
    225                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    226                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
    227                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
    228                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
    229                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
    230                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
    231                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
    232                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    233                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
     218  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
     219                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
     220                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
    234221}
    235222
    … …  
    241228# add InputArray InputArray
    242229  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    243   add InputArray ElectronEnergySmearing/electrons
     230  add InputArray ElectronMomentumSmearing/electrons
    244231  add InputArray MuonMomentumSmearing/muons
    245232  set OutputArray tracks
    … …  
    331318}
    332319
     320#################
     321# Electron filter
     322#################
     323
     324module PdgCodeFilter ElectronFilter {
     325  set InputArray Calorimeter/eflowTracks
     326  set OutputArray electrons
     327  set Invert true
     328  add PdgCode {11}
     329  add PdgCode {-11}
     330}
     331
    333332##########################
    334333# Track pile-up subtractor
    … …  
    338337# add InputArray InputArray OutputArray
    339338  add InputArray Calorimeter/eflowTracks eflowTracks
    340   add InputArray ElectronEnergySmearing/electrons electrons
     339  add InputArray ElectronFilter/electrons electrons
    341340  add InputArray MuonMomentumSmearing/muons muons
    342341

  • cards/delphes_card_FCC_basic.tcl

    rb3c0535 r934d037  
    1111
    1212  ChargedHadronMomentumSmearing
    13   ElectronEnergySmearing
     13  ElectronMomentumSmearing
    1414  MuonMomentumSmearing
    1515
    … …  
    2020  ECal
    2121  HCal
     22
     23  ElectronFilter
    2224
    2325  TowerMerger
    … …  
    136138}
    137139
    138 #################################
    139 # Energy resolution for electrons
    140 #################################
    141 
    142 module EnergySmearing ElectronEnergySmearing {
     140###################################
     141# Momentum resolution for electrons
     142###################################
     143
     144module MomentumSmearing ElectronMomentumSmearing {
    143145  set InputArray ElectronTrackingEfficiency/electrons
    144146  set OutputArray electrons
    … …  
    147149
    148150  # resolution formula for electrons
    149   set ResolutionFormula {                  (abs(eta) <= 4.0) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.007) +
    150                                            (abs(eta) <= 4.0) * (energy > 2.0e1)                    * sqrt(energy^2*0.005^2 + energy*0.02^2) +
    151                                            (abs(eta) > 4.0 && abs(eta) <= 6.0)                     * sqrt(energy^2*0.05^2 + energy*1.00^2)}
    152 
     151  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
     152                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
    153153}
    154154
    … …  
    177177# add InputArray InputArray
    178178  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    179   add InputArray ElectronEnergySmearing/electrons
     179  add InputArray ElectronMomentumSmearing/electrons
    180180  add InputArray MuonMomentumSmearing/muons
    181181  set OutputArray tracks
    … …  
    183183
    184184
    185 ################################
     185########################
    186186# Track angular smearing
    187 ################################
     187########################
    188188
    189189module AngularSmearing AngularSmearing {
    … …  
    200200}
    201201
    202 ################################
     202#################################
    203203# Track impact parameter smearing
    204 ################################
     204#################################
    205205
    206206module ImpactParameterSmearing ImpactParameterSmearing {
    … …  
    350350}
    351351
    352 ####################
     352#################
     353# Electron filter
     354#################
     355
     356module PdgCodeFilter ElectronFilter {
     357  set InputArray Calorimeter/eflowTracks
     358  set OutputArray electrons
     359  set Invert true
     360  add PdgCode {11}
     361  add PdgCode {-11}
     362}
     363
     364###################################################
    353365# Tower Merger (in case not using e-flow algorithm)
    354 ####################
     366###################################################
    355367
    356368module Merger TowerMerger {
    … …  
    395407}
    396408
    397 #####################
     409#################
    398410# Neutrino Filter
    399 #####################
     411#################
    400412
    401413module PdgCodeFilter NeutrinoFilter {
    … …  
    488500
    489501
    490 ##########################
     502#############
    491503# tau-tagging
    492 ##########################
     504#############
    493505
    494506
    … …  
    525537  add Branch ECal/eflowPhotons Photon Photon
    526538
    527   add Branch ElectronEnergySmearing/electrons Electron Electron
     539  add Branch ElectronFilter/electrons Electron Electron
    528540  add Branch MuonMomentumSmearing/muons Muon Muon
    529541  add Branch JetEnergyScale/jets Jet Jet

  • cards/delphes_card_ILD.tcl

    rb3c0535 r934d037  
    1313
    1414  ChargedHadronMomentumSmearing
    15   ElectronEnergySmearing
     15  ElectronMomentumSmearing
    1616  MuonMomentumSmearing
    1717
    … …  
    2222  ECal
    2323  HCal
     24
     25  ElectronFilter
    2426
    2527  TowerMerger
    … …  
    128130}
    129131
    130 #################################
    131 # Energy resolution for electrons
    132 #################################
    133 
    134 module EnergySmearing ElectronEnergySmearing {
     132###################################
     133# Momentum resolution for electrons
     134###################################
     135
     136module MomentumSmearing ElectronMomentumSmearing {
    135137  set InputArray ElectronTrackingEfficiency/electrons
    136138  set OutputArray electrons
    … …  
    167169# add InputArray InputArray
    168170  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    169   add InputArray ElectronEnergySmearing/electrons
     171  add InputArray ElectronMomentumSmearing/electrons
    170172  add InputArray MuonMomentumSmearing/muons
    171173  set OutputArray tracks
    … …  
    173175
    174176
    175 ################################
     177########################
    176178# Track angular smearing
    177 ################################
     179########################
    178180
    179181module AngularSmearing AngularSmearing {
    … …  
    190192}
    191193
    192 ################################
     194#################################
    193195# Track impact parameter smearing
    194 ################################
     196#################################
    195197
    196198module ImpactParameterSmearing ImpactParameterSmearing {
    … …  
    334336}
    335337
    336 ####################
     338#################
     339# Electron filter
     340#################
     341
     342module PdgCodeFilter ElectronFilter {
     343  set InputArray Calorimeter/eflowTracks
     344  set OutputArray electrons
     345  set Invert true
     346  add PdgCode {11}
     347  add PdgCode {-11}
     348}
     349
     350###################################################
    337351# Tower Merger (in case not using e-flow algorithm)
    338 ####################
     352###################################################
    339353
    340354module Merger TowerMerger {
    … …  
    379393}
    380394
    381 #####################
     395#################
    382396# Neutrino Filter
    383 #####################
     397#################
    384398
    385399module PdgCodeFilter NeutrinoFilter {
    … …  
    472486
    473487
    474 ##########################
     488#############
    475489# tau-tagging
    476 ##########################
     490#############
    477491
    478492
    … …  
    509523  add Branch ECal/eflowPhotons Photon Photon
    510524
    511   add Branch ElectronEnergySmearing/electrons Electron Electron
     525  add Branch ElectronFilter/electrons Electron Electron
    512526  add Branch MuonMomentumSmearing/muons Muon Muon
    513527  add Branch JetEnergyScale/jets Jet Jet

  • cards/delphes_card_LHCb.tcl

    rb3c0535 r934d037  
    1313
    1414  ChargedHadronMomentumSmearing
    15   ElectronEnergySmearing
     15  ElectronMomentumSmearing
    1616  MuonMomentumSmearing
    1717
    … …  
    111111}
    112112
    113 #################################
    114 # Energy resolution for electrons
    115 #################################
    116 
    117 module EnergySmearing ElectronEnergySmearing {
     113###################################
     114# Momentum resolution for electrons
     115###################################
     116
     117module MomentumSmearing ElectronMomentumSmearing {
    118118  set InputArray ParticlePropagator/electrons
    119119  set OutputArray electrons
    … …  
    122122
    123123  # resolution formula for electrons
    124   set ResolutionFormula { (eta > 2.0  && eta <= 5.0) * (energy > 0.1   && energy <= 8.0) * (energy*0.05) +
    125                           (eta > 2.0  && eta <= 5.0) * (energy > 8.0)                    *  sqrt(energy^2*0.015^2 + energy*0.10^2)}
    126   }
     124  set ResolutionFormula {(eta > 2.0  && eta <= 5.0)      * (pt > 0.5)* (0.005)}
     125}
    127126
    128127###############################
    … …  
    146145# add InputArray InputArray
    147146  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    148   add InputArray ElectronEnergySmearing/electrons
     147  add InputArray ElectronMomentumSmearing/electrons
    149148  add InputArray MuonMomentumSmearing/muons
    150149  set OutputArray tracks
Note: See TracChangeset for help on using the changeset viewer.