getMinMax.cxx

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <iomanip>
#include <cstring>
#include <string>
#include <vector>
#include <TFile.h>
#include <TTree.h>
#include <TString.h>
 
using namespace std;
 
// int main(int argc, char* argv[]){
//      char runlist[200];
//      char rawlist[200];
//      if(argc > 2){
//    sscanf(argv[1], "%s", runlist);
//    sscanf(argv[2], "%s", rawlist);
//      } else {
//    cout << "too few arguments" << endl;
//    return 0;
//      }
 
int main(int argc, char *argv[]){
  if(argc<2){
    cout << "please append a root file name " << endl;
    return 1;
  }
 
  TString str_root_file(argv[1]);
  str_root_file.Append(argv[2]);
  TFile f(str_root_file);
  TTree *t(0);
  if(argc>3) t = (TTree*)f.Get(argv[3]);
  else{
    TList *l = f.GetListOfKeys();
    if(l->GetEntries()<1){ 
      cout << "the root file " << argv[2] << " is empty, check it!" << endl;
      return 1;
    }
    else t = (TTree*)f.Get(l->At(0)->GetName());
  }
  Long_t m_totnum = t->GetEntries();
  double m_run(0), m_min(9999999), m_max(0), m_cur_run(0);
  t->SetBranchAddress("run", &m_run);
 
  for(Long_t i=0; i<m_totnum; i++){
    t->GetEntry(i);
    if(i%10000==0) cout << "i = " << i << "   run: " << m_run  << endl; 
    if(m_run != m_cur_run && m_run>1){
      if(m_min>m_run) m_min=m_run;
      if(m_max<m_run) m_max=m_run;
      m_cur_run = m_run;
    }
  }
  cout << "min " << (long)m_min << "  max " << (long)m_max << endl;
 
  return 0;
}