NCSelection.C

Go to the documentation of this file.

// NCSelection.C
//
// Author     : J. List
// Created    : 31/05/2005
// Last update: $Date: 2005/06/08 16:35:08 $
//          by: $Author: boehmej $
// 
// Comment    : 
//      This implements an NC Selection 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/h1/H1HatVariables.h"

#include <TROOT.h>
#include <TCanvas.h>
#include <TPostScript.h>
#include <TStyle.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 << "NCSelection 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);
 
//...declare histogram classes 
  Char_t OutputTextFile[200];

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

//---- hat selection (if needed to online preselect the sample) ...
  printf("Nb events before hat selection: %d \n",gH1Tree->GetEntries());
  H1Tree::Instance()->SelectHat("Q2e>500.");
  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 
  }
      
  //some style stuff
  gROOT->SetStyle("Plain");
  gStyle->SetOptTitle(1);
  gStyle->SetHistLineWidth(2);
  gStyle->SetCanvasBorderMode(0);
  gStyle->SetPadBorderMode(0);
  gStyle->SetPadColor(0);
  gStyle->SetCanvasColor(0);
  gStyle->SetTitleColor(1);
  gStyle->SetStatColor(0);
  gStyle->SetOptFit(1);
  gStyle->SetOptStat(1);

  TPostScript *ps  = new TPostScript("NC.ps", 111);
  TCanvas *c = new TCanvas("NC", "NCSelection", 600, 800);
  evtLoop.Draw(c, ps, "E0", false, 1, 3, false, false);
  ps->Close(); 
  cout << "Postscript written to NC.ps"  << endl;
 
  cout << "\n>>>>>>>>  End of events loop : " << event.num << " Events  <<<<<<<<<" << endl; 

  return 0;
 
}

---