CbmMcbm2019TimeWinEventBuilderAlgo.cxx

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *              GNU Lesser General Public Licence (LGPL) version 3,             *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
#include "CbmMcbm2019TimeWinEventBuilderAlgo.h"

/// CBM headers
#include "CbmEvent.h"
#include "TimesliceMetaData.h"
#include "CbmStsDigi.h"
#include "CbmMuchBeamTimeDigi.h"
#include "CbmTrdDigi.h"
#include "CbmTofDigi.h"
#include "CbmRichDigi.h"
#include "CbmPsdDigi.h"

#include "CbmDigiManager.h"

/// FAIRROOT headers
#include "FairLogger.h"
#include "FairRootManager.h"
#include "FairRunOnline.h"

/// FAIRSOFT headers (geant, boost, ...)
#include "TClonesArray.h"
#include "TH2.h"
#include "TH1.h"
#include "TCanvas.h"
#include "THttpServer.h"

/// C/C++ headers

// ---- Default constructor --------------------------------------------
CbmMcbm2019TimeWinEventBuilderAlgo::CbmMcbm2019TimeWinEventBuilderAlgo()
{
}

// ---- Destructor -----------------------------------------------------
CbmMcbm2019TimeWinEventBuilderAlgo::~CbmMcbm2019TimeWinEventBuilderAlgo()
{
}

// ---- Init -----------------------------------------------------------
Bool_t CbmMcbm2019TimeWinEventBuilderAlgo::InitAlgo()
{
   LOG(info) << "CbmMcbm2019TimeWinEventBuilderAlgo::InitAlgo => Starting sequence";

   // Get a handle from the IO manager
   FairRootManager* ioman = FairRootManager::Instance();

   // Get a pointer to the previous already existing data level
   fDigiMan = CbmDigiManager::Instance();
   fDigiMan->UseMuchBeamTimeDigi();
   fDigiMan->Init();

   // T0 is not included in DigiManager
   fT0DigiVec = ioman->InitObjectAs<std::vector<CbmTofDigi> const *>("T0Digi");
   if( ! fT0DigiVec )
   {
      LOG(info) << "No T0 digi input.";
      if( ECbmModuleId::kT0 == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fT0DigiVec )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kSts) )
   {
      LOG(info) << "No STS digi input.";
      if( ECbmModuleId::kSts == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kSts) )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kMuch) )
   {
      LOG(info) << "No MUCH digi input.";
      if( ECbmModuleId::kMuch == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kMuch) )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kTrd) )
   {
      LOG(info) << "No TRD digi input.";
      if( ECbmModuleId::kTrd == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kTrd) )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kTof) )
   {
      LOG(info) << "No TOF digi input.";
      if( ECbmModuleId::kTof == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kTof) )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kRich) )
   {
      LOG(info) << "No RICH digi input.";
      if( ECbmModuleId::kRich == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kRich) )

   if( ! fDigiMan->IsPresent(ECbmModuleId::kPsd) )
   {
      LOG(info) << "No PSD digi input.";
      if( ECbmModuleId::kPsd == fRefDet )
      {
         LOG(fatal) << "No digi input for reference detector, stopping there!";
      } // if( ECbmModuleId::kT0 == fRefDet )
   } // if( ! fDigiMan->IsPresent(ECbmModuleId::kPsd) )

   /// Access the TS metadata to know TS start tim
   fTimeSliceMetaDataArray = dynamic_cast<TClonesArray*>(ioman->GetObject("TimesliceMetaData"));
   if ( ! fTimeSliceMetaDataArray )
      LOG(fatal) << "No TS metadata input found";

   /// Store the time window for the reference detector
   switch( fRefDet )
   {
      case ECbmModuleId::kSts:
      {
         fdRefTimeWinBeg = fdStsTimeWinBeg;
         fdRefTimeWinEnd = fdStsTimeWinEnd;
         break;
      } // case ECbmModuleId::kSts:
      case ECbmModuleId::kMuch:
      {
         fdRefTimeWinBeg = fdMuchTimeWinBeg;
         fdRefTimeWinEnd = fdMuchTimeWinEnd;
         break;
      } // case ECbmModuleId::kMuch:
      case ECbmModuleId::kTrd:
      {
         fdRefTimeWinBeg = fdTrdTimeWinBeg;
         fdRefTimeWinEnd = fdTrdTimeWinEnd;
         break;
      } // case ECbmModuleId::kTrd:
      case ECbmModuleId::kTof:
      {
         fdRefTimeWinBeg = fdTofTimeWinBeg;
         fdRefTimeWinEnd = fdTofTimeWinEnd;
         break;
      } // case ECbmModuleId::kTof:
      case ECbmModuleId::kRich:
      {
         fdRefTimeWinBeg = fdRichTimeWinBeg;
         fdRefTimeWinEnd = fdRichTimeWinEnd;
         break;
      } // case ECbmModuleId::kRich:
      case ECbmModuleId::kPsd:
      {
         fdRefTimeWinBeg = fdPsdTimeWinBeg;
         fdRefTimeWinEnd = fdPsdTimeWinEnd;
         break;
      } // case ECbmModuleId::kPsd:
      case ECbmModuleId::kT0:
      {
         fdRefTimeWinBeg = fdT0TimeWinBeg;
         fdRefTimeWinEnd = fdT0TimeWinEnd;
         break;
      } // case ECbmModuleId::kT0:
      default:
      {
         LOG( fatal ) << "CbmMcbm2019TimeWinEventBuilderAlgo::Init => "
                      << "Trying to use unsupported detectore as reference: "
                      << fRefDet;
         break;
      } // default:
   } // switch( fRefDet )

   if( fbFillHistos )
   {
      CreateHistograms();
   } // if( fbFillHistos )

   LOG(info) << "CbmMcbm2019TimeWinEventBuilderAlgo::InitAlgo => Done";

   return kTRUE;
}

// ---- ProcessTs ------------------------------------------------------
void CbmMcbm2019TimeWinEventBuilderAlgo::ProcessTs()
{
   LOG_IF(info, fuNrTs%1000 == 0) <<"Begin of TS " << fuNrTs;

   InitTs();

   BuildEvents();