Merge pull request #27208 from Sam-Harper/ZeroPFClusEUncert_1060

zeroing PF clustering uncertainty : 106X
cms-sw · Jun 17, 2019 · d032c93 · d032c93
2 parents 6af6c5f + 91ba701
commit d032c93
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 4 deletions.
diff --git a/RecoParticleFlow/PFClusterProducer/interface/PFClusterEMEnergyCorrector.h b/RecoParticleFlow/PFClusterProducer/interface/PFClusterEMEnergyCorrector.h
@@ -53,6 +53,7 @@ class PFClusterEMEnergyCorrector {
   bool srfAwareCorrection_;
   bool applyCrackCorrections_;
   bool applyMVACorrections_;
+  bool setEnergyUncertainty_;
 
   bool autoDetectBunchSpacing_;
   int bunchSpacingManual_;

diff --git a/RecoParticleFlow/PFClusterProducer/plugins/CorrectedECALPFClusterProducer.cc b/RecoParticleFlow/PFClusterProducer/plugins/CorrectedECALPFClusterProducer.cc
@@ -132,6 +132,7 @@ void CorrectedECALPFClusterProducer::fillDescriptions(edm::ConfigurationDescript
     psd0.add<bool>("applyCrackCorrections",false);
     psd0.add<bool>("applyMVACorrections",false);
     psd0.add<bool>("srfAwareCorrection",false);    
+    psd0.add<bool>("setEnergyUncertainty",false);    
     psd0.add<bool>("autoDetectBunchSpacing",true);
     psd0.add<int>("bunchSpacing",25);
     psd0.add<double>("maxPtForMVAEvaluation",-99.);

diff --git a/RecoParticleFlow/PFClusterProducer/src/PFClusterEMEnergyCorrector.cc b/RecoParticleFlow/PFClusterProducer/src/PFClusterEMEnergyCorrector.cc
@@ -19,7 +19,7 @@ PFClusterEMEnergyCorrector::PFClusterEMEnergyCorrector(const edm::ParameterSet&
    applyCrackCorrections_ = conf.getParameter<bool>("applyCrackCorrections");
    applyMVACorrections_ = conf.getParameter<bool>("applyMVACorrections");
    srfAwareCorrection_ = conf.getParameter<bool>("srfAwareCorrection");
-
+   setEnergyUncertainty_ = conf.getParameter<bool>("setEnergyUncertainty");
    maxPtForMVAEvaluation_ = conf.getParameter<double>("maxPtForMVAEvaluation");
 
    if (applyMVACorrections_) {
@@ -250,8 +250,9 @@ void PFClusterEMEnergyCorrector::correctEnergies(const edm::Event &evt,
       double sigmacor = sigma*ecor;
 
       cluster.setCorrectedEnergy(ecor);
-      cluster.setCorrectedEnergyUncertainty(sigmacor);
-
+      if(setEnergyUncertainty_) cluster.setCorrectedEnergyUncertainty(sigmacor);
+      else cluster.setCorrectedEnergyUncertainty(0.);
+
     }
     return;
   }
@@ -414,7 +415,9 @@ void PFClusterEMEnergyCorrector::correctEnergies(const edm::Event &evt,
 					   << exp(mean) << " " << ecor;
 
     cluster.setCorrectedEnergy(ecor);
-    cluster.setCorrectedEnergyUncertainty(sigmacor);
+    if(setEnergyUncertainty_) cluster.setCorrectedEnergyUncertainty(sigmacor);
+    else cluster.setCorrectedEnergyUncertainty(0.);
+
   }
 
 }