bring to order N sigma [sgn, sidebands]

Author: lubynets

PWGHF/D2H/Macros/HFInvMassFitter.cxx

@@ -583,9 +583,9 @@ void HFInvMassFitter::drawFit(TVirtualPad* pad, const std::vector<std::string>&
   textFitMetrics->AddText(Form("S = %.0f #pm %.0f ", mRawYield, mRawYieldErr));
   textFitMetrics->AddText(Form("S_{count} = %.0f #pm %.0f ", mRawYieldCounted, mRawYieldCountedErr));
   if (mTypeOfBkgPdf != NoBkg) {
-    textFitMetrics->AddText(Form("B (%.1f#sigma) = %.0f #pm %.0f", mNSigmaForSidebands, mBkgYield, mBkgYieldErr));
-    textFitMetrics->AddText(Form("S/B (%.1f#sigma) = %.4g ", mNSigmaForSidebands, mRawYield / mBkgYield));
-    textFitMetrics->AddText(Form("Significance (%.1f#sigma) = %.1f #pm %.1f ", mNSigmaForSidebands, mSignificance, mSignificanceErr));
+    textFitMetrics->AddText(Form("B (%.1f#sigma) = %.0f #pm %.0f", mNSigmaForSgn, mBkgYield, mBkgYieldErr));
+    textFitMetrics->AddText(Form("S/B (%.1f#sigma) = %.4g ", mNSigmaForSgn, mRawYield / mBkgYield));
+    textFitMetrics->AddText(Form("Significance (%.1f#sigma) = %.1f #pm %.1f ", mNSigmaForSgn, mSignificance, mSignificanceErr));
     textFitMetrics->AddText(Form("#chi^{2} / ndf  =  %.3f", mChiSquareOverNdfTotal));
   }
   if (mReflPdf != nullptr) {

PWGHF/D2H/Macros/HFInvMassFitter.h

@@ -108,6 +108,7 @@ class HFInvMassFitter : public TNamed
   void setFixRatioToGausSigma(double sigmaFrac);
   void setFixSignalYield(double yield) { mFixedRawYield = yield; }
   void setNumberOfSigmaForSidebands(double numberOfSigma) { mNSigmaForSidebands = numberOfSigma; }
+  void setNumberOfSigmaForSignal(double numberOfSigma) { mNSigmaForSgn = numberOfSigma; }
   void setFixDscbAlphaL(double alphaL);
   void setFixDscbAlphaR(double alphaR);
   void setFixDscbNL(double nL);

PWGHF/D2H/Macros/runMassFitter.C

@@ -183,6 +183,8 @@ void runMassFitter(const std::string& configFileName)
   const bool drawBgPrefit = readJsonField<bool>(config, "DrawBgPrefit", true);
   const bool highlightPeakRegion = readJsonField<bool>(config, "HighlightPeakRegion", true);
   const int randomSeed = readJsonField<int>(config, "RandomSeed", -1);
+  const double nSigmaForSideband = readJsonField<double>(config, "NSigmaForSideband", 3.);
+  const double nSigmaForSignal = readJsonField<double>(config, "NSigmaForSignal", 3.);
 
   readJsonVector(dscbAlphaLInitial, config, "DscbAlphaLInitial");
   readJsonVector(dscbAlphaLLower, config, "DscbAlphaLLower");
@@ -475,6 +477,8 @@ void runMassFitter(const std::string& configFileName)
 
     HFInvMassFitter* massFitter = new HFInvMassFitter(hMass[iSliceVar], massMin[iSliceVar], massMax[iSliceVar], bkgFunc[iSliceVar], sgnFunc[iSliceVar], randomSeed);
     massFitter->setDrawBgPrefit(drawBgPrefit);
+    massFitter->setNumberOfSigmaForSidebands(nSigmaForSideband);
+    massFitter->setNumberOfSigmaForSignal(nSigmaForSignal);
     massFitter->setHighlightPeakRegion(highlightPeakRegion);
     massFitter->setInitialGaussianMean(massPDG);
     massFitter->setParticlePdgMass(massPDG);