JpsiMods.C

Go to the documentation of this file.

// JpsiMods.C
//
// Author     : J. List
// Created    : 31/05/2005
// Last update: $Date: 2005/06/03 12:32:22 $
//          by: $Author: boehmej $
// 
// Comment    : 
//      This implements the H1OO jpsimods example
//      in SFHMarana



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

#include "SFHMarana/SFHMarEvent.h"
#include "SFHMarana/MaranaEventLoop.h"
#include "jbltools/sfh/FloatFun.h"
#include "jbltools/sfh/h1examples/JPsiFuns.h"
#include "jbltools/sfh/h1/H1ModsJPsis.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 << "JpsiMods input_parameters " <<endl;
    exit(1);
  }

//...declare the run manager
  TRunHisto run(argv[1]);

//...declare the  event --> use "local" SFHMarEvent (derived from TMarEvent)
  SFHMarEvent event(&run);
 
  new JpsiElasticCut (mods_JPsis);
  new JpsiTTCut (mods_JPsis);
  new JpsiNoCosmicCut (mods_JPsis);
  new Jpsi2MuCut (mods_JPsis);
  new Jpsi2IronMuCut (mods_JPsis);
  
//...declare histogram classes 
  Char_t OutputTextFile[200];

//--- create object which manages histogram booking & filling 
//    (derived from THisto derived from THistoManager)
  MaranaEventLoop evtLoop(run,"JpsiHistos","JpsiEvents.txt");

//---- hat selection (if needed to online preselect the sample) ...
  printf("Nb events before hat selection: %d \n",gH1Tree->GetEntries());
  H1Tree::Instance()->SelectHat("NumJPsi>0");
  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; 

//--- disable other thight finders
  event.find_Jet=kFALSE;
  event.find_NP=kFALSE;
  
    
//---- Here start the main loop on events
  while (event.Next(run)) { 
    if (event.num%5000==0)cout << "--> process event "<< event.num<<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;
 
}

---