[PWGLF,PWGMM] Add filling scheme information to lumi stability task (#14480)

Author: fchinu

PWGMM/Lumi/Tasks/lumiStabilityPP.cxx

@@ -31,6 +31,7 @@
 
 #include <limits>
 #include <map>
+#include <memory>
 #include <string>
 #include <vector>
 
@@ -114,6 +115,7 @@ struct LumiStabilityPP {
   parameters::GRPLHCIFData* mLHCIFdata = nullptr;
   int runNumber{-1};
   ctpRateFetcher mRateFetcher;
+  std::string injectionScheme;
 
   HistogramRegistry registry{"registry"};
 
@@ -124,6 +126,7 @@ struct LumiStabilityPP {
   std::map<int, std::shared_ptr<TH1>> histTfPerMin;
   std::map<int, std::shared_ptr<TH1>> histBcHasFT0;
   std::map<int, std::shared_ptr<TH1>> histBcHasFDD;
+  std::map<int, std::shared_ptr<TH1>> histFillingScheme;
 
   static constexpr std::string_view NBCsVsTimeHistNames[NTriggerAliases][NBCCategories] =
     {{"AllBCs/BC_A/nBCsVsTime", "AllBCs/BC_B/nBCsVsTime", "AllBCs/BC_C/nBCsVsTime", "AllBCs/BC_E/nBCsVsTime", "AllBCs/BC_L/nBCsVsTime", "AllBCs/BC_SL/nBCsVsTime"},
@@ -137,7 +140,7 @@ struct LumiStabilityPP {
      {"FT0CE/BC_A/nBCsVsBCID", "FT0CE/BC_B/nBCsVsBCID", "FT0CE/BC_C/nBCsVsBCID", "FT0CE/BC_E/nBCsVsBCID", "FT0CE/BC_L/nBCsVsBCID", "FT0CE/BC_SL/nBCsVsBCID"},
      {"FDD/BC_A/nBCsVsBCID", "FDD/BC_B/nBCsVsBCID", "FDD/BC_C/nBCsVsBCID", "FDD/BC_E/nBCsVsBCID", "FDD/BC_L/nBCsVsBCID", "FDD/BC_SL/nBCsVsBCID"}};
 
-  const AxisSpec timeAxis{1440, 0., 1440., "#bf{t-t_{SOF} (min)}"}, bcIDAxis{nBCsPerOrbit, -0.5, static_cast<float>(nBCsPerOrbit) - 0.5, "#bf{BC ID in orbit}"};
+  const AxisSpec timeAxis{2880, 0., 2880., "#bf{t-t_{SOF} (min)}"}, bcIDAxis{nBCsPerOrbit, -0.5, static_cast<float>(nBCsPerOrbit) - 0.5, "#bf{BC ID in orbit}"};
 
   int64_t bcSOR;
   int nBCsPerTF;
@@ -154,6 +157,7 @@ struct LumiStabilityPP {
     histNBcsVsTime[runNumber] = registry.add<TH1>(Form("%d/FT0Vtx_EvSel/nBCsVsTime", runNumber), "Time of TVX triggered BCs since the start of fill;;#bf{#it{N}_{BC}}", HistType::kTH1D, {timeAxis});
     histNBcsVsBcId[runNumber] = registry.add<TH1>(Form("%d/nBCsVsBCID", runNumber), "Time of TVX triggered BCs since the start of fill;#bf{t-t_{SOF} (min)};#bf{#it{N}_{BC}}", HistType::kTH1D, {bcIDAxis});
     histTfPerMin[runNumber] = registry.add<TH1>(Form("%d/TFsPerMinute", runNumber), "TFs seen in this minute (to account for failed jobs);#bf{t-t_{SOF} (min)};#bf{#it{N}_{TFs}}", HistType::kTH1D, {timeAxis});
+    histFillingScheme[runNumber] = registry.add<TH1>(Form("%d/FillingScheme", runNumber), "Filling Scheme;Filling Scheme;", HistType::kTH1D, {{1, 0, 1}});
 
     histBcHasFT0[runNumber] = registry.add<TH2>(Form("%d/FITQA/BCHasFT0", runNumber), "Does the BC have FT0?;BC has FT0;TVX triggered according to CTP;#bf{#it{N}_{BC}}", HistType::kTH2D, {{2, -0.5, 1.5}, {2, -0.5, 1.5}});
     histBcHasFT0[runNumber]->GetYaxis()->SetBinLabel(1, "No CTP trigger");
@@ -195,6 +199,8 @@ struct LumiStabilityPP {
     LOG(info) << "LHCIF data fetched for run " << runNumber << " and timestamp " << timeStamp;
     createHistograms();
 
+    histFillingScheme[runNumber]->Fill(mLHCIFdata->getInjectionScheme().c_str(), 0);
+
     beamPatternA = mLHCIFdata->getBunchFilling().getBeamPattern(0);
     beamPatternC = mLHCIFdata->getBunchFilling().getBeamPattern(1);
     bcPatternA = beamPatternA & ~beamPatternC;
@@ -287,7 +293,7 @@ struct LumiStabilityPP {
 
       if (thisTFid != currentTFid) {
         currentTFid = thisTFid;
-        histNBcsVsTime[runNumber]->Fill(timeSinceSOF);
+        histTfPerMin[runNumber]->Fill(timeSinceSOF);
       }
 
       std::bitset<64> ctpInputMask(bc.inputMask());