PPbar.h

Go to the documentation of this file.

#ifndef Physics_Analysis_PPbar_H
#define Physics_Analysis_PPbar_H
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/ITHistSvc.h"
#include "GaudiKernel/NTuple.h"
 
#include "TH1.h"
 
class PPbar : public Algorithm {
 
public:
  PPbar( 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