TMeshAnalysis.cxx

Go to the documentation of this file.

// $Id$
//-----------------------------------------------------------------------
//       The GSI Online Offline Object Oriented (Go4) Project
//         Experiment Data Processing at EE department, GSI
//-----------------------------------------------------------------------
// Copyright (C) 2000- GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
//                     Planckstr. 1, 64291 Darmstadt, Germany
// Contact:            http://go4.gsi.de
//-----------------------------------------------------------------------
// This software can be used under the license agreements as stated
// in Go4License.txt file which is part of the distribution.
//-----------------------------------------------------------------------

#include "TMeshAnalysis.h"

#include "TH1.h"

#include "TGo4Version.h"
#include "TMeshParameter.h"
#include "TGo4StepFactory.h"
#include "TGo4AnalysisStep.h"

//***********************************************************
TMeshAnalysis::TMeshAnalysis() :
   TGo4Analysis()
{
   TGo4Log::Error("Wrong constructor TMeshAnalysis()!");
}

//***********************************************************
// this constructor is used
TMeshAnalysis::TMeshAnalysis(int argc, char **argv) :
   TGo4Analysis(argc, argv)
{
   if (!TGo4Version::CheckVersion(__GO4BUILDVERSION__)) {
      TGo4Log::Error("Go4 version mismatch");
      exit(-1);
   }

   TGo4Log::Info("Create TMeshAnalysis %s", GetName());

   SetStepChecking(kFALSE); // necessary for non-subsequent mesh analysis

   // first step definitions:
   TGo4StepFactory* factory1 = new TGo4StepFactory("UnpackFactory");
   factory1->DefOutputEvent("RawEvent", "TMeshRawEvent");
   factory1->DefEventProcessor("Unpacker","TMeshUnpackProc");
   TGo4EventSourceParameter *source1 = new TGo4MbsFileParameter(GetDefaultTestFileName());
   TGo4EventStoreParameter *store1   = new TGo4FileStoreParameter("MeshRawEvent");
   TGo4AnalysisStep *step1  = new TGo4AnalysisStep("Unpack",factory1,source1,store1,nullptr);
   AddAnalysisStep(step1);
   step1->SetSourceEnabled(kTRUE);
   step1->SetStoreEnabled(kFALSE);
   step1->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory2 = new TGo4StepFactory("RawProvider1Factory");
   //  factory2->DefInputEvent("RawEvent", "TMeshRawEvent"); // read full raw event without partial io
   //  factory2->DefEventProcessor("RawEvent_1","TMeshProviderProc"); // processorname must match name of input event + "_"
   factory2->DefInputEvent("RawEvent.fxSub1", "TMeshRawEvent"); // try partial io: eventname must match branch name!
   factory2->DefEventProcessor("RawEvent.fxSub1_1","TMeshProviderProc"); // processorname must match name of input event +"_"
   factory2->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step2 = new TGo4AnalysisStep("Input1Provider",factory2,nullptr,nullptr,nullptr);
   AddAnalysisStep(step2);
   step2->SetSourceEnabled(kFALSE);
   step2->SetStoreEnabled(kFALSE);
   step2->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory3 = new TGo4StepFactory("RawProvider2Factory");
   //  factory3->DefInputEvent("RawEvent", "TMeshRawEvent"); // read full raw event without partial io
   //  factory3->DefEventProcessor("RawEvent_2","TMeshProviderProc"); // name must match name of input event!
   factory3->DefInputEvent("RawEvent.fxSub2", "TMeshRawEvent"); // build full raw event, but fill only subbranch as defined in name
   factory3->DefEventProcessor("RawEvent.fxSub2_1","TMeshProviderProc");
   factory3->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step3 = new TGo4AnalysisStep("Input2Provider",factory3,nullptr,nullptr,nullptr);
   AddAnalysisStep(step3);
   step3->SetSourceEnabled(kFALSE);
   step3->SetStoreEnabled(kFALSE);
   step3->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory4 = new TGo4StepFactory("RawProvider3Factory");
   //  factory4->DefInputEvent("RawEvent", "TMeshRawEvent"); // read full raw event without partial io
   //  factory4->DefEventProcessor("RawEvent_3","TMeshProviderProc"); // name must match name of input event!
   factory4->DefInputEvent("RawEvent.fxSub3", "TMeshRawEvent"); // build full raw event, but fill only subbranch as defined in name
   factory4->DefEventProcessor("RawEvent.fxSub3_1","TMeshProviderProc");
   factory4->DefOutputEvent("Dummy", "TMeshDummyEvent");
   // provider factories without definition of output event and proc
   TGo4AnalysisStep *step4 = new TGo4AnalysisStep("Input3Provider",factory4,nullptr,nullptr,nullptr);
   AddAnalysisStep(step4);
   step4->SetSourceEnabled(kFALSE);
   step4->SetStoreEnabled(kFALSE);
   step4->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory5 = new TGo4StepFactory("Branch1Factory");
   factory5->DefOutputEvent("EventB1", "TMeshB1OutputEvent");
   factory5->DefEventProcessor("Branch1","TMeshB1AnlProc");
   TGo4EventStoreParameter *store5   = new TGo4FileStoreParameter("MeshB1");
   TGo4AnalysisStep *step5 = new TGo4AnalysisStep("Exec1",factory5,nullptr,store5,nullptr);
   // no own event source for execution step
   AddAnalysisStep(step5);
   step5->SetSourceEnabled(kFALSE);
   step5->SetStoreEnabled(kTRUE);
   step5->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory6 = new TGo4StepFactory("Branch2Factory");
   factory6->DefOutputEvent("EventB2", "TMeshB2OutputEvent");
   factory6->DefEventProcessor("Branch2","TMeshB2AnlProc");
   TGo4EventStoreParameter *store6   = new TGo4FileStoreParameter("MeshB2");
   TGo4AnalysisStep *step6  = new TGo4AnalysisStep("Exec2",factory6,nullptr,store6,nullptr);
   // no own event source for execution step
   AddAnalysisStep(step6);
   step6->SetSourceEnabled(kFALSE);
   step6->SetStoreEnabled(kTRUE);
   step6->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory7 = new TGo4StepFactory("Branch3Factory");
   factory7->DefOutputEvent("EventB3", "TMeshB3OutputEvent");
   factory7->DefEventProcessor("Branch3","TMeshB3AnlProc");
   TGo4EventStoreParameter *store7   = new TGo4FileStoreParameter("MeshB3");
   TGo4AnalysisStep *step7 = new TGo4AnalysisStep("Exec3",factory7,nullptr,store7,nullptr);
   // no own event source for execution step
   AddAnalysisStep(step7);
   step7->SetSourceEnabled(kFALSE);
   step7->SetStoreEnabled(kTRUE);
   step7->SetProcessEnabled(kTRUE);


   TGo4StepFactory* factory8 = new TGo4StepFactory("OutProvider1Factory");
   factory8->DefInputEvent("EventB1", "TMeshB1OutputEvent");
   factory8->DefEventProcessor("EventB1_1","TMeshProviderProc"); // name must match name of input event!
   factory8->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step8 = new TGo4AnalysisStep("Output1Provider",factory8,nullptr,nullptr,nullptr);
   AddAnalysisStep(step8);
   step8->SetSourceEnabled(kFALSE);
   step8->SetStoreEnabled(kFALSE);
   step8->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory9 = new TGo4StepFactory("OutProvider2Factory");
   factory9->DefInputEvent("EventB2", "TMeshB2OutputEvent");
   factory9->DefEventProcessor("EventB2_1","TMeshProviderProc"); // name must match name of input event!
   factory9->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step9 = new TGo4AnalysisStep("Output2Provider",factory9,nullptr,nullptr,nullptr);
   AddAnalysisStep(step9);
   step9->SetSourceEnabled(kFALSE);
   step9->SetStoreEnabled(kFALSE);
   step9->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory10 = new TGo4StepFactory("OutProvider3Factory");
   factory10->DefInputEvent("EventB3", "TMeshB3OutputEvent");
   factory10->DefEventProcessor("EventB3_1","TMeshProviderProc"); // name must match name of input event!
   factory10->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step10 = new TGo4AnalysisStep("Output3Provider",factory10,nullptr,nullptr,nullptr);
   AddAnalysisStep(step10);
   step10->SetSourceEnabled(kFALSE);
   step10->SetStoreEnabled(kFALSE);
   step10->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory11 = new TGo4StepFactory("Branch12Factory");
   factory11->DefOutputEvent("EventB12", "TMeshB12OutputEvent");
   factory11->DefEventProcessor("Branch12","TMeshB12AnlProc");
   TGo4EventStoreParameter *store11 = new TGo4FileStoreParameter("MeshB12");
   TGo4AnalysisStep *step11 = new TGo4AnalysisStep("Exec12",factory11,nullptr,store11,nullptr);
   // no own event source for execution step
   AddAnalysisStep(step11);
   step11->SetSourceEnabled(kFALSE);
   step11->SetStoreEnabled(kTRUE);
   step11->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory12 = new TGo4StepFactory("OutProvider12Factory");
   factory12->DefInputEvent("EventB12", "TMeshB12OutputEvent");
   factory12->DefEventProcessor("EventB12_1","TMeshProviderProc"); // name must match name of input event!
   factory12->DefOutputEvent("Dummy", "TMeshDummyEvent");
   TGo4AnalysisStep *step12 = new TGo4AnalysisStep("Output12Provider",factory12,nullptr,nullptr,nullptr);
   AddAnalysisStep(step12);
   step12->SetSourceEnabled(kFALSE);
   step12->SetStoreEnabled(kFALSE);
   step12->SetProcessEnabled(kTRUE);

   TGo4StepFactory* factory13 = new TGo4StepFactory("Branch13Factory");
   factory13->DefOutputEvent("FinalEvent", "TMeshFinalEvent");
   factory13->DefEventProcessor("Collector","TMeshFinalProc");
   TGo4EventStoreParameter *store13 = new TGo4FileStoreParameter("MeshFinal");
   TGo4AnalysisStep *step13 = new TGo4AnalysisStep("Final",factory13,nullptr,store13,nullptr);
   // no own event source for execution step
   AddAnalysisStep(step13);
   step13->SetSourceEnabled(kFALSE);
   step13->SetStoreEnabled(kTRUE);
   step13->SetProcessEnabled(kTRUE);

   // uncomment following line to define custom passwords for analysis server
   // DefineServerPasswords("Meshadmin", "Meshctrl", "Meshview");

   // At this point, autosave file has not yet been read!
   // Therefore parameter values set here will be overwritten
   // if an autosave file is there.
   fPar = new TMeshParameter("MeshSetup");
   AddParameter(fPar);
}

//***********************************************************
TMeshAnalysis::~TMeshAnalysis()
{
   TGo4Log::Info("TMeshAnalysis: Delete");
}

//-----------------------------------------------------------
Int_t TMeshAnalysis::UserPreLoop()
{
   TGo4Log::Info("TMeshAnalysis: PreLoop");
   // we update the pointers to the current event structures here:
   fMbsEvent = dynamic_cast<TGo4MbsEvent *>(GetInputEvent("Unpack"));   // of step "Unpack"
   fEvents = 0;
   fLastEvent = 0;

   // create histogram for UserEventFunc
   // At this point, the histogram has been restored from autosave file if any.
   fSize = (TH1D *)GetHistogram("Eventsize");
   if (!fSize) {
      // no autosave read, create new and register
      fSize = new TH1D("Eventsize", "Event size [b]", 160, 1, 160);
      AddHistogram(fSize);
   }
   // we use a fitter envelope parameters to exchange fit results:

   return 0;
}
//-----------------------------------------------------------
Int_t TMeshAnalysis::UserPostLoop()
{
   TGo4Log::Info("TMeshAnalysis: PostLoop");
   TGo4Log::Info("Last event: %d  Total events: %d", fLastEvent, fEvents);
   if(fMbsEvent)
   {
      // we can check some properties of last event here:
      //fMbsEvent->PrintEvent(); // header and data content

      // fileheader structure:
      fMbsEvent->PrintMbsFileHeader();

      // mbs buffer header structure:
      fMbsEvent->PrintMbsBufferHeader();

   }


   fMbsEvent = nullptr; // reset to avoid invalid pointer if analysis is changed in between
   fEvents = 0;
   return 0;
}

//-----------------------------------------------------------
Int_t TMeshAnalysis::UserEventFunc()
{
   Int_t value = 0;
   Int_t count = 0;
   if(fMbsEvent) value = fMbsEvent->GetDlen()/2+2; // total longwords
   fSize->Fill(value); // fill histogram
   fEvents++;
   if(fEvents == 1 || IsNewInputFile()) {
      if(fMbsEvent) {
         count = fMbsEvent->GetCount();
         TGo4Log::Info("First event #: %d", count);
         // mbs buffer header structure:
         fMbsEvent->PrintMbsBufferHeader();
      }
   }
   fLastEvent = count;
   return 0;
}