TRun.C

Go to the documentation of this file.

// Class TRun
//
// Author     : F. Cassol Brunner
// Created    : 1/10/2002
// Last update: 
//          by: 
// Comment: class to perform open/close operation of the data files

#ifndef __TRUN
#include "Marana/TRun.h"


TRun::TRun(string par_file)
{

 try {

   string in_file;
   

//...initialize the parameters

   debug=0;
   nMaxEvt=10000000;
   lumi_file="not_given";
   lumi_file_bad_runs="not_given";
   Type="Data.0";
   Year="9900";
   HfsCalib=1;
   ElCalib=0;
   Lumi=1;

//... read Parameter file 
  ifstream inParams(par_file.c_str());
   
  if(!inParams)    
  {
      throw string("TRun::TRun = Error in opening steering file!\n");


  } 


  cout << endl;
  cout << "<=========================================================>" << endl;
  cout << " Parameters file in use: " << par_file <<endl;
  cout << "<=========================================================>" << endl;
  cout << endl << " The read parameters are: "<<endl;

  char item[255];
  while(inParams >> item)
  {
   
    if (strcmp(item,"in_file_list:") == 0)
    { 
      inParams >> in_file ;
      cout << "in_file_list= " <<in_file << endl;
    }  

    if (strcmp(item,"nMaxEvt:") == 0)
    { 
      inParams >> nMaxEvt ;
      cout << "nMaxEvt= " <<nMaxEvt << endl;
    }
    if (strcmp(item,"debug:") == 0)
    { 
      inParams >> debug ;
      cout << "debug= " <<debug << endl;
    }
    if (strcmp(item,"runlist:") == 0) { 
      inParams >>  lumi_file;
      cout << "runlist= " <<lumi_file << endl;
    }

    if (strcmp(item,"lumi:") == 0) { 
      inParams >> Lumi;
      cout << "lumi= " << Lumi<< " "<<endl;
    }
    if (strcmp(item,"genCut:") == 0) {
       cout <<  "genCut= "<< endl;
        TMarCutMCList newList;
        inParams >> newList;
        cout  << newList << endl;
        genCut.push_back(newList);  
              
    }

    if (strcmp(item,"run_type:") == 0) {
      inParams >> Year  >> Type;
      cout << "run_type = "<<Year <<" "<<  Type <<endl;
     
    
    }
    if (strcmp(item,"calibration:") == 0) {
      inParams >>  HfsCalib;
      cout << "calibration (Hfs)= "<< HfsCalib <<endl;
    }    
    if (strcmp(item,"elcalibration:") == 0) {
      inParams >>  ElCalib;
      cout << "calibration (electrons)= "<< ElCalib <<endl;
    }    

  }
  cout <<endl;
   
//... close the Parameters file
  inParams.close(); 


#if USERTREE
  // Make user event
  UserEvent *myUserEvent = new UserEvent();
#endif
 
//...create  tree
  H1Tree::Instance();

//... open input file
  fstream in;
  in.open(in_file.c_str(),ios::in);     // open file 

  if (!in){
    ostringstream error;
    error << "Run::Run = Error in opening input list file! "
          <<in_file << " - stop" << ends;
    string mess=error.str();      
    throw mess;   
  }


//...read data files
  string filename="", type="";
  Int_t i=0; 
  while(1) { 
    in >> filename >> type; 
    
    if(filename[0]=='#')continue;
    
    if (!in.good()||in.eof()) break; 

    if (type=="LIST"){//---- to read eventlist files 
        gH1Tree->AddList(filename.c_str());
        string filename_short=filename.substr(filename.find_last_of("/")+1);
        list_files.push_back(filename_short);
    } else {//--- otherwise read mods/hat/usertree list
#if USERTREE
        if(type=="USER") gH1Tree->AddFile(filename.c_str(), "User", myUserEvent); 
        else gH1Tree->AddFile(filename.c_str(),type.c_str());  
#else
        gH1Tree->AddFile(filename.c_str(),type.c_str());
#endif
    
        string filename_short=filename.substr(filename.find_last_of("/")+1);
  
        if(type=="HAT") hat_files.push_back(filename_short);
        else if(type=="MODS") mods_files.push_back(filename_short);
        else if(type=="ODS") ods_files.push_back(filename_short);
        else if(type=="USER") user_files.push_back(filename_short);
    }
    i++; 
  }  
  in.close();

  cout <<"\n >>>>>> Total of "<<i<<" files opened  <<<<<<<\n"<<endl;   



//... Open fields in the Tree  
  gH1Tree->Open();  


 }
//... try 
 catch(string error)
 {
    cout << error <<endl;;
    exit(0); 
  
 }
 catch(...)
 {
     throw string("TRun::TRun = Undefined error! - stop");
     exit(0); 
 }     
  

};

TRun::~TRun(){


  gH1Tree->Close();

}
#endif

---