HarpoEvent.cxx

Go to the documentation of this file.

// 
// File HarpoEvent.cxx
//
#include "HarpoEvent.h"
#include "HarpoRecoEvent.h"
#include "HarpoDccEvent.h"
#include "HarpoSimEvent.h"
#include <cstdlib>
#include <cstring>
#include <iostream>
// for swab 
#include <unistd.h> 

ClassImp(HarpoEvent)

HarpoEvent::HarpoEvent()
{
  //  std::cout << "+++ New Event " << this << std::endl;
  if(gHarpoDebug>1)
    Info("HarpoEvent"," ");
  fEventHeader = NULL;
  fDetEvent = new HarpoDetEvent *[gkNDetectors];  
  UInt_t nd;
  for(nd=0;nd<gkNDetectors;nd++) {
    fDetEvent[nd] = (HarpoDetEvent *) NULL;
  }
  fRecoEvent = (HarpoRecoEvent *)NULL;//new HarpoRecoEvent();
}


HarpoEvent::HarpoEvent(const HarpoEvent &evt): TObject(evt),
  fEventHeader(evt.fEventHeader), fDetEvent(evt.fDetEvent), fRecoEvent(evt.fRecoEvent)
{
}

 HarpoEvent::~HarpoEvent()
{
  //std::cout << "=== Del Event " << this << std::endl;
  UInt_t nd;
  for(nd=0;nd<gkNDetectors;nd++) {
    if(gHDetSet->isExist(nd)) {
      if (fDetEvent[nd] != NULL ){
        if(gHarpoDebug>1)
          Info("~HarpoEvent","Delete event det %d (%p)",nd,fDetEvent[nd]);
        fDetEvent[nd]->Delete();
        //delete fDetEvent[nd];
        //      fDetEvent[nd]->ResetData();
      } else {
        std::cout << "---Det " << nd << " data == NULL" << std::endl;
      }
    }
  }
  delete [] fDetEvent;
  if(fRecoEvent!=NULL) 
    fRecoEvent->Delete();
  //delete fRecoEvent;
  if (fEventHeader != NULL) delete fEventHeader;
}

void          HarpoEvent::SetHeader(HarpoEventHeader *data)
{
  fEventHeader = data;
}

HarpoEventHeader *HarpoEvent::GetHeader()
{
 return fEventHeader;
}

void HarpoEvent::SetDetEvent(HarpoDetEvent *event,ULong_t plane) {
  if ( /* (plane < 0) || */ (plane >= gkNDetectors)) return; 
  if(gHarpoDebug>0) Info("SetDetEvent","%ld",plane);
  gHDetSet->Set(plane);   //Mark detectors as existing
  fDetEvent[plane] = event;
  fDetEvent[plane]->SetDet(plane);
}

HarpoDetEvent *HarpoEvent::GetDetEvent(Long_t plane) {
  //  if ((plane < 0) || (plane >= gkNDetectors)) return NULL; 
  if (gHDetSet->isExist(plane))
    {
       return fDetEvent[plane];
    } else {
    return NULL;
  }
}

// void          HarpoEvent::SetDccMap(HarpoDccMap *data,Long_t plane) {
//   if ((plane < 0) || (plane >= gkNDetectors)) return; 
//   fEventData[plane] = data;
// }

HarpoDccMap   *HarpoEvent::GetDccMap(Long_t plane) {
  if(gHarpoDebug>2) 
    Info("GetDccMap","%li",plane);
  // if(gHDetSet->isSim()  && fDetEvent[SIMDET] != NULL)
  //     return ((HarpoSimEvent*)fDetEvent[SIMDET])->GetMap(plane);
  if (gHDetSet->isExist(plane) && fDetEvent[plane] != NULL){
    if(plane == XDCC)
      return ((HarpoDccEvent*)fDetEvent[plane])->GetMap();
    if(plane == YDCC)
      return ((HarpoDccEvent*)fDetEvent[plane])->GetMap();
  }// else {
    return NULL;
    //}
}

void          HarpoEvent::print()
{
  //  std::cout << "---in Event print " << fEventHeader << std::endl;
  if ( fEventHeader != NULL ) fEventHeader->print();
  UInt_t nd;
  for(nd=0;nd<gkNDetectors;nd++) {
    if(gHDetSet->isExist(nd)) {
      if (fDetEvent[nd] != NULL ){
        fDetEvent[nd]->print();
      } else {
        std::cout << "---Det " << nd << " data == NULL" << std::endl;
      }
    }
  }
}

// Bool_t      HarpoEvent::IsZeroSuppressed(Long_t plane)
// {
//   if (gHDetSet->isExist(plane))
//     {
//       return fDetEvent[plane]->IsZeroSuppressed();
//     } else {
//     return false;
//   }
// }

ULong_t      HarpoEvent::GetTimeStamp(Long_t plane)
{

  if(gHDetSet->isSim())
    return 0;//fDetEvent[SIMDET]->GetTimeStamp();
  if (gHDetSet->isExist(plane))
    {
      if(gHarpoDebug>1) 
        std::cout << plane << " => " << fDetEvent[plane] << std::endl;
      return fDetEvent[plane]->GetTimeStamp();
    } else {
    return 0;
  }
}


void      HarpoEvent::ResetEvent()
{

  if (gHDetSet->isExist(SIMDET)) {
    fDetEvent[SIMDET]->Delete();
  }
  fRecoEvent->Delete();

}