beamtime/tof-tdctest/go4/ROOTCONV/TRootConverterProc.cxx (r4864/r4613)

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),
fGet4v1InputEvent(0),
fCrateInputEvent(0),
fTriglogInputEvent(0),
fEpicsInputEvent(0),
fVftxInputEvent(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);

   for( UInt_t uRpcIndex = 0; uRpcIndex < NB_RPC; uRpcIndex++ )
      fRpcInputEvent[uRpcIndex] = 0;
   for( UInt_t uPlasticIndex = 0; uPlasticIndex < NB_PLASTICS; uPlasticIndex++ )
      fPlasticsInputEvent[uPlasticIndex] = 0;

#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;

   fParAnalysis = (TGsiAug12Param*) GetParameter("GsiAug12Par");

   if( kTRUE == fParAnalysis->bWithGet4 )
   {
      fGet4v1Par = (TGet4v1Param*) GetParameter("Get4v1Par");
      if( 0 == fGet4v1Par )
         fGet4v1Par = (TGet4v1Param*) MakeParameter("Get4v1Par", "TGet4v1Param");
   } // if( kTRUE == fParAnalysis->bWithGet4 )

   // Initialize detector parameter pointers to 0
   for( UInt_t uRpcIndex = 0; uRpcIndex < NB_RPC; uRpcIndex++ )
      fParRpc[uRpcIndex] = 0;
   for( UInt_t uPlasticIndex = 0; uPlasticIndex < NB_PLASTICS; uPlasticIndex++ )
      fParPlastics[uPlasticIndex] = 0;

   // 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 && kTRUE == fParAnalysis->bWithEpics )
   {
      hTimeEpicsUpdates = MakeTH1('I', "RootConv/RC_TimeEpicsUpdate","Timing of epics updates; Time [s]",
            7200, 0, 7200);
   } // if( 1 == fRootConvPar->uSlowControlInTree && kTRUE == fParAnalysis->bWithEpics )
   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 && kTRUE == fParAnalysis->bWithEpics )
          {
             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*>(GetOutputEvent("Unpack"));
      if(btevent)
      {
         if( 1 == fRootConvPar->uGet4InTree && fGet4v1InputEvent==0 && kTRUE == fParAnalysis->bWithGet4 )
         {
            fGet4v1InputEvent=dynamic_cast<TGet4v1Event*>(btevent->GetSubEvent("ROCGET4V10"));

            if(fGet4v1InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Unpack Output event is/has not a TGet4v1Event!!! STOP GO4");
            }
         } // if( 1 == fRootConvPar->uGet4InTree && fGet4v1InputEvent==0 && kTRUE == fParAnalysis->bWithGet4 )
         if(( (1 == fRootConvPar->uCaenInTree && kTRUE == fParAnalysis->bWithMbs ) || 1 == fRootConvPar->uTriglogInputInTree)
             && fCrateInputEvent==0)
         {
            fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(btevent->GetSubEvent("MBSCRATE"));

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

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

            if(fEpicsInputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Unpack Output  event is/has not a TEpicsEvent!!! STOP GO4");
            }
         } // if( 0 == fEpicsInputEvent && 1 == fRootConvPar->uSlowControlInTree && kTRUE == fParAnalysis->bWithEpics )
      } // if(btevent)
         else if( ( 1 == fRootConvPar->uCaenInTree && kTRUE == fParAnalysis->bWithMbs ) ||
                   ( 1 == fRootConvPar->uTriglogInputInTree && kTRUE == fParAnalysis->bWithTriglog ) ||
                   ( 1 == fRootConvPar->uSlowControlInTree && kTRUE == fParAnalysis->bWithEpics ) )
            GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Unpack Output event not found!!! STOP GO4");

      btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetOutputEvent("Calibration"));
      if(btevent)
      {
         if( 0 ==fVftxInputEvent && 1 == fRootConvPar->uVftxInTree && kTRUE == fParAnalysis->bWithVftxCalib )
         {
               fVftxInputEvent=dynamic_cast<TVftxEvent*>(btevent->GetSubEvent("VFTX"));

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

      // Param need to be done!!!!
      btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetOutputEvent("Detector"));
      if(btevent)
      {
         if( 1 == fRootConvPar->uRpcInTree )
            for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )
            {
               if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
               {
                  if( 0 == fParRpc[uRpcIndex])
                  {
                     fParRpc[uRpcIndex] = (TRpcParam*) GetParameter( Form("RpcPar_%02d", uRpcIndex));
                     if( 0 == fParRpc[uRpcIndex])
                     {
                        fRootConvPar->uRpcEnabled[uRpcIndex] = 0;
                        cout<<"Parameter for Rpc"<<uRpcIndex<<" not found => Disabled in Event output"<<endl;
                        continue; // Jump this RPC as disabled anyway
                     } // if( 0 == fParRpc[uRpcIndex])
                     if( 0 == fParRpc[uRpcIndex]->uNbStrips )
                     {
                        fRootConvPar->uRpcEnabled[uRpcIndex] = 0;
                        cout<<"Rpc"<<uRpcIndex<<" empty (uNbStrips = 0) => Disabled in Event output"<<endl;
                        continue; // Jump this RPC as disabled anyway
                     } // if( 0 == fParRpc[uRpcIndex]->uNbStrips )
                  } // if( 0 == fParRpc[uRpcIndex])

                  if( 0 == fRpcInputEvent[uRpcIndex] )
                  {
                     fRpcInputEvent[uRpcIndex]  =dynamic_cast<TRpcEvent*>( btevent->GetSubEvent(
                           Form("RPC_%02d", uRpcIndex) ) );

                     if( fRpcInputEvent[uRpcIndex] ==0) {
                        GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Detector Output event is/has not a TRpcEvent for Rpc %02d!!! STOP GO4",
                              uRpcIndex);
                     }
                  } // if(  0 == fRpcInputEvent[uRpcIndex] )
               } // if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
            } // for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )

         /* Plastics */
         if( ( 1 == fRootConvPar->uPlasticsInTree  ) )
            for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )
            {
               if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
               {
                  if( 0 == fParPlastics[uPlasticIndex] )
                  {
                     fParPlastics[uPlasticIndex] = (TPlasticsParam*) GetParameter( Form("PlasticsPar_%02d", uPlasticIndex));
                     if( 0 == fParPlastics[uPlasticIndex])
                     {
                        fRootConvPar->uPlasticsEnabled[uPlasticIndex] = 0;
                        cout<<"Parameter for Plastics"<<uPlasticIndex<<" not found => Disabled in Event output"<<endl;
                        continue; // Jump this PLASTICS as disabled anyway
                     } // if( 0 == fParPlastics[uPlasticIndex])
                     if( 0 == fParPlastics[uPlasticIndex]->uNbPlastics )
                     {
                        fRootConvPar->uPlasticsEnabled[uPlasticIndex] = 0;
                        cout<<"Plastics"<<uPlasticIndex<<" empty (uNbPlastics = 0) => Disabled in Event output"<<endl;
                        continue; // Jump this PLASTICS as disabled anyway
                     } // if( 0 == fParPlastics[uPlasticIndex]->uNbPlastics )
                  } // if( 0 == fParPlastics[uPlasticIndex] )
                  if( 0 == fPlasticsInputEvent[uPlasticIndex] )
                  {
                     fPlasticsInputEvent[uPlasticIndex] =dynamic_cast<TPlasticsEvent*>(btevent->GetSubEvent(
                           Form("PLASTICS_%02d", uPlasticIndex) ) );

                     if( fPlasticsInputEvent[uPlasticIndex]==0) {
                        GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Detector Output event is/has not a TPlasticsEvent for Plastics %02d!!! STOP GO4",
                              uPlasticIndex);
                     }
                  } // if( 0 == fPlasticsInputEvent[uPlasticIndex] )
               } // if( ( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] && 0 == fPlasticsInputEvent[uPlasticIndex] ) )
            } // for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )

         /* 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: Detector Output event is/has not a TPadRpcEvent for Pad Rpc!!! STOP GO4");
            }
         } // if( ( 1 == fRootConvPar->uPadRpcInTree && fPadRpcInputEvent ==0 ) )
*/
         if( 0 == fScalOrMu1InputEvent && 1 == fRootConvPar->uScalOrMuInTree && kTRUE == fParAnalysis->bWithScalersRate )
         {
            fScalOrMu1InputEvent=dynamic_cast<TScalersEvent*>(btevent->GetSubEvent("SCALERS"));

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

            if(fScalOrMu2InputEvent==0) {
               GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Detector Output event is/has not a TScalersEvent!!! STOP GO4");
            }
         } // if( 0 == fScalOrMu2InputEvent && 1 == fRootConvPar->uScalOrMu2InTree )
      } // if(btevent)
         else if( 1 == fRootConvPar->uRpcInTree ||
                   1 == fRootConvPar->uPlasticsInTree ||
                   1 == fRootConvPar->uPadRpcInTree ||
                   1 == fRootConvPar->uScalOrMuInTree ||
                   1 == fRootConvPar->uScalOrMu2InTree )
            GO4_STOP_ANALYSIS_MESSAGE("**** TRootConverterProc: Fatal error: Detector Output event not found!!! STOP GO4");

      if( 1 == fRootConvPar->uMatchingInTree && kTRUE == fParAnalysis->bWithMatching)
      {
         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 && kTRUE == fParAnalysis->bWithMatching )

   } // if( kFALSE == fRocPar->noTrigger && ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
      else if( 1 == fRootConvPar->uSlowControlInTree && kTRUE == fParAnalysis->bWithEpics )
      {
//         TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent("Calibration"));
         TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetOutputEvent("Unpack"));
         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 && kTRUE == fParAnalysis->bWithEpics  )

   // 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
                  && kTRUE == fParAnalysis->bWithMbs )
            {
               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 && kTRUE == fParAnalysis->bWithMbs )
            if( -1 != fRootConvPar->iReferenceTdc && -1 != fRootConvPar->iReferenceChannel2 && 1 == fRootConvPar->uCaenInTree
                  && kTRUE == fParAnalysis->bWithMbs )
            {
               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 && kTRUE == fParAnalysis->bWithMbs )
            
            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 && kTRUE == fParAnalysis->bWithTriglog )
            {
               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 && kTRUE == fParAnalysis->bWithTriglog )
            if( 1 == fRootConvPar->uGet4InTree && kTRUE == fParAnalysis->bWithGet4 )
            {
               if( 1 == fRootConvPar->uMultiEventsAllowed)
               {
                  cleanTree->Branch("Get4v1Event",      "TGet4v1Event", &fGet4v1InputEvent);
               } // if( 1 == fRootConvPar->uMultiEventsAllowed)
                  else
                  {
                     for( Int_t iGet4Board = 0; iGet4Board<MAX_GET4; iGet4Board ++)
                        if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
                        {
                           // No multi event support !!!!!
                           /*
                           cleanTree->Branch( Form("Get4%d", iGet4Board), "Get4v1Tdc",
                                 &((fGet4v1InputEvent->fEvents[0]).fGet4Boards[iGet4Board] ));
                                 */
                           cleanTree->Branch( Form("Get4_%d", iGet4Board), "Get4v1Tdc", &( fGet4Board[iGet4Board] ) );
                        }
                  } // else of if( 1 == fRootConvPar->uMultiEventsAllowed)
            }
            if( 1 == fRootConvPar->uVftxInTree && kTRUE == fParAnalysis->bWithVftxCalib )
            {
               for( Int_t iVftxBoard = 0; iVftxBoard<MAX_FPGA_TDC; iVftxBoard ++)
                  if( 1 == fRootConvPar->uVftxEnabled[iVftxBoard])
                  {
                     if( 1 == fRootConvPar->uVftxFullBoardData )
                     {
                        cleanTree->Branch( Form("Vftx%02dData", iVftxBoard), "TVftxBoardData",
                                             &(fVftxInputEvent->fVftxBoards[iVftxBoard]) );
                     }
                        else
                        {
                           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 )  );
                        } // else of if( 1 == fRootConvPar->uVftxFullBoardData )
                  } // if( 1 == fRootConvPar->uVftxEnabled[iVftxBoard])
            } // if( 1 == fRootConvPar->uVftxInTree && kTRUE == fParAnalysis->bWithVftxCalib )
            if( 1 == fRootConvPar->uScalOrMuInTree && kTRUE == fParAnalysis->bWithScalersRate )
            {
               cleanTree->Branch("ScalOrMuScalers",  (fCrateInputEvent->fScalOrMu).uScaler,
                                 Form("ScalOrMuScalers[%u]/i", SCALORMU_NB_SCAL )  );
            } // if( 1 == fRootConvPar->uScalOrMuInTree && kTRUE == fParAnalysis->bWithScalersRate )
            if( 1 == fRootConvPar->uRpcInTree )
            {
               if( 1 == fRootConvPar->uMultiEventsAllowed)
               {
                  for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )
                     if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
                     {
                        cleanTree->Branch( Form("Rpc%02dEvent", uRpcIndex), "TRpcEvent",
                              &(fRpcInputEvent[uRpcIndex]) );
                     }
               } // if( 1 == fRootConvPar->uMultiEventsAllowed)
                  else for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )
                     if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
                     {
                        /*
                        cleanTree->Branch( Form("Rpc%02dEvent", uRpcIndex),  "Rpc_Event",
                              &(fRpcInputEvent[uRpcIndex]->fEvents[0]) );
                              */
                        cleanTree->Branch( Form("Rpc%02dEvent", uRpcIndex),  "Rpc_Event",
                              &( fRpcSingleEvent[uRpcIndex] ) );
                     } // if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
            } // if( 1 == fRootConvPar->uRpcInTree )
            if( 1 == fRootConvPar->uPlasticsInTree )
            {
               if( 1 == fRootConvPar->uMultiEventsAllowed)
               {
                  for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )
                     if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
                     {
                        cleanTree->Branch( Form("Plastics%02dEvent", uPlasticIndex), "TPlasticsEvent",
                              &(fPlasticsInputEvent[uPlasticIndex]) );
                     }
               } // if( 1 == fRootConvPar->uMultiEventsAllowed)
                  else for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )
                     if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
                     {
                        /*
                        cleanTree->Branch( Form("Plastics%02dEvent", uPlasticIndex),     "Plastics_Event",
                              &(fPlasticsInputEvent[uPlasticIndex]->fEvents[0]) );
                              */
                        cleanTree->Branch( Form("Plastics%02dEvent", uPlasticIndex),     "Plastics_Event",
                              &( fPlasticsSingleEvent[uPlasticIndex] ) );
                     } // if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
            } // if( 1 == fRootConvPar->uRpcInTree )
            /*
            if( 1 == fRootConvPar->uPadRpcInTree)
            {
               cleanTree->Branch("PadRpcEvent",     "TPadRpcEvent", &fPadRpcInputEvent);
            }
            */
            if( 1 == fRootConvPar->uScalOrMuInTree && kTRUE == fParAnalysis->bWithScalersRate )
            {
//               cleanTree->Branch("ScalOrMu1Event",  "TScalersEvent", &fScalOrMu1InputEvent);
               cleanTree->Branch("ScalOrMu1Event",  "Scalers_Event", &(fScalOrMu1InputEvent->fEvent) );
            } // if( 1 == fRootConvPar->uScalOrMuInTree && kTRUE == fParAnalysis->bWithScalersRate )
            if( 1 == fRootConvPar->uScalOrMu2InTree)
            {
//               cleanTree->Branch("ScalOrMu2Event",  "TScalersEvent", &fScalOrMu2InputEvent);
               cleanTree->Branch("ScalOrMu2Event",  "Scalers_Event", &(fScalOrMu2InputEvent->fEvent) );
            }

            if( 1 == fRootConvPar->uMatchingInTree && kTRUE == fParAnalysis->bWithMatching)
            {
               cleanTree->Branch("MatchingEvent",  "TMatchingEvent", &fMatchingInputEvent);
            } // if( 1 == fRootConvPar->uMatchingInTree && kTRUE == fParAnalysis->bWithMatching)
         } // 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 && kTRUE == fParAnalysis->bWithEpics )
         {
            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 && kTRUE == fParAnalysis->bWithEpics )

         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( -1 < fParAnalysis->iTriggerRejection && kTRUE == fParAnalysis->bWithTriglog)
      if( 1 == ((fTriglogInputEvent->fVulomTriggerPattern>>fParAnalysis->iTriggerRejection)&0x1 )  )
      {
          // The trigger Rejection is active and we have rejected trigger
         return;
      }

   if( (1 == fRootConvPar->uCaenInTree && kTRUE == fParAnalysis->bWithMbs ) || 1 == fRootConvPar->uTriglogInputInTree )
      if(!fCrateInputEvent->IsValid() )
         return; // skip incomplete mbs events from first step
   if( 1 == fRootConvPar->uTriglogInputInTree && kTRUE == fParAnalysis->bWithTriglog )
      if( !fTriglogInputEvent->IsValid() )
         return; // skip incomplete triglog events from first step
   if( 1 == fRootConvPar->uVftxInTree && kTRUE == fParAnalysis->bWithVftxCalib )
      if( !fVftxInputEvent->IsValid() )
         cout<<"Invalid VFTX event! event "<<uNbEvents<<endl;

   uNbEvents++;

   // Fill the tree
   if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )
   {
      if( 1 == fRootConvPar->uRpcInTree )
      {
         if( 0 == fRootConvPar->uMultiEventsAllowed )
         {
            for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )
               if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
               {
                  fRpcSingleEvent[uRpcIndex].Clear();
                  fRpcSingleEvent[uRpcIndex] = fRpcInputEvent[uRpcIndex]->fEvents[0];
               } // if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
         } // if( 0 == fRootConvPar->uMultiEventsAllowed )
      } // if( 1 == fRootConvPar->uRpcInTree )

      if( 1 == fRootConvPar->uPlasticsInTree )
      {
         if( 0 == fRootConvPar->uMultiEventsAllowed )
         {
            for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )
               if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
               {
                  fPlasticsSingleEvent[uPlasticIndex].Clear();
                  fPlasticsSingleEvent[uPlasticIndex] = fPlasticsInputEvent[uPlasticIndex]->fEvents[0];
               } // if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
         } // if( 0 == fRootConvPar->uMultiEventsAllowed )
      } // if( 1 == fRootConvPar->uPlasticsInTree )

      if( 1 == fRootConvPar->uGet4InTree && kTRUE == fParAnalysis->bWithGet4 )
      {
         if( 0 == fRootConvPar->uMultiEventsAllowed )
         {
            if( 0 < (fGet4v1InputEvent->fEvents).size())
            {
               for( Int_t iGet4Board = 0; iGet4Board<MAX_GET4; iGet4Board ++)
                  if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
                  {
                     fGet4Board[iGet4Board].Clear();
                     fGet4Board[iGet4Board] = (fGet4v1InputEvent->fEvents[0]).fGet4Boards[iGet4Board];
                  } // if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
            } // if( 0 < (fGet4v1InputEvent->fEvents).size())
            else for( Int_t iGet4Board = 0; iGet4Board<MAX_GET4; iGet4Board ++)
               if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
                  fGet4Board[iGet4Board].Clear();
         } // if( 0 == fRootConvPar->uMultiEventsAllowed )
      } // if( 1 == fRootConvPar->uGet4InTree && kTRUE == fParAnalysis->bWithGet4 )
#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
                  && kTRUE == fParAnalysis->bWithTriglog )
            {
               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 && kTRUE == fParAnalysis->bWithTriglog )
         } // if( 0 < fHodoInputEvent->NumHits() )
      } // if( 1 == fRootConvPar->uHodoInputInTree )
#endif
      if( 1 == fRootConvPar->uSlowControlInTree && kTRUE == fParAnalysis->bWithEpics )
      {
         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 && kTRUE == fParAnalysis->bWithEpics )

      cleanTree->Fill();

      // Case: RAW DATA mode and GET4 enabled => only GET4 data possible!
      // => Multi-event by definition
      if( kTRUE == fParAnalysis->bWithGet4 )
         if( kTRUE == fGet4v1Par->bRawDataMode )
         {
            for( UInt_t uGet4EventIndex = 1; uGet4EventIndex < (fGet4v1InputEvent->fEvents).size(); uGet4EventIndex++ )
            {
               uNbEvents++;
               for( Int_t iGet4Board = 0; iGet4Board<MAX_GET4; iGet4Board ++)
                  if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
                  {
                     fGet4Board[iGet4Board].Clear();
                     fGet4Board[iGet4Board] = (fGet4v1InputEvent->fEvents[uGet4EventIndex]).fGet4Boards[iGet4Board];
                  } // if( 1 == fRootConvPar->uGet4Enabled[iGet4Board])
               for( UInt_t uRpcIndex = 0; uRpcIndex < fParAnalysis->uNbRpc; uRpcIndex++ )
                  if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )
                  {
                     fRpcSingleEvent[uRpcIndex].Clear();
                     fRpcSingleEvent[uRpcIndex] = fRpcInputEvent[uRpcIndex]->fEvents[uGet4EventIndex];
                  } // if( 1 == fRootConvPar->uRpcEnabled[uRpcIndex] )

               for( UInt_t uPlasticIndex = 0; uPlasticIndex < fParAnalysis->uNbPlastics; uPlasticIndex++ )
                  if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )
                  {
                     fPlasticsSingleEvent[uPlasticIndex].Clear();
                     fPlasticsSingleEvent[uPlasticIndex] = fPlasticsInputEvent[uPlasticIndex]->fEvents[uGet4EventIndex];
                  } // if( 1 == fRootConvPar->uPlasticsEnabled[uPlasticIndex] )

               cleanTree->Fill();
            } // for( Int_t iGet4EventIndex = 1; iGet4EventIndex < (fGet4v1InputEvent->fEvents).size(); iGet4EventIndex++ )
         } // if( kTRUE == fParAnalysis->bWithGet4 ) && if( kTRUE == fGet4v1Par->bRawDataMode )
   } // if( ( 1 == fRootConvPar->uCleanOuputEnable || 1 == fRootConvPar->uRootDisplay) )

   fOutputEvent->SetValid(kTRUE);
}