beamtime/cosy-nov11/go4/ROOTFORMAT/TRootConverterProc.cxx (r4864/r3171)

Go to the documentation of this file.

/* Generated by Together */

#include "TRootConverterProc.h"
#include "TCBMBeamtimeEvent.h"
#include "TGo4Version.h"
#if __GO4BUILDVERSION__ > 40502
#include "go4iostream.h"
#else
#include "Riostream.h"
#endif
#include "TGo4UserException.h"
#include "TGo4MbsEvent.h"

#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TProfile.h"

#include "TGo4Analysis.h"

#include <algorithm>

//***********************************************************
TRootConverterProc::TRootConverterProc(const char* name) :
TCBMBeamtimeProc(name), 
fRpcInputEvent(0),
fCrateInputEvent(0),
fTriglogInputEvent(0),
fEpicsInputEvent(0),
fOutputEvent(0),
fBeamEvent(0),
cleanTree(0)
{
   TGo4Log::Info("**** TRootConverterProc: Create instance %s");
   
   fRocPar = (TRocParam*) GetParameter("RocPar");
   cout << "ROCPAR = " << fRocPar << endl;

   fRpcPar = (TRPCParam*) GetParameter("RPCPar");
   cout << "RPCPAR = " << fRpcPar << endl;

   fRootConvPar = (TRootConverterParam*) MakeParameter("RootConvPar", "TRootConverterParam");
   cout << "RootConvPar = " << fRootConvPar << endl;

   if (fRocPar->numRocs > 0) fNumRocs = (UInt_t) fRocPar->numRocs;

   // ROC general variables
   uNbEvents = 0;

   if(1 == fRootConvPar->uHodoDebugHisto)
   {
      TString path = GetName();
      path += "/";
      path += GetName();

      // Booking histograms
      hNbClustersPerEventHodoscope = MakeTH1('I', path+"_NbClustersPerEventHodoscope","Nb hosdoscope clusters per event",
            16,-0.5,15.5);

      hHitXvsMilliSecondsInSpill = MakeTH2('I', path+"_HitPosXVsTimeMilSecInSpill",
            "Hodoscope hit position versus time in Spill; Time in spill [ms]; X [mm]; Counts [mm]",
            8000,-50, 79950, 32,-14.5,15.5 );
      hHitYvsMilliSecondsInSpill = MakeTH2('I', path+"_HitPosYVsTimeMilSecInSpill",
            "Hodoscope hit position versus time in Spill; Time in spill [ms]; Y [mm]; Counts [mm]",
            8000,-50, 79950, 32,-14.5,15.5 );
      // No MakeTH3 function in go4, so do same operations by hand
      // => Need to check if loaded from autosave and create by ourselves
      path = GetName();
      path += "/";
      TString sName = GetName();
      hHit2DvsSeconds        = (TH3*)GetHistogram( path+ sName +"_HitPosVsTimeSec" );
      if( 0 == hHit2DvsSeconds )
      {
         hHit2DvsSeconds        = new TH3I( sName +"_HitPosVsTimeSec",
                                 "Hodoscope hit position versus time; Time [s]; X [mm]; Y [mm]",
                                 1000,-0.5, 1999.5, 32,-14.5,15.5, 32,-14.5,15.5 );
         AddHistogram( hHit2DvsSeconds, path );
      }
      hHit2DvsSecondsInSpill        = (TH3*)GetHistogram( path + sName +"_HitPosVsTimeSecInSpill" );
      if( 0 == hHit2DvsSecondsInSpill )
      {
         hHit2DvsSecondsInSpill = new TH3I( sName +"_HitPosVsTimeSecInSpill",
                                 "Hodoscope hit position versus time in Spill; Time in spill [s]; X [mm]; Y [mm]",
                                 80,-0.5, 79.5, 32,-14.5,15.5, 32,-14.5,15.5 );
         AddHistogram( hHit2DvsSecondsInSpill, path );
      }
      hHit2DvsMilliSecondsInSpill        = (TH3*)GetHistogram( path + sName +"_HitPosVsTimeMilSecInSpill" );
      if( 0 == hHit2DvsMilliSecondsInSpill )
      {
         hHit2DvsMilliSecondsInSpill = new TH3I( sName +"_HitPosVsTimeMilSecInSpill",
                                 "Hodoscope hit position versus time in Spill; Time in spill [ms]; X [mm]; Y [mm]",
                                 8000,-50, 79950, 32,-14.5,15.5, 32,-14.5,15.5 );
         AddHistogram( hHit2DvsMilliSecondsInSpill, path );
      }
      // For hodoscope hits map
      uFirstMbsTime = 0;
      uSpillStartMbsTime = 0;
      uLastMbsTime = 0;

      uFirstMbsTimeMs = 0;
      uSpillStartMbsTimeMs = 0;
   } // if(1 == fRootConvPar->uHodoDebugHisto)

   TString path = GetName();
   path += "/";
   path += GetName();
   hTimeEpicsUpdates = MakeTH1('I', path+"_TimeEpicsUpdate","Timing of epics updates; Time [s]",
         7200, 0, 7200);
   uEpicsTimeOrigin = 0;
      
   //printf("RPC Histograms created \n");
   fflush ( stdout);
}

TRootConverterProc::~TRootConverterProc()
{
   if( (1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay)
         && 0 != cleanTree )
      {
          if( 1 == fRootConvPar->uRootDisplay )
             TGo4Analysis::Instance()->RemoveTree(cleanTree);
          if( 1 == fRootConvPar->uCleanOuputEnable )
          {
             fRootOut->Write("");
             fRootOut->Close();
          } // if( 1 == fRootConvPar->uCleanOuputEnable )

          if( 0 == fRootConvPar->uCleanOuputEnable )
             delete cleanTree;

          if( 1 == fRootConvPar->uAddGet4DataInTree )
          {
             for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
             {
                (*vuHitNumber)[uChannelTemp]  = 0;
                (*vuTime)[uChannelTemp].clear();
                (*vuEpochTime)[uChannelTemp].clear();
                (*viTot)[uChannelTemp].clear();
             }

             delete vuHitNumber;
             delete vuTime;
             delete vuEpochTime;
             delete viTot;
          } // if( 1 == fRootConvPar->uAddGet4DataInTree )

          if( 1 == fRootConvPar->uSlowControlInTree )
          {
             vdPositiveHvI.clear();
             vdPositiveHvU.clear();
             vdNegativeHvI.clear();
             vdNegativeHvU.clear();
          }

//          if( 1 == fRootConvPar->uCleanOuputEnable )
//             delete fRootOut;
      }

   cout << "**** TRootConverterProc: Delete instance " << endl;
}
//***********************************************************

void TRootConverterProc::InitEvent(TGo4EventElement* outevnt)
{
   // first assign input event:
   // since input event object is never discarded within processor lifetime,
   // we just search for subevent by name once to speed up processing
   if(fRpcInputEvent==0 || (0==fHodoInputEvent && 1 == fRootConvPar->uHodoInputInTree) )
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent());
      //TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("CosyNov11DetectorEvent"));
      if(btevent)
      {
         fRpcInputEvent=dynamic_cast<TRPCEvent*>(btevent->GetSubEvent("RPC"));
         fBeamEvent=dynamic_cast<TBeamMonitorEvent*>(btevent->GetSubEvent("BEAM"));
         if( 1 == fRootConvPar->uHodoInputInTree )
            fHodoInputEvent=dynamic_cast<TFiberHodEvent*>(btevent->GetSubEvent("Hodo1"));
      }

      if(fRpcInputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event has no TRPCEvent!!! STOP GO4");
      }
      if(0==fHodoInputEvent && 1 == fRootConvPar->uHodoInputInTree) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event has no TFiberHodEvent!!! STOP GO4");
      }
   }

   // Optional: if triggered mode and tree, look for the VME crate data
   // since input event object is never discarded within processor lifetime,
   // we just search for subevent by name once to speed up processing
   if( kFALSE == fRocPar->noTrigger && ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Detector"));
      if(fCrateInputEvent==0)
      {
         if(btevent)
         {
            fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(btevent->GetSubEvent("MBSCRATE"));
         }

         if(fCrateInputEvent==0) {
            GO4_STOP_ANALYSIS_MESSAGE("**** TRPCProc: Fatal error: input event is/has not a TMbsCrateEvent!!! STOP GO4");
         }
      } // if(fCrateInputEvent==0)
      if( 1 == fRootConvPar->uTriglogInputInTree )
      {
         if(fTriglogInputEvent==0)
         {
            //TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("CosyNov11UnpackEvent"));
            if(btevent)
            {
               fTriglogInputEvent=dynamic_cast<TTriglogEvent*>(btevent->GetSubEvent("TRIGLOG"));
            }

            if(fTriglogInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRPCProc: Fatal error: input event is/has not a TTriglogEvent!!! STOP GO4");
            }
         } // if(fTriglogInputEvent==0)
      } // if( 1 == fRootConvPar->uTriglogInputInTree )
      if( 1 == fRootConvPar->uSlowControlInTree )
      {
         if(fEpicsInputEvent==0)
         {
            if(btevent)
            {
               fEpicsInputEvent=dynamic_cast<TEpicsEvent*>(btevent->GetSubEvent("EPICS"));
            }

            if(fEpicsInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRPCProc: Fatal error: input event is/has not a TEpicsEvent!!! STOP GO4");
            }
         } // if(fEpicsInputEvent==0)
      } // if( 1 == fRootConvPar->uSlowControlInTree )
   } // if( kFALSE == fRocPar->noTrigger && ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
      else if( 1 == fRootConvPar->uSlowControlInTree )
      {
         TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Detector"));
         if(fEpicsInputEvent==0)
         {
            if(btevent)
            {
               fEpicsInputEvent=dynamic_cast<TEpicsEvent*>(btevent->GetSubEvent("EPICS"));
            }

            if(fEpicsInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRPCProc: Fatal error: input event is/has not a TEpicsEvent!!! STOP GO4");
            }
         } // if(fEpicsInputEvent==0)
      } //  else if( 1 == fRootConvPar->uSlowControlInTree )

   // then assign output event
   // since output event object is never discarded within processor lifetime,
   // we just search for subevent by name once to speed up processing
   if(fOutputEvent==0)
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(outevnt);
      if(btevent)
      {
         fOutputEvent=dynamic_cast<TRootConverterEvent*>(btevent->GetSubEvent("ROOTCONV"));
      }
      else
      {

            fOutputEvent= dynamic_cast<TRootConverterEvent*>(outevnt);
      }
      if(fOutputEvent==0) {
      GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: output event is not a TRootConverterEvent!!! STOP GO4");
      }

   }

   // Prepare root tree and tree output file
   if( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay)
   {
      if( 0 == cleanTree)
      {
         TGo4Log::Info("**** TRootConverter: Create tree ");

         if( 1 == fRootConvPar->uCleanOuputEnable )
         {
            oldDir = gDirectory;

            sOutputFilename = fRootConvPar->sOutputFilename;
            fRootOut= new TFile(sOutputFilename, "RECREATE", "TTree file from GET4 and CAEN v1290", 9);
            sOutputFilename += ":/";
            if (fRootOut->IsZombie()) {
               cout << "Error opening file for TTree output in TRootConverter" << endl;
            }

            gDirectory->Cd(sOutputFilename);
         }

         if( 0 == TGo4Analysis::Instance()->GetTree( "CosyRpc" ) )
         {
            cleanTree = new TTree("CosyRpc", "RPC Caen and Get4");
            TGo4Log::Info("RPC Tree created in %s", gDirectory->GetPath());
         } // if( 0 == TGo4Analysis::Instance()->GetTree( "CosyRpc" ) )
            else
            {
               cleanTree = TGo4Analysis::Instance()->GetTree( "CosyRpc" );
               TGo4Log::Info("RPC Tree recovered from Go4");
            } // else of if( 0 == TGo4Analysis::Instance()->GetTree( "CosyRpc" ) )

         // To prevent having multiple trees when opening the file with Go4
         cleanTree->SetAutoFlush(0);

         // Data from CAEN v1290 & Lecroy QDC VME boards (triggered)
         if( kFALSE == fRocPar->noTrigger )
         {
            if( -1 != fRootConvPar->iPmtTdc && -1 != fRootConvPar->iPmtQdc )
            {
               cleanTree->Branch("PmtTdc_Leading",  (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).lead,
                                 Form("PmtTdc_Leading[%u]/I", T1290Data::NumChan )  );
               cleanTree->Branch("PmtTdc_Trailing", (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).trail,
                                 Form("PmtTdc_Trailing[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("PmtTdc_Leading_NbHit", (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).hit_lead,
                                    Form("PmtTdc_Leading_NbHit[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("PmtTdc_Trailing_NbHit", (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).hit_trail,
                                 Form("PmtTdc_Trailing_NbHit[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("PmtTdc_Leading_multi",  (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).lead_multi,
                                 Form("PmtTdc_Leading_multi[%u][%u]/I", T1290Data::NumChan, T1290Data::MaxHits )  );
               cleanTree->Branch("PmtTdc_Trailing_multi", (fCrateInputEvent->fMtdc[fRootConvPar->iPmtTdc]).trail_multi,
                                 Form("PmtTdc_Trailing_multi[%u][%u]/I", T1290Data::NumChan, T1290Data::MaxHits ) );
               cleanTree->Branch("PmtQdc", fCrateInputEvent->fData1182[fRootConvPar->iPmtQdc], Form("PmtQdc[%u]/s", NUM_1182_CH ) );
            } // if( -1 != fRootConvPar->iPmtTdc && -1 != fRootConvPar->iPmtQdc )
            if( -1 != fRootConvPar->iRpcTdcLeading && -1 != fRootConvPar->iRpcTdcTrailing )
            {
               cleanTree->Branch("CaenRpc_Leading",  (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcLeading]).lead,
                                 Form("CaenRpc_Leading[%u]/I", T1290Data::NumChan )  );
               cleanTree->Branch("CaenRpc_Trailing", (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcTrailing]).trail,
                                 Form("CaenRpc_Trailing[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("CaenRpc_Leading_NbHit", (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcLeading]).hit_lead,
                                    Form("CaenRpc_Leading_NbHit[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("CaenRpc_Trailing_NbHit", (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcTrailing]).hit_trail,
                                 Form("CaenRpc_Trailing_NbHit[%u]/I", T1290Data::NumChan ) );
               cleanTree->Branch("CaenRpc_Leading_multi",  (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcLeading]).lead_multi,
                                 Form("CaenRpc_Leading_multi[%u][%u]/I", T1290Data::NumChan, T1290Data::MaxHits )  );
               cleanTree->Branch("CaenRpc_Trailing_multi", (fCrateInputEvent->fMtdc[fRootConvPar->iRpcTdcTrailing]).trail_multi,
                                 Form("CaenRpc_Trailing_multi[%u][%u]/I", T1290Data::NumChan, T1290Data::MaxHits ) );
            } // if( -1 != fRootConvPar->iRpcTdcLeading && -1 != fRootConvPar->iRpcTdcTrailing )
            if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel1 )
            {
               cleanTree->Branch("CaenRef1_Leading",
                                 &((fCrateInputEvent->fMtdc[fRootConvPar->iReferenceTdc]).lead[fRootConvPar->iReferenceChannel1]),
                                 "CaenRef1_Leading/I" );
               cleanTree->Branch("CaenRef1_Trailing",
                                 &((fCrateInputEvent->fMtdc[fRootConvPar->iReferenceTdc]).trail[fRootConvPar->iReferenceChannel1]),
                                 "CaenRef1_Trailing/I" );
            } // if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel )
            if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel2 )
            {
               cleanTree->Branch("CaenRef2_Leading",
                                 &((fCrateInputEvent->fMtdc[fRootConvPar->iReferenceTdc]).lead[fRootConvPar->iReferenceChannel2]),
                                 "CaenRef2_Leading/I" );
               cleanTree->Branch("CaenRef2_Trailing",
                                 &((fCrateInputEvent->fMtdc[fRootConvPar->iReferenceTdc]).trail[fRootConvPar->iReferenceChannel2]),
                                 "CaenRef2_Trailing/I" );
            } // if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel )

            if( 1 == fRootConvPar->uTriglogInputInTree)
            {
               cleanTree->Branch("TriglogSyncNb",  &(fTriglogInputEvent->fVulomSyncNumber),
                                 "TriglogSyncNb/i"  );
               cleanTree->Branch("TriglogPattern",  &(fTriglogInputEvent->fVulomTriggerPattern),
                                 "TriglogPattern/i"  );
               cleanTree->Branch("TriglogTimeSec",  &(fTriglogInputEvent->fMbsTimeSecs),
                                 "TriglogTimeSec/i"  );
               cleanTree->Branch("TriglogTimeMilsec",  &(fTriglogInputEvent->fMbsTimeMillisec),
                                 "TriglogTimeMilsec/i"  );
               cleanTree->Branch("TriglogScalers1",  (fTriglogInputEvent->fVulomScaler[0]),
                                 Form("TriglogScalers1[%u]/i", N_SCALERS_CH )  );
               cleanTree->Branch("TriglogScalers2",  (fTriglogInputEvent->fVulomScaler[1]),
                                 Form("TriglogScalers2[%u]/i", N_SCALERS_CH )  );
               cleanTree->Branch("TriglogScalers3",  (fTriglogInputEvent->fVulomScaler[2]),
                                 Form("TriglogScalers3[%u]/i", N_SCALERS_CH )  );
            } // if( 1 == fRootConvPar->uTriglogInputInTree)
         } // if( kFALSE == fRocPar->noTrigger )
         if( 1 == fRootConvPar->uHodoInputInTree )
         {
            cleanTree->Branch("HodoscopeX",  &(hitHodoscope.X),"HodoscopeX/D");
            cleanTree->Branch("HodoscopeY",  &(hitHodoscope.Y),"HodoscopeY/D");
         } // if( 1 == fRootConvPar->uHodoInputInTree )

         if( 1 == fRootConvPar->uSlowControlInTree )
         {
            uSlowControlUpdateTime = 0;
            cleanTree->Branch( "TimeEpicsUpdate", &(uSlowControlUpdateTime),
                               "TimeEpicsUpdate/i"  );
            vdPositiveHvI.clear();
            vdPositiveHvU.clear();
            vdNegativeHvI.clear();
            vdNegativeHvU.clear();
            vdPositiveHvI.resize(fRootConvPar->uNbSlowControlChan);
            vdPositiveHvU.resize(fRootConvPar->uNbSlowControlChan);
            vdNegativeHvI.resize(fRootConvPar->uNbSlowControlChan);
            vdNegativeHvU.resize(fRootConvPar->uNbSlowControlChan);
            for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
            {
               cleanTree->Branch( fRootConvPar->sScOutputNamePosI[uChannelHV], &(vdPositiveHvI.at(uChannelHV)),
                                 fRootConvPar->sScOutputNamePosI[uChannelHV] + Form("/D" )  );
               cleanTree->Branch( fRootConvPar->sScOutputNamePosU[uChannelHV], &(vdPositiveHvU.at(uChannelHV)),
                                 fRootConvPar->sScOutputNamePosU[uChannelHV] + Form("/D" )  );
               cleanTree->Branch( fRootConvPar->sScOutputNameNegI[uChannelHV], &(vdNegativeHvI.at(uChannelHV)),
                                 fRootConvPar->sScOutputNameNegI[uChannelHV] + Form("/D" )  );
               cleanTree->Branch( fRootConvPar->sScOutputNameNegU[uChannelHV], &(vdNegativeHvU.at(uChannelHV)),
                                 fRootConvPar->sScOutputNameNegU[uChannelHV] + Form("/D" )  );
               vdPositiveHvI.at(uChannelHV) = 0;
               vdPositiveHvU.at(uChannelHV) = 0;
               vdNegativeHvI.at(uChannelHV) = 0;
               vdNegativeHvU.at(uChannelHV) = 0;
            } // for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
         } // if( 1 == fRootConvPar->uSlowControlInTree )

         if( 1 == fRootConvPar->uAddGet4DataInTree )
         {
            // Data from Get4-Roc
            vuHitNumber = new std::vector< UInt_t >;
            vuHitNumber->resize( fRocPar->uNbFeets*8 );

            vuTime      = new std::vector< std::vector< UInt_t > >;
            vuEpochTime = new std::vector< std::vector< UInt_t > >;
            viTot       = new std::vector< std::vector<  Int_t > >;

            vuTime->resize(      fRocPar->uNbFeets*8 );
            vuEpochTime->resize( fRocPar->uNbFeets*8 );
            viTot->resize(       fRocPar->uNbFeets*8 );

            for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
            {
               (*vuHitNumber)[uChannelTemp]  = 0;
               (*vuTime)[uChannelTemp].clear();
               (*vuEpochTime)[uChannelTemp].clear();
               (*viTot)[uChannelTemp].clear();
            }

            cleanTree->Bronch( "vuHitNumber",     "std::vector< UInt_t >",   &vuHitNumber );
            cleanTree->Bronch( "get4_time",       "std::vector< std::vector< UInt_t > >", &(vuTime) );
            cleanTree->Bronch( "get4_epoch_time", "std::vector< std::vector< UInt_t > >", &(vuEpochTime) );
            cleanTree->Bronch( "get4_tot",        "std::vector< std::vector<  Int_t > >", &(viTot) );
         }

         cleanTree->SetDirectory(gDirectory);

         if( 1 == fRootConvPar->uRootDisplay && 0 == TGo4Analysis::Instance()->GetTree( "CosyRpc" ))
         {
            TGo4Analysis::Instance()->AddTree(cleanTree );
            TGo4Log::Info("Root Tree registered in Go4: %p", TGo4Analysis::Instance()->GetTree( "CosyRpc" ) );
         }
      } // if( 0 == cleanTree)
   } // if( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay)
}


void TRootConverterProc::FinalizeEvent()
{
   if( 0 == fRootConvPar->uCleanOuputEnable && 0 == fRootConvPar->uRootDisplay)
         return;
        //if(!fRpcInputEvent->IsValid() &&  1 == fRootConvPar->uAddGet4DataInTree) return; // skip incomplete rpc events from second step
   if(!fCrateInputEvent->IsValid() ) return; // skip incomplete mbs events from first step
   if( 1 == fRootConvPar->uTriglogInputInTree )
      if( !fTriglogInputEvent->IsValid() )
         return; // skip incomplete triglog events from first step

// put additional RPC processing here 

   //loop over FEET rocs:
   Bool_t bTestAllValid = kTRUE;
   for(Int_t iFeetRoc=0; iFeetRoc < (Int_t)fRocPar->uNbRocsGet4; iFeetRoc ++)
   {
      if( kFALSE == FillHitsFromRoc( ) )
      {
         bTestAllValid = kFALSE;

         if( (1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) && 1 == fRootConvPar->uAddGet4DataInTree )
         {
            for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
            {
               (*vuHitNumber)[uChannelTemp]  = 0;
               (*vuTime)[uChannelTemp].clear();
               (*vuEpochTime)[uChannelTemp].clear();
               (*viTot)[uChannelTemp].clear();
            } // for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
         } //  if( (1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) && 1 == fRootConvPar->uAddGet4DataInTree )
      } // if( kFALSE == BuildHits( theRoc ) )
   } // for(Int_t iFeetRoc=0; iFeetRoc < (Int_t)fRocPar->uNbRocsGet4; iFeetRoc ++)

   uNbEvents++;

   // Fill the tree
   if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
   {
      if( 1 == fRootConvPar->uHodoInputInTree )
      {
         if( 1 == fRootConvPar->uHodoDebugHisto)
            hNbClustersPerEventHodoscope->Fill( fHodoInputEvent->NumHits() );

         if( 0 < fHodoInputEvent->NumHits() )
         {
            hitHodoscope = fHodoInputEvent->Hit(0);
            if( 1 == fRootConvPar->uTriglogInputInTree && 1 == fRootConvPar->uHodoDebugHisto )
            {
               if( 0 == uFirstMbsTime)
               {
                  uFirstMbsTime = fTriglogInputEvent->fMbsTimeSecs;
                  uSpillStartMbsTime = fTriglogInputEvent->fMbsTimeSecs- uFirstMbsTime;
                  uLastMbsTime = fTriglogInputEvent->fMbsTimeSecs;

                  uFirstMbsTimeMs = fTriglogInputEvent->fMbsTimeMillisec;
                  uSpillStartMbsTimeMs = fTriglogInputEvent->fMbsTimeMillisec- uFirstMbsTimeMs;
               }

               hHit2DvsSeconds->Fill( fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime,
                     fHodoInputEvent->Hit(0).X, fHodoInputEvent->Hit(0).Y);

               if( 20 < fTriglogInputEvent->fMbsTimeSecs - uLastMbsTime )
               {
                  uSpillStartMbsTime = fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime;
                  uSpillStartMbsTimeMs = fTriglogInputEvent->fMbsTimeMillisec - uFirstMbsTimeMs;
               }

               hHit2DvsSecondsInSpill->Fill(
                     fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime - uSpillStartMbsTime,
                     fHodoInputEvent->Hit(0).X, fHodoInputEvent->Hit(0).Y);

               hHit2DvsMilliSecondsInSpill->Fill(
                     (fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime - uSpillStartMbsTime)*1000 +
                     fTriglogInputEvent->fMbsTimeMillisec - uFirstMbsTimeMs - uSpillStartMbsTimeMs,
                     fHodoInputEvent->Hit(0).X, fHodoInputEvent->Hit(0).Y);

               hHitXvsMilliSecondsInSpill->Fill(
                     (fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime - uSpillStartMbsTime)*1000 +
                     fTriglogInputEvent->fMbsTimeMillisec - uFirstMbsTimeMs - uSpillStartMbsTimeMs,
                     fHodoInputEvent->Hit(0).X );
               hHitYvsMilliSecondsInSpill->Fill(
                     (fTriglogInputEvent->fMbsTimeSecs - uFirstMbsTime - uSpillStartMbsTime)*1000 +
                     fTriglogInputEvent->fMbsTimeMillisec - uFirstMbsTimeMs - uSpillStartMbsTimeMs,
                     fHodoInputEvent->Hit(0).Y );

               uLastMbsTime  = fTriglogInputEvent->fMbsTimeSecs;

            } // if( 1 == fRootConvPar->uTriglogInputInTree && 1 == fRootConvPar->uHodoDebugHisto )
         } // if( 0 < fHodoInputEvent->NumHits() )
      } // if( 1 == fRootConvPar->uHodoInputInTree )

      if( 1 == fRootConvPar->uSlowControlInTree )
      {
         if( kTRUE == fEpicsInputEvent->IsValid() )
         {
            TString sVariableName = "";
            // 100000 limit needed because of bug sometimes on Valid flag
//            if( fEpicsInputEvent->GetNumDoubles() <= 0 || 100000 <= fEpicsInputEvent->GetNumDoubles() )
//            {
//               cout<<"Test "<<kTRUE<<" "<<fEpicsInputEvent->IsValid()<<" "<<((fEpicsInputEvent->IsValid())&kTRUE);
//               cout<<" "<<fEpicsInputEvent->GetNumLongs();
//               cout<<" "<<fEpicsInputEvent->GetNumDoubles()<<" "<<fEpicsInputEvent<<endl;
//               uSlowControlUpdateTime = 0;
//               for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
//               {
//                  vdPositiveHvI.at(uChannelHV) = 0;
//                  vdPositiveHvU.at(uChannelHV) = 0;
//                  vdNegativeHvI.at(uChannelHV) = 0;
//                  vdNegativeHvU.at(uChannelHV) = 0;
//               }
//            }
//            else
//            {
//               cout<<"Test 1 "<<fEpicsInputEvent->GetNumDoubles()<<endl;
            if( 0 == uEpicsTimeOrigin )
               uEpicsTimeOrigin = fEpicsInputEvent->fUTimeSeconds;
            //hTimeEpicsUpdates->Fill(fEpicsInputEvent->fUTimeSeconds - uEpicsTimeOrigin);
            hTimeEpicsUpdates->Fill( (fEpicsInputEvent->fUTimeSeconds )% 1000);
            uSlowControlUpdateTime = fEpicsInputEvent->fUTimeSeconds;
            for(UInt_t uEpicsVariableIndex = 0; uEpicsVariableIndex < fEpicsInputEvent->GetNumDoubles(); uEpicsVariableIndex++ )
            {
               //cout<<"Test2 "<<uEpicsVariableIndex<<" "<<fEpicsInputEvent->GetNumDoubles()<<endl;
               sVariableName = fEpicsInputEvent->GetDoubleName( uEpicsVariableIndex );
               //cout<<"Test2 "<<sVariableName<<" "<<uEpicsVariableIndex<<" "<<fEpicsInputEvent->GetNumDoubles()<<endl;
               for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
               {
                  if( sVariableName == fRootConvPar->sScInputNamePosI[uChannelHV])
                  {
                     //cout<<"I Pos "<<uEpicsVariableIndex<<endl;
                     vdPositiveHvI.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     //cout<<"I Pos "<<uEpicsVariableIndex<<endl;
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNamePosU[uChannelHV])
                  {
                     //cout<<"U Pos "<<uEpicsVariableIndex<<endl;
                     vdPositiveHvU.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     //cout<<"U Pos "<<uEpicsVariableIndex<<endl;
                     //cout<<endl;
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNameNegI[uChannelHV])
                  {
                     // cout<<"I Neg "<<uEpicsVariableIndex<<endl;
                     vdNegativeHvI.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     //cout<<"I Neg "<<uEpicsVariableIndex<<endl;
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNameNegU[uChannelHV])
                  {
                     //cout<<"U Neg "<<uEpicsVariableIndex<<endl;
                     vdNegativeHvU.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     //cout<<"U Neg "<<uEpicsVariableIndex<<endl;
                     break;
                  }
               } //  for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
            } // for(UInt_t uEpicsVariableIndex = 0; uEpicsVariableIndex < fEpicsInputEvent->GetNumDoubles(); uEpicsVariableIndex++ )
//            } // else of if( fEpicsInputEvent->GetNumDoubles() <= 0 || 100000 <= fEpicsInputEvent->GetNumDoubles() )
         } // if( kTRUE == fEpicsInputEvent->IsValid() )
      } // if( 1 == fRootConvPar->uSlowControlInTree )

      cleanTree->Fill();
      if( 1 == fRootConvPar->uAddGet4DataInTree )
      {
         for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
         {
            (*vuHitNumber)[uChannelTemp]  = 0;
            (*vuTime)[uChannelTemp].clear();
            (*vuEpochTime)[uChannelTemp].clear();
            (*viTot)[uChannelTemp].clear();
         } // for( UInt_t uChannelTemp = 0; uChannelTemp < fRocPar->uNbFeets*8; uChannelTemp++)
      } // if( 1 == fRootConvPar->uAddGet4DataInTree )
   } // if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )

   fOutputEvent->SetValid(kTRUE);
}

Bool_t TRootConverterProc::FillHitsFromRoc()
{
   TGet4HitUnpacked hitTemp;

   for( Int_t iHitIndex = 0; iHitIndex < (Int_t)((fRpcInputEvent->fUnpackedHits).size()); iHitIndex ++)
   {
      hitTemp = (fRpcInputEvent->fUnpackedHits).at(iHitIndex);
      // Fill root tree
      if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) && 1 == fRootConvPar->uAddGet4DataInTree )
      {
         UInt_t uFullChan = hitTemp.GetFullChan();
         (*vuHitNumber)[uFullChan] ++;
         (*vuTime)[uFullChan].push_back( hitTemp.GetTsLeading() );
         (*vuEpochTime)[uFullChan].push_back( hitTemp.GetEpochLeading() );
         (*viTot)[uFullChan].push_back( hitTemp.GetRawTOT() );
      } // if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) && 1 == fRootConvPar->uAddGet4DataInTree )
   }

   return kTRUE;
}