reader.cxx

Go to the documentation of this file.

 
 
 
#include <stdlib.h>
 
#include <iostream>
 
#include "TChain.h"
#include "TFile.h"
#include "TTree.h"
 
#include "TrigInDetAnalysis/TIDAEvent.h"
 
#include "utils.h"
 
 
int usage(std::ostream& s, int argc, char** argv) {
  if ( argc>0 ) { 
    s << "Usage: " << argv[0] << " [OPTIONS] [-r] <input1.root input2.root ...>\n\n";
    s << "Configuration: \n";
    s << "Options: \n";
    s << "    -r,  --release  \t show the release data\n";
    s << "    -ro, --releaseonly \t show the release data and then quit\n";
    s << "    -h, --help     \t this help\n";
    //   s << "\nReport bugs to <" << PACKAGE_BUGREPORT << ">";
    s << std::endl;
  }
  return 0;
}
 
 
 
//extern bool TIDAprinttracks;
 
int main(int argc, char** argv) { 
 
  bool show_release = false;
  bool quit_after_release = false;
 
  std::vector<std::string> files;
  
  if ( argc<2 ) { 
    //    std::cerr << "usage: reader [-v] <filename>" << std::endl;
    usage( std::cerr, argc, argv );
    exit(-1);
  }
  else { 
    for ( int i=1 ; i<argc ; i++ ) { 
      std::string argi = std::string(argv[i]);
      if      ( argi == "-r"  || argi == "--release" ) show_release = true;
      if      ( argi == "-ro" || argi == "--releaseonly" ) show_release = quit_after_release = true;
      else if ( argi == "-h"  || argi == "--help" )    return usage( std::cout, argc, argv );
      else    files.push_back( std::move(argi) );
    }
  }
 
 
  for ( unsigned i=0 ; i<files.size() ; i++ ) {
    
    TFile* finput = TFile::Open( files[i].c_str() );
 
    if ( finput==0 || !finput->IsOpen() || finput->IsZombie() ) {
      std::cerr << "Error: could not open output file" << std::endl;
      exit(-1);
    }
 
    if ( show_release || quit_after_release ) {
      TTree*  dataTree = (TTree*)finput->Get("dataTree");
      TString* releaseData = new TString("");
      
      if ( dataTree ) { 
    dataTree->SetBranchAddress( "ReleaseMetaData", &releaseData);
    
    for (unsigned int i=0; i<dataTree->GetEntries() ; i++ ) {
      dataTree->GetEntry(i);      
      std::cout << "releaseData: " << *releaseData << std::endl;
    }
      }
      else { 
    std::cerr << "no release data" << std::endl;
      }
      if ( quit_after_release ) return 0;
    }      
 

    TTree* data = (TTree*)finput->Get("tree");
    
    TIDA::Event* track_ev = new TIDA::Event();
    data->SetBranchAddress("TIDA::Event",&track_ev);
    //    data->AddFile( argv[i] );
    
    for (unsigned int i=0; i<data->GetEntries() ; i++ ) {
      data->GetEntry(i);      
      std::cout << *track_ev << std::endl;
    }
    
    finput->Close();
    
  }
    
  return 0;
}