TestV.h

Go to the documentation of this file.

#ifndef Physics_Analysis_TestV_H
#define Physics_Analysis_TestV_H 
 
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/ITHistSvc.h"
 
#include "TH1.h"
 
class TestV : public Algorithm {
 
public:
  TestV(const std::string& name, ISvcLocator* pSvcLocator);
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();  
 
private:
  ITHistSvc *m_thistsvc;
 
  TH1F *m_pp_mass;      
  TH1F *m_pp_angle;     
  TH1F *m_pp_deltatof;  
  
  TH1F *m_pp_x_pp;      
  TH1F *m_pp_y_pp;      
  TH1F *m_pp_z_pp;      
  TH1F *m_pp_r_pp;      
  TH1F *m_pp_px_pp;     
  TH1F *m_pp_py_pp;     
  TH1F *m_pp_pz_pp;     
  TH1F *m_pp_p_pp;     
  TH1F *m_pp_cos_pp;    
  TH1F *m_pp_emc_pp;    
  TH1F *m_pp_pidchidedx_pp; 
  TH1F *m_pp_pidchitof1_pp; 
  TH1F *m_pp_pidchitof2_pp; 
 
  TH1F *m_pp_x_pb;      
  TH1F *m_pp_y_pb;      
  TH1F *m_pp_z_pb;      
  TH1F *m_pp_r_pb;      
  TH1F *m_pp_px_pb;     
  TH1F *m_pp_py_pb;     
  TH1F *m_pp_pz_pb;     
  TH1F *m_pp_p_pb;     
  TH1F *m_pp_cos_pb;    
  TH1F *m_pp_emc_pb;    
  TH1F *m_pp_pidchidedx_pb; 
  TH1F *m_pp_pidchitof1_pb; 
  TH1F *m_pp_pidchitof2_pb; 
  
 
  
  // Declare r0, z0 cut for charged tracks
  double m_vr0cut;
  double m_vz0cut;
  double m_coscut;
  bool   m_useKalTrk; 
  bool   m_tagPPbar; 
 
  //Declare energy, dphi, dthe cuts for fake gamma's
  double m_energyThreshold;
  double m_gammaPhiCut;
  double m_gammaThetaCut;
 
  // 
  int m_test4C;
  int m_test5C;
 
  // 
  int m_checkDedx;
  int m_checkTof;
  
 
  // define Ntuples here
 
  NTuple::Tuple*  m_tuple1;      // charged track vertex
 
  NTuple::Item<long>  m_run  ;
  NTuple::Item<long>  m_event;
  NTuple::Item<long>  m_nchrg;
  NTuple::Item<long>  m_nneu ;
  NTuple::Item<long>  m_ngch ;
 
  NTuple::Array<long> m_ch_truth ;
  NTuple::Array<long> m_ch_kal   ;
 
  NTuple::Array<long> m_charge ;
  NTuple::Array<double> m_vx0   ;
  NTuple::Array<double> m_vy0   ;
  NTuple::Array<double> m_vz0   ;
  NTuple::Array<double> m_vr0   ;
  NTuple::Array<double> m_px    ;
  NTuple::Array<double> m_py    ;
  NTuple::Array<double> m_pz    ;
  NTuple::Array<double> m_p     ;
  NTuple::Array<double> m_cos   ;
 
  NTuple::Array<double> m_bst_px    ;
  NTuple::Array<double> m_bst_py    ;
  NTuple::Array<double> m_bst_pz    ;
  NTuple::Array<double> m_bst_p     ;
  NTuple::Array<double> m_bst_cos   ;
 
  NTuple::Array<double> m_kal_vx0   ;
  NTuple::Array<double> m_kal_vy0   ;
  NTuple::Array<double> m_kal_vz0   ;
  NTuple::Array<double> m_kal_vr0   ;
  NTuple::Array<double> m_kal_px    ;
  NTuple::Array<double> m_kal_py    ;
  NTuple::Array<double> m_kal_pz    ;
  NTuple::Array<double> m_kal_p     ;
  NTuple::Array<double> m_kal_cos   ;
 
  NTuple::Array<double> m_bst_kal_px    ;
  NTuple::Array<double> m_bst_kal_py    ;
  NTuple::Array<double> m_bst_kal_pz    ;
  NTuple::Array<double> m_bst_kal_p     ;
  NTuple::Array<double> m_bst_kal_cos   ;
 
  NTuple::Array<double> m_vtx_px    ;
  NTuple::Array<double> m_vtx_py    ;
  NTuple::Array<double> m_vtx_pz    ;
  NTuple::Array<double> m_vtx_p     ;
  NTuple::Array<double> m_vtx_cos   ;
 
  NTuple::Array<double> m_probPH;
  NTuple::Array<double> m_normPH;
  NTuple::Array<double> m_chie  ;
  NTuple::Array<double> m_chimu ;
  NTuple::Array<double> m_chipi ;
  NTuple::Array<double> m_chik  ;
  NTuple::Array<double> m_chip  ;
  NTuple::Array<double> m_ghit  ;
  NTuple::Array<double> m_thit  ;
 
  NTuple::Array<double> m_e_emc ;
 
  NTuple::Array<double> m_tof_b1;
  NTuple::Array<double> m_beta_b1;
  NTuple::Array<double> m_tof_b2;
  NTuple::Array<double> m_beta_b2;
  NTuple::Array<double> m_tof_e;
  NTuple::Array<double> m_beta_e;
 
  NTuple::Array<double> m_clus_tof     ;
  NTuple::Array<double> m_clus_texp_e  ;
  NTuple::Array<double> m_clus_texp_mu ;
  NTuple::Array<double> m_clus_texp_pi ;
  NTuple::Array<double> m_clus_texp_k  ;
  NTuple::Array<double> m_clus_texp_p  ;
  NTuple::Array<double> m_clus_toff_e  ;
  NTuple::Array<double> m_clus_toff_mu ;
  NTuple::Array<double> m_clus_toff_pi ;
  NTuple::Array<double> m_clus_toff_k  ;
  NTuple::Array<double> m_clus_toff_p  ;
  NTuple::Array<double> m_clus_tsig_e  ;
  NTuple::Array<double> m_clus_tsig_mu ;
  NTuple::Array<double> m_clus_tsig_pi ;
  NTuple::Array<double> m_clus_tsig_k  ;
  NTuple::Array<double> m_clus_tsig_p  ;
  NTuple::Array<double> m_clus_path    ;
  NTuple::Array<double> m_clus_zrhit   ;
  NTuple::Array<double> m_clus_ph      ;
  NTuple::Array<double> m_clus_beta    ;
  NTuple::Array<double> m_clus_qual    ;
  NTuple::Array<double> m_clus_t0      ;
 
  NTuple::Array<double>   m_cntr_etof  ;
  NTuple::Array<double>   m_ptot_etof  ;
  NTuple::Array<double>   m_ph_etof    ;
  NTuple::Array<double>   m_rhit_etof  ;
  NTuple::Array<double>   m_qual_etof  ;
  NTuple::Array<double>   m_tof_etof   ;
  NTuple::Array<double>   m_te_etof    ;
  NTuple::Array<double>   m_tmu_etof   ;
  NTuple::Array<double>   m_tpi_etof   ;
  NTuple::Array<double>   m_tk_etof    ;
  NTuple::Array<double>   m_tp_etof    ;
 
  NTuple::Array<double>   m_cntr_btof1 ;
  NTuple::Array<double>   m_ptot_btof1 ;
  NTuple::Array<double>   m_ph_btof1   ;
  NTuple::Array<double>   m_zhit_btof1 ;
  NTuple::Array<double>   m_qual_btof1 ;
  NTuple::Array<double>   m_tof_btof1  ;
  NTuple::Array<double>   m_te_btof1   ;
  NTuple::Array<double>   m_tmu_btof1  ;
  NTuple::Array<double>   m_tpi_btof1  ;
  NTuple::Array<double>   m_tk_btof1   ;
  NTuple::Array<double>   m_tp_btof1   ;
 
  NTuple::Array<double>   m_cntr_btof2 ;
  NTuple::Array<double>   m_ptot_btof2 ;
  NTuple::Array<double>   m_ph_btof2   ;
  NTuple::Array<double>   m_zhit_btof2 ;
  NTuple::Array<double>   m_qual_btof2 ;
  NTuple::Array<double>   m_tof_btof2  ;
  NTuple::Array<double>   m_te_btof2   ;
  NTuple::Array<double>   m_tmu_btof2  ;
  NTuple::Array<double>   m_tpi_btof2  ;
  NTuple::Array<double>   m_tk_btof2   ;
  NTuple::Array<double>   m_tp_btof2   ;
 
  NTuple::Array<double>    m_ptrk_pid ;
  NTuple::Array<double>    m_cost_pid ;
  NTuple::Array<double>    m_dedx_pid ;
  NTuple::Array<double>    m_tof1_pid ;
  NTuple::Array<double>    m_tof2_pid ;
  NTuple::Array<double>    m_prob_pi  ;
  NTuple::Array<double>    m_prob_k   ;
  NTuple::Array<double>    m_prob_p   ;
 
  NTuple::Item<long>   m_np ;
  NTuple::Item<long>   m_npb;
 
  NTuple::Item<double>   m_m2p   ;
  NTuple::Item<double>   m_angle ;
  NTuple::Item<double>   m_deltatof ;
 
  NTuple::Item<double>   m_kal_m2p   ;
  NTuple::Item<double>   m_kal_angle ;
  
  NTuple::Item<double>   m_vtx_m2p   ;
  NTuple::Item<double>   m_vtx_angle ;
  
//  NTuple::Item<double>   m_Jpx   ;
//  NTuple::Item<double>   m_Jpy   ;
//  NTuple::Item<double>   m_Jpz   ;
//  NTuple::Item<double>   m_Jpp   ;
//  NTuple::Item<double>   m_m2p  ;
//                   
//  NTuple::Item<double>   m_p1pp  ;
//  NTuple::Item<double>   m_p1px  ;
//  NTuple::Item<double>   m_p1py  ;
//  NTuple::Item<double>   m_p1pz  ;
//  NTuple::Item<double>   m_p1cos ;
//
//  NTuple::Item<double>   m_p2pp ;
//  NTuple::Item<double>   m_p2px ;
//  NTuple::Item<double>   m_p2py ;
//  NTuple::Item<double>   m_p2pz ;
//  NTuple::Item<double>   m_p2cos;
//
//  NTuple::Item<double>   m_angle ;
 
};
 
#endif