beamtime/cern-oct11/go4/BUCHAREST/TTRDBucharestProc.cxx (r4864/r2480)

Go to the documentation of this file.

#include "TTRDBucharestProc.h"

#include "TMbsCrateParam.h"

TTRDBucharestProc::TTRDBucharestProc(const char* name) : TCBMBeamtimeProc(name),fCrateInputEvent(0), fOutputEvent(0),fBeamEvent(0)
{
   cout << "**** TTRDBucharestProc: Create instance" << name << endl;

   fPar = (TTRDBucharestParam*) MakeParameter("TRDBucharestPar", "TTRDBucharestParam");

   TMbsCrateParam* mbspar = dynamic_cast<TMbsCrateParam*> (GetParameter("MbsCratePar"));

   TString obname, obtitle;   

   if (mbspar && !mbspar->fBucharestSetup)

   /*********** MADCs *************/
   for (int mad = 0; mad < 2; mad++) {


      for (int ch = 0; ch < N_MADC_CHA; ch++) {

         obname.Form("BUCHAREST/MADC%d/MADC%d_ch%d_pion", mad,mad,ch);
         obtitle.Form("MADC %d Channel %d pion", mad, ch);
         fMadc_pion[mad][ch] = MakeTH1('I', obname.Data(), obtitle.Data(), 0x3FF, 0, 0x3FF);

         obname.Form("BUCHAREST/MADC%d/MADC%d_ch%d_electron",mad,mad,ch);
         obtitle.Form("MADC %d Channel %d electron", mad, ch);
         fMadc_electron[mad][ch] = MakeTH1('I', obname.Data(), obtitle.Data(), 0x3FF, 0, 0x3FF);

      }


      obname.Form("MADC%d_pion_Buch", mad);
      if (GetPicture(obname.Data()) == 0) {
         obtitle.Form("Pion channels overview MADC %2d", mad);
         TGo4Picture* pic = new TGo4Picture(obname.Data(), obtitle.Data());
         pic->SetDivision(8, 4);
         for(int y=0; y<8;++y)
            for(int x=0;x<4;++x) {
               pic->Pic(y, x)->AddObject(fMadc_pion[mad][4*y+x]);
               pic->Pic(y, x)->SetFillAtt(5, 3001);
               pic->Pic(y, x)->SetRangeX(0, 200);
               pic->Pic(y, x)->SetLogScale(1, 1);
            }
         AddPicture(pic);
      }


      obname.Form("MADC%d_electron_Buch", mad);
      if (GetPicture(obname.Data()) == 0) {
         obtitle.Form("Electron channels overview MADC %2d", mad);
         TGo4Picture* pic = new TGo4Picture(obname.Data(), obtitle.Data());
         pic->SetDivision(8, 4);
         for(int y=0; y<8;++y)
            for(int x=0;x<4;++x) {
               pic->Pic(y, x)->AddObject(fMadc_electron[mad][4*y+x]);
               pic->Pic(y, x)->SetFillAtt(5, 3001);
               pic->Pic(y, x)->SetRangeX(0, 200);
               pic->Pic(y, x)->SetLogScale(1, 1);
            }
         AddPicture(pic);
      }
   }
   printf("TTRDBucharestProc Histograms created \n");
   fflush ( stdout);
}


TTRDBucharestProc::~TTRDBucharestProc()
{
}

void TTRDBucharestProc::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"));
      }
      else
      {
         fCrateInputEvent=dynamic_cast<TMbsCrateEvent*>(GetInputEvent());

      }
      if(fCrateInputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TTRDBucharestProc: Fatal error: input event is not a TMbsCrateEvent!!! STOP GO4");
      }
   }
   // 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<TTRDBucharestEvent*>(btevent->GetSubEvent("TRDBUCHAREST"));
         // since output event object is never discarded within processor lifetime,
         // we just search for subevent by name once to speed up processing
         if(fBeamEvent==0)
            fBeamEvent=dynamic_cast<TBeamMonitorEvent*>(btevent->GetSubEvent("BEAM"));
      }
      else
      {

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

void TTRDBucharestProc::FinalizeEvent()
{
   // correlate with beam only if we are in combined mode
   if (fBeamEvent && fCrateInputEvent->IsValid())
   {

      for (int mad = 0; mad < 2; mad++) {
         for (int ch = 0; ch < N_MADC_CHA; ch++) {
            UShort_t adc=fCrateInputEvent->fMadc[mad][ch];
            if(adc==0) continue; // skip not filled channels

            if (fBeamEvent->IsPulser()) {
               // to be done ...
            }


            if (fBeamEvent->fIsPion)
            {
               fMadc_pion[mad][ch]->Fill(adc);
               fOutputEvent->fMadc_Pion[mad][ch] = adc;
            }
            if (fBeamEvent->fIsElectron)
            {
               fMadc_electron[mad][ch]->Fill(adc);
               fOutputEvent->fMadc_Electron[mad][ch] = adc;
            }
         } // for ch
      }// for mad
   } // if BeamEvent
}