File indexing completed on 2023-03-17 11:26:10
0001 #ifndef TtSemiLepJetCombEval_h
0002 #define TtSemiLepJetCombEval_h
0003
0004 #include "PhysicsTools/MVAComputer/interface/Variable.h"
0005 #include "TopQuarkAnalysis/TopJetCombination/interface/TtSemiLepJetComb.h"
0006
0007
0008
0009
0010
0011
0012
0013 inline void evaluateTtSemiLepJetComb(PhysicsTools::Variable::ValueList& values, const TtSemiLepJetComb& jetComb) {
0014
0015
0016
0017
0018 values.add("massHadTop", jetComb.topVar(JetComb::kHad, JetComb::kMass));
0019 values.add("ptHadTop", jetComb.topVar(JetComb::kHad, JetComb::kPt));
0020 values.add("etaHadTop", jetComb.topVar(JetComb::kHad, JetComb::kEta));
0021 values.add("phiHadTop", jetComb.topVar(JetComb::kHad, JetComb::kPhi));
0022 values.add("thetaHadTop", jetComb.topVar(JetComb::kHad, JetComb::kTheta));
0023
0024 values.add("massLepTop", jetComb.topVar(JetComb::kLep, JetComb::kMass));
0025 values.add("ptLepTop", jetComb.topVar(JetComb::kLep, JetComb::kPt));
0026 values.add("etaLepTop", jetComb.topVar(JetComb::kLep, JetComb::kEta));
0027 values.add("phiLepTop", jetComb.topVar(JetComb::kLep, JetComb::kPhi));
0028 values.add("thetaLepTop", jetComb.topVar(JetComb::kLep, JetComb::kTheta));
0029
0030 values.add("massHadW", jetComb.wBosonVar(JetComb::kHad, JetComb::kMass));
0031 values.add("ptHadW", jetComb.wBosonVar(JetComb::kHad, JetComb::kPt));
0032 values.add("etaHadW", jetComb.wBosonVar(JetComb::kHad, JetComb::kEta));
0033 values.add("phiHadW", jetComb.wBosonVar(JetComb::kHad, JetComb::kPhi));
0034 values.add("thetaHadW", jetComb.wBosonVar(JetComb::kHad, JetComb::kTheta));
0035
0036 values.add("ptHadB", jetComb.bQuarkVar(JetComb::kHad, JetComb::kPt));
0037 values.add("etaHadB", jetComb.bQuarkVar(JetComb::kHad, JetComb::kEta));
0038 values.add("phiHadB", jetComb.bQuarkVar(JetComb::kHad, JetComb::kPhi));
0039 values.add("thetaHadB", jetComb.bQuarkVar(JetComb::kHad, JetComb::kTheta));
0040
0041 values.add("ptLepB", jetComb.bQuarkVar(JetComb::kLep, JetComb::kPt));
0042 values.add("etaLepB", jetComb.bQuarkVar(JetComb::kLep, JetComb::kEta));
0043 values.add("phiLepB", jetComb.bQuarkVar(JetComb::kLep, JetComb::kPhi));
0044 values.add("thetaLepB", jetComb.bQuarkVar(JetComb::kLep, JetComb::kTheta));
0045
0046 values.add("ptLightQ", jetComb.lightQVar(false, JetComb::kPt));
0047 values.add("etaLightQ", jetComb.lightQVar(false, JetComb::kEta));
0048 values.add("phiLightQ", jetComb.lightQVar(false, JetComb::kPhi));
0049 values.add("thetaLightQ", jetComb.lightQVar(false, JetComb::kTheta));
0050
0051 values.add("ptLightQBar", jetComb.lightQVar(true, JetComb::kPt));
0052 values.add("etaLightQBar", jetComb.lightQVar(true, JetComb::kEta));
0053 values.add("phiLightQBar", jetComb.lightQVar(true, JetComb::kPhi));
0054 values.add("thetaLightQBar", jetComb.lightQVar(true, JetComb::kTheta));
0055
0056
0057
0058
0059 values.add("deltaMHadTopLepTop", jetComb.compareHadTopLepTop(JetComb::kDeltaM));
0060 values.add("deltaRHadTopLepTop", jetComb.compareHadTopLepTop(JetComb::kDeltaR));
0061 values.add("deltaPhiHadTopLepTop", jetComb.compareHadTopLepTop(JetComb::kDeltaPhi));
0062 values.add("deltaThetaHadTopLepTop", jetComb.compareHadTopLepTop(JetComb::kDeltaTheta));
0063
0064 values.add("deltaMHadWLepW", jetComb.compareHadWLepW(JetComb::kDeltaM));
0065 values.add("deltaRHadWLepW", jetComb.compareHadWLepW(JetComb::kDeltaR));
0066 values.add("deltaPhiHadWLepW", jetComb.compareHadWLepW(JetComb::kDeltaPhi));
0067 values.add("deltaThetaHadWLepW", jetComb.compareHadWLepW(JetComb::kDeltaTheta));
0068
0069 values.add("deltaRHadBLepB", jetComb.compareHadBLepB(JetComb::kDeltaR));
0070 values.add("deltaPhiHadBLepB", jetComb.compareHadBLepB(JetComb::kDeltaPhi));
0071 values.add("deltaThetaHadBLepB", jetComb.compareHadBLepB(JetComb::kDeltaTheta));
0072
0073 values.add("deltaRHadQHadQBar", jetComb.compareLightQuarks(JetComb::kDeltaR));
0074 values.add("deltaPhiHadQHadQBar", jetComb.compareLightQuarks(JetComb::kDeltaPhi));
0075 values.add("deltaThetaHadQHadQBar", jetComb.compareLightQuarks(JetComb::kDeltaTheta));
0076
0077 values.add("deltaMHadTopHadW", jetComb.compareTopW(JetComb::kHad, JetComb::kHad, JetComb::kDeltaM));
0078 values.add("deltaRHadTopHadW", jetComb.compareTopW(JetComb::kHad, JetComb::kHad, JetComb::kDeltaR));
0079 values.add("deltaPhiHadTopHadW", jetComb.compareTopW(JetComb::kHad, JetComb::kHad, JetComb::kDeltaPhi));
0080 values.add("deltaThetaHadTopHadW", jetComb.compareTopW(JetComb::kHad, JetComb::kHad, JetComb::kDeltaTheta));
0081
0082 values.add("deltaMLepTopLepW", jetComb.compareTopW(JetComb::kLep, JetComb::kLep, JetComb::kDeltaM));
0083 values.add("deltaRLepTopLepW", jetComb.compareTopW(JetComb::kLep, JetComb::kLep, JetComb::kDeltaR));
0084 values.add("deltaPhiLepTopLepW", jetComb.compareTopW(JetComb::kLep, JetComb::kLep, JetComb::kDeltaPhi));
0085 values.add("deltaThetaLepTopLepW", jetComb.compareTopW(JetComb::kLep, JetComb::kLep, JetComb::kDeltaTheta));
0086
0087 values.add("deltaMHadTopLepW", jetComb.compareTopW(JetComb::kHad, JetComb::kLep, JetComb::kDeltaM));
0088 values.add("deltaRHadTopLepW", jetComb.compareTopW(JetComb::kHad, JetComb::kLep, JetComb::kDeltaR));
0089 values.add("deltaPhiHadTopLepW", jetComb.compareTopW(JetComb::kHad, JetComb::kLep, JetComb::kDeltaPhi));
0090 values.add("deltaThetaHadTopLepW", jetComb.compareTopW(JetComb::kHad, JetComb::kLep, JetComb::kDeltaTheta));
0091
0092 values.add("deltaMLepTopHadW", jetComb.compareTopW(JetComb::kLep, JetComb::kHad, JetComb::kDeltaM));
0093 values.add("deltaRLepTopHadW", jetComb.compareTopW(JetComb::kLep, JetComb::kHad, JetComb::kDeltaR));
0094 values.add("deltaPhiLepTopHadW", jetComb.compareTopW(JetComb::kLep, JetComb::kHad, JetComb::kDeltaPhi));
0095 values.add("deltaThetaLepTopHadW", jetComb.compareTopW(JetComb::kLep, JetComb::kHad, JetComb::kDeltaTheta));
0096
0097 values.add("deltaRHadTopHadB", jetComb.compareTopB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaR));
0098 values.add("deltaPhiHadTopHadB", jetComb.compareTopB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaPhi));
0099 values.add("deltaThetaHadTopHadB", jetComb.compareTopB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaTheta));
0100
0101 values.add("deltaRLepTopLepB", jetComb.compareTopB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaR));
0102 values.add("deltaPhiLepTopLepB", jetComb.compareTopB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaPhi));
0103 values.add("deltaThetaLepTopLepB", jetComb.compareTopB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaTheta));
0104
0105 values.add("deltaRHadTopLepB", jetComb.compareTopB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaR));
0106 values.add("deltaPhiHadTopLepB", jetComb.compareTopB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaPhi));
0107 values.add("deltaThetaHadTopLepB", jetComb.compareTopB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaTheta));
0108
0109 values.add("deltaRLepTopHadB", jetComb.compareTopB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaR));
0110 values.add("deltaPhiLepTopHadB", jetComb.compareTopB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaPhi));
0111 values.add("deltaThetaLepTopHadB", jetComb.compareTopB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaTheta));
0112
0113 values.add("deltaRHadWHadB", jetComb.compareWB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaR));
0114 values.add("deltaPhiHadWHadB", jetComb.compareWB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaPhi));
0115 values.add("deltaThetaHadWHadB", jetComb.compareWB(JetComb::kHad, JetComb::kHad, JetComb::kDeltaTheta));
0116
0117 values.add("deltaRLepWLepB", jetComb.compareWB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaR));
0118 values.add("deltaPhiLepWLepB", jetComb.compareWB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaPhi));
0119 values.add("deltaThetaLepWLepB", jetComb.compareWB(JetComb::kLep, JetComb::kLep, JetComb::kDeltaTheta));
0120
0121 values.add("deltaRHadWLepB", jetComb.compareWB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaR));
0122 values.add("deltaPhiHadWLepB", jetComb.compareWB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaPhi));
0123 values.add("deltaThetaHadWLepB", jetComb.compareWB(JetComb::kHad, JetComb::kLep, JetComb::kDeltaTheta));
0124
0125 values.add("deltaRLepWHadB", jetComb.compareWB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaR));
0126 values.add("deltaPhiLepWHadB", jetComb.compareWB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaPhi));
0127 values.add("deltaThetaLepWHadB", jetComb.compareWB(JetComb::kLep, JetComb::kHad, JetComb::kDeltaTheta));
0128
0129 values.add("deltaRHadTopLepton", jetComb.compareTopLepton(JetComb::kHad, JetComb::kDeltaR));
0130 values.add("deltaPhiHadTopLepton", jetComb.compareTopLepton(JetComb::kHad, JetComb::kDeltaPhi));
0131 values.add("deltaThetaHadTopLepton", jetComb.compareTopLepton(JetComb::kHad, JetComb::kDeltaTheta));
0132
0133 values.add("deltaRLepTopLepton", jetComb.compareTopLepton(JetComb::kLep, JetComb::kDeltaR));
0134 values.add("deltaPhiLepTopLepton", jetComb.compareTopLepton(JetComb::kLep, JetComb::kDeltaPhi));
0135 values.add("deltaThetaLepTopLepton", jetComb.compareTopLepton(JetComb::kLep, JetComb::kDeltaTheta));
0136
0137 values.add("deltaRHadTopNeutrino", jetComb.compareTopNeutrino(JetComb::kHad, JetComb::kDeltaR));
0138 values.add("deltaPhiHadTopNeutrino", jetComb.compareTopNeutrino(JetComb::kHad, JetComb::kDeltaPhi));
0139 values.add("deltaThetaHadTopNeutrino", jetComb.compareTopNeutrino(JetComb::kHad, JetComb::kDeltaTheta));
0140
0141 values.add("deltaRLepTopNeutrino", jetComb.compareTopNeutrino(JetComb::kLep, JetComb::kDeltaR));
0142 values.add("deltaPhiLepTopNeutrino", jetComb.compareTopNeutrino(JetComb::kLep, JetComb::kDeltaPhi));
0143 values.add("deltaThetaLepTopNeutrino", jetComb.compareTopNeutrino(JetComb::kLep, JetComb::kDeltaTheta));
0144
0145 values.add("deltaRHadWLepton", jetComb.compareWLepton(JetComb::kHad, JetComb::kDeltaR));
0146 values.add("deltaPhiHadWLepton", jetComb.compareWLepton(JetComb::kHad, JetComb::kDeltaPhi));
0147 values.add("deltaThetaHadWLepton", jetComb.compareWLepton(JetComb::kHad, JetComb::kDeltaTheta));
0148
0149 values.add("deltaRHadWNeutrino", jetComb.compareWNeutrino(JetComb::kHad, JetComb::kDeltaR));
0150 values.add("deltaPhiHadWNeutrino", jetComb.compareWNeutrino(JetComb::kHad, JetComb::kDeltaPhi));
0151 values.add("deltaThetaHadWNeutrino", jetComb.compareWNeutrino(JetComb::kHad, JetComb::kDeltaTheta));
0152
0153 values.add("deltaRHadBLepton", jetComb.compareBLepton(JetComb::kHad, JetComb::kDeltaR));
0154 values.add("deltaPhiHadBLepton", jetComb.compareBLepton(JetComb::kHad, JetComb::kDeltaPhi));
0155 values.add("deltaThetaHadBLepton", jetComb.compareBLepton(JetComb::kHad, JetComb::kDeltaTheta));
0156
0157 values.add("deltaRLepBLepton", jetComb.compareBLepton(JetComb::kLep, JetComb::kDeltaR));
0158 values.add("deltaPhiLepBLepton", jetComb.compareBLepton(JetComb::kLep, JetComb::kDeltaPhi));
0159 values.add("deltaThetaLepBLepton", jetComb.compareBLepton(JetComb::kLep, JetComb::kDeltaTheta));
0160
0161 values.add("deltaRHadBNeutrino", jetComb.compareBNeutrino(JetComb::kHad, JetComb::kDeltaR));
0162 values.add("deltaPhiHadBNeutrino", jetComb.compareBNeutrino(JetComb::kHad, JetComb::kDeltaPhi));
0163 values.add("deltaThetaHadBNeutrino", jetComb.compareBNeutrino(JetComb::kHad, JetComb::kDeltaTheta));
0164
0165 values.add("deltaRLepBNeutrino", jetComb.compareBNeutrino(JetComb::kLep, JetComb::kDeltaR));
0166 values.add("deltaPhiLepBNeutrino", jetComb.compareBNeutrino(JetComb::kLep, JetComb::kDeltaPhi));
0167 values.add("deltaThetaLepBNeutrino", jetComb.compareBNeutrino(JetComb::kLep, JetComb::kDeltaTheta));
0168
0169
0170
0171 values.add("relativePtHadronicTop", jetComb.relativePtHadronicTop());
0172 values.add("bOverLightQPt", jetComb.bOverLightQPt());
0173
0174
0175
0176
0177 values.add("bTagHadBTrkCntHighEff", jetComb.bTag(JetComb::kHad, JetComb::kTrackCountHighEff));
0178 values.add("bTagHadBTrkCntHighPur", jetComb.bTag(JetComb::kHad, JetComb::kTrackCountHighPur));
0179 values.add("bTagHadBSoftMuon", jetComb.bTag(JetComb::kHad, JetComb::kSoftMuon));
0180 values.add("bTagHadBSimpSecondVtx", jetComb.bTag(JetComb::kHad, JetComb::kSimpleSecondVtx));
0181 values.add("bTagHadBCombSecondVtx", jetComb.bTag(JetComb::kHad, JetComb::kCombSecondVtx));
0182 values.add("bTagHadBImpactParaMVA", jetComb.bTag(JetComb::kHad, JetComb::kCombSecondVtxMVA));
0183
0184 values.add("bTagLepBTrkCntHighEff", jetComb.bTag(JetComb::kLep, JetComb::kTrackCountHighEff));
0185 values.add("bTagLepBTrkCntHighPur", jetComb.bTag(JetComb::kLep, JetComb::kTrackCountHighPur));
0186 values.add("bTagLepBSoftMuon", jetComb.bTag(JetComb::kLep, JetComb::kSoftMuon));
0187 values.add("bTagLepBSimpSecondVtx", jetComb.bTag(JetComb::kLep, JetComb::kSimpleSecondVtx));
0188 values.add("bTagLepBCombSecondVtx", jetComb.bTag(JetComb::kLep, JetComb::kCombSecondVtx));
0189 values.add("bTagLepBImpactParaMVA", jetComb.bTag(JetComb::kLep, JetComb::kCombSecondVtxMVA));
0190
0191 values.add("bTagSumTrkCntHighEff", jetComb.combinedBTags(JetComb::kTrackCountHighEff, JetComb::kAdd));
0192 values.add("bTagSumTrkCntHighPur", jetComb.combinedBTags(JetComb::kTrackCountHighPur, JetComb::kAdd));
0193 values.add("bTagSumSoftMuon", jetComb.combinedBTags(JetComb::kSoftMuon, JetComb::kAdd));
0194 values.add("bTagSumSimpSecondVtx", jetComb.combinedBTags(JetComb::kSimpleSecondVtx, JetComb::kAdd));
0195 values.add("bTagSumCombSecondVtx", jetComb.combinedBTags(JetComb::kCombSecondVtx, JetComb::kAdd));
0196 values.add("bTagSumImpactParaMVA", jetComb.combinedBTags(JetComb::kCombSecondVtxMVA, JetComb::kAdd));
0197
0198 values.add("bTagProdImpactParaMVA", jetComb.combinedBTags(JetComb::kCombSecondVtxMVA, JetComb::kMult));
0199
0200 values.add("bTagSumHadQHadQBarTrkCntHighEff",
0201 jetComb.combinedBTagsForLightQuarks(JetComb::kTrackCountHighEff, JetComb::kAdd));
0202 values.add("bTagSumHadQHadQBarTrkCntHighPur",
0203 jetComb.combinedBTagsForLightQuarks(JetComb::kTrackCountHighPur, JetComb::kAdd));
0204 values.add("bTagSumHadQHadQBarSoftMuon", jetComb.combinedBTagsForLightQuarks(JetComb::kSoftMuon, JetComb::kAdd));
0205 values.add("bTagSumHadQHadQBarSimpSecondVtx",
0206 jetComb.combinedBTagsForLightQuarks(JetComb::kSimpleSecondVtx, JetComb::kAdd));
0207 values.add("bTagSumHadQHadQBarCombSecondVtx",
0208 jetComb.combinedBTagsForLightQuarks(JetComb::kCombSecondVtx, JetComb::kAdd));
0209 values.add("bTagSumHadQHadQBarImpactParaMVA",
0210 jetComb.combinedBTagsForLightQuarks(JetComb::kCombSecondVtxMVA, JetComb::kAdd));
0211
0212 values.add("bTagProdHadQHadQBarImpactParaMVA",
0213 jetComb.combinedBTagsForLightQuarks(JetComb::kCombSecondVtxMVA, JetComb::kMult));
0214 }
0215
0216 #endif