beamtime/cern-gem12/go4/Gem1/TGem1Proc.cxx (r4864/r3571)

Go to the documentation of this file.

/* Generated by Together */

#include "TGem1Proc.h"

int xGSI[2][128] = {   {8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,8,6,4,2,9,11,13,15,
                13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6,13,15,9,11,4,2,8,6},
                {7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,7,5,3,1,10,12,14,16,
                14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5,14,16,10,12,3,1,7,5}  };
int yGSI[2][128] = {   {16,16,16,16,16,16,16,16,15,15,15,15,15,15,15,15,14,14,14,14,14,14,14,14,13,13,13,13,13,13,13,13,12,12,12,12,12,12,12,12,11,11,11,11,11,11,11,11,
                10,10,10,10,10,10,10,10,9,9,9,9,9,9,9,9,8,8,8,8,8,8,8,8,7,7,7,7,7,7,7,7,6,6,6,6,6,6,6,6,5,5,5,5,5,5,5,5,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,
                1,1,1,1,1,1,1,1},
                {16,16,16,16,16,16,16,16,15,15,15,15,15,15,15,15,14,14,14,14,14,14,14,14,13,13,13,13,13,13,13,13,12,12,12,12,12,12,12,12,11,11,11,11,11,11,11,11,
                10,10,10,10,10,10,10,10,9,9,9,9,9,9,9,9,8,8,8,8,8,8,8,8,7,7,7,7,7,7,7,7,6,6,6,6,6,6,6,6,5,5,5,5,5,5,5,5,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,
                1,1,1,1,1,1,1,1}  };



TGem1Proc::TGem1Proc(const char* name) : TCBMBeamtimeProc(name),
                                                                         fOutputEvent(0)
{
  cout << "**** TGem1ProcTTRDSpadicProc: Create instance " << name << endl;

  fPar = (TGem1Param*) MakeParameter("Gem1Par", "TGem1Param", "set_Gem1Par.C");

  fMappingGSI_GEM1 = MakeTH2('I', "Gem1/MappingGSI_GEM1", "MAPPING GSI for station 1 ", 20, 0., 19., 20, 0., 19., "Mapping_X", "Mapping_Y");
  fMappingGSI_GEM2 = MakeTH2('I', "Gem1/MappingGSI_GEM2", "MAPPING GSI for station 2 ", 20, 0., 19., 20, 0., 19., "Mapping_X", "Mapping_Y");
  fMappingGSI_GEM3 = MakeTH2('I', "Gem1/MappingGSI_GEM3", "MAPPING GSI for station 3 ", 20, 0., 19., 20, 0., 19., "Mapping_X", "Mapping_Y");

  fRocInputEvent = 0;
  fOutputEvent = 0;

}

TGem1Proc::~TGem1Proc()
{
}

void TGem1Proc::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(fRocInputEvent==0)
   {
      TCBMBeamtimeEvent* btevent=dynamic_cast<TCBMBeamtimeEvent*>(GetInputEvent());
      if(btevent)
      {
         fRocInputEvent = dynamic_cast<TRocEvent*>(btevent->GetSubEvent("ROC"));
      }
      else
      {
         fRocInputEvent = dynamic_cast<TRocEvent*>(GetInputEvent());
      }

      if (fRocInputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TRICHProc: Fatal error: input event is not a TRocEvent!!! 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<TGem1Event*>(btevent->GetSubEvent("Gem1"));
      }
      else
      {
         fOutputEvent= dynamic_cast<TGem1Event*>(outevnt);
      }
      if(fOutputEvent==0) {
         GO4_STOP_ANALYSIS_MESSAGE("**** TGem1Proc: Fatal error: output event is not a TGem1Event!!! STOP GO4");
      }
   }
}

void TGem1Proc::FinalizeEvent()
{

   for (int nn=0; nn<3; nn++) {
     // cout << "  Processing ROC "<< RICHROC << endl;
      // Loop over all ROCs which deliver RICH data

      int rocid(0);
      if (nn == 0) rocid = fPar->fRoc1; else
      if (nn == 1) rocid = fPar->fRoc2; else
      if (nn == 2) rocid = fPar->fRoc3;

      TRocData* RocData=dynamic_cast<TRocData*>(fRocInputEvent->getEventElement(rocid));
      if (RocData==0) continue;  //no data from this ROC available;

      for (unsigned MessageNr=0; MessageNr<RocData->fExtMessages.size(); MessageNr++) {
         //cout << "  Processing message " << MessageNr << "from ROC " << RICHROC << endl;
         //cout << "Msg.size() = " << RocData->fExtMessages.size() << endl;
         // Loop over all messages from particular RICH-ROC

         TRocMessageExtended msg = RocData->fExtMessages.at(MessageNr);
         if (msg.GetMessageType() != roc::MSG_HIT) continue;  // only interrested in HIT messages
         //cout << "roc::MSG_HIT" << endl;

         // extraxt message information
         UInt_t dataROC=msg.GetRocNumber();
         UInt_t nxid = msg.GetNxNumber();  // 0 or 2
         UInt_t nxch=msg.GetNxChNum();   // 0..127
         Double_t dataHitTime=msg.GetTriggerDeltaT();
         Int_t dataHitADCraw=msg.GetNxADC();
         Int_t dataHitADCcor=msg.GetCorrectedNxADC();

         int z;

         if(rocid == 2)
         {
                 if(nxid ==0)
                 {
                         z = 1;
                 }
                 else if(nxid ==2)
                 {
                         z =0;
                 }
                 fMappingGSI_GEM1->Fill((17-xGSI[z][nxch]),yGSI[z][nxch]);
         }
         if(rocid == 3)
         {
                 if(nxid ==0)
                 {
                         z = 1;
                 }
                 else if(nxid ==2)
                 {
                         z =0;
                 }
                 fMappingGSI_GEM2->Fill((17-xGSI[z][nxch]),yGSI[z][nxch]);
         }
         if(rocid == 4)
         {
                 if(nxid ==0)
                 {
                         z = 1;
                 }
                 else if(nxid ==2)
                 {
                         z =0;
                 }
                 fMappingGSI_GEM3->Fill((17-xGSI[z][nxch]),yGSI[z][nxch]);
         }

      }
   }

}