HarpoReaderSim.cxx

Go to the documentation of this file.

// 
// File HarpoReaderSim.cxx
//
#include "HarpoReaderSim.h"
#include "HarpoDetReader.h"
#include "HarpoConfig.h"
#include "HarpoRunHeader.h"
#include "HarpoEvent.h"
#include "HarpoTimeStamp.h"

#include "HarpoDccReader.h"
#include "HarpoSimReader.h"
#include "HarpoFeminosHeader.h"

#include <climits>
#include <cstdio>
#include <arpa/inet.h>
#include <cstdlib>
//#include <cstring>
#include <iostream>


ClassImp(HarpoReaderSim)


HarpoReaderSim::HarpoReaderSim() : HarpoReaderBase()
{
  if(gHarpoDebug>0) Info("HarpoReaderSim()"," ");

  fMode = hSyncTime;
}


HarpoReaderSim::HarpoReaderSim(HarpoConfig *cfg) : HarpoReaderBase(cfg)
 {
  if(gHarpoDebug>0) Info("HarpoReaderSim(cfg)"," ");

  fMode = hSyncTime;
  fRunHeader = new HarpoRunHeader(cfg);

}

HarpoReaderSim::~HarpoReaderSim()
{
  if(gHarpoDebug>0) Info("~HarpoReaderSim()"," ");
  // ??? Delete Det Constructors
  // if ( fNewFile != NULL) {
  //   fNewFile->Close();
  // }
}

void HarpoReaderSim::SetDetFile(TString * /*name */,Int_t /* plane */)
{
  //  fInFile[plane] = name; //???
}

 void HarpoReaderSim::SetRootFile(TString *name)
 {
   //if(gHarpoDebug>0) 
   Info("Dummy SetRootFile","%s",name->Data());
 }
  
Bool_t HarpoReaderSim::NextEvent()
{
  if(gHarpoDebug>1) Info("NextEvent"," ");
  Bool_t res = true;
  //ULong_t ndet;
  //Long_t deltat;
  //ULong_t ts[2],tso[2];
 
  res = res && fDetReader[SIMDET]->NextEvent();
  if ( ! res ) return res; // No more event or errors

  fNevts++;
  GetEvent();

  return res;
}

HarpoEvent *HarpoReaderSim::GetEvent()
{
  if(gHarpoDebug>0) Info("GetEvent","event %ld",fNevts);
  HarpoDetEvent *revt;
  ULong_t ndet;

  /* if(fEvent==NULL) */ fEvent = new HarpoEvent();
  /* if(fEventHeader==NULL) */ fEventHeader = new HarpoEventHeader();

  fEventHeader->SetEvtNo(fNevts);  // Event header
  fEvent->SetHeader(fEventHeader);

  for(ndet=0;ndet<gkNDetectors;ndet++) {
    if (fDetSet->isExist(ndet)) {
      revt = ((HarpoSimReader*)fDetReader[SIMDET])->GetEvent(ndet);
      fEvent->SetDetEvent(revt,ndet);
    }
  }
  
  return fEvent;
}

HarpoRunHeader *HarpoReaderSim::GetRunHeader()
{
  if(gHarpoDebug>0) Info("GetRunHeader"," ");
  return fRunHeader;
}


Bool_t HarpoReaderSim::Init(hReadMode mode)
{
  if(gHarpoDebug>0) Info("Init","Readout mode %d",mode);
  Bool_t res = true;
  fMode = mode;
  fDetSet->print();
  
  std::cout << "Det is SIM " << fDetSet->isExist(SIMDET) << std::endl;
  fDetReader[SIMDET] = new HarpoSimReader();
  fDetReader[SIMDET]->SetDet(SIMDET);
  fDetReader[SIMDET]->SetRunNo(fRunNo);
  fDetReader[SIMDET]->SetDetFile(gHConfig->GetInFile(SIMDET));
  res = res && fDetReader[SIMDET]->Init();

  HarpoDetHeader *rawhdr;
  // ULong_t ndet;

  // for(ndet=0;ndet<gkNDetectors;ndet++) {
  //   if (fDetSet->isSim()) {
  //     rawhdr = fDetReader[SIMDET]->GetDetHeader();
  //     fRunHeader->SetDetHeader(rawhdr,ndet);
  //     continue;
  //   }
  // }

  // Save det mask
  fDetSet->Set(SIMDET);  // always set for simulation
  fRunHeader->SetSimDataFlag(true);
  fRunHeader->SetDetMask(fDetSet->GetMask());

  for(UInt_t ndet=0;ndet<gkNDetectors;ndet++) {
    //  std::cout << "DET : " << ndet << "/" << gkNDetectors << " => " << fDetSet->isExist(ndet)<< std::endl;
    if (!fDetSet->isExist(ndet)) continue;
    rawhdr = ((HarpoSimReader*)fDetReader[SIMDET])->GetDetHeader(ndet);
    if(gHarpoDebug>0) {
      std::cout << "Add Det Header "<< ndet << " " << rawhdr << std::endl;
      rawhdr->print();
    }
    fRunHeader->SetDetHeader(rawhdr,ndet);
  }
  //  std::cout << res << std::endl;
  return res;
}

Long_t      HarpoReaderSim::Loop(Long_t maxevents)
{
  Long_t cnt=0;
  if(gHarpoDebug>0) Info("Loop", "maxevents = %ld",maxevents);

  while (NextEvent()) {
    if(gHarpoDebug>1) 
      Info("Loop", "Get Event %ld",fNevts);
    HarpoEvent* evt = GetEvent();
    if(gHarpoDebug>1) 
      Info("Loop", "Process event %p",evt);

    // Pre
    Int_t st = hAnaOK;
    fEvent->GetHeader()->SetEvAnaStatus(st);

    for(Int_t i = 0; i<fNanalyses; i++)
      {
        fAnalyse[i]->ProcessEvent(fEvent);
        st = fEvent->GetHeader()->GetEvAnaStatus();
        if ( st < hAnaNextEvent ) break;
      }

    if(gHarpoDebug>1) 
      Info("Loop", "Event processed");
   // if (gHDetSet->isExist(SIMDET)){
   //   HarpoSimEvent *simEvt = (HarpoSimEvent *)(evt->GetDetEvent(SIMDET));
   //   Int_t nIonTr = simEvt->GetNtracks();
   // }
   evt->Delete();
   if(gHarpoDebug>1) 
     Info("Loop", "Event Deleted");

   if (st == hAnaNextFile )
     {
       //how to go to next file ?
       break;
     }
   else if (st == hAnaStop)
     {
       break;
     }
   else if (st <= hAnaAbort)
     {
       Error(__FUNCTION__,
             "Fatal Analysis Error %d aborting prgram.\n",st);
       std::abort();
     }

   if (++cnt >= maxevents) break;
  }
  return cnt;  
}

Long_t      HarpoReaderSim::Find(Long_t /* eventNo */ )
{
  //todo Find
  return true;
}

void      HarpoReaderSim::SkipEvent()
{
  //todo Skip
}