beamtime/cosy-jan12/go4/TBeamMonitorProc.cxx (r4864/r3159)

Go to the documentation of this file.

#include "TBeamMonitorProc.h"
#include "TGo4Version.h"
#if __GO4BUILDVERSION__ > 40502
#include "go4iostream.h"
#else
#include "Riostream.h"
#endif
#include "TTimeStamp.h"
#include "TSystem.h"
#include "TROOT.h"
#include "TGo4Log.h"

#include "TGo4UserException.h"
#include "TGo4MbsEvent.h"
#include "roc/Message.h"
#include "roc/Board.h"


TBeamMonitorProc::TBeamMonitorProc(const char* name) :
TCBMBeamtimeProc(name),
fCrateInputEvent(0),
fRocInputEvent(0),
fEpicsInputEvent(0),
fOutputEvent(0)
{
   cout << "**** TBeamMonitorProc: Create instance " << name << endl;
   TString obname, obtitle;

   fPar = (TBeamMonitorParam*) MakeParameter("BeamMonitorPar", "TBeamMonitorParam");

   TString setupmacro = "set_BeamMonitorPar.C";
   if (!gSystem->AccessPathName(setupmacro.Data())) {
      TGo4Log::Info("Executing beam monitor condition setup script %s", setupmacro.Data());
      gROOT->ProcessLine(Form(".x %s", setupmacro.Data()));

   } else {
      TGo4Log::Info("NO beam monitor  setup script %s. Use previous values!", setupmacro.Data());
   }

   // fCh_Pb = MakeTH2('I', "BeamMonitor/Ch1_Pb", "Cherenkov - Pb", 512, 0., 4096., 512, 0., 4096., "Cherenkov1", "Pb");

   fPmtDiff = MakeTH1('I', "BeamMonitor/PmtDiff", "Time difference between TDC[5][2] and TDC[5][3]", 4096, -2048, +2048, "time");

   // following is for slow control display:
   //   fMotorPos = MakeTH2('I', Form("EPICS/MotorPosition"), Form("GEM Motor"),
   //         60, 0, 300, 60, 0, 300, "X (mm)", "Y (mm)");
   //   fRoluPos = MakeTH2('I', Form("EPICS/ROLUPosition"), Form("Rolu positions"),
   //         35, 0, 35, 35, 0, 35, "X (mm)", "Y (mm)");

   printf("TBeamMonitorProc Histograms created \n");
   fflush ( stdout);
}


TBeamMonitorProc::~TBeamMonitorProc()
{
}



void TBeamMonitorProc::InitEvent(TGo4EventElement* outevnt)
{
   // first assign input event:
   // since input event object is never discarded within processor lifetime,
   // we just search for subevent by name once to speed up processing
   if(fCrateInputEvent==0)
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent());
      if(btevent)
      {
         fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(btevent->GetSubEvent("MBSCRATE"));
         fRocInputEvent=dynamic_cast<TRocEvent*>(btevent->GetSubEvent("ROC"));
         fEpicsInputEvent=dynamic_cast<TEpicsEvent*>(btevent->GetSubEvent("EPICS"));
      }
      else
      {
         fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(GetInputEvent());

      }
      if(fCrateInputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TBeamMonitorProc: Fatal error: no input event TMbsCrateEvent!!! STOP GO4");
      }
      if(fRocInputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TBeamMonitorProc: Fatal error: no input event TRocEvent!!! STOP GO4");
      }
      if(fEpicsInputEvent==0) {
         Message(2,"**** TBeamMonitorProc: could not find TEpicsEvent!"); // but do not stop...
      }

   }
   // 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<TBeamMonitorEvent*>(btevent->GetSubEvent("BEAM"));
      }
      else
      {
         fOutputEvent= dynamic_cast<TBeamMonitorEvent*>(outevnt);
      }
      if(fOutputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TBeamMonitorProc: Fatal error: output event is not a TBeamMonitorEvent!!! STOP GO4");
      }
   }
}



void TBeamMonitorProc::FinalizeEvent()
{
   if(fCrateInputEvent && fCrateInputEvent->IsValid())
   {
      fPmtDiff->Fill(fCrateInputEvent->fMtdc[5].lead[2] - fCrateInputEvent->fMtdc[5].lead[3]);

      // fCh_Pb->Fill(fCrateInputEvent->fData1182[1][0], fCrateInputEvent->fData1182[1][5]);

      // fOutputEvent->fIsPion = fPionCond->Test(fCrateInputEvent->fData1182[1][0], fCrateInputEvent->fData1182[1][5]);
   }

   if( 0 && fEpicsInputEvent && fEpicsInputEvent->IsValid()) // check if new update has been done
   {
      // TODO: put reasonable display of interesting slow control variables here



      // the following is example from cosy december 2010 beamtime
      // note access to EPICS process variables by full name here
      /*     int motorx = fEpicsInputEvent->GetDouble("CBM:MOTOR:GETPOSX");
     int motory = fEpicsInputEvent->GetDouble("CBM:MOTOR:GETPOSY");
     int rolu_left = fEpicsInputEvent->GetDouble("CBM:rolu:step_to_mm_L");
     int rolu_right = fEpicsInputEvent->GetDouble("CBM:rolu:step_to_mm_R");
     int rolu_up = fEpicsInputEvent->GetDouble("CBM:rolu:step_to_mm_O");
     int rolu_down = fEpicsInputEvent->GetDouble("CBM:rolu:step_to_mm_U");

     fMotorPos->SetTitle(
           Form("Motor position, EPICS Time:%s -",
                 fEpicsInputEvent->GetUpdateTimeString()));
     fMotorPos->Fill(motorx, motory);

     fRoluPos->SetTitle(
           Form("ROLU positions, EPICS Time:%s -",
                 fEpicsInputEvent->GetUpdateTimeString()));
     Int_t roludim = 35;
     Int_t middle = roludim / 2;
     Int_t displaywidth = 6;
     for (int i = 0; i < 2 * displaywidth; ++i) {
        fRoluPos->Fill(middle - rolu_left, middle - displaywidth + i);
        fRoluPos->Fill(middle + rolu_right, middle - displaywidth + i);
        fRoluPos->Fill(middle - displaywidth + i, middle - rolu_down);
        fRoluPos->Fill(middle - displaywidth + i, rolu_up + middle);
     }
       */

   } // end EPICS block
}