beamtime/cern-oct11/go4/DUBNA/TTRDDubnaProc.cxx (r4864/r2480)

Go to the documentation of this file.

#include "TTRDDubnaProc.h"

#include "TMbsCrateParam.h"

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

   fPar = (TTRDDubnaParam*) MakeParameter("TRDDubnaPar", "TTRDDubnaParam");

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

   TString obname, obtitle;   

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

      int amad = mad - 2;

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

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

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

      obname.Form("MADC%d_pion_Dubna", 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[amad][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_Dubna", 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[amad][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("TTRDDubnaProc Histograms created \n");
   fflush ( stdout);
}


TTRDDubnaProc::~TTRDDubnaProc()
{
}

void TTRDDubnaProc::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("**** TTRDDubnaProc: 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<TTRDDubnaEvent*>(btevent->GetSubEvent("TRDDUBNA"));
         // 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<TTRDDubnaEvent*>(outevnt);
      }
      if(fOutputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TTRDDubnaProc: Fatal error: output event is not a TTRDDubnaEvent!!! STOP GO4");
      }
   }
}

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

     for (int mad = 2; mad < MAX_MADC; mad++) {
         int amad = mad - 2;
         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[amad][ch]->Fill(adc);
               fOutputEvent->fMadc_Pion[amad][ch] = adc;
            }
            if (fBeamEvent->fIsElectron)
            {
               fMadc_electron[amad][ch]->Fill(adc);
               fOutputEvent->fMadc_Electron[amad][ch] = adc;
            }
         } // for ch
      }// for mad
   } // if BeamEvent
}