TRunHisto.C

Go to the documentation of this file.

// Class TRun
//
// Author     : F. Cassol Brunner
// Created    : 24/10/2002
// Last update: 
//          by: 
// Comment: class to perform job steering in programs that fill histograms

#ifndef __TRUNHISTO
#include "Marana/TRunHisto.h"


TRunHisto::TRunHisto(string par_file): TRun(par_file)
{
   
try{

  
  out_file_mods_hat="";
  
  string out_file;

   
//... read specific parameters from the Parameter file 
  ifstream inParams(par_file.c_str());
  if(!inParams)    
  {
      throw string("TRunHisto::TRunHisto = Error in opening steering file!");

  } 

  char item[255];
  while(inParams >> item)
  {
    
    if (strcmp(item,"in_histo_par:") == 0)    
    { 
      inParams >> histo_par_file ;
      cout << "histo_par_file= " <<histo_par_file<< endl;
    }  
    if (strcmp(item,"out_file_histo:") == 0)
    { 
      inParams >> out_file ;
      cout << "out_file_histo= " <<out_file << endl;
    }   
// //---- to write out H1eventlist of selected events
//     if (strcmp(item,"out_file_list:") == 0)
//     { 
//       inParams >> out_file_list ;
//       cout << "out_file_list= " <<out_file_list << endl;
//     }   

    if (strcmp(item,"write_mods_hat:") == 0)
    { 
      string tag="";
      Int_t num=0;  
      cout << "write_mods_hat= ";
      inParams >> out_file_mods_hat >> num;
      cout <<out_file_mods_hat << " " <<num;

      for (Int_t i=0;i<num;i++){
       inParams >> tag;
       TagToBeWritten.push_back(tag); 
       cout << " " << tag;
      }
      cout <<  endl;
    }   
    if (strcmp(item,"systematics:") == 0){
       cout << "systematics= ";
       inParams >>  Syst;
       cout << Syst ;
    }

  }
  cout <<endl;
    
//... close the Parameters file
  inParams.close(); 
 
 
// //... Open H1Eventlist 
// //    cout << " avant le event list " << endl;
//   output_file_list= new TFile(out_file_list.c_str(),"RECREATE","Event list file");
//   list= new H1EventList("mylist");
// 
//   if (output_file_list==NULL)
//   {
//     ostrstream error;
//     error << "TRunHisto::TRunHisto = Error in opening H1Eventlist file! "
//           << out_file_list << " - stop" << ends;
//     string mess=error.str();      
//     throw mess;     
//   }
 

//... open output histo file 
  output_file_root= new TFile(out_file.c_str(),"RECREATE","Histos from ModsHat",2);

  if (output_file_root==NULL)
  {
    ostringstream error;
    error << "TRunHisto::TRunHisto = Error in opening output file! "
          << out_file<< " - stop" << ends;
    string mess=error.str();      
    throw mess;     
  }

//... make output directories
  output_file_root->cd();

  
  for(list<string>::iterator i=Syst.ScalesToBeDone.begin();i!=Syst.ScalesToBeDone.end();i++){
     string dir=(*i);  
     if(dir=="/")continue; 
     TDirectory* DirPtr=output_file_root->mkdir(dir.c_str());
      
     if(DirPtr==0) throw string("TRunHisto::TRunHisto = problem in mkdir");
    
  }

 

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

  return;   

};


TRunHisto::~TRunHisto()
{
 gDebug=0;
 
// gDirectory->ls();
// output_file_root->ls();

// output_file_root->Print();

 output_file_root->Close(); 
 

 delete output_file_root;
 
}


// void TRunHisto::AddToList()
// {
//   list->Enter();
// }
// 
// 
// void TRunHisto::CloseList()
// {
// 
//  output_file_list->cd();
// 
//  list->Write();
// 
//  list->Print();
// 
// //...Write event list to file 
// 
//  output_file_list->Write(); 
//  output_file_list->Close(); 
// }



#endif

---