hrootwriter.cxx

Go to the documentation of this file.

//
// Demo code for harpo root tree writing
//

#include <stdlib.h>
#include <stdio.h>

#include "TDirectory.h"
#include "HarpoConfig.h"
#include "HarpoReader.h"
//#include "HarpoReaderRoot.h"
#include "HarpoRunHeader.h"

#include "HarpoAnalyseWriter.h"

int main(int argc, char **argv)
{
  Long_t processed = -1;  // Processed events counter , negative on error
  
  HarpoRunHeader *hdr;
  //HarpoEvent *hevent;

  HarpoConfig *cfg = new HarpoConfig(argc,argv);
  cfg->Init();

  cfg->print();

  // printf ("Before Reader Init  ...\n");
  // gDirectory->pwd();

  HarpoReader *rdr = new HarpoReader( cfg );
  //HarpoReaderRoot *rdr = new HarpoReaderRoot( cfg );
  // hNoCMPTime,  // Read without timeshump comparision
  // hFirstTime, // return only plane which smaller time stamp
  // hSyncTime   // skip plane event with smaller timeshatp, return only

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

  HarpoAnalyseWriter *hana = new HarpoAnalyseWriter();

  // printf ("Before Set Header  ...\n");
  // gDirectory->pwd();

  hana->SetRunHeader(hdr);
  // printf ("Before Reader Init  ...\n");
  // gDirectory->pwd();
  hana->Init(); // if errors ???
  // printf ("After Reader Init  ...\n");
  // gDirectory->pwd();
  

  rdr->SetAnalyseFunction(hana);

  //if ( rdr->NextEvent() ) printf("\n Next Event OK\n");
  processed = rdr->LoopAll();
  
  hana->Save();
  delete hana;
  delete rdr;
  delete cfg;
  printf("\n Done processed Events %li\n",processed);
  return 0;
}