Fork me on GitHub

Changes in / [1f34dac:0e2f49b] in git


Ignore:
Location:
cards/CMS_PhaseII
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl

    r1f34dac r0e2f49b  
    2121
    2222  ChargedHadronMomentumSmearing
    23   ElectronEnergySmearing
     23  ElectronMomentumSmearing
    2424  MuonMomentumSmearing
    2525
     
    168168}
    169169
    170 #################################
    171 # Energy resolution for electrons
    172 #################################
    173 
    174 module EnergySmearing ElectronEnergySmearing {
    175   set InputArray ElectronTrackingEfficiency/electrons
     170########################################
     171# Momentum resolution for electrons
     172########################################
     173
     174module MomentumSmearing ElectronMomentumSmearing {
     175  ## hadrons after having applied the tracking efficiency
     176  set InputArray  ElectronTrackingEfficiency/electrons
    176177  set OutputArray electrons
    177 
    178   # set ResolutionFormula {resolution formula as a function of eta and energy}
    179 
    180178  # resolution formula for electrons
    181  
    182   # taking something flat in energy for now, ECAL will take over at high energy anyway.
    183   # inferred from hep-ex/1306.2016 and 1502.02701
    184   set ResolutionFormula { 
    185  
    186                         (abs(eta) <= 1.5)  * (1+0.64*abs(eta)^2)*(energy*0.028) +
    187     (abs(eta) > 1.5  && abs(eta) <= 1.75)  * (energy*0.037) +
    188     (abs(eta) > 1.75  && abs(eta) <= 2.15) * (energy*0.038) +
    189     (abs(eta) > 2.15  && abs(eta) <= 3.00) * (energy*0.044) +
    190     (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}   
    191 
    192 }
     179  # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
     180  source trackMomentumResolution.tcl 
     181
     182}
     183
    193184
    194185###############################
     
    215206# add InputArray InputArray
    216207  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    217   add InputArray ElectronEnergySmearing/electrons
     208  add InputArray ElectronMomentumSmearing/electrons
    218209  add InputArray MuonMomentumSmearing/muons
    219210  set OutputArray tracks
     
    310301  # set ResolutionFormula {resolution formula as a function of eta and energy}
    311302 
    312   # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
     303  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
    313304  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    314305
    315   set ResolutionFormula {  (abs(eta) <= 1.50)                    * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
     306  set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
    316307                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    317308                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    318309                           (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
    319                            (abs(eta) >= 3.0 && abs(eta) <= 5.0)  * sqrt(energy^2*0.08^2 + energy*1.98^2)}
     310                           (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
     311
    320312}
    321313
  • cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl

    r1f34dac r0e2f49b  
    2424
    2525  ChargedHadronMomentumSmearing
    26   ElectronEnergySmearing
     26  ElectronMomentumSmearing
    2727  MuonMomentumSmearing
    2828
     
    218218}
    219219
    220 #################################
    221 # Energy resolution for electrons
    222 #################################
    223 
    224 module EnergySmearing ElectronEnergySmearing {
    225   set InputArray ElectronTrackingEfficiency/electrons
     220########################################
     221# Momentum resolution for electrons
     222########################################
     223
     224module MomentumSmearing ElectronMomentumSmearing {
     225  ## hadrons after having applied the tracking efficiency
     226  set InputArray  ElectronTrackingEfficiency/electrons
    226227  set OutputArray electrons
    227 
    228   # set ResolutionFormula {resolution formula as a function of eta and energy}
    229 
    230228  # resolution formula for electrons
    231  
    232   # taking something flat in energy for now, ECAL will take over at high energy anyway.
    233   # inferred from hep-ex/1306.2016 and 1502.02701
    234   set ResolutionFormula { 
    235  
    236                         (abs(eta) <= 1.5)  * (1+0.64*abs(eta)^2)*(energy*0.028) +
    237     (abs(eta) > 1.5  && abs(eta) <= 1.75)  * (energy*0.037) +
    238     (abs(eta) > 1.75  && abs(eta) <= 2.15) * (energy*0.038) +
    239     (abs(eta) > 2.15  && abs(eta) <= 3.00) * (energy*0.044) +
    240     (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}   
    241 
    242 }
     229  # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
     230  source trackMomentumResolution.tcl 
     231
     232}
     233
    243234
    244235###############################
     
    267258# add InputArray InputArray
    268259  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    269   add InputArray ElectronEnergySmearing/electrons
     260  add InputArray ElectronMomentumSmearing/electrons
    270261  add InputArray MuonMomentumSmearing/muons
    271262  set OutputArray tracks
     
    362353  # set ResolutionFormula {resolution formula as a function of eta and energy}
    363354 
    364   # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
     355  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
    365356  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    366357
    367   set ResolutionFormula {  (abs(eta) <= 1.50)                    * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
     358  set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
    368359                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    369360                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    370361                           (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
    371                            (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
     362                           (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
    372363
    373364}
Note: See TracChangeset for help on using the changeset viewer.