Fork me on GitHub

Changes in / [20ca0cd:9327245] in git


Ignore:
Files:
2 added
5 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl

    r20ca0cd r9327245  
    3333  PhotonEnergySmearing
    3434  ElectronFilter
     35
    3536  TrackPileUpSubtractor
     37  RecoPuFilter
    3638
    3739  TowerMerger
    3840  NeutralEFlowMerger
    39   EFlowMergerAllTracks
     41
    4042  EFlowMerger
     43  EFlowMergerCHS
     44  Rho
    4145
    4246  LeptonFilterNoLep
    … …  
    4852
    4953  PhotonIsolation
     54  PhotonIsolationCHS
    5055  PhotonEfficiency
     56  PhotonEfficiencyCHS
    5157
    5258  ElectronIsolation
     59  ElectronIsolationCHS
     60
    5361  ElectronEfficiency
     62  ElectronEfficiencyCHS
    5463
    5564  MuonIsolation
     65  MuonIsolationCHS
     66
    5667  MuonLooseIdEfficiency
    5768  MuonTightIdEfficiency
     69
     70  MuonLooseIdEfficiencyCHS
     71  MuonTightIdEfficiencyCHS
    5872
    5973  NeutrinoFilter
    … …  
    6882  FastJetFinder
    6983  FastJetFinderAK8
     84  JetPileUpSubtractor
     85  JetPileUpSubtractorAK8
    7086  FastJetFinderPUPPI
    7187  FastJetFinderPUPPIAK8
    … …  
    121137
    122138  # maximum spread in the beam direction in m
    123   set ZVertexSpread 0
     139  set ZVertexSpread 0.25
    124140
    125141  # maximum spread in time in s
    126   set TVertexSpread 0
     142  set TVertexSpread 800E-12
    127143
    128144  # 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))}
    129   set VertexDistributionFormula 1
     145  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    130146
    131147}
    … …  
    273289  source muonMomentumResolution.tcl
    274290}
    275 
    276 
    277 
    278291
    279292##############
    … …  
    527540}
    528541
     542########################
     543# Reco PU filter
     544########################
     545
     546module RecoPuFilter RecoPuFilter {
     547  set InputArray HCal/eflowTracks
     548  set OutputArray eflowTracks
     549}
    529550
    530551###################################################
    … …  
    539560}
    540561
    541 
    542562####################
    543563# Neutral eflow erger
    … …  
    551571}
    552572
    553 
    554 ####################
     573#####################
    555574# Energy flow merger
    556 ####################
     575#####################
    557576
    558577module Merger EFlowMerger {
    … …  
    564583}
    565584
    566 ##################################
    567 # Energy flow merger (all tracks)
    568 ##################################
    569 
    570 module Merger EFlowMergerAllTracks {
     585############################
     586# Energy flow merger no PU
     587############################
     588
     589module Merger EFlowMergerCHS {
    571590# add InputArray InputArray
    572   add InputArray TrackMerger/tracks
     591  add InputArray RecoPuFilter/eflowTracks
    573592  add InputArray PhotonEnergySmearing/eflowPhotons
    574593  add InputArray HCal/eflowNeutralHadrons
    … …  
    695714}
    696715
    697 
    698716#####################
    699717# MC truth jet finder
    … …  
    735753}
    736754
    737 ############
     755
     756#############
     757# Rho pile-up
     758#############
     759
     760module 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##############
    738786# Jet finder
    739 ############
     787##############
    740788
    741789module FastJetFinder FastJetFinder {
    742790#  set InputArray TowerMerger/towers
    743   set InputArray EFlowMerger/eflow
     791  set InputArray EFlowMergerCHS/eflow
    744792
    745793  set OutputArray jets
     794
     795  set AreaAlgorithm 5
    746796
    747797  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    755805module FastJetFinder FastJetFinderAK8 {
    756806#  set InputArray TowerMerger/towers
    757   set InputArray EFlowMerger/eflow
     807  set InputArray EFlowMergerCHS/eflow
    758808
    759809  set OutputArray jets
     810
     811  set AreaAlgorithm 5
    760812
    761813  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    784836}
    785837
     838###########################
     839# Jet Pile-Up Subtraction
     840###########################
     841
     842module 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
     855module 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
    786864module FastJetFinder FastJetFinderPUPPI {
    787865#  set InputArray TowerMerger/towers
    … …  
    833911
    834912module EnergyScale JetEnergyScale {
    835   set InputArray FastJetFinder/jets
     913  set InputArray JetPileUpSubtractor/jets
    836914  set OutputArray jets
    837915
    … …  
    841919
    842920module EnergyScale JetEnergyScaleAK8 {
    843   set InputArray FastJetFinderAK8/jets
     921  set InputArray JetPileUpSubtractorAK8/jets
    844922  set OutputArray jets
    845923
    … …  
    907985}
    908986
     987
     988########################
     989# Photon isolation CHS #
     990########################
     991
     992module 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}
    9091013
    9101014
    … …  
    9291033
    9301034
     1035#####################
     1036# Photon efficiency #
     1037#####################
     1038
     1039module 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
    9311054######################
    9321055# Electron isolation #
    … …  
    9491072}
    9501073
     1074
     1075##########################
     1076# Electron isolation CHS #
     1077##########################
     1078
     1079module 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}
    9511093
    9521094
    … …  
    9951137}
    9961138
     1139###########################
     1140# Electron efficiency CHS #
     1141###########################
     1142
     1143module 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
    9971184##################
    9981185# Muon isolation #
    … …  
    10131200}
    10141201
    1015 
    1016 ##################
    1017 # Muon Loose Id  #
    1018 ##################
     1202######################
     1203# Muon isolation CHS #
     1204######################
     1205
     1206module 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#####################
    10191224
    10201225module Efficiency MuonLooseIdEfficiency {
    … …  
    10381243}
    10391244
     1245
     1246#####################
     1247# Muon Loose Id CHS #
     1248#####################
     1249
     1250module 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
     1263module Efficiency MuonTightIdEfficiencyCHS {
     1264    set InputArray MuonIsolationCHS/muons
     1265    set OutputArray muons
     1266    # tracking + TightID efficiency formula for muons
     1267    source muonTightId.tcl
     1268}
    10401269
    10411270
    … …  
    12501479module StatusPidFilter GenParticleFilter {
    12511480
    1252     set InputArray  Delphes/allParticles
     1481    set InputArray Delphes/allParticles
    12531482    set OutputArray filteredParticles
    12541483    set PTMin 5.0
    … …  
    12781507  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    12791508  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
    12801514
    12811515  add Branch JetEnergyScale/jets Jet Jet
    … …  
    12891523  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    12901524  add Branch ScalarHT/energy ScalarHT ScalarHT
    1291 }
     1525
     1526}

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl

    r20ca0cd r9327245  
    3333  PhotonEnergySmearing
    3434  ElectronFilter
     35
    3536  TrackPileUpSubtractor
     37  RecoPuFilter
    3638
    3739  TowerMerger
    3840  NeutralEFlowMerger
    39   EFlowMergerAllTracks
     41
    4042  EFlowMerger
     43  EFlowMergerCHS
     44  Rho
    4145
    4246  LeptonFilterNoLep
    … …  
    4852
    4953  PhotonIsolation
     54  PhotonIsolationCHS
    5055  PhotonEfficiency
     56  PhotonEfficiencyCHS
    5157
    5258  ElectronIsolation
     59  ElectronIsolationCHS
     60
    5361  ElectronEfficiency
     62  ElectronEfficiencyCHS
    5463
    5564  MuonIsolation
     65  MuonIsolationCHS
     66
    5667  MuonLooseIdEfficiency
    5768  MuonTightIdEfficiency
     69
     70  MuonLooseIdEfficiencyCHS
     71  MuonTightIdEfficiencyCHS
    5872
    5973  NeutrinoFilter
    … …  
    6882  FastJetFinder
    6983  FastJetFinderAK8
     84  JetPileUpSubtractor
     85  JetPileUpSubtractorAK8
    7086  FastJetFinderPUPPI
    7187  FastJetFinderPUPPIAK8
    … …  
    126142  set TVertexSpread 800E-12
    127143
    128   # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
     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))}
    129145  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    130146
    … …  
    138154module ParticlePropagator ParticlePropagator {
    139155  set InputArray PileUpMerger/stableParticles
     156  #set InputArray Delphes/stableParticles
    140157
    141158  set OutputArray stableParticles
    … …  
    272289  source muonMomentumResolution.tcl
    273290}
    274 
    275 
    276 
    277291
    278292##############
    … …  
    526540}
    527541
     542########################
     543# Reco PU filter
     544########################
     545
     546module RecoPuFilter RecoPuFilter {
     547  set InputArray HCal/eflowTracks
     548  set OutputArray eflowTracks
     549}
    528550
    529551###################################################
    … …  
    538560}
    539561
    540 
    541562####################
    542563# Neutral eflow erger
    … …  
    550571}
    551572
    552 
    553 ####################
     573#####################
    554574# Energy flow merger
    555 ####################
     575#####################
    556576
    557577module Merger EFlowMerger {
    … …  
    563583}
    564584
    565 ##################################
    566 # Energy flow merger (all tracks)
    567 ##################################
    568 
    569 module Merger EFlowMergerAllTracks {
     585############################
     586# Energy flow merger no PU
     587############################
     588
     589module Merger EFlowMergerCHS {
    570590# add InputArray InputArray
    571   add InputArray TrackMerger/tracks
     591  add InputArray RecoPuFilter/eflowTracks
    572592  add InputArray PhotonEnergySmearing/eflowPhotons
    573593  add InputArray HCal/eflowNeutralHadrons
    … …  
    694714}
    695715
    696 
    697716#####################
    698717# MC truth jet finder
    … …  
    735754
    736755
    737 
    738 ############
     756#############
     757# Rho pile-up
     758#############
     759
     760module 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##############
    739786# Jet finder
    740 ############
     787##############
    741788
    742789module FastJetFinder FastJetFinder {
    743790#  set InputArray TowerMerger/towers
    744   set InputArray EFlowMerger/eflow
     791  set InputArray EFlowMergerCHS/eflow
    745792
    746793  set OutputArray jets
     794
     795  set AreaAlgorithm 5
    747796
    748797  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    756805module FastJetFinder FastJetFinderAK8 {
    757806#  set InputArray TowerMerger/towers
    758   set InputArray EFlowMerger/eflow
     807  set InputArray EFlowMergerCHS/eflow
    759808
    760809  set OutputArray jets
     810
     811  set AreaAlgorithm 5
    761812
    762813  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    785836}
    786837
     838###########################
     839# Jet Pile-Up Subtraction
     840###########################
     841
     842module 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
     855module 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
    787864module FastJetFinder FastJetFinderPUPPI {
    788865#  set InputArray TowerMerger/towers
    … …  
    834911
    835912module EnergyScale JetEnergyScale {
    836   set InputArray FastJetFinder/jets
     913  set InputArray JetPileUpSubtractor/jets
    837914  set OutputArray jets
    838915
    … …  
    842919
    843920module EnergyScale JetEnergyScaleAK8 {
    844   set InputArray FastJetFinderAK8/jets
     921  set InputArray JetPileUpSubtractorAK8/jets
    845922  set OutputArray jets
    846923
    … …  
    908985}
    909986
     987
     988########################
     989# Photon isolation CHS #
     990########################
     991
     992module 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}
    9101013
    9111014
    … …  
    9301033
    9311034
     1035#####################
     1036# Photon efficiency #
     1037#####################
     1038
     1039module 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
    9321054######################
    9331055# Electron isolation #
    … …  
    9501072}
    9511073
     1074
     1075##########################
     1076# Electron isolation CHS #
     1077##########################
     1078
     1079module 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}
    9521093
    9531094
    … …  
    9961137}
    9971138
     1139###########################
     1140# Electron efficiency CHS #
     1141###########################
     1142
     1143module 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
    9981184##################
    9991185# Muon isolation #
    … …  
    10141200}
    10151201
    1016 
    1017 ##################
    1018 # Muon Loose Id  #
    1019 ##################
     1202######################
     1203# Muon isolation CHS #
     1204######################
     1205
     1206module 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#####################
    10201224
    10211225module Efficiency MuonLooseIdEfficiency {
    … …  
    10391243}
    10401244
     1245
     1246#####################
     1247# Muon Loose Id CHS #
     1248#####################
     1249
     1250module 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
     1263module Efficiency MuonTightIdEfficiencyCHS {
     1264    set InputArray MuonIsolationCHS/muons
     1265    set OutputArray muons
     1266    # tracking + TightID efficiency formula for muons
     1267    source muonTightId.tcl
     1268}
    10411269
    10421270
    … …  
    12511479module StatusPidFilter GenParticleFilter {
    12521480
    1253     set InputArray  Delphes/allParticles
     1481    set InputArray Delphes/allParticles
    12541482    set OutputArray filteredParticles
    12551483    set PTMin 5.0
    … …  
    12791507  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    12801508  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
    12811514
    12821515  add Branch JetEnergyScale/jets Jet Jet
    … …  
    12901523  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    12911524  add Branch ScalarHT/energy ScalarHT ScalarHT
    1292 }
     1525
     1526}

  • modules/ModulesLinkDef.h

    r20ca0cd r9327245  
    5858#include "modules/StatusPidFilter.h"
    5959#include "modules/PdgCodeFilter.h"
     60#include "modules/RecoPuFilter.h"
    6061#include "modules/Cloner.h"
    6162#include "modules/Weighter.h"
    … …  
    103104#pragma link C++ class StatusPidFilter+;
    104105#pragma link C++ class PdgCodeFilter+;
     106#pragma link C++ class RecoPuFilter+;
    105107#pragma link C++ class Cloner+;
    106108#pragma link C++ class Weighter+;

  • modules/PhotonConversions.cc

    r20ca0cd r9327245  
    5959  fItInputArray(0), fConversionMap(0), fDecayXsec(0)
    6060{
    61   fDecayXsec = new TF1;
     61  fDecayXsec = new TF1("decayXsec","1.0 - 4.0/3.0 * x * (1.0 - x)", 0.0, 1.0);
    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);
    9386
    9487  // import array with output from filter/classifier module

  • readers/DelphesCMSFWLite.cpp

    r20ca0cd r9327245  
    6565  ExRootTreeBranch *branchEvent, ExRootTreeBranch *branchRwgt,
    6666  DelphesFactory *factory, TObjArray *allParticleOutputArray,
    67   TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray)
     67  TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray, Bool_t firstEvent)
    6868{
     69
    6970  fwlite::Handle< GenEventInfoProduct > handleGenEventInfo;
    70 
    7171  fwlite::Handle< LHEEventProduct > handleLHEEvent;
    72 
    7372  fwlite::Handle< vector< reco::GenParticle > > handleParticle;
     73
    7474  vector< reco::GenParticle >::const_iterator itParticle;
    7575
    … …  
    7878
    7979  handleGenEventInfo.getByLabel(event, "generator");
    80   handleLHEEvent.getByLabel(event, "externalLHEProducer");
    81   handleParticle.getByLabel(event, "genParticles");
     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());
    82109
    83110  HepMCEvent *element;
    … …  
    91118  Double_t px, py, pz, e, mass;
    92119  Double_t x, y, z;
    93 
    94   const vector< gen::WeightsInfo > &vectorWeightsInfo = handleLHEEvent->weights();
    95   vector< gen::WeightsInfo >::const_iterator itWeightsInfo;
    96120
    97121  element = static_cast<HepMCEvent *>(branchEvent->NewEntry());
    … …  
    117141  element->ProcTime = 0.0;
    118142
    119   for(itWeightsInfo = vectorWeightsInfo.begin(); itWeightsInfo != vectorWeightsInfo.end(); ++itWeightsInfo)
    120   {
    121     weight = static_cast<Weight *>(branchRwgt->NewEntry());
    122     weight->Weight = itWeightsInfo->wgt;
     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    } 
    123154  }
    124155
    … …  
    207238  Int_t i;
    208239  Long64_t eventCounter, numberOfEvents;
     240  Bool_t firstEvent = kTRUE;
    209241
    210242  if(argc < 4)
    … …  
    281313      modularDelphes->Clear();
    282314      treeWriter->Clear();
     315
    283316      for(event.toBegin(); !event.atEnd() && !interrupted; ++event)
    284317      {
    285318        ConvertInput(event, eventCounter, branchEvent, branchRwgt, factory,
    286           allParticleOutputArray, stableParticleOutputArray, partonOutputArray);
     319          allParticleOutputArray, stableParticleOutputArray, partonOutputArray, firstEvent);
    287320        modularDelphes->ProcessTask();
     321         
     322        firstEvent = kFALSE;
    288323
    289324        treeWriter->Fill();
Note: See TracChangeset for help on using the changeset viewer.