Monitoring for STS included in new unpacker scheme.

macro/run/run_unpack_tsa.C

@@ -26,6 +26,7 @@
 
 std::shared_ptr<CbmTrdUnpackMonitor> GetTrdMonitor(std::string treefilename);
 std::shared_ptr<CbmTrdSpadic> GetTrdSpadic(bool useAvgBaseline = false);
+std::shared_ptr<CbmStsUnpackMonitor> GetStsMonitor(std::string treefilename);
 
 void run_unpack_tsa(std::string infile = "test.tsa", UInt_t runid = 0, const char* setupName = "mcbm_beam_2021_03",
                     std::int32_t nevents = -1, std::string outpath = "")
@@ -119,6 +120,7 @@ void run_unpack_tsa(std::string infile = "test.tsa", UInt_t runid = 0, const cha
     stsconfig->SetDoWriteOutput();
     std::string parfilesbasepathSts = Form("%s/macro/beamtime/mcbm2021/", srcDir.Data());
     stsconfig->SetParFilesBasePath(parfilesbasepathSts);
+    stsconfig->SetMonitor(GetStsMonitor(outfilename));
     stsconfig->SetSystemTimeOffset(-2221);  // [ns] value to be updated
   }
   // -------------
@@ -301,6 +303,7 @@ std::shared_ptr<CbmTrdUnpackMonitor> GetTrdMonitor(std::string treefilename)
   return monitor;
 }
 
+
 /**
  * @brief Get the Trd Spadic
  * @return std::shared_ptr<CbmTrdSpadic>
@@ -313,3 +316,40 @@ std::shared_ptr<CbmTrdSpadic> GetTrdSpadic(bool useAvgBaseline)
 
   return spadic;
 }
+
+/**
+ * @brief Get the Sts Monitor. Extra function to keep default macro part more silent.
+ * @return std::shared_ptr<CbmStsUnpackMonitor>
+*/
+std::shared_ptr<CbmStsUnpackMonitor> GetStsMonitor(std::string treefilename)
+{
+  // -----   Output filename and path   -------------------------------------
+  std::string outpath  = "";
+  std::string filename = "";
+  auto filenamepos     = treefilename.find_last_of("/");
+  if (filenamepos != treefilename.npos) {
+    outpath  = treefilename.substr(0, filenamepos);
+    filename = treefilename.substr(filenamepos++);
+  }
+  if (outpath.empty()) outpath = gSystem->GetWorkingDirectory();
+  std::string mydir = "/qa";
+  outpath += mydir;
+
+  auto currentdir = gSystem->GetWorkingDirectory();
+
+  if (!gSystem->cd(outpath.data())) gSystem->MakeDirectory(outpath.data());
+  else
+    gSystem->cd(currentdir.data());
+
+  std::string outfilename = outpath + filename;
+  auto filetypepos        = outfilename.find(".digi.root");
+  if (filetypepos != outfilename.npos) outfilename.replace(filetypepos, 10, ".mon.sts.root");
+  else
+    outfilename += ".mon.sts.root";
+  // ------------------------------------------------------------------------
+
+  auto monitor = std::make_shared<CbmStsUnpackMonitor>();
+  monitor->SetHistoFileNameFull(outfilename);
+  //monitor->SetDebugMode(true);
+  return monitor;
+}

reco/detectors/sts/unpack/CbmStsUnpackAlgo.cxx

@@ -242,7 +242,6 @@ Bool_t CbmStsUnpackAlgo::initParSet(CbmMcbm2018StsPar* parset)
   LOG(debug) << "Unpacking data in bin sorter FW mode";
   initInternalStatus(parset);
 
-
   return kTRUE;
 }
 
@@ -696,6 +695,7 @@ bool CbmStsUnpackAlgo::unpack(const fles::Timeslice* ts, std::uint16_t icomp, UI
       const double dMsTime = (1e-9) * static_cast<double>(fMsStartTime);
       if (icomp < fMonitor->GetMaxNbFlibLinks()) {
         if (fdStartTimeMsSz < 0) fdStartTimeMsSz = dMsTime;
+        fMonitor->CreateMsComponentSizeHistos(icomp);
         fMonitor->FillMsSize(icomp, uSize);
         fMonitor->FillMsSizeTime(icomp, dMsTime - fdStartTimeMsSz, uSize);
       }
@@ -749,7 +749,7 @@ bool CbmStsUnpackAlgo::unpack(const fles::Timeslice* ts, std::uint16_t icomp, UI
 
   //Output debugging info
   /** @todo @experts this is printout debugging which depends on monitor settings. Not sure whether this is a good way. Please check and decide. (It was this way already before to be clear) */
-  if (fMonitor)
+  if (fMonitor) {
     if (fMonitor->GetDebugMode()) {
       if (18967040000 == fMsStartTime || 18968320000 == fMsStartTime) { LOG(debug) << sMessPatt; }
       if (static_cast<uint16_t>(fles::MicrosliceFlags::CrcValid) != msDescriptor.flags) {
@@ -764,8 +764,13 @@ bool CbmStsUnpackAlgo::unpack(const fles::Timeslice* ts, std::uint16_t icomp, UI
                    << (!bError && 400 != vNbMessType[1]);
       }
     }
+    for (auto itHit = fOutputVec.begin(); itHit != fOutputVec.end(); ++itHit) {
+      fMonitor->FillDigisTimeInRun(itHit->GetTime());
+    }
+    fMonitor->FillVectorSize(ts->index(), fOutputVec.size());
+    //fMonitor->DrawCanvases();
+  }
   return true;
 }
 
-
 ClassImp(CbmStsUnpackAlgo)

reco/detectors/sts/unpack/CbmStsUnpackAlgo.h

@@ -89,7 +89,7 @@ protected:
   /** @brief Finish function for this algorithm base clase */
   void finish()
   {
-    // if (fMonitor) fMonitor->Finish();
+    if (fMonitor) fMonitor->Finish();
     return;
   }
 
@@ -267,7 +267,6 @@ protected:
   /** @brief Vector of channel masks, [ NbFeb ][ NbCHanInFeb ], used only if fbUseChannelMask is true */
   std::vector<std::vector<bool>> fvvbMaskedChannels = {};  //!
 
-
   /** @brief  FEB type, [ NbDpb ][ NbCrobPerDpb ][ NbFebsPerCrob ], 0 = A, 1 = B, -1 if inactive */
   std::vector<std::vector<std::vector<int32_t>>> fviFebType = {};  //!
 

reco/detectors/sts/unpack/CbmStsUnpackConfig.cxx

@@ -54,7 +54,10 @@ void CbmStsUnpackConfig::InitUnpacker()
   // Now we have all information required to initialise the algorithm
   algo->Init();
 
-  if (fMonitor) algo->SetMonitor(fMonitor);
+  if (fMonitor) {
+    fMonitor->Init(static_cast<CbmMcbm2018StsPar*>(reqparvec->at(0).second.get()));
+    algo->SetMonitor(fMonitor);
+  }
 
   // Pass the algo to its member in the base class
   fAlgo = algo;

reco/detectors/sts/unpack/CbmStsUnpackMonitor.cxx

@@ -292,13 +292,13 @@ Bool_t CbmStsUnpackMonitor::CreateMsComponentSizeHistos(UInt_t component)
     TString sMsSizeTitle = Form("Size of MS for nDPB of link %02u; Ms Size [bytes]", component);
     fvhMsSize[component] = new TH1F(sMsSizeName.Data(), sMsSizeTitle.Data(), 30000, 0., 30000.);
     fvhMsSize[component]->SetCanExtend(TH2::kAllAxes);
-
-    sMsSizeName              = Form("MsSizeTime_link_%02u", component);
-    sMsSizeTitle             = Form("Size of MS vs time for gDPB of link %02u; Time[s] ; Ms Size [bytes]", component);
+    AddHistoToVector(fvhMsSize[component], "perComponent");
+  }
+  if (nullptr == fvhMsSizeTime[component]) {
+    TString sMsSizeName      = Form("MsSizeTime_link_%02u", component);
+    TString sMsSizeTitle     = Form("Size of MS vs time for gDPB of link %02u; Time[s] ; Ms Size [bytes]", component);
     fvhMsSizeTime[component] = new TProfile(sMsSizeName.Data(), sMsSizeTitle.Data(), 15000, 0., 300.);
     fvhMsSizeTime[component]->SetCanExtend(TH2::kAllAxes);
-
-    AddHistoToVector(fvhMsSize[component], "perComponent");
     AddHistoToVector(fvhMsSizeTime[component], "perComponent");
   }
   return kTRUE;
@@ -306,10 +306,8 @@ Bool_t CbmStsUnpackMonitor::CreateMsComponentSizeHistos(UInt_t component)
 
 Bool_t CbmStsUnpackMonitor::ResetMsComponentSizeHistos(UInt_t component)
 {
-  if (nullptr == fvhMsSize[component]) {
-    fvhMsSize[component]->Reset();
-    fvhMsSizeTime[component]->Reset();
-  }
+  if (nullptr != fvhMsSize[component]) { fvhMsSize[component]->Reset(); }
+  if (nullptr != fvhMsSizeTime[component]) { fvhMsSizeTime[component]->Reset(); }
   return kTRUE;
 }
 
@@ -678,5 +676,58 @@ Bool_t CbmStsUnpackMonitor::Init(CbmMcbm2018StsPar* parset)
   return kTRUE;
 }
 
+// ---- Finish ----
+void CbmStsUnpackMonitor::Finish()
+{
+  DrawCanvases();
+
+  /// Obtain vector of pointers on each histo (+ optionally desired folder)
+  std::vector<std::pair<TNamed*, std::string>> vHistos = GetHistoVector();
+
+  /// Obtain vector of pointers on each canvas (+ optionally desired folder)
+  std::vector<std::pair<CbmQaCanvas*, std::string>> vCanvases = GetCanvasVector();
+
+  /// Save old global file and folder pointer to avoid messing with FairRoot
+  TFile* oldFile     = gFile;
+  TDirectory* oldDir = gDirectory;
+  TFile* histoFile   = nullptr;
+
+  // open separate histo file in recreate mode
+  histoFile = new TFile(fHistoFileName, "RECREATE");
+  histoFile->cd();
+
+  /// Write histos to output file
+  for (UInt_t uHisto = 0; uHisto < vHistos.size(); ++uHisto) {
+    /// Make sure we end up in chosen folder
+    TString sFolder = vHistos[uHisto].second.data();
+    if (nullptr == gDirectory->Get(sFolder)) gDirectory->mkdir(sFolder);
+    gDirectory->cd(sFolder);
+
+    /// Write plot
+    vHistos[uHisto].first->Write();
+
+    histoFile->cd();
+  }
+
+  /// Write canvases to output file
+  for (UInt_t uCanvas = 0; uCanvas < vCanvases.size(); ++uCanvas) {
+    /// Make sure we end up in chosen folder
+    TString sFolder = vCanvases[uCanvas].second.data();
+    if (nullptr == gDirectory->Get(sFolder)) gDirectory->mkdir(sFolder);
+    gDirectory->cd(sFolder);
+
+    /// Write canvas
+    vCanvases[uCanvas].first->Write();
+
+    histoFile->cd();
+  }
+
+  /// Restore old global file and folder pointer to avoid messing with FairRoot
+  gFile      = oldFile;
+  gDirectory = oldDir;
+
+  histoFile->Close();
+}
+
 
 ClassImp(CbmStsUnpackMonitor)

reco/detectors/sts/unpack/CbmStsUnpackMonitor.h

@@ -41,6 +41,12 @@ public:
   Bool_t CreateMsComponentSizeHistos(UInt_t component);
   Bool_t ResetMsComponentSizeHistos(UInt_t component);
 
+  /** @brief Write all histograms and canvases to file */
+  void Finish();
+
+  void SetHistoFileName(TString nameIn) { fHistoFileName = "data/" + nameIn + ".root"; }
+  void SetHistoFileNameFull(TString nameIn) { fHistoFileName = nameIn; }
+
   void DrawCanvases();
 
   void AddHistoToVector(TNamed* pointer, std::string sFolder = "")
@@ -209,6 +215,8 @@ public:
   void SetDebugMode(bool value) { fDebugMode = value; }
 
 private:
+  TString fHistoFileName = "data/HistosUnpackerSts.root";
+
   /// Rate evolution histos
   //Bool_t fbLongHistoEnable;
   UInt_t fuLongHistoNbSeconds  = 3600;