Fix a bug in stochastic JER smearing

Resolution in pt is relative, so should not divide by pt.
andrey-popov · Aug 3, 2016 · 7967a1d · 7967a1d
1 parent da86372
commit 7967a1d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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;
         }