random seed: setter, read from json

Author: lubynets

PWGHF/D2H/Macros/HFInvMassFitter.cxx

@@ -1013,7 +1013,7 @@ double HFInvMassFitter::randomizeInitialParameter(const ParameterRanges& paramet
     ++nIter;
     if (nIter > MaximalNumberOfIterations) {
       printf("randomizeInitialFitParameter() - long while loop with lower = %f upper = %f initial = %f sigma = %f\n", parameterRanges.lower, parameterRanges.upper, parameterRanges.initial, sigma);
-      throw std::runtime_error("");
+      throw;
     }
   } while (result < parameterRanges.lower || result > parameterRanges.upper);
 

PWGHF/D2H/Macros/HFInvMassFitter.h

@@ -144,6 +144,7 @@ class HFInvMassFitter : public TNamed
   void drawResidual(TVirtualPad* c);
   void drawRatio(TVirtualPad* c);
   void drawReflection(TVirtualPad* c);
+  void setRandomSeed(int seed) { mRandomSeed = seed; }
 
  private:
   HFInvMassFitter(const HFInvMassFitter& source);

PWGHF/D2H/Macros/config_massfitter.json

@@ -158,5 +158,6 @@
         "3 for GausSec"
     ],
     "DrawBgPrefit": true,
-    "HighlightPeakRegion": true
+    "HighlightPeakRegion": true,
+    "RandomSeed": -1
 }

PWGHF/D2H/Macros/runMassFitter.C

@@ -154,6 +154,7 @@ void runMassFitter(const std::string& configFileName)
   const bool enableRefl = readJsonField<bool>(config, "EnableRefl", false);
   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 int nSliceVarBins = static_cast<int>(sliceVarMin.size());
   std::vector<double> sliceVarLimits(nSliceVarBins + 1);
@@ -381,6 +382,7 @@ void runMassFitter(const std::string& configFileName)
     } else {
       massFitter->setUseChi2Fit();
     }
+    massFitter->setRandomSeed(randomSeed);
 
     auto setFixedValue = [&iSliceVar](bool const& isFix, std::vector<double> const& fixManual, const TH1* histToFix, std::function<void(double)> setFunc, std::string const& var) -> void {
       if (isFix) {