TestSFHMarana.C

Go to the documentation of this file.

// TestSFHMarana.C
//
// Author     : J. List
// Created    : 31/05/2005
// Last update: $Date: 2005/07/28 09:14:27 $
//          by: $Author: boehmej $
// 
// Comment    : 
//      This is a test executable for the SFHMarana package.
//      It is based of ExempleAnalysis from MarExemple by E.Sauvan



#include <cstdio>
#include <cstdlib> 
#include <iostream> 
#include <fstream> 

#include "SFHMarana/SFHMarEvent.h"
#include "SFHMarana/MaranaEventLoop.h"
#include "SFHMarana/Funs.h"
#include "Marana/TMarEvent.h"
#include "Marana/TTools.h"

// ----- Init ROOT environment ---------------- 
TROOT the_root("EXEMPLE","Exemple analysis"); 


using namespace std;

int main(int argc, char *argv[]) 
{ 

  if (argc <2) { 
    cout << "Usage:" <<endl;
    cout << "TestSFHMarana input_parameters " <<endl;
    exit(1);
  }

//...declare the run manager
  cout << "Creating TRunHisto" <<endl;
  TRunHisto run(argv[1]);

//...declare the  event --> use "local" SFHMarEvent (derived from TMarEvent)
  cout << "Creating SFHMarEvent" <<endl;
  SFHMarEvent event(&run);
 
//...declare histogram classes 
  Char_t OutputTextFile[200];

//--- create object which manages histogram booking & filling 
//    (derived from THisto derived from THistoManager)
  cout << "Creating MaranaEventLoop" <<endl;
  MaranaEventLoop evtLoop(run,"TestHistos","TestEvents.txt");

//---- hat selection (if needed to online preselect the sample) ...
  printf("Nb events before hat selection: %d \n",gH1Tree->GetEntries());
//H1Tree::Instance()->SelectHat("(NumEmParts>0&&NumMuons>0)||NumMuons>=2");
  printf("Nb events after hat selection: %d \n",gH1Tree->GetEntries());

//...loop over events  
  cout << "\n>>> start loop \n"<<endl;
  
  event.find_ElePho=kTRUE;
  event.PtMinEle=0;
  event.ThetaMinEle=5;
  event.ThetaMaxEle=175;

  event.find_Muon=kTRUE;
  event.PtMinMu=2;
  event.ThetaMinMu=10;
  event.ThetaMaxMu=160;
  event.DTrTrMinMu=0.5; 

  event.find_Jet=kTRUE;    // parameters from 1999 LQ paper
  event.PtMinJet=7.;
  event.ThetaMinJet=7;  
  event.ThetaMaxJet=145;
  event.RMinJet=0.1;      // default in TMarEvent
  event.EmFracJet=0.95;   // max!
 
    
//--- disable other thight finders
  event.find_NP=kFALSE;

  FloatFun& ffNSelJets = *new NSelJetsFun(event.jet, "NSelJets");
  FloatFun& ffPtSelJet1 = *new PtSelJet1Fun(event.jet, "PtSelJet1");
  FloatFun& ffThSelJet1 = *new ThSelJet1Fun(event.jet, "ThSelJet1");

   
//---- Here start the main loop on events
  while (event.Next(run)) { 
    //if (event.num%5000==0)cout << "--> process event "<< event.num<<endl; 
    if (event.num%1000==0) cout << "--> process event " << event.num 
                                << ", NSelJets = " << ffNSelJets() 
                                << ", PtSelJet1 = " << ffPtSelJet1() 
                                << ", ThSelJet1 = " << ffThSelJet1() << endl; 

    if (!(event.MCFlag)) { //--- checks for data quality runs
      static H1ShortPtr RunQuality("RunQ");
      if (*RunQuality==3&&!event.IsData2003()) {
        printf("Poor Quality run rejected \n");
        continue;
      }
    }      
           
    if(!event.IsHV()||event.IsCosmic()) continue;       //--- HV, cosmic base cut

    if((fabs(event.Vertex[2])>35)||event.Vertex[2]==0) continue;        //--- HV, cosmic base cut

//---- Check Timing
    if(event.IsTimingOk() == kFALSE) continue;
    

//---- Filling of histograms
    evtLoop.Fill();        // trigger self-filling process 
  }
      
 
 
 cout << "\n>>>>>>>>  End of events loop : " << event.num << " Events  <<<<<<<<<" << endl; 

  return 0;
 
}

---