def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [
calculables.HT(),
calculables.rho(),
calculables.jecFactor(*pars["jec"]),
calculables.Filtered(pids=[-5, 5],
label="b",
key="Particle",
status=[3]),
calculables.Filtered(pids=[-15, 15],
label="tau",
key="Particle",
status=[3]),
calculables.Filtered(label="",
ptMin=20.0,
absEtaMax=2.4,
key="Jet"),
calculables.bTagged("Jets", mask=0x1),
calculables.tauTagged("Jets"),
calculables.tauTagged("bTagged_Jets_mask1"),
calculables.Duplicates(key1="bTagged_Jets_mask1",
key2="tauTagged_Jets",
minDR=0.2),
calculables.DeltaR("bTagged_Jets_mask1"),
calculables.DeltaR("tauTagged_Jets"),
# mbb
calculables.SumP4("bTagged_Jets_mask1"),
calculables.JetsFixedMass("bTagged_Jets_mask1",
m=mb,
label="bTagged"),
calculables.SumP4("JetsFixedMass_bTagged"),
calculables.JetsFixedMass("bTagged_Jets_mask1",
m=mb,
correctPt=True,
label="bTagged"),
calculables.SumP4("JetsFixedMass_bTagged_Corrected"),
# mtt
calculables.SumP4("tauTagged_Jets"),
calculables.JetsFixedMass("tauTagged_Jets",
m=mtau,
label="tauTagged"),
calculables.SumP4("JetsFixedMass_tauTagged"),
calculables.JetsFixedMass("tauTagged_Jets",
m=mtau,
correctPt=True,
label="tauTagged"),
calculables.SumP4("JetsFixedMass_tauTagged_Corrected"),
# mbbtt
calculables.LvSumP4(label="bbtt",
keys=[
"JetsFixedMass_bTagged_Corrected_SumP4",
"JetsFixedMass_tauTagged_Corrected_SumP4"
]),
]
return listOfCalculables
def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
if "checkJec" in pars:
listOfCalculables.append(calculables.jecFactor(*pars["checkJec"]))
listOfCalculables += [calculables.HT(),
calculables.rho(),
calculables.Filtered(pids=[23], label="Z", key="Particle", status=[3]),
calculables.Filtered(pids=[25], label="h", key="Particle", status=[3], indices=True),
calculables.Filtered(pids=[21], label="g", key="Particle", status=[3], indices=True, ptMin=30.0),
calculables.Following(key="Particle",
motherIndices="hParticle_Indices",
daughterIndices="gParticle_Indices",
nMax=2, shortName=True),
calculables.SumP4("gParticles"),
calculables.JetMatchedTo(sourceKey="gParticles"),
calculables.JetMatchedTo(sourceKey="gParticles", maxDR=0.3),
supy.calculables.other.values("JetMatchedTo_gParticles"),
calculables.SumP4("JetMatchedTo_gParticles.values"),
]
listOfCalculables += [# b
calculables.Filtered(pids=[-5, 5], label="b", status=[3], key="Particle", ptMin=0.0),
calculables.minPt("bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", maxDR=0.3),
calculables.SumP4("bParticles"),
supy.calculables.other.values("JetMatchedTo_bParticles"),
calculables.DeltaR("JetMatchedTo_bParticles.values"),
# mass 1
calculables.SumP4("JetMatchedTo_bParticles.values"),
# mass 2
calculables.JetsFixedMass("JetMatchedTo_bParticles.values", m=mb, label="bMatched"),
calculables.SumP4("JetsFixedMass_bMatched"),
calculables.JetsFixedMass("JetMatchedTo_bParticles.values", m=mb, correctPt=True, label="bMatched"),
calculables.SumP4("JetsFixedMass_bMatched_Corrected"),
]
listOfCalculables += [# tau
calculables.Filtered(pids=[-15, 15], label="tau", key="Particle", status=[3]),
calculables.Duplicates(key1="bParticles", key2="tauParticles", minDR=0.5),
calculables.JetMatchedTo(sourceKey="tauParticles"),
calculables.JetMatchedTo(sourceKey="tauParticles", maxDR=0.3),
calculables.SumP4("tauParticles"),
supy.calculables.other.values("JetMatchedTo_tauParticles"),
calculables.DeltaR("JetMatchedTo_tauParticles.values"),
# mass 1
calculables.SumP4("JetMatchedTo_tauParticles.values"),
# mass 2
calculables.JetsFixedMass("JetMatchedTo_tauParticles.values",
m=mtau, label="tauMatched"),
calculables.SumP4("JetsFixedMass_tauMatched"),
calculables.JetsFixedMass("JetMatchedTo_tauParticles.values",
m=mtau, correctPt=True, label="tauMatched"),
calculables.SumP4("JetsFixedMass_tauMatched_Corrected"),
]
return listOfCalculables
def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [supy.calculables.other.fixedValue("one", 1.0),
calculables.Filtered(pids=[-3, -2, -1, 1, 2, 3], label="uds", status=[3]),
calculables.Filtered(pids=[-4, 4], label="c", status=[3]),
calculables.Filtered(pids=[-5, 5], label="b", status=[3]),
calculables.Filtered(pids=[-6, 6], label="t", status=[3]),
calculables.Filtered(pids=[-11, 11], label="e", status=[3]),
calculables.Filtered(pids=[-13, 13], label="mu", status=[3]),
calculables.Filtered(pids=[-15, 15], label="tau", status=[3]),
calculables.Filtered(pids=[-16, -14, -12, 12, 14, 16], label="nu", status=[3]),
calculables.Filtered(pids=[21], label="g", status=[3]),
calculables.Filtered(pids=[25], label="h", status=[3]),
calculables.Filtered(pids=[-15, -5, 5, 15], label="btau", status=[3]),
calculables.Filtered(pids=[-5, 5], label="bKine", ptMin=30.0, absEtaMax=2.4, status=[3]),
calculables.Filtered(pids=[-15, 15], label="tauKine", ptMin=30.0, absEtaMax=2.4, status=[3]),
calculables.Filtered(pids=[-15, -5, 5, 15], label="btauKine", ptMin=30.0, absEtaMax=2.4, status=[3]),
calculables.Filtered(label="10GeV", ptMin=10.0, key="Track", ptSort=True),
calculables.Filtered(label="15GeV", ptMin=15.0, key="Tower", ptSort=True, pt="ET"),
#calculables.JetMatchedTo(sourceKey="udsParticles"),
calculables.JetMatchedTo(sourceKey="cParticles"),
calculables.JetMatchedTo(sourceKey="bParticles"),
calculables.JetMatchedTo(sourceKey="tauParticles"),
#calculables.JetMatchedTo(sourceKey="gParticles"),
calculables.bTaggedJets(),
calculables.HT(),
calculables.rho(),
calculables.SumP4("bParticles"),
calculables.SumP4("tauParticles"),
calculables.SumP4("btauParticles"),
]
return listOfCalculables
def listOfCalculables(self, pars):
maxDR = 0.3
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [calculables.HT(),
calculables.rho(),
calculables.JetMatchedTo(sourceKey="GenJet"),
calculables.JetMatchedTo(sourceKey="GenJet", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_GenJet"),
supy.calculables.other.keys("JetMatchedTo_GenJet"),
calculables.Filtered(label="70", key="GenJet", ptMin=70.0),
calculables.JetMatchedTo(sourceKey="70GenJets"),
calculables.JetMatchedTo(sourceKey="70GenJets", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_70GenJets"),
supy.calculables.other.keys("JetMatchedTo_70GenJets"),
calculables.Filtered(pids=[-5, 5], label="b", status=[3], key="Particle", ptMin=0.0),
calculables.SumP4("bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_bParticles"),
supy.calculables.other.keys("JetMatchedTo_bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", jetKey="GenJet"),
calculables.JetMatchedTo(sourceKey="bParticles", jetKey="GenJet", maxDR=maxDR),
supy.calculables.other.values("GenJetMatchedTo_bParticles"),
supy.calculables.other.keys("GenJetMatchedTo_bParticles"),
calculables.PtRatio("GenJetMatchedTo_bParticles"),
calculables.minItem("GenJetMatchedTo_bParticles_PtRatio"),
]
return listOfCalculables
def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [calculables.HT(),
calculables.rho(),
calculables.Filtered(label="", ptMin=20.0, absEtaMax=2.4, key="Jet"),
calculables.bTagged("Jets", mask=0x1),
calculables.tauTagged("Jets"),
]
return listOfCalculables
def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [
calculables.HT(),
calculables.rho(),
calculables.Filtered(label="",
ptMin=20.0,
absEtaMax=2.4,
key="Jet"),
calculables.bTagged("Jets", mask=0x1),
calculables.tauTagged("Jets"),
]
return listOfCalculables
def listOfCalculables(self, pars):
maxDR = 0.3
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
listOfCalculables += [
calculables.HT(),
calculables.rho(),
calculables.JetMatchedTo(sourceKey="GenJet"),
calculables.JetMatchedTo(sourceKey="GenJet", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_GenJet"),
supy.calculables.other.keys("JetMatchedTo_GenJet"),
calculables.Filtered(label="70", key="GenJet", ptMin=70.0),
calculables.JetMatchedTo(sourceKey="70GenJets"),
calculables.JetMatchedTo(sourceKey="70GenJets", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_70GenJets"),
supy.calculables.other.keys("JetMatchedTo_70GenJets"),
calculables.Filtered(pids=[-5, 5],
label="b",
status=[3],
key="Particle",
ptMin=0.0),
calculables.SumP4("bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", maxDR=maxDR),
supy.calculables.other.values("JetMatchedTo_bParticles"),
supy.calculables.other.keys("JetMatchedTo_bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", jetKey="GenJet"),
calculables.JetMatchedTo(sourceKey="bParticles",
jetKey="GenJet",
maxDR=maxDR),
supy.calculables.other.values("GenJetMatchedTo_bParticles"),
supy.calculables.other.keys("GenJetMatchedTo_bParticles"),
calculables.PtRatio("GenJetMatchedTo_bParticles"),
calculables.minItem("GenJetMatchedTo_bParticles_PtRatio"),
]
return listOfCalculables
def listOfCalculables(self, pars):
listOfCalculables = supy.calculables.zeroArgs(supy.calculables)
listOfCalculables += supy.calculables.zeroArgs(calculables)
if "checkJec" in pars:
listOfCalculables.append(calculables.jecFactor(*pars["checkJec"]))
listOfCalculables += [
calculables.HT(),
calculables.rho(),
calculables.Filtered(pids=[23],
label="Z",
key="Particle",
status=[3]),
calculables.Filtered(pids=[25],
label="h",
key="Particle",
status=[3],
indices=True),
calculables.Filtered(pids=[21],
label="g",
key="Particle",
status=[3],
indices=True,
ptMin=30.0),
calculables.Following(key="Particle",
motherIndices="hParticle_Indices",
daughterIndices="gParticle_Indices",
nMax=2,
shortName=True),
calculables.SumP4("gParticles"),
calculables.JetMatchedTo(sourceKey="gParticles"),
calculables.JetMatchedTo(sourceKey="gParticles", maxDR=0.3),
supy.calculables.other.values("JetMatchedTo_gParticles"),
calculables.SumP4("JetMatchedTo_gParticles.values"),
]
listOfCalculables += [ # b
calculables.Filtered(pids=[-5, 5],
label="b",
status=[3],
key="Particle",
ptMin=0.0),
calculables.minPt("bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles"),
calculables.JetMatchedTo(sourceKey="bParticles", maxDR=0.3),
calculables.SumP4("bParticles"),
supy.calculables.other.values("JetMatchedTo_bParticles"),
calculables.DeltaR("JetMatchedTo_bParticles.values"),
# mass 1
calculables.SumP4("JetMatchedTo_bParticles.values"),
# mass 2
calculables.JetsFixedMass("JetMatchedTo_bParticles.values",
m=mb,
label="bMatched"),
calculables.SumP4("JetsFixedMass_bMatched"),
calculables.JetsFixedMass("JetMatchedTo_bParticles.values",
m=mb,
correctPt=True,
label="bMatched"),
calculables.SumP4("JetsFixedMass_bMatched_Corrected"),
]
listOfCalculables += [ # tau
calculables.Filtered(pids=[-15, 15],
label="tau",
key="Particle",
status=[3]),
calculables.Duplicates(key1="bParticles",
key2="tauParticles",
minDR=0.5),
calculables.JetMatchedTo(sourceKey="tauParticles"),
calculables.JetMatchedTo(sourceKey="tauParticles", maxDR=0.3),
calculables.SumP4("tauParticles"),
supy.calculables.other.values("JetMatchedTo_tauParticles"),
calculables.DeltaR("JetMatchedTo_tauParticles.values"),
# mass 1
calculables.SumP4("JetMatchedTo_tauParticles.values"),
# mass 2
calculables.JetsFixedMass("JetMatchedTo_tauParticles.values",
m=mtau,
label="tauMatched"),
calculables.SumP4("JetsFixedMass_tauMatched"),
calculables.JetsFixedMass("JetMatchedTo_tauParticles.values",
m=mtau,
correctPt=True,
label="tauMatched"),
calculables.SumP4("JetsFixedMass_tauMatched_Corrected"),
]
return listOfCalculables