00001 #ifndef DiJetAnalysis_H
00002 #define DiJetAnalysis_H
00003
00004 #include "HepMC/GenEvent.h"
00005
00006 #include <iostream>
00007 #include <stdio.h>
00008 #include "HepMC/IO_GenEvent.h"
00009 #include "HepMC/GenEvent.h"
00010 #include "HepMC/IO_AsciiParticles.h"
00011 #include "HepMC/SimpleVector.h"
00012
00013 #include "fastjet/JetDefinition.hh"
00014
00015 #include "TH1.h"
00016 #include "baseAnalysis.h"
00017
00018
00019 using namespace std;
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 class DiJetAnalysis: public baseAnalysis
00032 {
00033 public:
00034 DiJetAnalysis();
00035 ~DiJetAnalysis();
00036
00037 int Init(double tr_max_eta, double tr_min_pt);
00038 int Process(HepMC::GenEvent* hepmcevt);
00039
00040 private:
00041 TH1D *m_jet_count;
00042 TH1D *m_jet_pt;
00043 TH1D *m_jet_pt_log;
00044 TH1D *m_leadingjet_pt;
00045 TH1D *m_leadingjet_pt_log;
00046 TH1D *m_secondleadingjet_pt;
00047 TH1D *m_secondleadingjet_pt_log;
00048
00049
00050 TH1D *m_jet_phi;
00051 TH1D *m_jet_eta;
00052
00053 TH1D *m_charged_particle_multiplicity;
00054 TH1D *m_charged_particle_mean_pt;
00055 TH1D *m_charged_particle_temp_pt;
00056 TH1D *m_charged_particle_pt;
00057 TH1D *m_charged_particle_rms_pt;
00058 TH1D *m_charged_particle_pdgID;
00059
00060 TH1D *m_jet_Deta_leadingJet_secondJet;
00061 TH1D *m_jet_Dphi_leadingJet_secondJet;
00062 TH1D *m_jet_DR_leadingJet_secondJet;
00063
00064 using baseAnalysis::m_histVector;
00065 using baseAnalysis::m_outputFileName;
00066 using baseAnalysis::m_outputRootDir;
00067 };
00068 #endif