QA for digi events. Checks the digi time distribution within the event.

This simple QA in the task format checks the (STS) digi time distribution within each constructed event. This distribution should be peaked, demonstrating that real events have been triggered on and were built, and not just some noise is sampled, which would result in a flat distribution.

The corresponding histogram is registered to the THttpServer of FairRunOnline, but I do not yet know how to get it on screen. So, for the time being, it is also written to the output file. Example plots are attached below.

N.b.: I notice that filling the (single) histogram (TH1F) is rather costly: the tasks consumes like 70 ms per timeslice (about 6.5 millions digis / calls to TH1F::Fill()). The task also calculates the first and second moment of the distribution numerically; without histogram filling, the time consumation is only about 10 ms per timeslice.

macro/reco/reco_steer.C

@@ -24,6 +24,7 @@ using std::string;
  ** @param tsaFile    Name of input file (w/o extension .tsa)
  ** @param outFile    Name of output file (w/o extension .digi.root)
  ** @param numTs      Number of timeslices to process. If not specified, all available will be used.
+ ** @param port       Port of http server. If 0, server will not be activated.
  **
  ** Reconstruction from timeslice level, making use of the steering class CbmReco.
  ** Currently included stages:
@@ -36,7 +37,8 @@ using std::string;
  ** the extension .tsa by .digi.root
  **/
 
-void reco_steer(TString tsaFile = "", TString outFile = "", int32_t numTs = std::numeric_limits<int32_t>::max())
+void reco_steer(TString tsaFile = "", TString outFile = "", int32_t numTs = std::numeric_limits<int32_t>::max(),
+                uint32_t port = 8080)
 {
 
   // ========================================================================
@@ -98,7 +100,7 @@ void reco_steer(TString tsaFile = "", TString outFile = "", int32_t numTs = std:
   // -----   Run reconstruction   -------------------------------------------
   TStopwatch timer;
   timer.Start();
-  CbmReco run(inFile.Data(), outFile.Data(), numTs, config);
+  CbmReco run(inFile.Data(), outFile.Data(), numTs, config, port);
   run.Run();
   timer.Stop();
   // ------------------------------------------------------------------------