Fork me on GitHub

Changes in / [9327245:20ca0cd] in git


Ignore:
Files:
2 deleted
5 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl

    r9327245 r20ca0cd  
    3333  PhotonEnergySmearing
    3434  ElectronFilter
    35 
    3635  TrackPileUpSubtractor
    37   RecoPuFilter
    3836
    3937  TowerMerger
    4038  NeutralEFlowMerger
    41 
     39  EFlowMergerAllTracks
    4240  EFlowMerger
    43   EFlowMergerCHS
    44   Rho
    4541
    4642  LeptonFilterNoLep
    … …  
    5248
    5349  PhotonIsolation
    54   PhotonIsolationCHS
    5550  PhotonEfficiency
    56   PhotonEfficiencyCHS
    5751
    5852  ElectronIsolation
    59   ElectronIsolationCHS
    60 
    6153  ElectronEfficiency
    62   ElectronEfficiencyCHS
    6354
    6455  MuonIsolation
    65   MuonIsolationCHS
    66 
    6756  MuonLooseIdEfficiency
    6857  MuonTightIdEfficiency
    69 
    70   MuonLooseIdEfficiencyCHS
    71   MuonTightIdEfficiencyCHS
    7258
    7359  NeutrinoFilter
    … …  
    8268  FastJetFinder
    8369  FastJetFinderAK8
    84   JetPileUpSubtractor
    85   JetPileUpSubtractorAK8
    8670  FastJetFinderPUPPI
    8771  FastJetFinderPUPPIAK8
    … …  
    137121
    138122  # maximum spread in the beam direction in m
    139   set ZVertexSpread 0.25
     123  set ZVertexSpread 0
    140124
    141125  # maximum spread in time in s
    142   set TVertexSpread 800E-12
     126  set TVertexSpread 0
    143127
    144128  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    145   set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
     129  set VertexDistributionFormula 1
    146130
    147131}
    … …  
    289273  source muonMomentumResolution.tcl
    290274}
     275
     276
     277
    291278
    292279##############
    … …  
    540527}
    541528
    542 ########################
    543 # Reco PU filter
    544 ########################
    545 
    546 module RecoPuFilter RecoPuFilter {
    547   set InputArray HCal/eflowTracks
    548   set OutputArray eflowTracks
    549 }
    550529
    551530###################################################
    … …  
    560539}
    561540
     541
    562542####################
    563543# Neutral eflow erger
    … …  
    571551}
    572552
    573 #####################
     553
     554####################
    574555# Energy flow merger
    575 #####################
     556####################
    576557
    577558module Merger EFlowMerger {
    … …  
    583564}
    584565
    585 ############################
    586 # Energy flow merger no PU
    587 ############################
    588 
    589 module Merger EFlowMergerCHS {
     566##################################
     567# Energy flow merger (all tracks)
     568##################################
     569
     570module Merger EFlowMergerAllTracks {
    590571# add InputArray InputArray
    591   add InputArray RecoPuFilter/eflowTracks
     572  add InputArray TrackMerger/tracks
    592573  add InputArray PhotonEnergySmearing/eflowPhotons
    593574  add InputArray HCal/eflowNeutralHadrons
    … …  
    714695}
    715696
     697
    716698#####################
    717699# MC truth jet finder
    … …  
    753735}
    754736
    755 
    756 #############
    757 # Rho pile-up
    758 #############
    759 
    760 module FastJetFinder Rho {
    761 #  set InputArray Calorimeter/towers
    762   set InputArray EFlowMergerCHS/eflow
    763 
    764   set ComputeRho true
    765   set RhoOutputArray rho
    766 
    767   # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
    768   set AreaAlgorithm 5
    769 
    770   # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    771   set JetAlgorithm 4
    772   set ParameterR 0.4
    773   set GhostEtaMax 5.0
    774 
    775   add RhoEtaRange -5.0 -4.0
    776   add RhoEtaRange -4.0 -1.5
    777   add RhoEtaRange -1.5 1.5
    778   add RhoEtaRange 1.5 4.0
    779   add RhoEtaRange 4.0 5.0
    780 
    781   set JetPTMin 0.0
    782 }
    783 
    784 
    785 ##############
     737############
    786738# Jet finder
    787 ##############
     739############
    788740
    789741module FastJetFinder FastJetFinder {
    790742#  set InputArray TowerMerger/towers
    791   set InputArray EFlowMergerCHS/eflow
     743  set InputArray EFlowMerger/eflow
    792744
    793745  set OutputArray jets
    794 
    795   set AreaAlgorithm 5
    796746
    797747  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    805755module FastJetFinder FastJetFinderAK8 {
    806756#  set InputArray TowerMerger/towers
    807   set InputArray EFlowMergerCHS/eflow
     757  set InputArray EFlowMerger/eflow
    808758
    809759  set OutputArray jets
    810 
    811   set AreaAlgorithm 5
    812760
    813761  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    836784}
    837785
    838 ###########################
    839 # Jet Pile-Up Subtraction
    840 ###########################
    841 
    842 module JetPileUpSubtractor JetPileUpSubtractor {
    843   set JetInputArray FastJetFinder/jets
    844   set RhoInputArray Rho/rho
    845 
    846   set OutputArray jets
    847 
    848   set JetPTMin 15.0
    849 }
    850 
    851 ##############################
    852 # Jet Pile-Up Subtraction AK8
    853 ##############################
    854 
    855 module JetPileUpSubtractor JetPileUpSubtractorAK8 {
    856   set JetInputArray FastJetFinderAK8/jets
    857   set RhoInputArray Rho/rho
    858 
    859   set OutputArray jets
    860 
    861   set JetPTMin 15.0
    862 }
    863 
    864786module FastJetFinder FastJetFinderPUPPI {
    865787#  set InputArray TowerMerger/towers
    … …  
    911833
    912834module EnergyScale JetEnergyScale {
    913   set InputArray JetPileUpSubtractor/jets
     835  set InputArray FastJetFinder/jets
    914836  set OutputArray jets
    915837
    … …  
    919841
    920842module EnergyScale JetEnergyScaleAK8 {
    921   set InputArray JetPileUpSubtractorAK8/jets
     843  set InputArray FastJetFinderAK8/jets
    922844  set OutputArray jets
    923845
    … …  
    985907}
    986908
    987 
    988 ########################
    989 # Photon isolation CHS #
    990 ########################
    991 
    992 module Isolation PhotonIsolationCHS {
    993 
    994   # particle for which calculate the isolation
    995   set CandidateInputArray PhotonFilter/photons
    996 
    997   # isolation collection
    998   set IsolationInputArray EFlowMerger/eflow
    999 
    1000   # output array
    1001   set OutputArray photons
    1002 
    1003   # isolation cone
    1004   set DeltaRMax 0.3
    1005 
    1006   # minimum pT
    1007   set PTMin     1.0
    1008 
    1009   # iso ratio to cut
    1010   set PTRatioMax 9999.
    1011 
    1012 }
    1013909
    1014910
    … …  
    1033929
    1034930
    1035 #####################
    1036 # Photon efficiency #
    1037 #####################
    1038 
    1039 module Efficiency PhotonEfficiencyCHS {
    1040 
    1041   ## input particles
    1042   set InputArray PhotonIsolationCHS/photons
    1043   ## output particles
    1044   set OutputArray photons
    1045   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    1046   # efficiency formula for photons
    1047   set EfficiencyFormula {                      (pt <= 10.0) * (0.00) + \
    1048                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.9635) + \
    1049          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0)  * (0.9624) + \
    1050          (abs(eta) > 4.0)                                   * (0.00)}
    1051 
    1052 }
    1053 
    1054931######################
    1055932# Electron isolation #
    … …  
    1072949}
    1073950
    1074 
    1075 ##########################
    1076 # Electron isolation CHS #
    1077 ##########################
    1078 
    1079 module Isolation ElectronIsolationCHS {
    1080 
    1081   set CandidateInputArray ElectronFilter/electrons
    1082 
    1083   # isolation collection
    1084   set IsolationInputArray EFlowMerger/eflow
    1085 
    1086   set OutputArray electrons
    1087 
    1088   set DeltaRMax 0.3
    1089   set PTMin 1.0
    1090   set PTRatioMax 9999.
    1091 
    1092 }
    1093951
    1094952
    … …  
    1137995}
    1138996
    1139 ###########################
    1140 # Electron efficiency CHS #
    1141 ###########################
    1142 
    1143 module Efficiency ElectronEfficiencyCHS {
    1144 
    1145   set InputArray ElectronIsolationCHS/electrons
    1146   set OutputArray electrons
    1147 
    1148   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    1149   # efficiency formula for electrons
    1150   set EfficiencyFormula {
    1151                                       (pt <= 4.0)  * (0.00) + \
    1152                          (abs(eta) <= 1.45 ) * (pt >  4.0 && pt <= 6.0)   * (0.50) + \
    1153                          (abs(eta) <= 1.45 ) * (pt >  6.0 && pt <= 8.0)   * (0.70) + \
    1154                          (abs(eta) <= 1.45 ) * (pt >  8.0 && pt <= 10.0)  * (0.85) + \
    1155                          (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \
    1156                          (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \
    1157                          (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \
    1158                          (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \
    1159                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt >  4.0 && pt <= 10.0)   * (0.35) + \
    1160                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \
    1161                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \
    1162                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \
    1163                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt >  4.0 && pt <= 10.0)  * (0.75) + \
    1164                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \
    1165                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \
    1166                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \
    1167                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \
    1168                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt >  4.0 && pt <= 10.0)  * (0.65) + \
    1169                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \
    1170                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \
    1171                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \
    1172                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \
    1173                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
    1174                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
    1175                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
    1176                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
    1177                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
    1178                                           (abs(eta) > 4.0) * (0.00)
    1179 
    1180   }
    1181 }
    1182 
    1183 
    1184997##################
    1185998# Muon isolation #
    … …  
    12001013}
    12011014
    1202 ######################
    1203 # Muon isolation CHS #
    1204 ######################
    1205 
    1206 module Isolation MuonIsolationCHS {
    1207   set CandidateInputArray MuonMomentumSmearing/muons
    1208 
    1209   # isolation collection
    1210   set IsolationInputArray EFlowMerger/eflow
    1211 
    1212   set OutputArray muons
    1213 
    1214   set DeltaRMax 0.3
    1215   set PTMin 1.0
    1216   set PTRatioMax 9999.
    1217 
    1218 }
    1219 
    1220 
    1221 #####################
    1222 # Muon Loose Id     #
    1223 #####################
     1015
     1016##################
     1017# Muon Loose Id  #
     1018##################
    12241019
    12251020module Efficiency MuonLooseIdEfficiency {
    … …  
    12431038}
    12441039
    1245 
    1246 #####################
    1247 # Muon Loose Id CHS #
    1248 #####################
    1249 
    1250 module Efficiency MuonLooseIdEfficiencyCHS {
    1251     set InputArray MuonIsolationCHS/muons
    1252     set OutputArray muons
    1253     # tracking + LooseID efficiency formula for muons
    1254     source muonLooseId.tcl
    1255 
    1256 }
    1257 
    1258 
    1259 ######################
    1260 # Muon Tight Id  CHS #
    1261 ######################
    1262 
    1263 module Efficiency MuonTightIdEfficiencyCHS {
    1264     set InputArray MuonIsolationCHS/muons
    1265     set OutputArray muons
    1266     # tracking + TightID efficiency formula for muons
    1267     source muonTightId.tcl
    1268 }
    12691040
    12701041
    … …  
    14791250module StatusPidFilter GenParticleFilter {
    14801251
    1481     set InputArray Delphes/allParticles
     1252    set InputArray  Delphes/allParticles
    14821253    set OutputArray filteredParticles
    14831254    set PTMin 5.0
    … …  
    15071278  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    15081279  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    1509 
    1510   add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
    1511   add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
    1512   add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
    1513   add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
    15141280
    15151281  add Branch JetEnergyScale/jets Jet Jet
    … …  
    15231289  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    15241290  add Branch ScalarHT/energy ScalarHT ScalarHT
    1525 
    1526 }
     1291}

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl

    r9327245 r20ca0cd  
    3333  PhotonEnergySmearing
    3434  ElectronFilter
    35 
    3635  TrackPileUpSubtractor
    37   RecoPuFilter
    3836
    3937  TowerMerger
    4038  NeutralEFlowMerger
    41 
     39  EFlowMergerAllTracks
    4240  EFlowMerger
    43   EFlowMergerCHS
    44   Rho
    4541
    4642  LeptonFilterNoLep
    … …  
    5248
    5349  PhotonIsolation
    54   PhotonIsolationCHS
    5550  PhotonEfficiency
    56   PhotonEfficiencyCHS
    5751
    5852  ElectronIsolation
    59   ElectronIsolationCHS
    60 
    6153  ElectronEfficiency
    62   ElectronEfficiencyCHS
    6354
    6455  MuonIsolation
    65   MuonIsolationCHS
    66 
    6756  MuonLooseIdEfficiency
    6857  MuonTightIdEfficiency
    69 
    70   MuonLooseIdEfficiencyCHS
    71   MuonTightIdEfficiencyCHS
    7258
    7359  NeutrinoFilter
    … …  
    8268  FastJetFinder
    8369  FastJetFinderAK8
    84   JetPileUpSubtractor
    85   JetPileUpSubtractorAK8
    8670  FastJetFinderPUPPI
    8771  FastJetFinderPUPPIAK8
    … …  
    142126  set TVertexSpread 800E-12
    143127
    144   # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
     128  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
    145129  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    146130
    … …  
    154138module ParticlePropagator ParticlePropagator {
    155139  set InputArray PileUpMerger/stableParticles
    156   #set InputArray Delphes/stableParticles
    157140
    158141  set OutputArray stableParticles
    … …  
    289272  source muonMomentumResolution.tcl
    290273}
     274
     275
     276
    291277
    292278##############
    … …  
    540526}
    541527
    542 ########################
    543 # Reco PU filter
    544 ########################
    545 
    546 module RecoPuFilter RecoPuFilter {
    547   set InputArray HCal/eflowTracks
    548   set OutputArray eflowTracks
    549 }
    550528
    551529###################################################
    … …  
    560538}
    561539
     540
    562541####################
    563542# Neutral eflow erger
    … …  
    571550}
    572551
    573 #####################
     552
     553####################
    574554# Energy flow merger
    575 #####################
     555####################
    576556
    577557module Merger EFlowMerger {
    … …  
    583563}
    584564
    585 ############################
    586 # Energy flow merger no PU
    587 ############################
    588 
    589 module Merger EFlowMergerCHS {
     565##################################
     566# Energy flow merger (all tracks)
     567##################################
     568
     569module Merger EFlowMergerAllTracks {
    590570# add InputArray InputArray
    591   add InputArray RecoPuFilter/eflowTracks
     571  add InputArray TrackMerger/tracks
    592572  add InputArray PhotonEnergySmearing/eflowPhotons
    593573  add InputArray HCal/eflowNeutralHadrons
    … …  
    714694}
    715695
     696
    716697#####################
    717698# MC truth jet finder
    … …  
    754735
    755736
    756 #############
    757 # Rho pile-up
    758 #############
    759 
    760 module FastJetFinder Rho {
    761 #  set InputArray Calorimeter/towers
    762   set InputArray EFlowMergerCHS/eflow
    763 
    764   set ComputeRho true
    765   set RhoOutputArray rho
    766 
    767   # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
    768   set AreaAlgorithm 5
    769 
    770   # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    771   set JetAlgorithm 4
    772   set ParameterR 0.4
    773   set GhostEtaMax 5.0
    774 
    775   add RhoEtaRange -5.0 -4.0
    776   add RhoEtaRange -4.0 -1.5
    777   add RhoEtaRange -1.5 1.5
    778   add RhoEtaRange 1.5 4.0
    779   add RhoEtaRange 4.0 5.0
    780 
    781   set JetPTMin 0.0
    782 }
    783 
    784 
    785 ##############
     737
     738############
    786739# Jet finder
    787 ##############
     740############
    788741
    789742module FastJetFinder FastJetFinder {
    790743#  set InputArray TowerMerger/towers
    791   set InputArray EFlowMergerCHS/eflow
     744  set InputArray EFlowMerger/eflow
    792745
    793746  set OutputArray jets
    794 
    795   set AreaAlgorithm 5
    796747
    797748  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    805756module FastJetFinder FastJetFinderAK8 {
    806757#  set InputArray TowerMerger/towers
    807   set InputArray EFlowMergerCHS/eflow
     758  set InputArray EFlowMerger/eflow
    808759
    809760  set OutputArray jets
    810 
    811   set AreaAlgorithm 5
    812761
    813762  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    836785}
    837786
    838 ###########################
    839 # Jet Pile-Up Subtraction
    840 ###########################
    841 
    842 module JetPileUpSubtractor JetPileUpSubtractor {
    843   set JetInputArray FastJetFinder/jets
    844   set RhoInputArray Rho/rho
    845 
    846   set OutputArray jets
    847 
    848   set JetPTMin 15.0
    849 }
    850 
    851 ##############################
    852 # Jet Pile-Up Subtraction AK8
    853 ##############################
    854 
    855 module JetPileUpSubtractor JetPileUpSubtractorAK8 {
    856   set JetInputArray FastJetFinderAK8/jets
    857   set RhoInputArray Rho/rho
    858 
    859   set OutputArray jets
    860 
    861   set JetPTMin 15.0
    862 }
    863 
    864787module FastJetFinder FastJetFinderPUPPI {
    865788#  set InputArray TowerMerger/towers
    … …  
    911834
    912835module EnergyScale JetEnergyScale {
    913   set InputArray JetPileUpSubtractor/jets
     836  set InputArray FastJetFinder/jets
    914837  set OutputArray jets
    915838
    … …  
    919842
    920843module EnergyScale JetEnergyScaleAK8 {
    921   set InputArray JetPileUpSubtractorAK8/jets
     844  set InputArray FastJetFinderAK8/jets
    922845  set OutputArray jets
    923846
    … …  
    985908}
    986909
    987 
    988 ########################
    989 # Photon isolation CHS #
    990 ########################
    991 
    992 module Isolation PhotonIsolationCHS {
    993 
    994   # particle for which calculate the isolation
    995   set CandidateInputArray PhotonFilter/photons
    996 
    997   # isolation collection
    998   set IsolationInputArray EFlowMerger/eflow
    999 
    1000   # output array
    1001   set OutputArray photons
    1002 
    1003   # isolation cone
    1004   set DeltaRMax 0.3
    1005 
    1006   # minimum pT
    1007   set PTMin     1.0
    1008 
    1009   # iso ratio to cut
    1010   set PTRatioMax 9999.
    1011 
    1012 }
    1013910
    1014911
    … …  
    1033930
    1034931
    1035 #####################
    1036 # Photon efficiency #
    1037 #####################
    1038 
    1039 module Efficiency PhotonEfficiencyCHS {
    1040 
    1041   ## input particles
    1042   set InputArray PhotonIsolationCHS/photons
    1043   ## output particles
    1044   set OutputArray photons
    1045   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    1046   # efficiency formula for photons
    1047   set EfficiencyFormula {                      (pt <= 10.0) * (0.00) + \
    1048                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.9635) + \
    1049          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0)  * (0.9624) + \
    1050          (abs(eta) > 4.0)                                   * (0.00)}
    1051 
    1052 }
    1053 
    1054932######################
    1055933# Electron isolation #
    … …  
    1072950}
    1073951
    1074 
    1075 ##########################
    1076 # Electron isolation CHS #
    1077 ##########################
    1078 
    1079 module Isolation ElectronIsolationCHS {
    1080 
    1081   set CandidateInputArray ElectronFilter/electrons
    1082 
    1083   # isolation collection
    1084   set IsolationInputArray EFlowMerger/eflow
    1085 
    1086   set OutputArray electrons
    1087 
    1088   set DeltaRMax 0.3
    1089   set PTMin 1.0
    1090   set PTRatioMax 9999.
    1091 
    1092 }
    1093952
    1094953
    … …  
    1137996}
    1138997
    1139 ###########################
    1140 # Electron efficiency CHS #
    1141 ###########################
    1142 
    1143 module Efficiency ElectronEfficiencyCHS {
    1144 
    1145   set InputArray ElectronIsolationCHS/electrons
    1146   set OutputArray electrons
    1147 
    1148   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    1149   # efficiency formula for electrons
    1150   set EfficiencyFormula {
    1151                                       (pt <= 4.0)  * (0.00) + \
    1152                          (abs(eta) <= 1.45 ) * (pt >  4.0 && pt <= 6.0)   * (0.50) + \
    1153                          (abs(eta) <= 1.45 ) * (pt >  6.0 && pt <= 8.0)   * (0.70) + \
    1154                          (abs(eta) <= 1.45 ) * (pt >  8.0 && pt <= 10.0)  * (0.85) + \
    1155                          (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \
    1156                          (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \
    1157                          (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \
    1158                          (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \
    1159                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt >  4.0 && pt <= 10.0)   * (0.35) + \
    1160                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \
    1161                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \
    1162                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \
    1163                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt >  4.0 && pt <= 10.0)  * (0.75) + \
    1164                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \
    1165                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \
    1166                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \
    1167                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \
    1168                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt >  4.0 && pt <= 10.0)  * (0.65) + \
    1169                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \
    1170                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \
    1171                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \
    1172                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \
    1173                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
    1174                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
    1175                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
    1176                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
    1177                                           (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
    1178                                           (abs(eta) > 4.0) * (0.00)
    1179 
    1180   }
    1181 }
    1182 
    1183 
    1184998##################
    1185999# Muon isolation #
    … …  
    12001014}
    12011015
    1202 ######################
    1203 # Muon isolation CHS #
    1204 ######################
    1205 
    1206 module Isolation MuonIsolationCHS {
    1207   set CandidateInputArray MuonMomentumSmearing/muons
    1208 
    1209   # isolation collection
    1210   set IsolationInputArray EFlowMerger/eflow
    1211 
    1212   set OutputArray muons
    1213 
    1214   set DeltaRMax 0.3
    1215   set PTMin 1.0
    1216   set PTRatioMax 9999.
    1217 
    1218 }
    1219 
    1220 
    1221 #####################
    1222 # Muon Loose Id     #
    1223 #####################
     1016
     1017##################
     1018# Muon Loose Id  #
     1019##################
    12241020
    12251021module Efficiency MuonLooseIdEfficiency {
    … …  
    12431039}
    12441040
    1245 
    1246 #####################
    1247 # Muon Loose Id CHS #
    1248 #####################
    1249 
    1250 module Efficiency MuonLooseIdEfficiencyCHS {
    1251     set InputArray MuonIsolationCHS/muons
    1252     set OutputArray muons
    1253     # tracking + LooseID efficiency formula for muons
    1254     source muonLooseId.tcl
    1255 
    1256 }
    1257 
    1258 
    1259 ######################
    1260 # Muon Tight Id  CHS #
    1261 ######################
    1262 
    1263 module Efficiency MuonTightIdEfficiencyCHS {
    1264     set InputArray MuonIsolationCHS/muons
    1265     set OutputArray muons
    1266     # tracking + TightID efficiency formula for muons
    1267     source muonTightId.tcl
    1268 }
    12691041
    12701042
    … …  
    14791251module StatusPidFilter GenParticleFilter {
    14801252
    1481     set InputArray Delphes/allParticles
     1253    set InputArray  Delphes/allParticles
    14821254    set OutputArray filteredParticles
    14831255    set PTMin 5.0
    … …  
    15071279  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    15081280  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    1509 
    1510   add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
    1511   add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
    1512   add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
    1513   add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
    15141281
    15151282  add Branch JetEnergyScale/jets Jet Jet
    … …  
    15231290  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    15241291  add Branch ScalarHT/energy ScalarHT ScalarHT
    1525 
    1526 }
     1292}

  • modules/ModulesLinkDef.h

    r9327245 r20ca0cd  
    5858#include "modules/StatusPidFilter.h"
    5959#include "modules/PdgCodeFilter.h"
    60 #include "modules/RecoPuFilter.h"
    6160#include "modules/Cloner.h"
    6261#include "modules/Weighter.h"
    … …  
    104103#pragma link C++ class StatusPidFilter+;
    105104#pragma link C++ class PdgCodeFilter+;
    106 #pragma link C++ class RecoPuFilter+;
    107105#pragma link C++ class Cloner+;
    108106#pragma link C++ class Weighter+;

  • modules/PhotonConversions.cc

    r9327245 r20ca0cd  
    5959  fItInputArray(0), fConversionMap(0), fDecayXsec(0)
    6060{
    61   fDecayXsec = new TF1("decayXsec","1.0 - 4.0/3.0 * x * (1.0 - x)", 0.0, 1.0);
     61  fDecayXsec = new TF1;
    6262  fConversionMap = new DelphesCylindricalFormula;
    6363}
    … …  
    8484
    8585  fConversionMap->Compile(GetString("ConversionMap", "0.0"));
     86
     87#if  ROOT_VERSION_CODE < ROOT_VERSION(6,04,00)
     88  fDecayXsec->Compile("1.0 - 4.0/3.0 * x * (1.0 - x)");
     89#else
     90  fDecayXsec->GetFormula()->Compile("1.0 - 4.0/3.0 * x * (1.0 - x)");
     91#endif
     92  fDecayXsec->SetRange(0.0, 1.0);
    8693
    8794  // import array with output from filter/classifier module

  • readers/DelphesCMSFWLite.cpp

    r9327245 r20ca0cd  
    6565  ExRootTreeBranch *branchEvent, ExRootTreeBranch *branchRwgt,
    6666  DelphesFactory *factory, TObjArray *allParticleOutputArray,
    67   TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray, Bool_t firstEvent)
     67  TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray)
    6868{
    69 
    7069  fwlite::Handle< GenEventInfoProduct > handleGenEventInfo;
     70
    7171  fwlite::Handle< LHEEventProduct > handleLHEEvent;
     72
    7273  fwlite::Handle< vector< reco::GenParticle > > handleParticle;
    73 
    7474  vector< reco::GenParticle >::const_iterator itParticle;
    7575
    … …  
    7878
    7979  handleGenEventInfo.getByLabel(event, "generator");
    80 
    81   if (!((handleLHEEvent.getBranchNameFor(event, "source")).empty()))
    82   {
    83     handleLHEEvent.getByLabel(event, "source");
    84   }
    85   else if (!((handleLHEEvent.getBranchNameFor(event, "externalLHEProducer")).empty()))
    86   {
    87     handleLHEEvent.getByLabel(event, "externalLHEProducer");
    88   }
    89   else if (firstEvent)
    90   {
    91     std::cout<<"Wrong LHEEvent Label! Please, check the input file."<<std::endl;
    92   }
    93 
    94   if (!((handleParticle.getBranchNameFor(event, "genParticles")).empty()))
    95   {
    96     handleParticle.getByLabel(event, "genParticles");
    97   }
    98   else if (!((handleParticle.getBranchNameFor(event, "prunedGenParticles")).empty()))
    99   {
    100     handleParticle.getByLabel(event, "prunedGenParticles");
    101   }
    102   else
    103   {
    104     std::cout<<"Wrong GenParticle Label! Please, check the input file."<<std::endl;
    105     exit(-1);
    106   }
    107 
    108   Bool_t foundLHE = !((handleLHEEvent.getBranchNameFor(event, "source")).empty()) || !((handleLHEEvent.getBranchNameFor(event, "externalLHEProducer")).empty());
     80  handleLHEEvent.getByLabel(event, "externalLHEProducer");
     81  handleParticle.getByLabel(event, "genParticles");
    10982
    11083  HepMCEvent *element;
    … …  
    11891  Double_t px, py, pz, e, mass;
    11992  Double_t x, y, z;
     93
     94  const vector< gen::WeightsInfo > &vectorWeightsInfo = handleLHEEvent->weights();
     95  vector< gen::WeightsInfo >::const_iterator itWeightsInfo;
    12096
    12197  element = static_cast<HepMCEvent *>(branchEvent->NewEntry());
    … …  
    141117  element->ProcTime = 0.0;
    142118
    143 
    144   if(foundLHE)
    145   {
    146     const vector< gen::WeightsInfo > &vectorWeightsInfo = handleLHEEvent->weights();
    147     vector< gen::WeightsInfo >::const_iterator itWeightsInfo;
    148    
    149     for(itWeightsInfo = vectorWeightsInfo.begin(); itWeightsInfo != vectorWeightsInfo.end(); ++itWeightsInfo)
    150     {
    151       weight = static_cast<Weight *>(branchRwgt->NewEntry());
    152       weight->Weight = itWeightsInfo->wgt;
    153     } 
     119  for(itWeightsInfo = vectorWeightsInfo.begin(); itWeightsInfo != vectorWeightsInfo.end(); ++itWeightsInfo)
     120  {
     121    weight = static_cast<Weight *>(branchRwgt->NewEntry());
     122    weight->Weight = itWeightsInfo->wgt;
    154123  }
    155124
    … …  
    238207  Int_t i;
    239208  Long64_t eventCounter, numberOfEvents;
    240   Bool_t firstEvent = kTRUE;
    241209
    242210  if(argc < 4)
    … …  
    313281      modularDelphes->Clear();
    314282      treeWriter->Clear();
    315 
    316283      for(event.toBegin(); !event.atEnd() && !interrupted; ++event)
    317284      {
    318285        ConvertInput(event, eventCounter, branchEvent, branchRwgt, factory,
    319           allParticleOutputArray, stableParticleOutputArray, partonOutputArray, firstEvent);
     286          allParticleOutputArray, stableParticleOutputArray, partonOutputArray);
    320287        modularDelphes->ProcessTask();
    321          
    322         firstEvent = kFALSE;
    323288
    324289        treeWriter->Fill();
Note: See TracChangeset for help on using the changeset viewer.