HepMCAnalysis: src/UserAnalysis.cc Source File

00001 /**
00002   
00003 */
00004 
00005 #include <iostream>
00006 
00007 #include "HepMC/GenEvent.h"
00008 #include "HepMC/IO_GenEvent.h"
00009 #include "HepMC/GenParticle.h"
00010 #include "HepMC/GenVertex.h"
00011 #include "HepMC/IO_AsciiParticles.h"
00012 #include "HepMC/SimpleVector.h"
00013 #include "HepPDT/ParticleData.hh"
00014 #include "CLHEP/Vector/LorentzVector.h"
00015 
00016 #include "TH1.h"
00017 #include "TMath.h"
00018 #include "TLorentzVector.h"
00019 
00020 #include "fastjet/PseudoJet.hh"
00021 #include "fastjet/ClusterSequence.hh"
00022 #include "fastjet/JetDefinition.hh"
00023 #include "fastjet/SISConePlugin.hh"
00024 
00025 #include "../include/UserAnalysis.h"
00026 
00027 // ---------------------------------------------------------------------- 
00028 UserAnalysis::UserAnalysis()
00029 {}
00030 
00031 // ---------------------------------------------------------------------- 
00032 UserAnalysis::~UserAnalysis()
00033 {}
00034 
00035 // ---------------------------------------------------------------------- 
00036 int UserAnalysis::Init( double maxEta, double minPt )
00037 {
00038   m_max_eta = maxEta;
00039   m_min_pt  = minPt;
00040   
00041   m_outputFileName  = "UserAnalysis.root";
00042   m_outputRootDir   = "User";
00043  
00044   m_h_njets = initHist( "njets", "N_jets", "Number of Jets", 15, 0.0, 15.0 );
00045 
00046   return true; 
00047 }
00048 
00049 // ---------------------------------------------------------------------- 
00050 int UserAnalysis::Process( HepMC::GenEvent* /*hepmcevt*/ )
00051 { 
00052   m_h_njets->Fill( m_inclusive_jets.size() );
00053   
00054   return true;
00055 }
00056