00001 #ifndef Physics_Analysis_Rhopi_H 00002 #define Physics_Analysis_Rhopi_H 00003 00004 #include "GaudiKernel/AlgFactory.h" 00005 #include "GaudiKernel/Algorithm.h" 00006 #include "GaudiKernel/NTuple.h" 00007 //#include "VertexFit/ReadBeamParFromDb.h" 00008 00009 00010 class Rhopi : public Algorithm { 00011 00012 public: 00013 Rhopi(const std::string& name, ISvcLocator* pSvcLocator); 00014 StatusCode initialize(); 00015 StatusCode execute(); 00016 StatusCode finalize(); 00017 00018 private: 00019 00020 //ReadBeamParFromDb m_reader; 00021 // Declare r0, z0 cut for charged tracks 00022 double m_vr0cut; 00023 double m_vz0cut; 00024 00025 //Declare energy, dphi, dthe cuts for fake gamma's 00026 double m_energyThreshold; 00027 double m_gammaPhiCut; 00028 double m_gammaThetaCut; 00029 double m_gammaAngleCut; 00030 00031 // 00032 int m_test4C; 00033 int m_test5C; 00034 00035 // 00036 int m_checkDedx; 00037 int m_checkTof; 00038 00039 // define Ntuples here 00040 00041 NTuple::Tuple* m_tuple1; // charged track vertex 00042 NTuple::Item<double> m_vx0; 00043 NTuple::Item<double> m_vy0; 00044 NTuple::Item<double> m_vz0; 00045 NTuple::Item<double> m_vr0; 00046 NTuple::Item<double> m_rvxy0; 00047 NTuple::Item<double> m_rvz0; 00048 NTuple::Item<double> m_rvphi0; 00049 00050 NTuple::Tuple* m_tuple2; // fake photon 00051 NTuple::Item<double> m_dthe; 00052 NTuple::Item<double> m_dphi; 00053 NTuple::Item<double> m_dang; 00054 NTuple::Item<double> m_eraw; 00055 00056 NTuple::Tuple* m_tuple3; // rhopi: raw mgg, etot 00057 NTuple::Item<double> m_m2gg; 00058 NTuple::Item<double> m_etot; 00059 00060 NTuple::Tuple* m_tuple4; // rhopi 4C 00061 NTuple::Item<double> m_chi1; 00062 NTuple::Item<double> m_mpi0; 00063 00064 NTuple::Tuple* m_tuple5; // rhopi 5C 00065 NTuple::Item<double> m_chi2; 00066 NTuple::Item<double> m_mrh0; 00067 NTuple::Item<double> m_mrhp; 00068 NTuple::Item<double> m_mrhm; 00069 00070 NTuple::Tuple* m_tuple6; // photons 00071 NTuple::Item<double> m_fcos; 00072 NTuple::Item<double> m_elow; 00073 00074 NTuple::Tuple* m_tuple7; // dE/dx 00075 NTuple::Item<double> m_ptrk; 00076 NTuple::Item<double> m_chie; 00077 NTuple::Item<double> m_chimu; 00078 NTuple::Item<double> m_chipi; 00079 NTuple::Item<double> m_chik; 00080 NTuple::Item<double> m_chip; 00081 NTuple::Item<double> m_probPH; 00082 NTuple::Item<double> m_normPH; 00083 NTuple::Item<double> m_ghit; 00084 NTuple::Item<double> m_thit; 00085 00086 NTuple::Tuple* m_tuple8; // endcap tof 00087 NTuple::Item<double> m_ptot_etof; 00088 NTuple::Item<double> m_cntr_etof; 00089 NTuple::Item<double> m_te_etof; 00090 NTuple::Item<double> m_tmu_etof; 00091 NTuple::Item<double> m_tpi_etof; 00092 NTuple::Item<double> m_tk_etof; 00093 NTuple::Item<double> m_tp_etof; 00094 NTuple::Item<double> m_ph_etof; 00095 NTuple::Item<double> m_rhit_etof; 00096 NTuple::Item<double> m_qual_etof; 00097 00098 NTuple::Tuple* m_tuple9; // barrel inner tof 00099 NTuple::Item<double> m_ptot_btof1; 00100 NTuple::Item<double> m_cntr_btof1; 00101 NTuple::Item<double> m_te_btof1; 00102 NTuple::Item<double> m_tmu_btof1; 00103 NTuple::Item<double> m_tpi_btof1; 00104 NTuple::Item<double> m_tk_btof1; 00105 NTuple::Item<double> m_tp_btof1; 00106 NTuple::Item<double> m_ph_btof1; 00107 NTuple::Item<double> m_zhit_btof1; 00108 NTuple::Item<double> m_qual_btof1; 00109 00110 NTuple::Tuple* m_tuple10; // barrel outer tof 00111 NTuple::Item<double> m_ptot_btof2; 00112 NTuple::Item<double> m_cntr_btof2; 00113 NTuple::Item<double> m_te_btof2; 00114 NTuple::Item<double> m_tmu_btof2; 00115 NTuple::Item<double> m_tpi_btof2; 00116 NTuple::Item<double> m_tk_btof2; 00117 NTuple::Item<double> m_tp_btof2; 00118 NTuple::Item<double> m_ph_btof2; 00119 NTuple::Item<double> m_zhit_btof2; 00120 NTuple::Item<double> m_qual_btof2; 00121 00122 NTuple::Tuple* m_tuple11; // Particle ID info. 00123 NTuple::Item<double> m_ptrk_pid; 00124 NTuple::Item<double> m_cost_pid; 00125 NTuple::Item<double> m_dedx_pid; 00126 NTuple::Item<double> m_tof1_pid; 00127 NTuple::Item<double> m_tof2_pid; 00128 NTuple::Item<double> m_prob_pid; 00129 00130 }; 00131 00132 #endif