rhopi.h

Go to the documentation of this file.

#ifndef Physics_Analysis_rhopi_H
#define Physics_Analysis_rhopi_H 
 
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
 
#include "GaudiKernel/ITHistSvc.h"
#include "TH1.h"
 
class rhopi : public Algorithm {
 
public:
  rhopi(const std::string& name, ISvcLocator* pSvcLocator);
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();  
 
private:
  ITHistSvc *m_thistsvc;
     TH1F *m_pi_vx;
     TH1F *m_pi_vy;
     TH1F *m_pi_vz;
     TH1F *m_pip_mom;
     TH1F *m_pim_mom;
     TH1F *m_pi0_mass;
     TH1F *m_rho0_mass;
     TH1F *m_rhop_mass;
     TH1F *m_rhom_mass;
     TH1F *m_chisq_4c;
     TH1F *m_chisq_5c;
     TH1F *m_chipi_dedx;
     TH1F *m_chie_dedx;
     TH1F *m_chimu_dedx;
     TH1F *m_chik_dedx;
     TH1F *m_chip_dedx;
     TH1F *m_cos_pip;
     TH1F *m_cos_pim;
  // Declare r0, z0 cut for charged tracks
  double m_vr0cut;
  double m_vz0cut;
 
  //Declare energy, dphi, dthe cuts for fake gamma's
  double m_energyThreshold;
  double m_gammaPhiCut;
  double m_gammaThetaCut;
 
  int m_pass[12];
  // 
 
  NTuple::Tuple*  m_tuple4;     // 5C
  NTuple::Item<double>  m_ngoodneu;
  NTuple::Item<double>  m_ngoodch;
  NTuple::Item<double>  m_npip;
  NTuple::Item<double>  m_npim;
  NTuple::Item<double>  m_chisq4c;
  NTuple::Item<double>  m_ppi0;
  NTuple::Item<double>  m_mpi0;
  NTuple::Item<double>  m_chisq5c;
  NTuple::Item<double>  m_ppi0fit;
  NTuple::Item<double>  m_mpi0fit;
  NTuple::Item<double>  m_g1inpi0the;
  NTuple::Item<double>  m_g2inpi0the;
  NTuple::Item<double>  m_theta2pi;
  NTuple::Item<double>  m_ppip;
  NTuple::Item<double>  m_ppim;
  NTuple::Item<double>  m_p2pi;
  NTuple::Item<double>  m_m2pi;
  NTuple::Item<double>  m_ppip0;
  NTuple::Item<double>  m_mpip0;
  NTuple::Item<double>  m_ppim0;
  NTuple::Item<double>  m_mpim0;
  NTuple::Item<double>  m_eneumiss;
  NTuple::Item<double>  m_pneumiss;
  NTuple::Item<double>  m_mneumiss;
 
  NTuple::Array<double> m_kal_cos;
  
  NTuple::Tuple*  m_tuple5;     // 5C
  NTuple::Item<double>  m_m2graw;
  NTuple::Item<double>  m_emiss;
  NTuple::Item<double>  m_pmiss;
  NTuple::Item<double>  m_mmiss;
  NTuple::Item<double>  m_prho0;
  NTuple::Item<double>  m_mrho0;
  NTuple::Item<double>  m_pmrho0;
  NTuple::Item<double>  m_mmrho0;
  NTuple::Item<double>  m_prhop;
  NTuple::Item<double>  m_mrhop;
  NTuple::Item<double>  m_pmrhom;
  NTuple::Item<double>  m_mmrhom;
  NTuple::Item<double>  m_prhom;
  NTuple::Item<double>  m_ppipraw;
  NTuple::Item<double>  m_mrhom;
};
 
#endif