/** @file run_digi.C
 ** @author Volker Friese <v.friese@gsi.de>
 ** @since 8 June 2018
 **/


// Includes needed for IDE
#if !defined(__CLING__)
#include "FairSystemInfo.h"
#include "CbmDigitization.h"
#endif


/** @brief run_digi
 ** @author Volker Friese <v.friese@gsi.de>
 ** @since  8 June 2018
 ** @param nEvents    Number of events to process
 ** @param dataSet    Name of data set (input and output)
 ** @param eventRate  Interaction rate [1/s]
 ** @param tsLength   Length of time slice [ns]
 **
 ** Macro for detector response simulation with CBM, for regular tests.
 **
 ** The detector response produces a raw data file from the transport
 ** data, which serves as input for reconstruction. Raw data will
 ** be delivered in time-slice format (one tree entry per time-slice).
 **
 ** For the file names, the following conventions are applied:
 ** Input file:     [dataSet].tra.root
 ** Parameter file: [dataSet].par.root
 ** Output file:    [dataSet].raw.root
 ** If different names are wanted, they have to be specified
 ** explicitly in the macro.
 **
 ** The arguments to the macro call represent the minimal
 ** information to be passed to the digitization run. In addition, 
 ** the mode (time-based as default, or event-by-event) can be specified.
 ** By default, already existing output files will not be overwritten
 ** to prevent accidental data loss.
 **
 ** If the time-slice length is negative, all data will be in one
 ** time-slice.
 **
 ** For further options to modify the run settings, consult
 ** the documentation of the CbmDigitization class.
 **/
void run_digi(
	Int_t nEvents = 2,                // Number of events to process
	TString dataSet = "test",         // Dataset for file names
	Double_t eventRate = 1.e7,        // Interaction rate [1/s]
	Double_t timeSliceLength = 1.e4,  // Length of time-slice [ns]
	Bool_t eventMode = kFALSE         // Event-by-event mode
)
{

  // --- Logger settings ----------------------------------------------------
  FairLogger::GetLogger()->SetLogScreenLevel("INFO");
  FairLogger::GetLogger()->SetLogVerbosityLevel("LOW");
  // ------------------------------------------------------------------------


  // -----   Allow overwriting of output file   -----------------------------
  Bool_t overwrite = kTRUE;
  // ------------------------------------------------------------------------
 

  // -----   File names   ---------------------------------------------------
  TString inFile  = dataSet + ".tra.root";
  TString parFile = dataSet + ".par.root";
  TString outFile = dataSet + ".raw.root";
  TString monFile = dataSet + ".raw.moni.root";
  if ( eventMode ) {
    outFile = dataSet + ".event.raw.root";
    monFile = dataSet + ".event.raw.moni.root";
  }

   // -----   Timer   --------------------------------------------------------
  TStopwatch timer;
  timer.Start();
  // ------------------------------------------------------------------------
 
 
  // -----   Digitization run   ---------------------------------------------
  CbmDigitization run;
  
  run.AddInput(inFile, eventRate);
  run.SetOutputFile(outFile, overwrite);
  run.SetMonitorFile(monFile);
  run.SetParameterRootFile(parFile);
  run.SetTimeSliceLength(timeSliceLength);
  run.SetEventMode(eventMode);
  
  run.Run(nEvents);
  // ------------------------------------------------------------------------


  // -----   Finish   -------------------------------------------------------
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  std::cout << std::endl;
  std::cout << "Macro finished successfully." << std::endl;
  std::cout << "Output file is " << outFile << std::endl;
  std::cout << "Parameter file is " << parFile << std::endl;
  std::cout << "Real time " << rtime << " s, CPU time " << ctime
            << " s" << std::endl << std::endl;
  // ------------------------------------------------------------------------


  // -----   CTest resource monitoring   ------------------------------------
  FairSystemInfo sysInfo;
  Float_t maxMemory=sysInfo.GetMaxMemory();
  std::cout << "<DartMeasurement name=\"MaxMemory\" type=\"numeric/double\">";
  std::cout << maxMemory;
  std::cout << "</DartMeasurement>" << std::endl;
  std::cout << "<DartMeasurement name=\"WallTime\" type=\"numeric/double\">";
  std::cout << rtime;
  std::cout << "</DartMeasurement>" << std::endl;
  Float_t cpuUsage=ctime/rtime;
  std::cout << "<DartMeasurement name=\"CpuLoad\" type=\"numeric/double\">";
  std::cout << cpuUsage;
  std::cout << "</DartMeasurement>" << std::endl;
  // ------------------------------------------------------------------------


  // -----   Finish   -------------------------------------------------------
  std::cout << " Test passed" << std::endl;
  std::cout << " All ok " << std::endl;
  // ------------------------------------------------------------------------


} // End of macro