Sing.cxx

Go to the documentation of this file.

//
//  Sing.cxx is the code to take care of some variables for anti-D0 tags for each single D0 tag
//  mode. It transfers these variables into a Class or Memmory, from which one can get these
//  variables for doing double tag analysis for study of the D0 meson decays. Sing.cxx was
//  transfered from the Fortran routine "Sing.f" which was orignally used for study of  the
//  D0D0-bar production and D0 meson decays at the BES-II experiment during the time period
//  from 2002 to 2008.
//
//  The orignal Fortran routine "Sing.f" used at the BES-II experiment was coded by G. Rong in
//  2001.
//
//  Sing.cxx was transfered by G. Rong and J. Liu in December, 2005.
//
//  Since 2008, G. Rong and L.L. Jiang have been working on developing this code to analyze of
//  the data taken at 3.773 GeV with the BES-III detector at the BEPC-II collider.
//
//  During developing this code, many People made significant contributions to this code. These
//  are
//          G. Rong, L.L. Jiang, J. Liu, H.L. Ma, J.C. Chen, D.H. Zhang,
//          M.G. Zhao, B. Zheng, L. Li, Y. Fang, Z.Y. Yi, H.H. Liu, Z.Q. Liu et al.
//
//                                       By G. Rong and L.L. Jiang
//                                       March, 2009
//
//  ==========================================================================================
//
#include "SD0TagAlg/Sing.h"
#include "SD0TagAlg/K0kk.h"
#include "SD0TagAlg/K0kpi.h"
#include "SD0TagAlg/K0pi0.h"
#include "SD0TagAlg/K0pipi.h"
#include "SD0TagAlg/K0pipipi0.h"
#include "SD0TagAlg/K3pi.h"
#include "SD0TagAlg/K3pipi0.h"
#include "SD0TagAlg/Kk.h"
#include "SD0TagAlg/Kkpi0.h"
#include "SD0TagAlg/Kkpipi.h"
#include "SD0TagAlg/Kpi.h"
#include "SD0TagAlg/Kpipi0.h"
#include "SD0TagAlg/Kpipi0pi0.h"
#include "SD0TagAlg/Pipi.h"
#include "SD0TagAlg/Pipipi0.h"
//  --------------------------------------------
 
#include "SD0TagAlg/SingleBase.h"
 
Sing::Sing() {}
 
Sing::~Sing() {}
 
void Sing::Mdset( double event, SmartDataPtr<EvtRecTrackCol> evtRecTrkCol, Vint iGood,
                  Vint iGam, int mdset, double Ebeam, int PID_flag, int Charge_candidate_D ) {
  oktg = false;
 
  if ( mdset == 1 )
  {
    //   ---- mode 11 ----- D0-->KPi ---------
    Kpi kpi;
    kpi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kpimd = kpi.Getkpimd();
 
    if ( kpimd )
    {
      tagmd_temp    = kpi.Gettagmd();
      mass_bc_temp  = kpi.Getmass_bc();
      cqtm_temp     = kpi.GetCQtm();
      delE_tag_temp = kpi.GetdelE_tag();
      iGoodtag_temp = kpi.Gettagtrk1();
      iGamtag_temp  = kpi.GettagGam1();
      ptag_temp     = kpi.Gettagp1();
 
      oktg = true;
    }
  }
 
  if ( mdset == 2 )
  {
    //   ---- mode 12 ----- D0 --> KPiPi0---------
    Kpipi0 kpipi0;
    kpipi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kpipi0md = kpipi0.Getkpipi0md();
    if ( kpipi0md )
    {
      tagmd_temp    = kpipi0.Gettagmd();
      iGoodtag_temp = kpipi0.Gettagtrk1();
      iGamtag_temp  = kpipi0.GettagGam1();
      ptag_temp     = kpipi0.Gettagp1();
      delE_tag_temp = kpipi0.GetdelE_tag();
      mass_bc_temp  = kpipi0.Getmass_bc();
      cqtm_temp     = kpipi0.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 4 )
  {
    //   ---- mode 13 ----- D0-->K3Pi ---------
    K3pi k3pi;
    k3pi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k3pimd = k3pi.Getk3pimd();
    if ( k3pimd )
    {
      tagmd_temp    = k3pi.Gettagmd();
      iGoodtag_temp = k3pi.Gettagtrk1();
      iGamtag_temp  = k3pi.GettagGam1();
      ptag_temp     = k3pi.Gettagp1();
      delE_tag_temp = k3pi.GetdelE_tag();
      mass_bc_temp  = k3pi.Getmass_bc();
      cqtm_temp     = k3pi.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 8 )
  {
    //   ---- mode 14-----D0-->K0PiPi ---------
    K0pipi k0pipi;
    k0pipi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k0pipimd = k0pipi.Getk0pipimd();
    if ( k0pipimd )
    {
      tagmd_temp    = k0pipi.Gettagmd();
      iGoodtag_temp = k0pipi.Gettagtrk1();
      iGamtag_temp  = k0pipi.GettagGam1();
      ptag_temp     = k0pipi.Gettagp1();
      delE_tag_temp = k0pipi.GetdelE_tag();
      mass_bc_temp  = k0pipi.Getmass_bc();
      cqtm_temp     = k0pipi.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 16 )
  {
    //   ---- mode 15 ----- D0-->KK ---------
    Kk kk;
    kk.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kkmd = kk.Getkkmd();
    if ( kkmd )
    {
      tagmd_temp    = kk.Gettagmd();
      iGoodtag_temp = kk.Gettagtrk1();
      iGamtag_temp  = kk.GettagGam1();
      ptag_temp     = kk.Gettagp1();
      delE_tag_temp = kk.GetdelE_tag();
      mass_bc_temp  = kk.Getmass_bc();
      cqtm_temp     = kk.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 32 )
  {
    //   ---- mode 16 ----- D0-->PiPiPi0 ---------
    Pipipi0 pipipi0;
    pipipi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool pipipi0md = pipipi0.GetPipipi0md();
    if ( pipipi0md )
    {
      tagmd_temp    = pipipi0.Gettagmd();
      iGoodtag_temp = pipipi0.Gettagtrk1();
      iGamtag_temp  = pipipi0.GettagGam1();
      ptag_temp     = pipipi0.Gettagp1();
      delE_tag_temp = pipipi0.GetdelE_tag();
      mass_bc_temp  = pipipi0.Getmass_bc();
      cqtm_temp     = pipipi0.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 64 )
  {
    //   ---- mode 17 ----- D0-->K0PiPiPi0 ---------
    K0pipipi0 k0pipipi0;
    k0pipipi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k0pipipi0md = k0pipipi0.GetK0pipipi0md();
    if ( k0pipipi0md )
    {
      tagmd_temp    = k0pipipi0.Gettagmd();
      iGoodtag_temp = k0pipipi0.Gettagtrk1();
      iGamtag_temp  = k0pipipi0.GettagGam1();
      ptag_temp     = k0pipipi0.Gettagp1();
      delE_tag_temp = k0pipipi0.GetdelE_tag();
      mass_bc_temp  = k0pipipi0.Getmass_bc();
      cqtm_temp     = k0pipipi0.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 128 )
  {
    //   ---- mode 18 ----- D0-->K0Pi0  ---------
    K0pi0 k0pi0;
    k0pi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k0pi0md = k0pi0.GetK0pi0md();
    if ( k0pi0md )
    {
      tagmd_temp    = k0pi0.Gettagmd();
      iGoodtag_temp = k0pi0.Gettagtrk1();
      iGamtag_temp  = k0pi0.GettagGam1();
      ptag_temp     = k0pi0.Gettagp1();
      delE_tag_temp = k0pi0.GetdelE_tag();
      mass_bc_temp  = k0pi0.Getmass_bc();
      cqtm_temp     = k0pi0.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 256 )
  {
    //   ---- mode 19 ----- D0-->Pi+Pi- mode -----
    Pipi pipi;
    pipi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool pipimd = pipi.GetPipimd();
    if ( pipimd )
    {
      tagmd_temp    = pipi.Gettagmd();
      iGoodtag_temp = pipi.Gettagtrk1();
      iGamtag_temp  = pipi.GettagGam1();
      ptag_temp     = pipi.Gettagp1();
      delE_tag_temp = pipi.GetdelE_tag();
      mass_bc_temp  = pipi.Getmass_bc();
      cqtm_temp     = pipi.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 512 )
  {
    //   ---- mode 20 ----- D0-->K+K-Pi+Pi- mode -----
    Kkpipi kkpipi;
    kkpipi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kkpipimd = kkpipi.Getkkpipimd();
    if ( kkpipimd )
    {
      tagmd_temp    = kkpipi.Gettagmd();
      iGoodtag_temp = kkpipi.Gettagtrk1();
      iGamtag_temp  = kkpipi.GettagGam1();
      ptag_temp     = kkpipi.Gettagp1();
      delE_tag_temp = kkpipi.GetdelE_tag();
      mass_bc_temp  = kkpipi.Getmass_bc();
      cqtm_temp     = kkpipi.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 1024 )
  {
    //   ---- mode 21 ----- D0-->K0K+Pi- mode -----
    K0kpi k0kpi;
    k0kpi.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k0kpimd = k0kpi.Getk0kpimd();
    if ( k0kpimd )
    {
      tagmd_temp    = k0kpi.Gettagmd();
      iGoodtag_temp = k0kpi.Gettagtrk1();
      iGamtag_temp  = k0kpi.GettagGam1();
      ptag_temp     = k0kpi.Gettagp1();
      delE_tag_temp = k0kpi.GetdelE_tag();
      mass_bc_temp  = k0kpi.Getmass_bc();
      cqtm_temp     = k0kpi.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 2048 )
  {
    //   ---- mode 22 ----- D0-->K0K+K- mode -----
    K0kk k0kk;
    k0kk.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k0kkmd = k0kk.Getk0kkmd();
    if ( k0kkmd )
    {
      tagmd_temp    = k0kk.Gettagmd();
      iGoodtag_temp = k0kk.Gettagtrk1();
      iGamtag_temp  = k0kk.GettagGam1();
      ptag_temp     = k0kk.Gettagp1();
      delE_tag_temp = k0kk.GetdelE_tag();
      mass_bc_temp  = k0kk.Getmass_bc();
      cqtm_temp     = k0kk.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 4096 )
  {
    //   ---- mode 23 ----- D0-->K-K+Pi0 mode  -----
    Kkpi0 kkpi0;
    kkpi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kkpi0md = kkpi0.Getkkpi0md();
    if ( kkpi0md )
    {
      tagmd_temp    = kkpi0.Gettagmd();
      iGoodtag_temp = kkpi0.Gettagtrk1();
      iGamtag_temp  = kkpi0.GettagGam1();
      ptag_temp     = kkpi0.Gettagp1();
      delE_tag_temp = kkpi0.GetdelE_tag();
      mass_bc_temp  = kkpi0.Getmass_bc();
      cqtm_temp     = kkpi0.GetCQtm();
 
      oktg = true;
    }
  }
 
  if ( mdset == 8192 )
  {
    //  ---- mode 24 ----- D0-->K-PiPiPiPi0 mode -----
    K3pipi0 k3pipi0;
    k3pipi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool k3pipi0md = k3pipi0.Getk3pipi0md();
    if ( k3pipi0md )
    {
      tagmd_temp    = k3pipi0.Gettagmd();
      iGoodtag_temp = k3pipi0.Gettagtrk1();
      iGamtag_temp  = k3pipi0.GettagGam1();
      ptag_temp     = k3pipi0.Gettagp1();
      delE_tag_temp = k3pipi0.GetdelE_tag();
      mass_bc_temp  = k3pipi0.Getmass_bc();
      cqtm_temp     = k3pipi0.GetCQtm();
      oktg          = true;
    }
  }
 
  if ( mdset == 16384 )
  {
    //   ----mode 25 ----- D0-->K-PiPi0Pi0 mode -----
    Kpipi0pi0 kpipi0pi0;
    kpipi0pi0.MTotal( event, evtRecTrkCol, iGood, iGam, Ebeam, PID_flag, Charge_candidate_D );
    bool kpipi0pi0md = kpipi0pi0.Getkpipi0pi0md();
    if ( kpipi0pi0md )
    {
      tagmd_temp    = kpipi0pi0.Gettagmd();
      iGoodtag_temp = kpipi0pi0.Gettagtrk1();
      iGamtag_temp  = kpipi0pi0.GettagGam1();
      ptag_temp     = kpipi0pi0.Gettagp1();
      delE_tag_temp = kpipi0pi0.GetdelE_tag();
      mass_bc_temp  = kpipi0pi0.Getmass_bc();
      cqtm_temp     = kpipi0pi0.GetCQtm();
      oktg          = true;
    }
  }
 
  return;
}