LocalPionSelector Class Reference

#include <LocalPionSelector.h>

Inheritance diagram for LocalPionSelector:

Public Member Functions
	LocalPionSelector (const std::string &type, const std::string &name, const IInterface *parent)
bool	operator() (CDChargedPion &aPion) override
void	setpidtype (int type) override

Detailed Description

Definition at line 9 of file LocalPionSelector.h.

Constructor & Destructor Documentation

LocalPionSelector()

LocalPionSelector::LocalPionSelector	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 16 of file LocalPionSelector.cxx.

    : base_class( type, name, parent ) {
  m_pidtype = 0;
 
  // Declare the properties
  declareProperty( "RxyCut", m_VrCut = 1.0 );
  declareProperty( "Vz0Cut", m_VzCut = 10.0 );
  declareProperty( "CosThetaCut", m_CosThetaCut = 0.93 );
 
  declareProperty( "PID", m_PID = true );
  declareProperty( "useSimplePID", m_useSimplePID = true );
 
  declareProperty( "PidUseDedx", m_useDedx = true );
  declareProperty( "PidUseTof1", m_useTof1 = true );
  declareProperty( "PidUseTof2", m_useTof2 = true );
  declareProperty( "PidUseTofE", m_useTofE = false );
  declareProperty( "PidUseTofQ", m_useTofQ = false );
  declareProperty( "PidUseEmc", m_useEmc = false );
  declareProperty( "PidUseMuc", m_useMuc = false );
 
  declareProperty( "PidwithProbability", m_probability = true );
  declareProperty( "PidProbCut", m_PidProbCut = 0.001 );
  declareProperty( "RejectKaon", m_rejectKaon = true );
  declareProperty( "RejectProton", m_rejectProton = false );
 
  declareProperty( "PidwithLikelihood", m_likelihood = false );
  declareProperty( "PidwithNeuronNetwork", m_neuronNetwork = false );
  declareProperty( "NeuronValue", m_neuronValCut = { 1.5, 2.5 } );
}

Referenced by LocalPionSelector().

Member Function Documentation

operator()()

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

override

Definition at line 47 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;
  StatusCode    sc = serviceLocator()->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;
    StatusCode     sc = serviceLocator()->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;
}

setpidtype()

void LocalPionSelector::setpidtype ( int type )

inlineoverride

Definition at line 15 of file LocalPionSelector.h.

{ m_pidtype = type; }

---

The documentation for this class was generated from the following files:

• LocalPionSelector.h
• LocalPionSelector.cxx