heventpr.cxx

Go to the documentation of this file.

/* 
 * Look for Event by number and print event map if found
 */
#include <cstdlib>
#include <cstdio>
#include <iostream>

#include "HarpoConfig.h"
#include "HarpoReader.h"
#include "HarpoRunHeader.h"
#include "HarpoEvent.h"
#include "Pmm2Event.h"

int main(int argc, char **argv)
{
  std::cout << "Starting" << std::endl;

  Bool_t found  = false; 
  HarpoRunHeader *hdr;

  HarpoConfig *cfg = new HarpoConfig(argc,argv);
  cfg->Init(); // Last call before Get.., Reading Config file not hire

  cfg->print();

  //exit(0);

  HarpoReader *rdr = new HarpoReader( cfg );

  if ( rdr->Init() ) {
    hdr = rdr->GetRunHeader();
    hdr->print();
  }
  else
    {
      printf ("Int ERROR exiting ...\n");
      return 1;
    }

  //if ( rdr->NextEvent() ) printf("\n Next Event OK\n");
  std::cout << "Finding event " << cfg->GetNEvent() << std::endl;
  found = rdr->Find(cfg->GetNEvent());
  if (found) {
    std::cout << "Event " << cfg->GetNEvent() << " found." << std::endl;
    HarpoEvent *evt = rdr->GetEvent();
    evt->print();
    
    ULong_t ndet;
    for(ndet=0;ndet<gkNDetectors;ndet++) {
      if (gHDetSet->isExist(ndet)) {
        std::cout << "Found " << HarpoDet(ndet).GetDescription()
                  << " data for event " << cfg->GetNEvent() << std::endl;
        evt->GetDetEvent(ndet)->print();
        if ((ndet == XDCC) || ( ndet == YDCC )) { 
        HarpoDccMap * m = evt->GetDccMap(ndet);
        if ( m != NULL ) {
          m->print();
        } else {
          std::cout << "No Data Map for plane "<< ndet << std::endl;
        }
        } else if ( ndet == PMM2 ) {
          // Pmm2Event *anaEvt =static_cast<Pmm2Event *>(evt->GetDetEvent(ndet));
          // anaEvt->print();
        } else {
          // std::cout << "Found " << HarpoDet(ndet).GetDescription()
          //        << " data for event " << cfg->GetNEvent() << std::endl;
        }

      }
    }
  }
  else
    {
      std::cout << "Event " << cfg->GetNEvent() << " no found." << std::endl;
    }
  
  std::cout << " Done after event " << rdr->GetReadEvts()
            << " ( EOF ) : " << rdr->isDone()
            << " Err " << rdr->isError() << std::endl; 
  return 0;
}