Fork me on GitHub

Changeset 29b722a in git


Ignore:
Timestamp:
Dec 9, 2021, 8:42:12 PM (3 years ago)
Author:
christinaw97 <christina.wang@…>
Children:
a09b75f
Parents:
d612dec
Message:

Added DelphesCscClusterFormula

Files:
2 added
8 edited

Legend:

Unmodified
Added
Removed
  • Makefile

    rd612dec r29b722a  
    522522        classes/DelphesFactory.h \
    523523        classes/SortableObject.h
     524tmp/classes/DelphesCscClusterFormula.$(ObjSuf): \
     525                classes/DelphesCscClusterFormula.$(SrcSuf) \
     526                classes/DelphesCscClusterFormula.h
    524527tmp/classes/DelphesCylindricalFormula.$(ObjSuf): \
    525528        classes/DelphesCylindricalFormula.$(SrcSuf) \
     
    751754        classes/DelphesClasses.h \
    752755        classes/DelphesFactory.h \
    753         classes/DelphesFormula.h \
     756        classes/DelphesCscClusterFormula.h \
    754757        external/ExRootAnalysis/ExRootClassifier.h \
    755758        external/ExRootAnalysis/ExRootFilter.h \
     
    759762        modules/CscClusterId.h \
    760763        classes/DelphesClasses.h \
    761         classes/DelphesFactory.h \
    762         classes/DelphesFormula.h \
     764  classes/DelphesFactory.h \
     765        classes/DelphesCscClusterFormula.h \
    763766        external/ExRootAnalysis/ExRootClassifier.h \
    764767        external/ExRootAnalysis/ExRootFilter.h \
     
    11941197DELPHES_OBJ +=  \
    11951198        tmp/classes/DelphesClasses.$(ObjSuf) \
     1199        tmp/classes/DelphesCscClusterFormula.$(ObjSuf) \
    11961200        tmp/classes/DelphesCylindricalFormula.$(ObjSuf) \
    11971201        tmp/classes/DelphesFactory.$(ObjSuf) \
     
    18011805ifeq ($(HAS_PYTHIA8),true)
    18021806FASTJET_OBJ +=  \
    1803        
     1807
    18041808endif
    18051809
     
    18481852ifeq ($(HAS_PYTHIA8),true)
    18491853DISPLAY_OBJ +=  \
    1850        
     1854
    18511855endif
    18521856
     
    25942598
    25952599###
    2596 
    2597 
  • cards/delphes_card_CMS_CSCCluster.tcl

    rd612dec r29b722a  
    6565  TreeWriter
    6666}
     67
    6768
    6869#################################
     
    408409#################
    409410
    410 module PdgCodeFilter HiggsFilter {
    411 set InputArray Delphes/allParticles
    412   set OutputArray higgs
    413   set Invert true
    414   add PdgCode {25}
    415 }
    416 
    417 
    418 
    419411
    420412# filter out LLPs that decay in CSC
     
    426418  # DecayRegion = 1: select LLP that decays in CSC volume
    427419  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
    428   set RequireStatus true
    429   set Status 2
    430   add PdgCode {9000005}
     420  set RequireStatus false
     421  add PdgCode {1500001}
     422
    431423
    432424}
     
    436428  set OutputArray LLP
    437429  set DecayRegion 0
    438   set RequireStatus true
    439   set Status 2
     430  set RequireStatus false
    440431
    441432  # DecayRegion = 0: no cuts on decay region
    442433  # DecayRegion = 1: select LLP that decays in CSC volume
    443434  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
    444   add PdgCode {9000005}
     435  add PdgCode {1500001}
    445436
    446437}
     
    455446  set EfficiencyFormula {
    456447
    457     (pt > 3900 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0049) +
    458     (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0130) +
    459     (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0346) +
    460     (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0623) +
    461     (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0919) +
    462     (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.1086) +
    463     (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.1292) +
    464     (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.1106) +
    465     (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0249) +
    466     (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0285) +
    467     (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0501) +
    468     (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0841) +
    469     (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.1021) +
    470     (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.1129) +
    471     (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.1141) +
    472     (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.1370) +
    473     (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0282) +
    474     (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0445) +
    475     (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0643) +
    476     (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0903) +
    477     (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0998) +
    478     (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.1420) +
    479     (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.1429) +
    480     (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0882) +
    481     (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0594) +
    482     (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0521) +
    483     (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0605) +
    484     (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0791) +
    485     (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.1117) +
    486     (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0862) +
    487     (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0698) +
    488     (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0500) +
    489     (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0758) +
    490     (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0414) +
    491     (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0755) +
    492     (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.1027) +
    493     (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0440) +
    494     (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0811) +
    495     (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.1538) +
    496     (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0833) +
    497     (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0300) +
    498     (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0609) +
    499     (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0745) +
    500     (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0610) +
    501     (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.1224) +
    502     (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.1667) +
    503     (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0000) +
    504     (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0000) +
    505     (energy >= 150.0 && 0.0 == phi)*(0.0282) +
    506     (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0809) +
    507     (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0352) +
    508     (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0984) +
    509     (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0968) +
    510     (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.1282) +
    511     (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.2105) +
    512     (energy >= 150.0 && phi >= 200.0)*(0.0769)) +
    513     (eta > 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0184) +
    514     (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0772) +
    515     (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.2086) +
    516     (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.3091) +
    517     (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.3867) +
    518     (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.4500) +
    519     (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.4746) +
    520     (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.4906) +
    521     (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0955) +
    522     (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.1461) +
    523     (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.2594) +
    524     (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.3556) +
    525     (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.4165) +
    526     (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.4693) +
    527     (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.5054) +
    528     (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.5219) +
    529     (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.1472) +
    530     (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.1970) +
    531     (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.2974) +
    532     (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.3783) +
    533     (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.4335) +
    534     (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.4736) +
    535     (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.4937) +
    536     (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.5077) +
    537     (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.2053) +
    538     (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.2314) +
    539     (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.3114) +
    540     (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.3799) +
    541     (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.4420) +
    542     (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.4502) +
    543     (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.5348) +
    544     (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.5115) +
    545     (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.2198) +
    546     (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.2404) +
    547     (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.3295) +
    548     (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.3932) +
    549     (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.4327) +
    550     (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.4377) +
    551     (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.5175) +
    552     (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.6087) +
    553     (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.2147) +
    554     (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.2605) +
    555     (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.3442) +
    556     (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.3622) +
    557     (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.4407) +
    558     (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.5168) +
    559     (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.5056) +
    560     (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.4559) +
    561     (energy >= 150.0 && 0.0 == phi)*(0.2824) +
    562     (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.2447) +
    563     (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.3519) +
    564     (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.3772) +
    565     (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.4447) +
    566     (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.4703) +
    567     (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.4460) +
    568     (energy >= 150.0 && phi >= 200.0)*(0.4400)) +
    569     (pt < 2700 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0002) +
    570     (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
    571     (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0006) +
    572     (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0014) +
    573     (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0025) +
    574     (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.0046) +
    575     (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.0060) +
    576     (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.0136) +
    577     (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0000) +
    578     (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
    579     (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0006) +
    580     (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0015) +
    581     (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.0033) +
    582     (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.0051) +
    583     (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.0098) +
    584     (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.0146) +
    585     (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0000) +
    586     (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
    587     (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0003) +
    588     (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0015) +
    589     (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0038) +
    590     (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.0052) +
    591     (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.0114) +
    592     (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0181) +
    593     (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0000) +
    594     (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
    595     (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0005) +
    596     (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0022) +
    597     (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.0067) +
    598     (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0047) +
    599     (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0113) +
    600     (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0145) +
    601     (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0000) +
    602     (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
    603     (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0003) +
    604     (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.0016) +
    605     (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0110) +
    606     (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0029) +
    607     (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.0138) +
    608     (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0000) +
    609     (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0000) +
    610     (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
    611     (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) +
    612     (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0026) +
    613     (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.0047) +
    614     (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.0085) +
    615     (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0152) +
    616     (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0164) +
    617     (energy >= 150.0 && 0.0 == phi)*(0.0000) +
    618     (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
    619     (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) +
    620     (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0000) +
    621     (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0000) +
    622     (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.0000) +
    623     (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.0080) +
    624     (energy >= 150.0 && phi >= 200.0)*(0.0143))
    625  }
     448      (decayR > 3900 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0049) +
     449      (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0130) +
     450      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0346) +
     451      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0623) +
     452      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0919) +
     453      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1086) +
     454      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1292) +
     455      (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.1106) +
     456      (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0249) +
     457      (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0285) +
     458      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0501) +
     459      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0841) +
     460      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1021) +
     461      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1129) +
     462      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1141) +
     463      (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.1370) +
     464      (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0282) +
     465      (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0445) +
     466      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0643) +
     467      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0903) +
     468      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0998) +
     469      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1420) +
     470      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1429) +
     471      (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0882) +
     472      (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0594) +
     473      (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0521) +
     474      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0605) +
     475      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0791) +
     476      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1117) +
     477      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0862) +
     478      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0698) +
     479      (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0500) +
     480      (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0758) +
     481      (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0414) +
     482      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0755) +
     483      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.1027) +
     484      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0440) +
     485      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0811) +
     486      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1538) +
     487      (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0833) +
     488      (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0300) +
     489      (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0609) +
     490      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0745) +
     491      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0610) +
     492      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1224) +
     493      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1667) +
     494      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0000) +
     495      (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0000) +
     496      (Eem >= 150.0 && 0.0 == Ehad)*(0.0282) +
     497      (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0809) +
     498      (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0352) +
     499      (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0984) +
     500      (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0968) +
     501      (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1282) +
     502      (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.2105) +
     503      (Eem >= 150.0 && Ehad >= 200.0)*(0.0769)) +
     504      (decayZ > 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0184) +
     505      (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0772) +
     506      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2086) +
     507      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3091) +
     508      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.3867) +
     509      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4500) +
     510      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4746) +
     511      (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.4906) +
     512      (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0955) +
     513      (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1461) +
     514      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2594) +
     515      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3556) +
     516      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4165) +
     517      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4693) +
     518      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5054) +
     519      (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.5219) +
     520      (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.1472) +
     521      (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1970) +
     522      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2974) +
     523      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3783) +
     524      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4335) +
     525      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4736) +
     526      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4937) +
     527      (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.5077) +
     528      (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.2053) +
     529      (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2314) +
     530      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3114) +
     531      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3799) +
     532      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4420) +
     533      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4502) +
     534      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5348) +
     535      (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.5115) +
     536      (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.2198) +
     537      (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2404) +
     538      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3295) +
     539      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3932) +
     540      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4327) +
     541      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4377) +
     542      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5175) +
     543      (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.6087) +
     544      (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.2147) +
     545      (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2605) +
     546      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3442) +
     547      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3622) +
     548      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4407) +
     549      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.5168) +
     550      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5056) +
     551      (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.4559) +
     552      (Eem >= 150.0 && 0.0 == Ehad)*(0.2824) +
     553      (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2447) +
     554      (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3519) +
     555      (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3772) +
     556      (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4447) +
     557      (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4703) +
     558      (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4460) +
     559      (Eem >= 150.0 && Ehad >= 200.0)*(0.4400)) +
     560      (decayR < 2700 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0002) +
     561      (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
     562      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) +
     563      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0014) +
     564      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0025) +
     565      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0046) +
     566      (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0060) +
     567      (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.0136) +
     568      (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0000) +
     569      (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
     570      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) +
     571      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) +
     572      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0033) +
     573      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0051) +
     574      (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0098) +
     575      (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.0146) +
     576      (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0000) +
     577      (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
     578      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) +
     579      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) +
     580      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0038) +
     581      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0052) +
     582      (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0114) +
     583      (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0181) +
     584      (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0000) +
     585      (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
     586      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0005) +
     587      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0022) +
     588      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0067) +
     589      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0047) +
     590      (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0113) +
     591      (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0145) +
     592      (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0000) +
     593      (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
     594      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) +
     595      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0016) +
     596      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0110) +
     597      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0029) +
     598      (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0138) +
     599      (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0000) +
     600      (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0000) +
     601      (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
     602      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) +
     603      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0026) +
     604      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0047) +
     605      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0085) +
     606      (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0152) +
     607      (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0164) +
     608      (Eem >= 150.0 && 0.0 == Ehad)*(0.0000) +
     609      (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
     610      (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) +
     611      (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0000) +
     612      (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0000) +
     613      (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0000) +
     614      (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0080) +
     615      (Eem >= 150.0 && Ehad >= 200.0)*(0.0143)) }
    626616}
    627617
     
    632622  # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
    633623  set EfficiencyFormula {
    634     (pt > 3900 && eta < 6710) * ((0.0 == phi)*(0.0656) +
    635     (phi > 0.0 && phi < 25.0)*(0.0777) +
    636     (phi >= 25.0 && phi < 50.0)*(0.1607) +
    637     (phi >= 50.0 && phi < 75.0)*(0.2294) +
    638     (phi >= 75.0 && phi < 100.0)*(0.3146) +
    639     (phi >= 100.0 && phi < 150.0)*(0.3107) +
    640     (phi >= 150.0 && phi < 200.0)*(0.3177) +
    641     (phi >= 200.0)*(0.3229)) +
    642     (eta > 6710) * ((0.0 == phi)*(0.2987) +
    643     (phi > 0.0 && phi < 25.0)*(0.3100) +
    644     (phi >= 25.0 && phi < 50.0)*(0.4476) +
    645     (phi >= 50.0 && phi < 75.0)*(0.5335) +
    646     (phi >= 75.0 && phi < 100.0)*(0.5961) +
    647     (phi >= 100.0 && phi < 150.0)*(0.6368) +
    648     (phi >= 150.0 && phi < 200.0)*(0.6814) +
    649     (phi >= 200.0)*(0.6998)) +
    650     (pt < 2700 && eta < 6710) * ((0.0 == phi)*(0.8604) +
    651     (phi > 0.0 && phi < 25.0)*(0.3335) +
    652     (phi >= 25.0 && phi < 50.0)*(0.2457) +
    653     (phi >= 50.0 && phi < 75.0)*(0.1831) +
    654     (phi >= 75.0 && phi < 100.0)*(0.2100) +
    655     (phi >= 100.0 && phi < 150.0)*(0.2443) +
    656     (phi >= 150.0 && phi < 200.0)*(0.2532) +
    657     (phi >= 200.0)*(0.2404))
     624    (decayR > 3900 && decayZ < 6710) * ((0.0 == Ehad)*(0.0656) +
     625    (Ehad > 0.0 && Ehad < 25.0)*(0.0777) +
     626    (Ehad >= 25.0 && Ehad < 50.0)*(0.1607) +
     627    (Ehad >= 50.0 && Ehad < 75.0)*(0.2294) +
     628    (Ehad >= 75.0 && Ehad < 100.0)*(0.3146) +
     629    (Ehad >= 100.0 && Ehad < 150.0)*(0.3107) +
     630    (Ehad >= 150.0 && Ehad < 200.0)*(0.3177) +
     631    (Ehad >= 200.0)*(0.3229)) +
     632    (decayZ > 6710) * ((0.0 == Ehad)*(0.2987) +
     633    (Ehad > 0.0 && Ehad < 25.0)*(0.3100) +
     634    (Ehad >= 25.0 && Ehad < 50.0)*(0.4476) +
     635    (Ehad >= 50.0 && Ehad < 75.0)*(0.5335) +
     636    (Ehad >= 75.0 && Ehad < 100.0)*(0.5961) +
     637    (Ehad >= 100.0 && Ehad < 150.0)*(0.6368) +
     638    (Ehad >= 150.0 && Ehad < 200.0)*(0.6814) +
     639    (Ehad >= 200.0)*(0.6998)) +
     640    (decayR < 2700 && decayZ < 6710) * ((0.0 == Ehad)*(0.8604) +
     641    (Ehad > 0.0 && Ehad < 25.0)*(0.3335) +
     642    (Ehad >= 25.0 && Ehad < 50.0)*(0.2457) +
     643    (Ehad >= 50.0 && Ehad < 75.0)*(0.1831) +
     644    (Ehad >= 75.0 && Ehad < 100.0)*(0.2100) +
     645    (Ehad >= 100.0 && Ehad < 150.0)*(0.2443) +
     646    (Ehad >= 150.0 && Ehad < 200.0)*(0.2532) +
     647    (Ehad >= 200.0)*(0.2404))
    658648  }
    659649}
     
    10551045  #add Branch GenJetFinder/jets GenJet Jet
    10561046  #add Branch GenMissingET/momentum GenMissingET MissingET
    1057   add Branch HiggsFilter/higgs higgs GenParticle
    10581047
    10591048
     
    10711060  add Branch CSCFilter/LLP Cscllp CscCluster
    10721061
    1073   add Branch ClusterEfficiency/cluster CscCluster130 CscCluster
    1074 }
     1062  add Branch ClusterEfficiency/cluster CscCluster CscCluster
     1063}
  • classes/DelphesHepMC2Reader.cc

    rd612dec r29b722a  
    438438{
    439439  Candidate *candidate;
     440  Candidate *candidateDaughter;
    440441  map<int, pair<int, int> >::iterator itMotherMap;
    441442  map<int, pair<int, int> >::iterator itDaughterMap;
     
    446447    candidate = static_cast<Candidate *>(allParticleOutputArray->At(i));
    447448
     449
    448450    if(candidate->M1 > 0)
    449451    {
     
    477479        candidate->D1 = -1;
    478480        candidate->D2 = -1;
     481        const TLorentzVector &decayPosition = candidate->Position;
     482        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
    479483      }
    480484      else
     
    482486        candidate->D1 = itDaughterMap->second.first;
    483487        candidate->D2 = itDaughterMap->second.second;
    484       }
    485     }
    486   }
    487 }
    488 
    489 //---------------------------------------------------------------------------
     488        candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1));
     489        const TLorentzVector &decayPosition = candidateDaughter->Position;
     490        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
     491
     492      }
     493    }
     494  }
     495}
     496
     497//---------------------------------------------------------------------------
  • classes/DelphesHepMC3Reader.cc

    rd612dec r29b722a  
    458458  TObjArray *array;
    459459  Candidate *candidate;
     460  Candidate *candidateDaughter;
    460461  TParticlePDG *pdgParticle;
    461462  int pdgCode;
     
    573574        candidate->D1 = -1;
    574575        candidate->D2 = -1;
     576        const TLorentzVector &decayPosition = candidate->Position;
     577        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
    575578      }
    576579      else
     
    578581        candidate->D1 = itDaughterMap->second.first;
    579582        candidate->D2 = itDaughterMap->second.second;
    580       }
    581     }
    582   }
    583 }
    584 
    585 //---------------------------------------------------------------------------
     583        candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1));
     584        const TLorentzVector &decayPosition = candidateDaughter->Position;
     585        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
     586      }
     587    }
     588  }
     589}
     590
     591//---------------------------------------------------------------------------
  • modules/CscClusterEfficiency.cc

    rd612dec r29b722a  
    2929#include "classes/DelphesClasses.h"
    3030#include "classes/DelphesFactory.h"
    31 #include "classes/DelphesFormula.h"
     31#include "classes/DelphesCscClusterFormula.h"
    3232
    3333#include "ExRootAnalysis/ExRootClassifier.h"
     
    5555  fFormula(0), fItInputArray(0)
    5656{
    57   fFormula = new DelphesFormula;
     57  fFormula = new DelphesCscClusterFormula;
    5858}
    5959
     
    7070{
    7171  // read CscClusterEfficiency formula
    72 
     72  cout<<"initialize"<<endl;
    7373  fFormula->Compile(GetString("EfficiencyFormula", "1.0"));
    7474
     
    106106    Eem = candidate->Eem;
    107107    // apply an efficency formula
    108     //pt, eta, phi, energy
    109108    if(gRandom->Uniform() > fFormula->Eval(decayR, decayZ, Ehad, Eem)) continue;
     109
    110110
    111111    fOutputArray->Add(candidate);
  • modules/CscClusterEfficiency.h

    rd612dec r29b722a  
    3232class TIterator;
    3333class TObjArray;
    34 class DelphesFormula;
     34class DelphesCscClusterFormula;
    3535
    3636class CscClusterEfficiency: public DelphesModule
     
    4545
    4646private:
    47   DelphesFormula *fFormula; //!
     47  DelphesCscClusterFormula *fFormula; //!
    4848
    4949  TIterator *fItInputArray; //!
  • modules/CscClusterId.cc

    rd612dec r29b722a  
    2929#include "classes/DelphesClasses.h"
    3030#include "classes/DelphesFactory.h"
    31 #include "classes/DelphesFormula.h"
     31#include "classes/DelphesCscClusterFormula.h"
    3232
    3333#include "ExRootAnalysis/ExRootClassifier.h"
     
    5555  fFormula(0), fItInputArray(0)
    5656{
    57   fFormula = new DelphesFormula;
     57  fFormula = new DelphesCscClusterFormula;
    5858}
    5959
     
    113113    eta = (cosTheta == 1.0 ? signPz * 999.9 : momentum.Eta());
    114114
    115     NStationEff = fFormula->Eval(decayR, decayZ, Ehad); //pt is used as argument in DelphesFormula
     115    NStationEff = fFormula->Eval(decayR, decayZ, Ehad); //pt is used as argument in DelphesCscClusterFormula
    116116
    117117    // assign average station for the cluster
  • modules/CscClusterId.h

    rd612dec r29b722a  
    3232class TIterator;
    3333class TObjArray;
    34 class DelphesFormula;
     34class DelphesCscClusterFormula;
    3535
    3636class CscClusterId: public DelphesModule
     
    4545
    4646private:
    47   DelphesFormula *fFormula; //!
     47  DelphesCscClusterFormula *fFormula; //!
    4848
    4949  TIterator *fItInputArray; //!
Note: See TracChangeset for help on using the changeset viewer.