pmm2print.cxx

Go to the documentation of this file.

/* 
 * Print PMM2 DetEvents from PMM2v2DAQ raw file.
 * The option --from --to allow select printing region. 
 */
#include <cstdlib>
#include <cstdio>
#include <iostream>

#include "HarpoConfig.h"
#include "Pmm2Reader.h"
#include "Pmm2Event.h"

using namespace std;

int main(int argc, char **argv)
{
  Bool_t found  = false; 
  Pmm2Header *hdr;
  int rdformat = 1;

  HarpoConfig::gHReaderType = hDummyReader; //no high level reader

  HarpoConfig *cfg = new HarpoConfig(argc,argv,PMM2);
  cfg->Init();  // read configuration file
  int verbose = cfg->GetVerbose();
  if(verbose) cfg->print();
  if(verbose) gHDetSet->print();
  rdformat = cfg->GetDetRdrType(PMM2);
  if(verbose) cout << endl << "Pmm2 reader " << rdformat  << endl;

  Pmm2Reader *rdr = new Pmm2Reader( cfg );
  //  cout << "cfg Det " << cfg->GetDet() << " rdr Det " << rdr->GetDet() << endl;
  //rdr->SetDet(cfg->GetDet()); // FIXED

  rdr->SetReadOutMode(rdformat); 

  if ( rdr->Init() ) {
    hdr = rdr->GetDetHeader();
    if(verbose) hdr->print();
  }
  else
    {
      printf ("Int ERROR exiting ...\n");
      return 1;
    }
  if(verbose) printf("\nLooking for event %lld ...\n",cfg->GetFromEvent());
  
  found = rdr->Find(cfg->GetFromEvent());

  if (found) {
    Pmm2Status *st = hdr->getStatus();
    if (st != NULL) {
      if (verbose) st->print();
      if(verbose > 1) st->Dump();
    }
    do {
      Pmm2Event *evt = rdr->GetEvent();
      if(verbose)evt->print();
      if(verbose > 1)cout << " ** Start Event " << endl;

      Pmm2MesVect * mlist = evt->GetMesurements();

      if (mlist != NULL ) {
        // Print all 212 mesurements
        mlist->print();
      } else {
        std::cout << "--- Decoding error--" << std::endl;
      }
      if(verbose>1)std::cout << "** Done Event " << rdr->GetReadEvts() << std::endl; 
    } while (rdr->NextEvent() && ((rdr->GetReadEvts() -1) <= cfg->GetToEvent()));
  }
  else
    {
      std::cout << "Event " << cfg->GetFromEvent() << " no found." << std::endl;
    }
  
  if(verbose) std::cout << " Done after event " << rdr->GetReadEvts()
            << " ( EOF ) : " << rdr->isDone()
            << " Err " << rdr->isError() << std::endl;

  // valgrind happy with this "delete" ? Yes.  
  delete rdr;
  delete cfg;
  return 0;
}