Commit c040d9a9 authored by Oliver Majersky's avatar Oliver Majersky Committed by Nils Erik Krumnack
Browse files

AnalysisTop: Update large-R jet tagging uncertainties

parent 29e97740
......@@ -265,7 +265,7 @@ namespace top {
"CategoryJES_FullJER_FullJMS");
registerParameter("LargeRJetUncertainties_JMR_NPModel",
"Full (10 NP; aimed at the most precise jet-dependent measurements), Simple (1 NP; flat 20 percent uncertainty, as it was recommended in the past) - for LCTopo large-R jet mass resolution uncertainties",
"FullJMR");
"FullJMR_COMB");
registerParameter("AdvancedUsage_LargeRJetUncertaintiesConfigDir",
"Path to directory containing large-R jet uncertainties config",
"rel21/Winter2021");
......
......@@ -556,6 +556,10 @@ namespace top {
const SG::AuxElement::Accessor< char > accRange("passedRangeCheck_" + fullName);
const std::string sfNameNominal = sfNames.at(fullName);
const SG::AuxElement::Accessor< float > accSF(sfNameNominal);
// accessor to retrieve nominal efficiency decoration from BoostedJetTaggers
std::string taggerName = sfNameNominal;
taggerName.erase(taggerName.length()-3);
const SG::AuxElement::Accessor<float> accEff(taggerName + "_efficiency");
for(const CP::SystematicSet& sys : m_tagSFUncorrelatedSystematics[fullName]) {
......@@ -579,6 +583,13 @@ namespace top {
top::check(tool->applyCorrection(*shallowJet), "Failed to applyCorrection");
float sf = accSF.isAvailable(*shallowJet) ? accSF(*shallowJet) : -999.;
jet->auxdecor<float>(sfNameShifted.c_str()) = sf;
// decorate efficiencies for inefficiency SF variations
if (sys.name().find("TagEffUnc") != std::string::npos) {
if (accEff.isAvailable(*shallowJet)) {
jet->auxdecor<float>(fullName + "_" + sys.name() + "_efficiency") = accEff(*shallowJet);
}
}
}
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment