THisto.C

Go to the documentation of this file.

// Class THisto
//
// Author     : F. Cassol Brunner
// Created    : 6/11/2002
// Last update: 
//          by: 


#include "Marana/THisto.h"




ClassImp(THisto)

THisto::THisto()
{
//--- empty constructor

}

THisto::THisto(const TRunHisto& run,string Tag,string textfile)
{
  debug=run.debug;
  
  NumSelEve=0;
  NumEve=0;
  NumSelEvePerRun=0;
  NumSelEvePerRunOld=0;
  fp=NULL;
  hatCopy=NULL; 
  modsCopy=NULL; 

//... Open file to write selected events
  if(textfile!=""){
    
    if ((fp  = fopen(textfile.c_str(),"w+")) == NULL) {
       cout << "THisto::THisto : error in opening "<< textfile <<endl;
    }
  }
//...read the parameters and book the histograms in the / dir
  ReadParameters(run.histo_par_file,Tag,0);

//...read the parameters and book the histograms in each sub-directory of run.output_file_root
  TList* dir=run.output_file_root->GetListOfKeys();
  for(int i=0;i<dir->GetSize();i++){
    cout << "--> Create histo in  dir "<< dir->At(i)->GetName() << endl;
    run.output_file_root->cd(dir->At(i)->GetName());
    
    ReadParameters(run.histo_par_file,Tag,i+1);
  }
  run.output_file_root->cd();

  ntuple=NULL;

//... open hat/mods file to write selected events
  
  if(run.out_file_mods_hat!=""){
    

    for(list<string>::const_iterator itag=run.TagToBeWritten.begin();itag!=run.TagToBeWritten.end();itag++)
    {
        

        if((*itag)==Tag){ 
          string dir="";
          string file="";
          unsigned int pos=run.out_file_mods_hat.rfind("/");
          if(pos!=string::npos){
             dir=run.out_file_mods_hat.substr(0,pos+1);   
             file=run.out_file_mods_hat.substr(pos+1) ;          
             
          }  else {
             file=run.out_file_mods_hat; 
          
          }   

            
          if ( gH1Tree->GetChain(3)) {
            string hatname;

            hatname=dir+"HAT."+Tag+"."+file;
            cout << "\n--> output hat: "<< hatname<<endl;      
        
            hatCopy = gH1Tree->CopyTree("HAT", hatname.c_str() );
            if (!hatCopy) throw string("THisto::THisto= could not copy HAT-tree!");
            
        
          }

          if ( gH1Tree->GetChain(2) ) {
            string modsname;
            modsname=dir+"MODS."+Tag+"."+file;
            cout << "\n--> output mods: "<< modsname<<endl;    

            modsCopy = gH1Tree->CopyTree("MODS", modsname.c_str());
            if (!modsCopy)throw string("THisto::THisto= could not copy MODS-tree!");
          }
        }

     }  
  }
  

}
THisto::~THisto()
{
// Base destructor
  if ( hatCopy )  hatCopy->BuildIndex("RunNumber","EventNumber");
  if ( modsCopy ) modsCopy->BuildIndex("RunNumber","EventNumber");

  
  if (fp) fclose(fp);
  
}
void THisto::Print(const TMarEvent &event){
  if (fp){
  // Write selected events of channel in text file 
    Char_t OutString[50];
    sprintf(OutString,"%d %d\n",event.RunNumber,event.EventNumber);
    fputs(OutString,fp);
  }
  else {
    cout << "THisto::Print = error, no output file opened!"<<endl;
  }
}


Bool_t THisto::Cut(TMarEvent &event){

//...Fill Lumi info before the cut
  if(event.Syst_Type==0){
    FillLumi(event);
    NumEve++;
  }
  
  if(event.IsCosmic()) return kFALSE;

  if(fabs(event.Vertex[2])>50) return kFALSE;

  if(event.Syst_Type==0){  
    NumSelEve++;  
    NumSelEvePerRun++;
  }

  return kTRUE; 

} 


void THisto::FillLumi(TMarEvent &event)
{

//... fill lumi info if new run
  if (event.lumi==NULL)return;

  if (!event.NewRun)
  {
     NumSelEvePerRunOld=NumSelEvePerRun;
     return;
  }  


  if(event.num>1){
//...here one has to fill the lumi info objects      

      cout <<  "  - lumi : "<<event.lumi->TotalLumiOld << " " << event.lumi->LumiPerRunOld <<" "<< NumSelEve<< 
           " "<<NumSelEvePerRunOld<< endl; 
   
//     FillHisto(1,event.lumi->TotalLumiOld,NumSelEve,event.MCWeight);
//     FillHisto(2,event.lumi->TotalLumiOld,NumSelEvePerRunOld/event.lumi->LumiPerRunOld,event.MCWeight);
  
//     if(nt_lumi!=NULL)nt_lumi->Fill(event.lumi->TotalLumiOld,event.lumi->LumiPerRunOld,(Float_t)NumSelEve,(Float_t)NumSelEvePerRunOld);

     NumSelEvePerRun=0;  
    
 
  }  


}

---