From 58300274d2c924f7074854d4d3393a75ac56ba30 Mon Sep 17 00:00:00 2001 From: John Alison Date: Fri, 11 Dec 2015 15:18:14 -0600 Subject: [PATCH 1/2] Protection for FatAssJet Uncertianties --- Root/JetCalibrator.cxx | 13 ++++++++++++- xAODAnaHelpers/JetCalibrator.h | 1 + 2 files changed, 13 insertions(+), 1 deletion(-) diff --git a/Root/JetCalibrator.cxx b/Root/JetCalibrator.cxx index 28ce2397db..c4a6a20258 100644 --- a/Root/JetCalibrator.cxx +++ b/Root/JetCalibrator.cxx @@ -91,6 +91,7 @@ JetCalibrator :: JetCalibrator (std::string className) : m_saveAllCleanDecisions = false; m_jetCleanUgly = false; m_cleanParent = false; + m_applyFatJetPreSel = false; //recalculate JVT using calibrated jets m_redoJVT = false; @@ -146,6 +147,7 @@ EL::StatusCode JetCalibrator :: configure () m_jetCleanUgly = config->GetValue("JetCleanUgly", m_jetCleanUgly ); m_saveAllCleanDecisions = config->GetValue("SaveAllCleanDecisions", m_saveAllCleanDecisions); m_cleanParent = config->GetValue("CleanParent", m_cleanParent); + m_applyFatJetPreSel = config->GetValue("ApplyFatJetPreSel", m_applyFatJetPreSel); m_redoJVT = config->GetValue("RedoJVT", m_redoJVT); @@ -537,8 +539,17 @@ EL::StatusCode JetCalibrator :: execute () Error("execute()", "Cannot configure JetUncertaintiesTool for systematic %s", m_systName.c_str()); return EL::StatusCode::FAILURE; } - for ( auto jet_itr : *(calibJetsSC.first) ) { + + for ( auto jet_itr : *(calibJetsSC.first) ) { if ( m_runSysts ) { + + if (m_applyFatJetPreSel) { + bool validForJES = (jet_itr->pt() >= 150e3 && jet_itr->pt() < 3000e3); + validForJES &= (jet_itr->m()/jet_itr->pt() >= 0 && jet_itr->m()/jet_itr->pt() < 1); + validForJES &= (fabs(jet_itr->eta()) < 2); + if (!validForJES) continue; + } + if ( m_JESUncertTool->applyCorrection( *jet_itr ) == CP::CorrectionCode::Error ) { Error("execute()", "JetUncertaintiesTool reported a CP::CorrectionCode::Error"); Error("execute()", "%s", m_name.c_str()); diff --git a/xAODAnaHelpers/JetCalibrator.h b/xAODAnaHelpers/JetCalibrator.h index 57ceb37a35..b0b90eab4e 100644 --- a/xAODAnaHelpers/JetCalibrator.h +++ b/xAODAnaHelpers/JetCalibrator.h @@ -58,6 +58,7 @@ class JetCalibrator : public xAH::Algorithm bool m_sort; //Apply jet cleaning to parent jet bool m_cleanParent; + bool m_applyFatJetPreSel; // systematics bool m_runSysts; From 1e63c37f32d21a86f16f02b85626a7634955d0d8 Mon Sep 17 00:00:00 2001 From: John Alison Date: Mon, 14 Dec 2015 08:28:54 -0600 Subject: [PATCH 2/2] Update ttbar allhad x-section --- data/metadata/ttbar_allhad_crosssections_13TeV.txt | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/data/metadata/ttbar_allhad_crosssections_13TeV.txt b/data/metadata/ttbar_allhad_crosssections_13TeV.txt index 0381685104..3d6e0460d3 100644 --- a/data/metadata/ttbar_allhad_crosssections_13TeV.txt +++ b/data/metadata/ttbar_allhad_crosssections_13TeV.txt @@ -9,5 +9,5 @@ id/I:name/C:xsec/F:kfac/F:eff/F:relunc/F 303725 PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_allhad_mtt_4 696.24 1. 0.00043823 1. 303726 PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_allhad_mtt_5 696.24 1. 0.00025914 1. -410007 PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_allhad 696.21 1. 0.45618 1. +410007 PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_allhad 831.76 1. 0.45618 1.