UEAnalysis.h

Go to the documentation of this file.

#ifndef UEAnalysis_H
#define UEAnalysis_H

#include "HepMC/GenEvent.h"

#include <iostream>
#include <stdio.h>
#include "HepMC/IO_GenEvent.h"
#include "HepMC/GenEvent.h"
#include "HepMC/IO_AsciiParticles.h"
#include "HepMC/SimpleVector.h"

#include "fastjet/JetDefinition.hh"

#include "TH1.h"
#include "TH2.h"
#include "baseAnalysis.h"

using namespace std;
/**
@class UEAnalysis.h
@brief Class to study underlying events. Currently the UE is studied for 
       events with at least one  jet, e.g. QCD 2->2, pp->W+jet, pp->Z+jet.
       Extension to DrellYan events will follow.
       Examples to generate e.g events are given in the specific generator
       directory, e.g. examples/pythia6 ; the histogramm output is stored
       in the UE folder of the root file.


@author Sebastian Johnert, Judith Katzy, Zhonghua Qin Dec 2008 */





class UEAnalysis: public baseAnalysis
{
  public:
  UEAnalysis();
  ~UEAnalysis();

  int Init(double tr_max_eta, double tr_min_pt);
  int Process(HepMC::GenEvent* hepmcevt);
  int finalize(TFile* output);

  private:
  
  TH1D *m_NchargedToward[26]; 
  TH1D *m_NchargedTransverse[26]; 
  TH1D *m_NchargedAway[26]; 
  
  TH1D *m_NchargedMeanToward; 
  TH1D *m_NchargedMeanTransverse; 
  TH1D *m_NchargedMeanAway; 

  TH1D *m_PtsumToward[26]; 
  TH1D *m_PtsumTransverse[26]; 
  TH1D *m_PtsumAway[26]; 
  
  TH1D *m_PtsumMeanToward; 
  TH1D *m_PtsumMeanTransverse; 
  TH1D *m_PtsumMeanAway; 


  TH1D *m_Njet; 
  TH1D *m_Ptjet;
  TH1D *m_Ptjet_log;
  TH1D *m_Ptleadingjet;
  TH1D *m_Ptleadingjet_log;

  TH1D *m_charged_particle_pt;

  /** definition for the variable bin size for the <Ncharged> vs pT calculation*/
  int m_nbin_pT;
  //size of m_nbinRange has to be nbin + 1 due to 0 at the beginning
  Double_t m_nbinRange_pT[26];

  using  baseAnalysis::m_histVector;
  using  baseAnalysis::m_outputFileName;
  using  baseAnalysis::m_outputRootDir;
};

#endif

---