PipiJpsi.h

Go to the documentation of this file.

// psi'--> J/psi pion pion, J/psi --> di-leptons
// Kai Zhu (zhuk@ihep.ac.cn)
#ifndef Physics_Analysis_PipiJpsi_H
#define Physics_Analysis_PipiJpsi_H
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "trkInfo.h"
 
class PipiJpsi : public Algorithm {
 
public:
  PipiJpsi( const std::string& name, ISvcLocator* pSvcLocator );
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
private:
  // Declare r0, z0 cut for charged tracks
  double m_vr0cut;
  double m_vz0cut;
 
  //
  bool m_checkDedx;
  bool m_checkTof;
  bool m_eventrate;
  int  m_chan_det;
  // declare the track angle cut and track's transverse momentum
  double m_cosThetaCut;
 
  // declare the cut angle between two pions to drop gamma conversion
  double m_pipi_dang_cut;
 
  // declare whether pick up sub-sample
  bool m_subsample_flag, m_trigger_flag;
 
  // declare energy/moentum to distinguish e and muon
  double m_distin_emuon;
 
  // define Ntuples here
  NTuple::Tuple*       m_tuple1; // charged track vertex
  NTuple::Item<double> m_vx0;
  NTuple::Item<double> m_vy0;
  NTuple::Item<double> m_vz0;
  NTuple::Item<double> m_vr0;
 
  NTuple::Tuple*       m_tuple2; // fake photon
  NTuple::Item<double> m_dthe;
  NTuple::Item<double> m_dphi;
  NTuple::Item<double> m_dang;
  NTuple::Item<double> m_eraw;
  NTuple::Item<long>   m_nGam;
 
  NTuple::Tuple*       m_tuple3; // dE/dx
  NTuple::Item<double> m_ptrk;
  NTuple::Item<double> m_chie;
  NTuple::Item<double> m_chimu;
  NTuple::Item<double> m_chipi;
  NTuple::Item<double> m_chik;
  NTuple::Item<double> m_chip;
  NTuple::Item<double> m_probPH;
  NTuple::Item<double> m_normPH;
  NTuple::Item<double> m_ghit;
  NTuple::Item<double> m_thit;
 
  NTuple::Tuple*       m_tuple4; // endcap tof
  NTuple::Item<double> m_ptot_etof;
  NTuple::Item<double> m_cntr_etof;
  NTuple::Item<double> m_path_etof;
  NTuple::Item<double> m_tof_etof;
  NTuple::Item<double> m_te_etof;
  NTuple::Item<double> m_tmu_etof;
  NTuple::Item<double> m_tpi_etof;
  NTuple::Item<double> m_tk_etof;
  NTuple::Item<double> m_tp_etof;
  NTuple::Item<double> m_ph_etof;
  NTuple::Item<double> m_rhit_etof;
  NTuple::Item<double> m_qual_etof;
 
  NTuple::Tuple*       m_tuple5; // barrel inner tof
  NTuple::Item<double> m_ptot_btof1;
  NTuple::Item<double> m_cntr_btof1;
  NTuple::Item<double> m_path_btof1;
  NTuple::Item<double> m_tof_btof1;
  NTuple::Item<double> m_te_btof1;
  NTuple::Item<double> m_tmu_btof1;
  NTuple::Item<double> m_tpi_btof1;
  NTuple::Item<double> m_tk_btof1;
  NTuple::Item<double> m_tp_btof1;
  NTuple::Item<double> m_ph_btof1;
  NTuple::Item<double> m_zhit_btof1;
  NTuple::Item<double> m_qual_btof1;
 
  NTuple::Tuple*       m_tuple6; // barrel outer tof
  NTuple::Item<double> m_ptot_btof2;
  NTuple::Item<double> m_cntr_btof2;
  NTuple::Item<double> m_path_btof2;
  NTuple::Item<double> m_tof_btof2;
  NTuple::Item<double> m_te_btof2;
  NTuple::Item<double> m_tmu_btof2;
  NTuple::Item<double> m_tpi_btof2;
  NTuple::Item<double> m_tk_btof2;
  NTuple::Item<double> m_tp_btof2;
  NTuple::Item<double> m_ph_btof2;
  NTuple::Item<double> m_zhit_btof2;
  NTuple::Item<double> m_qual_btof2;
 
  // with the method of momentum selection
  NTuple::Tuple*       m_tuple8;
  NTuple::Item<double> m_mom_lepm;
  NTuple::Item<double> m_mom_lepp;
  NTuple::Item<double> m_mom_pionp;
  NTuple::Item<double> m_mom_pionm;
  NTuple::Item<double> m_pipi_dang;
  NTuple::Item<double> m_cms_lepp;
  NTuple::Item<double> m_cms_lepm;
  NTuple::Item<double> m_mass_twopi;
  NTuple::Item<double> m_mass_jpsi;
  NTuple::Item<double> m_mass_recoil;
  NTuple::Item<double> m_inv_mass;
  NTuple::Item<double> m_tot_e;
  NTuple::Item<double> m_tot_px;
  NTuple::Item<double> m_tot_py;
  NTuple::Item<double> m_tot_pz;
  NTuple::Item<double> m_ep_ratio;
  NTuple::Item<long>   m_event_flag; // 3 or 4 tracks, 4=>4 tracks, 0=> miss pi+, 1=> miss pi-,
                                     // 2=> miss lepton+, 3=> miss lepton-
  NTuple::Item<double>   m_trans_ratio_lepm;
  NTuple::Item<double>   m_trans_ratio_lepp;
  NTuple::Item<double>   m_trans_ratio_pionp;
  NTuple::Item<double>   m_trans_ratio_pionm;
  NTuple::Item<long>     m_run;
  NTuple::Item<long>     m_event;
  NTuple::Item<long>     m_index;
  NTuple::Array<double>  m_cos_theta;
  NTuple::Array<double>  m_phi;
  NTuple::Matrix<double> m_four_mom;
  // we don't need smear, actually, assume one track of pi pi is missing, the recoil mass is
  // actually the invariant of pair-lepton
  NTuple::Item<long> m_pion_matched;
  NTuple::Item<long> m_lep_matched;
  // book MCtruth
  NTuple::Item<long>   m_idxmc;
  NTuple::Array<long>  m_pdgid;
  NTuple::Array<long>  m_motheridx;
  NTuple::Item<double> m_true_pionp; // pions' momentum from MC truth
  NTuple::Item<double> m_true_pionm;
};
 
#endif