beamtime/gsi-aug12/prod/go4/ROOTCONV/TRootConverterProc.cxx (r4864/r3870)

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),
fCrateInputEvent(0),
fTriglogInputEvent(0),
fEpicsInputEvent(0),
fVftxInputEvent(0),
fRpcHdInputEvent(0),
fRpcBucInputEvent(0),
fRpcBucMod1InputEvent(0),
fRpcBucMod2InputEvent(0),
fRpcBucMod3InputEvent(0),
fRpcBucMod4InputEvent(0),
fPlasticsInputEvent(0),
fPadRpcInputEvent(0),
fScalOrMu1InputEvent(0),
fScalOrMu2InputEvent(0),
fMatchingInputEvent(0),
#ifdef ROOTCONV_WITHHODO
fHodoInputEvent(0),
#endif
fOutputEvent(0),
cleanTree(0)
{
   TGo4Log::Info("**** TRootConverterProc: Create instance %s", name);

#if defined(ROOTCONV_WITHHODO)
   fRocPar = (TRocParam*) GetParameter("RocPar");
   cout << "ROCPAR = " << fRocPar << endl;
   
   if (fRocPar->numRocs > 0) fNumRocs = (UInt_t) fRocPar->numRocs;
#endif


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


   // ROC general variables
   uNbEvents = 0;
   
#ifdef ROOTCONV_WITHHODO
   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)
#endif

   //~ TString path = GetName();
   //~ path += "/";
   //~ path += GetName();
   if( 1 == fRootConvPar->uSlowControlInTree )
   {
      hTimeEpicsUpdates = MakeTH1('I', "RootConv/RC_TimeEpicsUpdate","Timing of epics updates; Time [s]",
            7200, 0, 7200);
   } // if( 1 == fRootConvPar->uSlowControlInTree )
   uEpicsTimeOrigin = 0;
      
   //printf("RPC Histograms created \n");
   fflush ( stdout);
}

TRootConverterProc::~TRootConverterProc()
{
//   TGo4Log::Info("Proc  %p", this);
   if( (1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay)
         && 0 != cleanTree )
      {
          if( 1 == fRootConvPar->uRootDisplay )
             TGo4Analysis::Instance()->RemoveTree(cleanTree);
          if( 1 == fRootConvPar->uCleanOuputEnable )
          {
             fRootOut->Write("",TObject::kOverwrite);
             //~ gDirectory->Cd(sOutputFilename);
//~ 
             //~ cleanTree->SetDirectory(gDirectory);
             //~ cleanTree->Write(0,TObject::kOverwrite); // to avoid double trees in file: cf http://root.cern.ch/root/roottalk/roottalk01/0842.html
             //~ gDirectory->Cd(oldDir->GetPath());
             fRootOut->Close();
          } // if( 1 == fRootConvPar->uCleanOuputEnable )

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

          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
#ifdef ROOTCONV_WITHHODO
   if( 0==fHodoInputEvent && 1 == fRootConvPar->uHodoInputInTree )
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Detector"));
      if(btevent)
      {
         fHodoInputEvent=dynamic_cast<TFiberHodEvent*>(btevent->GetSubEvent("Hodo1"));
      }
      if( 0==fHodoInputEvent ) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event has no TFiberHodEvent!!! STOP GO4");
      }
   }
#endif

   // 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 defined(ROOTCONV_WITHHODO)
   if( kFALSE == fRocPar->noTrigger && ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
#else
   if( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay )
#endif
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Calibration"));
      if(btevent)
      {
         if(( 1 == fRootConvPar->uCaenInTree || 1 == fRootConvPar->uTriglogInputInTree)
             && fCrateInputEvent==0)
         {
            fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(btevent->GetSubEvent("MBSCRATE"));

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

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

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

      btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Detector"));
      if(btevent)
      {
         if( 0 ==fVftxInputEvent && 1 == fRootConvPar->uVftxInTree)
         {
               fVftxInputEvent=dynamic_cast<TVftxEvent*>(btevent->GetSubEvent("VFTX"));

            if(fVftxInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TVftxEvent!!! STOP GO4");
            }
         } // if( 0 == fVftxInputEvent && 1 == fRootConvPar->uVftxInTree)
      } // if(btevent)
         else if( 1 == fRootConvPar->uVftxInTree )
            GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Detector input event not found!!! STOP GO4");

      // Param need to be done!!!!
      btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Tracking"));
      if(btevent)
      {
         if( 1 == fRootConvPar->uRpcHdInTree  && fRpcHdInputEvent  ==0 )
         {
            fRpcHdInputEvent  =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_HD"));

            if( fRpcHdInputEvent ==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for HD!!! STOP GO4");
            }
         } // if( 1 == fRootConvPar->uRpcHdInTree  && fRpcHdInputEvent  ==0 )

         if( ( 1 == fRootConvPar->uRpcBucInTree && fRpcBucInputEvent ==0 ) )
         {
            fRpcBucInputEvent =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_BUC_REF"));

            if( fRpcBucInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for Buc!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uRpcBucInTree && fRpcBucInputEvent ==0 ) )

         if( ( 1 == fRootConvPar->uRpcBucM1InTree && fRpcBucMod1InputEvent ==0 ) )
         {
            fRpcBucMod1InputEvent =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_BUC_MOD1"));

            if( fRpcBucMod1InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for Buc Mod1!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uRpcBucM1InTree && fRpcBucMod1InputEvent ==0 ) )

         if( ( 1 == fRootConvPar->uRpcBucM2InTree && fRpcBucMod2InputEvent ==0 ) )
         {
            fRpcBucMod2InputEvent =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_BUC_MOD2"));

            if( fRpcBucMod2InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for Buc Mod2!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uRpcBucM2InTree && fRpcBucMod2InputEvent ==0 ) )

         if( ( 1 == fRootConvPar->uRpcBucM3InTree && fRpcBucMod3InputEvent ==0 ) )
         {
            fRpcBucMod3InputEvent =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_BUC_MOD3"));

            if( fRpcBucMod3InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for Buc Mod3!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uRpcBucM3InTree && fRpcBucMod3InputEvent ==0 ) )

         if( ( 1 == fRootConvPar->uRpcBucM4InTree && fRpcBucMod4InputEvent ==0 ) )
         {
            fRpcBucMod4InputEvent =dynamic_cast<TRpcEvent*>(btevent->GetSubEvent("RPC_BUC_MOD4"));

            if( fRpcBucMod4InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TRpcEvent for Buc Mod4!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uRpcBucM4InTree && fRpcBucMod4InputEvent ==0 ) )

         /* Plastics */
         if( ( 1 == fRootConvPar->uPlasticsInTree && fPlasticsInputEvent ==0 ) )
         {
            fPlasticsInputEvent =dynamic_cast<TPlasticsEvent*>(btevent->GetSubEvent("PLASTICS"));

            if( fPlasticsInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TPlasticsEvent for Plastics!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uPlasticsInTree && fPlasticsInputEvent ==0 ) )

         /* Pad Rpc */
         if( ( 1 == fRootConvPar->uPadRpcInTree && fPadRpcInputEvent ==0 ) )
         {
            fPadRpcInputEvent =dynamic_cast<TPadRpcEvent*>(btevent->GetSubEvent("PADRPC"));

            if( fPadRpcInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TPadRpcEvent for Pad Rpc!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uPadRpcInTree && fPadRpcInputEvent ==0 ) )

         if( 0 == fScalOrMu1InputEvent && 1 == fRootConvPar->uScalOrMuInTree )
         {
            fScalOrMu1InputEvent=dynamic_cast<TScalersEvent*>(btevent->GetSubEvent("SCALERS"));

            if(fScalOrMu1InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TScalersEvent!!! STOP GO4");
            }
         } // if( 1 == fRootConvPar->uTriglogInputInTree )
         if( 0 == fScalOrMu2InputEvent && 1 == fRootConvPar->uScalOrMu2InTree )
         {
            fScalOrMu2InputEvent=dynamic_cast<TScalersEvent*>(btevent->GetSubEvent("SCALERS2"));

            if(fScalOrMu2InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TScalersEvent!!! STOP GO4");
            }
         } // if( 1 == fRootConvPar->uTriglogInputInTree )
      } // if(btevent)
         else if(  1 == fRootConvPar->uRpcHdInTree ||
                   1 == fRootConvPar->uRpcBucInTree ||
                   1 == fRootConvPar->uRpcBucM1InTree ||
                   1 == fRootConvPar->uRpcBucM2InTree ||
                   1 == fRootConvPar->uRpcBucM3InTree ||
                   1 == fRootConvPar->uRpcBucM4InTree ||
                   1 == fRootConvPar->uPlasticsInTree ||
                   1 == fRootConvPar->uPadRpcInTree ||
                   1 == fRootConvPar->uScalOrMuInTree ||
                   1 == fRootConvPar->uScalOrMu2InTree )
            GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Tracking input event not found!!! STOP GO4");

      if( 1 == fRootConvPar->uMatchingInTree )
      {
         if(fMatchingInputEvent==0)
         {
            TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent());
            if(btevent)
            {
               fMatchingInputEvent=dynamic_cast<TMatchingEvent*>(btevent->GetSubEvent("MATCHING"));
            }

            if(fMatchingInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: input event is/has not a TMatchingEvent!!! STOP GO4");
            }
         } // if(fMatchingInputEvent==0)
      } // if( 1 == fRootConvPar->uMatchingInTree )

   } // if( kFALSE == fRocPar->noTrigger && ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
      else if( 1 == fRootConvPar->uSlowControlInTree )
      {
         TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Calibration"));
         if(fEpicsInputEvent==0)
         {
            if(btevent)
            {
               fEpicsInputEvent=dynamic_cast<TEpicsEvent*>(btevent->GetSubEvent("EPICS"));
            }

            if(fEpicsInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: 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( "GsiAug12" ) )
         {
            cleanTree = new TTree("GsiAug12", "GSI August 2012 data");
            TGo4Log::Info("Output Tree created in %s", gDirectory->GetPath());
         } // if( 0 == TGo4Analysis::Instance()->GetTree( "GsiAug12" ) )
            else
            {
               cleanTree = TGo4Analysis::Instance()->GetTree( "GsiAug12" );
               TGo4Log::Info("Output Tree recovered from Go4");
            } // else of if( 0 == TGo4Analysis::Instance()->GetTree( "GsiAug12" ) )

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

         // Data from CAEN v1290 & Lecroy QDC VME boards (triggered)
#if defined(ROOTCONV_WITHHODO)
         if( kFALSE == fRocPar->noTrigger )
#endif
         {
            if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel1 && 1 == fRootConvPar->uCaenInTree)
            {
               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 && 1 == fRootConvPar->uCaenInTree )
            if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel2 && 1 == fRootConvPar->uCaenInTree )
            {
               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 && 1 == fRootConvPar->uCaenInTree )
            
            if( 1 == fRootConvPar->uCaenMode )
            {
               for( Int_t iCaenBoard = 0; iCaenBoard<MAX_1290; iCaenBoard ++)
                  if( 1 == fRootConvPar->uGeneralCaenTdcEnabled[iCaenBoard])
                  {
                     cleanTree->Branch(Form("Caen%d_TriggerTime", iCaenBoard),  
                                       &( (fCrateInputEvent->fMtdc[iCaenBoard]).extended_trigger_time ),
                                       Form("Caen%d_TriggerTime/I", iCaenBoard )  );
                                       
                     // The multiple hits array for CAEN v1290 now also include first hit
                     // => either add to tree the first hit or add to tree only the multiple hit array
                     if( 1 == fRootConvPar->uCaenMulti )
                     {
               cleanTree->Branch(Form("Caen%d_Leading_NbHit", iCaenBoard), 
                                 (fCrateInputEvent->fMtdc[iCaenBoard]).hit_lead,
                                 Form("Caen%d_Leading_NbHit[%u]/I", iCaenBoard, T1290Data::NumChan ) );
               cleanTree->Branch(Form("Caen%d_Trailing_NbHit", iCaenBoard), 
                                 (fCrateInputEvent->fMtdc[iCaenBoard]).hit_trail,
                                 Form("Caen%d_Trailing_NbHit[%u]/I", iCaenBoard, T1290Data::NumChan ) );
               cleanTree->Branch(Form("Caen%d_Leading_multi", iCaenBoard),  
                                 (fCrateInputEvent->fMtdc[iCaenBoard]).lead_multi,
                                 Form("Caen%d_Leading_multi[%u][%u]/I", iCaenBoard, T1290Data::NumChan, T1290Data::MaxHits )  );
               cleanTree->Branch(Form("Caen%d_Trailing_multi", iCaenBoard), 
                                 (fCrateInputEvent->fMtdc[iCaenBoard]).trail_multi,
                                 Form("Caen%d_Trailing_multi[%u][%u]/I", iCaenBoard, T1290Data::NumChan, T1290Data::MaxHits ) );
                     } // if( 1 == fRootConvPar->uCaenMulti )
                        else
                        {
                     cleanTree->Branch(Form("Caen%d_Leading", iCaenBoard),  
                                       (fCrateInputEvent->fMtdc[iCaenBoard]).lead,
                                       Form("Caen%d_Leading[%u]/I", iCaenBoard, T1290Data::NumChan )  );
                     cleanTree->Branch(Form("Caen%d_Trailing", iCaenBoard), 
                                       (fCrateInputEvent->fMtdc[iCaenBoard]).trail,
                                       Form("Caen%d_Trailing[%u]/I", iCaenBoard, T1290Data::NumChan ) );
                        } // else of if( 1 == fRootConvPar->uCaenMulti )
                  } // if( 1 == fRootConvPar->uGeneralCaenTdcEnabled[iCaenBoard])
            } // if( 1 == fRootConvPar->uCaenMode )
            
            if( -1 != fRootConvPar->iPmtTdc && 0 == fRootConvPar->uCaenMode)
            {
               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 ) );
            } // if( -1 != fRootConvPar->iPmtTdc )
            if( -1 != fRootConvPar->iRpcTdcLeading && -1 != fRootConvPar->iRpcTdcTrailing && 0 == fRootConvPar->uCaenMode )
            {
               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 && 0 == fRootConvPar->uCaenMode )
            
            if( -1 != fRootConvPar->iPmtQdc )
            {
               cleanTree->Branch("PmtQdc", fCrateInputEvent->fData1182[fRootConvPar->iPmtQdc], Form("PmtQdc[%u]/s", NUM_1182_CH ) );
            } // if( -1 != fRootConvPar->iPmtQdc )

            if( 1 == fRootConvPar->uTriglogInputInTree)
            {
               cleanTree->Branch("TriglogSyncNb",  &(fTriglogInputEvent->fVulomSyncNumber),
                                 "TriglogSyncNb/i"  );
               cleanTree->Branch("TriglogPattern",  &(fTriglogInputEvent->fVulomTriggerPattern),
                                 "TriglogPattern/i"  );
               cleanTree->Branch("TriglogInputPattern",  &(fCrateInputEvent->fTriglogInputPattern),
                                 "TriglogInputPattern/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( 1 == fRootConvPar->uVftxInTree)
            {
               for( Int_t iVftxBoard = 0; iVftxBoard<MAX_FPGA_TDC; iVftxBoard ++)
                  if( 1 == fRootConvPar->uVftxEnabled[iVftxBoard])
                  {
                     TString sBrNameTrigg  = Form("Vftx%dTriggerTime", iVftxBoard);
                     TString sBrNameMulti  = Form("Vftx%dMultiplicity", iVftxBoard);
                     TString sBrNameFineT  = Form("Vftx%dFineTime", iVftxBoard);
                     TString sBrNameCoarT  = Form("Vftx%dCoarseTime", iVftxBoard);
                     TString sBrNameFutuB  = Form("Vftx%dFutureBit", iVftxBoard);
                     TString sBrNameTimeC  = Form("Vftx%dTimeCorr", iVftxBoard);
                     TString sBrNameFtCorr = Form("Vftx%dFineTimeCorr", iVftxBoard);

               cleanTree->Branch( sBrNameTrigg,  &((fVftxInputEvent->fVftxBoards[iVftxBoard]).iTriggerTime),
                                  sBrNameTrigg + "/I" );
               cleanTree->Branch( sBrNameMulti,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).iMultiplicity),
                                  sBrNameMulti + Form("[%u]/I", FPGA_TDC_NBCHAN )  );
                     if( 1 == fRootConvPar->uVftxFullDataMode )
                     {
               cleanTree->Branch( sBrNameFineT,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).iFineTime),
                                  sBrNameFineT + Form("[%u][%u]/I", FPGA_TDC_NBCHAN, TVftxBoardData::MaxMult )  );
               cleanTree->Branch( sBrNameCoarT,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).iCoarseTime),
                                  sBrNameCoarT + Form("[%u][%u]/I", FPGA_TDC_NBCHAN, TVftxBoardData::MaxMult )  );
               cleanTree->Branch( sBrNameFutuB,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).uFutureBit),
                                  sBrNameFutuB + Form("[%u][%u]/i", FPGA_TDC_NBCHAN, TVftxBoardData::MaxMult )  );
               cleanTree->Branch( sBrNameFtCorr,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).dFineTimeCorr),
                                  sBrNameFtCorr + Form("[%u][%u]/D", FPGA_TDC_NBCHAN, TVftxBoardData::MaxMult )  );
                     } //if( 1 == fRootConvPar->uVftxFullDataMode )

               cleanTree->Branch( sBrNameTimeC,  ((fVftxInputEvent->fVftxBoards[iVftxBoard]).dTimeCorr),
                                  sBrNameTimeC + Form("[%u][%u]/D", FPGA_TDC_NBCHAN, TVftxBoardData::MaxMult )  );
                  } // if( 1 == fRootConvPar->uVftxEnabled[iVftxBoard])
            } // if( 1 == fRootConvPar->uVftxInTree)
            if( 1 == fRootConvPar->uScalOrMuInTree)
            {
               cleanTree->Branch("ScalOrMuScalers",  (fCrateInputEvent->fScalOrMu).uScaler,
                                 Form("ScalOrMuScalers[%u]/i", SCALORMU_NB_SCAL )  );
            } // if( 1 == fRootConvPar->uVftxInTree)

            if( 1 == fRootConvPar->uRpcHdInTree)
            {
               cleanTree->Branch("RpcHdEvent",      "TRpcEvent", &fRpcHdInputEvent);
            }
            if( 1 == fRootConvPar->uRpcBucInTree)
            {
               cleanTree->Branch("RpcBucEvent",     "TRpcEvent", &fRpcBucInputEvent);
            }
            if( 1 == fRootConvPar->uRpcBucM1InTree)
            {
               cleanTree->Branch("RpcBucMod1Event", "TRpcEvent", &fRpcBucMod1InputEvent);
            }
            if( 1 == fRootConvPar->uRpcBucM2InTree)
            {
               cleanTree->Branch("RpcBucMod2Event", "TRpcEvent", &fRpcBucMod2InputEvent);
            }
            if( 1 == fRootConvPar->uRpcBucM3InTree)
            {
               cleanTree->Branch("RpcBucMod3Event", "TRpcEvent", &fRpcBucMod3InputEvent);
            }
            if( 1 == fRootConvPar->uRpcBucM4InTree)
            {
               cleanTree->Branch("RpcBucMod4Event", "TRpcEvent", &fRpcBucMod4InputEvent);
            }
            if( 1 == fRootConvPar->uPlasticsInTree)
            {
               cleanTree->Branch("PlasticsEvent",   "TPlasticsEvent", &fPlasticsInputEvent);
            }
            if( 1 == fRootConvPar->uPadRpcInTree)
            {
               cleanTree->Branch("PadRpcEvent",     "TPadRpcEvent", &fPadRpcInputEvent);
            }
            if( 1 == fRootConvPar->uScalOrMuInTree)
            {
               cleanTree->Branch("ScalOrMu1Event",  "TScalersEvent", &fScalOrMu1InputEvent);
            }
            if( 1 == fRootConvPar->uScalOrMu2InTree)
            {
               cleanTree->Branch("ScalOrMu2Event",  "TScalersEvent", &fScalOrMu2InputEvent);
            }

            if( 1 == fRootConvPar->uMatchingInTree)
            {
               cleanTree->Branch("MatchingEvent",  "TMatchingEvent", &fMatchingInputEvent);
            }
         } // if( kFALSE == fRocPar->noTrigger )
         
#ifdef ROOTCONV_WITHHODO
         if( 1 == fRootConvPar->uHodoInputInTree )
         {
            cleanTree->Branch("HodoscopeX",  &(hitHodoscope.X),"HodoscopeX/D");
            cleanTree->Branch("HodoscopeY",  &(hitHodoscope.Y),"HodoscopeY/D");
         } // if( 1 == fRootConvPar->uHodoInputInTree )
#endif

         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 )

         cleanTree->SetDirectory(gDirectory);

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


void TRootConverterProc::FinalizeEvent()
{
   if( 0 == fRootConvPar->uCleanOuputEnable && 0 == fRootConvPar->uRootDisplay)
         return;
   if( 0 == fRootConvPar->uCleanOuputEnable && 1 == fRootConvPar->uRootDisplay &&
       fRootConvPar->uEventsUntilReset <= cleanTree->GetEntries())
      cleanTree->Reset("");

   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
   if( 1 == fRootConvPar->uVftxInTree )
      if( !fVftxInputEvent->IsValid() )
         cout<<"Invalid VFTX event! event "<<uNbEvents<<endl;

   uNbEvents++;

   // Fill the tree
   if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
   {
#ifdef ROOTCONV_WITHHODO
      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 )
#endif
      if( 1 == fRootConvPar->uSlowControlInTree )
      {
         if( kTRUE == fEpicsInputEvent->IsValid() )
         {
            TString sVariableName = "";
            if( 0 == uEpicsTimeOrigin )
               uEpicsTimeOrigin = fEpicsInputEvent->fUTimeSeconds;
            hTimeEpicsUpdates->Fill( (fEpicsInputEvent->fUTimeSeconds )% 1000);
            uSlowControlUpdateTime = fEpicsInputEvent->fUTimeSeconds;
            for(UInt_t uEpicsVariableIndex = 0; uEpicsVariableIndex < fEpicsInputEvent->GetNumDoubles(); uEpicsVariableIndex++ )
            {
               sVariableName = fEpicsInputEvent->GetDoubleName( uEpicsVariableIndex );
               for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
               {
                  if( sVariableName == fRootConvPar->sScInputNamePosI[uChannelHV])
                  {
                     vdPositiveHvI.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNamePosU[uChannelHV])
                  {
                     vdPositiveHvU.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNameNegI[uChannelHV])
                  {
                     vdNegativeHvI.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     break;
                  }
                  if( sVariableName == fRootConvPar->sScInputNameNegU[uChannelHV])
                  {
                     vdNegativeHvU.at(uChannelHV) = fEpicsInputEvent->GetDouble(uEpicsVariableIndex);
                     break;
                  }
               } //  for( UInt_t uChannelHV = 0; uChannelHV < fRootConvPar->uNbSlowControlChan; uChannelHV++)
            } // for(UInt_t uEpicsVariableIndex = 0; uEpicsVariableIndex < fEpicsInputEvent->GetNumDoubles(); uEpicsVariableIndex++ )
         } // if( kTRUE == fEpicsInputEvent->IsValid() )
      } // if( 1 == fRootConvPar->uSlowControlInTree )

      cleanTree->Fill();
   } // if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )

   fOutputEvent->SetValid(kTRUE);
}