From 7967a1d98fccb21fed20c85644dc7081bb55f2d8 Mon Sep 17 00:00:00 2001 From: Andrey Popov Date: Wed, 3 Aug 2016 17:08:25 +0200 Subject: [PATCH] Fix a bug in stochastic JER smearing Resolution in pt is relative, so should not divide by pt. --- modules/extensions/src/JetCorrectorService.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/modules/extensions/src/JetCorrectorService.cpp b/modules/extensions/src/JetCorrectorService.cpp index b641bf6..57286bc 100644 --- a/modules/extensions/src/JetCorrectorService.cpp +++ b/modules/extensions/src/JetCorrectorService.cpp @@ -153,7 +153,7 @@ double JetCorrectorService::Eval(Jet const &jet, double rho, SystType syst /*= S jerProvider->getResolution({{JME::Binning::JetPt, jecCorrPt}, {JME::Binning::JetEta, jet.Eta()}, {JME::Binning::Rho, rho}}); double const jerFactor = 1. + rGen->Gaus(0., ptResolution) * - std::sqrt(std::max(std::pow(jerSF, 2) - 1., 0.)) / jecCorrPt; + std::sqrt(std::max(std::pow(jerSF, 2) - 1., 0.)); corrFactor *= jerFactor; }