Fork me on GitHub

Changeset df7986a in git for cards


Ignore:
Timestamp:
May 16, 2016, 8:09:01 PM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
db7f5e4
Parents:
4bf6061
Message:

fixed elecrtron resolution

Location:
cards/CMS_PhaseII
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl

    r4bf6061 rdf7986a  
    2121
    2222  ChargedHadronMomentumSmearing
    23   ElectronMomentumSmearing
     23  ElectronEnergySmearing
    2424  MuonMomentumSmearing
    2525
    … …  
    168168}
    169169
    170 ########################################
    171 # Momentum resolution for electrons
    172 ########################################
    173 
    174 module MomentumSmearing ElectronMomentumSmearing {
    175   ## hadrons after having applied the tracking efficiency
    176   set InputArray  ElectronTrackingEfficiency/electrons
     170#################################
     171# Energy resolution for electrons
     172#################################
     173
     174module EnergySmearing ElectronEnergySmearing {
     175  set InputArray ElectronTrackingEfficiency/electrons
    177176  set OutputArray electrons
     177
     178  # set ResolutionFormula {resolution formula as a function of eta and energy}
     179
    178180  # resolution formula for electrons
    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 
     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}
    184193
    185194###############################
    … …  
    206215# add InputArray InputArray
    207216  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    208   add InputArray ElectronMomentumSmearing/electrons
     217  add InputArray ElectronEnergySmearing/electrons
    209218  add InputArray MuonMomentumSmearing/muons
    210219  set OutputArray tracks
    … …  
    301310  # set ResolutionFormula {resolution formula as a function of eta and energy}
    302311 
    303   # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
     312  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
    304313  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    305314
    306   set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
    307                            (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    308                            (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    309                            (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
    310                            (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
    311 
     315  set ResolutionFormula {  (abs(eta) <= 1.50)                    * 1.00 * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.50^2) +
     316                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * 1.00 * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
     317                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * 1.00 * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
     318                           (abs(eta) > 2.15 && abs(eta) <= 3.00) * 1.00 * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
     319                           (abs(eta) >= 3.0 && abs(eta) <= 5.0)  * 1.00 * sqrt(energy^2*0.08^2 + energy*1.98^2)}
    312320}
    313321

  • cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl

    r4bf6061 rdf7986a  
    2424
    2525  ChargedHadronMomentumSmearing
    26   ElectronMomentumSmearing
     26  ElectronEnergySmearing
    2727  MuonMomentumSmearing
    2828
    … …  
    218218}
    219219
    220 ########################################
    221 # Momentum resolution for electrons
    222 ########################################
    223 
    224 module MomentumSmearing ElectronMomentumSmearing {
    225   ## hadrons after having applied the tracking efficiency
    226   set InputArray  ElectronTrackingEfficiency/electrons
     220#################################
     221# Energy resolution for electrons
     222#################################
     223
     224module EnergySmearing ElectronEnergySmearing {
     225  set InputArray ElectronTrackingEfficiency/electrons
    227226  set OutputArray electrons
     227
     228  # set ResolutionFormula {resolution formula as a function of eta and energy}
     229
    228230  # resolution formula for electrons
    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 
     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}
    234243
    235244###############################
    … …  
    258267# add InputArray InputArray
    259268  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    260   add InputArray ElectronMomentumSmearing/electrons
     269  add InputArray ElectronEnergySmearing/electrons
    261270  add InputArray MuonMomentumSmearing/muons
    262271  set OutputArray tracks
    … …  
    353362  # set ResolutionFormula {resolution formula as a function of eta and energy}
    354363 
    355   # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
     364  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
    356365  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    357366
    358   set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
    359                            (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    360                            (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    361                            (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
    362                            (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
     367  set ResolutionFormula {  (abs(eta) <= 1.50)                    * 1.00 * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.50^2) +
     368                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * 1.00 * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
     369                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * 1.00 * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
     370                           (abs(eta) > 2.15 && abs(eta) <= 3.00) * 1.00 * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
     371                           (abs(eta) >= 3.0 && abs(eta) <= 5.0)  * 1.00 * sqrt(energy^2*0.08^2 + energy*1.98^2)}
    363372
    364373}
Note: See TracChangeset for help on using the changeset viewer.