UserAnalysis.cc

Go to the documentation of this file.

// header for this analysis class
#include "../include/UserAnalysis.h"

#include <iostream>

#include "HepMC/GenEvent.h"
#include "HepMC/IO_GenEvent.h"
#include "HepMC/GenParticle.h"
#include "HepMC/GenVertex.h"
#include "HepMC/IO_AsciiParticles.h"
#include "HepMC/SimpleVector.h"
#include "HepPDT/ParticleData.hh"
#include "CLHEP/Vector/LorentzVector.h"

// ROOT headers
#include "TH1.h"
#include "TMath.h"
#include "TLorentzVector.h"

// fastjet headers
#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"
#include "fastjet/JetDefinition.hh"
#include "fastjet/SISConePlugin.hh"

//**********************

/// empty default constructor
UserAnalysis::UserAnalysis()
{
}

/// empty default destructor
UserAnalysis::~UserAnalysis()
{
}


int UserAnalysis::Init(double tr_max_eta, double tr_min_pt)
{
  // specify default eta cut
  m_max_eta=tr_max_eta;
  
  // specify default pt cut
  m_min_pt=tr_min_pt;
  
  // default Output file name
  m_outputFileName="UserAnalysis.root";
  m_outputRootDir="User";
 
  m_h_njets   = initHist("njets", "N_jets", "Number of Jets", 15, 0., 15.);

  return true; 
}


int UserAnalysis::Process(HepMC::GenEvent* hepmcevt){ 

//   // print number of events
//   if (hepmcevt->event_number()%1 == 0) {
//     std::cout << "Processing event " << hepmcevt->event_number() << std::endl;
//   }
  
  // fill number of jets
  m_h_njets->Fill(m_inclusive_jets.size());
  
  return true;

} //end of Process

---