#include <LocalPionSelector.h>

Inheritance diagram for LocalPionSelector:

Public Member Functions
	LocalPionSelector ()

bool	operator() (CDChargedPion &aPion)

void	setpidtype (int type)

	LocalPionSelector ()

bool	operator() (CDChargedPion &aPion)

void	setpidtype (int type)

Public Member Functions inherited from DCSelectionFunction< CDChargedPion >
	DCSelectionFunction ()

	DCSelectionFunction ()

virtual	~DCSelectionFunction ()

virtual	~DCSelectionFunction ()

virtual bool	operator() (CDChargedPion &)=0

bool	operator() (CDChargedPion &iArg) const

virtual bool	operator() (CDChargedPion &)=0

bool	operator() (CDChargedPion &iArg) const

Detailed Description

Definition at line 7 of file InstallArea/include/DTagAlg/DTagAlg/LocalPionSelector.h.

Constructor & Destructor Documentation

◆ LocalPionSelector() [1/2]

LocalPionSelector::LocalPionSelector ( )

Definition at line 20 of file LocalPionSelector.cxx.

{
  m_pidtype=0;
 
  IJobOptionsSvc* jobSvc;
  Gaudi::svcLocator()->service("JobOptionsSvc", jobSvc);
  
  PropertyMgr m_propMgr;
  
  //Declare the properties
  m_propMgr.declareProperty("RxyCut",               m_VrCut           = 1.0);
  m_propMgr.declareProperty("Vz0Cut",               m_VzCut           = 10.0);
  m_propMgr.declareProperty("CosThetaCut",          m_CosThetaCut     = 0.93);
 
    m_propMgr.declareProperty("PID",                  m_PID     = true);
    m_propMgr.declareProperty("useSimplePID",         m_useSimplePID    = true);
 
  m_propMgr.declareProperty("PidUseDedx",           m_useDedx         = true);
  m_propMgr.declareProperty("PidUseTof1",           m_useTof1         = true);
  m_propMgr.declareProperty("PidUseTof2",           m_useTof2         = true);
  m_propMgr.declareProperty("PidUseTofE",           m_useTofE         = false);
  m_propMgr.declareProperty("PidUseTofQ",           m_useTofQ         = false);
  m_propMgr.declareProperty("PidUseEmc",            m_useEmc          = false);
  m_propMgr.declareProperty("PidUseMuc",            m_useMuc          = false);
  
  m_propMgr.declareProperty("PidwithProbability",   m_probability     = true);
  m_propMgr.declareProperty("PidProbCut",           m_PidProbCut      = 0.001);
  m_propMgr.declareProperty("RejectKaon",           m_rejectKaon      = true);
  m_propMgr.declareProperty("RejectProton",         m_rejectProton    = false);
  
  m_propMgr.declareProperty("PidwithLikelihood",    m_likelihood      = false);
  m_propMgr.declareProperty("PidwithNeuronNetwork", m_neuronNetwork   = false);
  m_neuronValCut.clear();
  m_neuronValCut.push_back(1.5);
  m_neuronValCut.push_back(2.5);
  m_propMgr.declareProperty("NeuronValue",      m_neuronValCut);
  
  jobSvc->setMyProperties("LocalPionSelector", &m_propMgr);
}

◆ LocalPionSelector() [2/2]

LocalPionSelector::LocalPionSelector ( )

Member Function Documentation

◆ operator()() [1/2]

bool LocalPionSelector::operator() ( CDChargedPion & aPion )

virtual

Implements DCSelectionFunction< CDChargedPion >.

Definition at line 60 of file LocalPionSelector.cxx.

                                                        {
 
  aPion.setUserTag(1);
  EvtRecTrack* recTrk = const_cast<EvtRecTrack*>( aPion.track() );
 
 
  // MDC track selection
 
  RecMdcKalTrack::setPidType(RecMdcKalTrack::pion);
  if ( !recTrk->isMdcKalTrackValid() ) return false;
  RecMdcKalTrack* mdcKalTrk = recTrk->mdcKalTrack();
  if ( mdcKalTrk->charge()==0 ) return false;
 
  Hep3Vector xorigin(0,0,0);
  IVertexDbSvc*  vtxsvc;
  Gaudi::svcLocator()->service("VertexDbSvc", vtxsvc);
  if(vtxsvc->isVertexValid()){
    double* dbv = vtxsvc->PrimaryVertex();
    double*  vv = vtxsvc->SigmaPrimaryVertex();
    xorigin.setX(dbv[0]);
    xorigin.setY(dbv[1]);
    xorigin.setZ(dbv[2]);
  }
  HepVector a = mdcKalTrk->getZHelix();
  HepSymMatrix Ea = mdcKalTrk->getZError();
  HepPoint3D point0(0.,0.,0.);
  HepPoint3D IP(xorigin[0],xorigin[1],xorigin[2]);
  VFHelix helixip3(point0,a,Ea);
  helixip3.pivot(IP);
  HepVector  vecipa = helixip3.a();
 
  double dr=fabs(vecipa[0]);
  double dz=fabs(vecipa[3]);
  double costheta=cos(mdcKalTrk->theta());
  if (  dr>= m_VrCut ) return false;
  if (  dz>= m_VzCut ) return false;
  if ( fabs(costheta) >= m_CosThetaCut ) return false;
 
    if(!m_PID) return true;
 
  if(m_useSimplePID){
    // PID
    ISimplePIDSvc*         m_simplePIDSvc; 
    Gaudi::svcLocator()->service("SimplePIDSvc", m_simplePIDSvc);
    m_simplePIDSvc->preparePID(recTrk);
    if( !m_simplePIDSvc->ispion() ){
      aPion.setUserTag(2);
    }
  }
  else{
    //particleID pacakge
    ParticleID *pid = ParticleID::instance();
    pid->init();
    pid->setMethod(pid->methodProbability());
    pid->setChiMinCut(4);
    pid->setRecTrack(recTrk);
    pid->usePidSys(pid->useDedx() | pid->useTofCorr() );
    pid->identify(pid->onlyPion() | pid->onlyKaon());
    pid->calculate();
 
    if( !(pid->probPion()>0 && pid->probPion()>pid->probKaon()) )
      aPion.setUserTag(2);
 
  }
  
  if(m_pidtype==0) return true;
  
  if(aPion.userTag()==1)
    return true;
  else
    return false;
  
}