#include <LocalPionSelector.h>

Inheritance diagram for LocalPionSelector:

Public Member Functions
	LocalPionSelector ()

bool	operator() (CDChargedPion &aPion)

void	setpidtype (int type)

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

virtual	~DCSelectionFunction ()

virtual bool	operator() (CDChargedPion &)=0

bool	operator() (CDChargedPion &iArg) const

Detailed Description

Definition at line 7 of file LocalPionSelector.h.

Constructor & Destructor Documentation

◆ LocalPionSelector()

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("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);
}

Member Function Documentation

◆ operator()()

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

virtual

Implements DCSelectionFunction< CDChargedPion >.

Definition at line 57 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;
 
 
  // PID
  ISimplePIDSvc*         m_simplePIDSvc; 
  Gaudi::svcLocator()->service("SimplePIDSvc", m_simplePIDSvc);
  m_simplePIDSvc->preparePID(recTrk);
  if( !m_simplePIDSvc->ispion() ){
    aPion.setUserTag(2);
  }
  
  if(m_pidtype==0) return true;
  
  if(aPion.userTag()==1)
    return true;
  else
    return false;
  
}

◆ setpidtype()

void LocalPionSelector::setpidtype ( int type )

inline

Definition at line 14 of file LocalPionSelector.h.

14{m_pidtype=type;}

Referenced by ChargedDReconstruction::execute(), DsReconstruction::execute(), and NeutralDReconstruction::execute().