MvaVertexSelectionAlgorithm class. More...

#include <MvaVertexSelectionAlgorithm.h>

Inheritance diagram for lar_content::MvaVertexSelectionAlgorithm< T >:

Public Member Functions
	MvaVertexSelectionAlgorithm ()
	Default constructor. More...

Public Member Functions inherited from lar_content::TrainedVertexSelectionAlgorithm
	TrainedVertexSelectionAlgorithm ()
	Default constructor. More...

Public Member Functions inherited from lar_content::VertexSelectionBaseAlgorithm
	VertexSelectionBaseAlgorithm ()
	Default constructor. More...

Protected Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Protected Member Functions inherited from lar_content::TrainedVertexSelectionAlgorithm
void	CalculateShowerClusterList (const pandora::ClusterList &inputClusterList, ShowerClusterList &showerClusterList) const
	Calculate the shower cluster map for a cluster list. More...

void	GetShowerLikeClusterEndPoints (const pandora::ClusterList &clusterList, pandora::ClusterList &showerLikeClusters, ClusterEndPointsMap &clusterEndPointsMap) const
	Add the endpoints of any shower-like clusters to the map. More...

void	PopulateKdTree (const pandora::ClusterList &clusterList, HitKDTree2D &kdTree, HitToClusterMap &hitToClusterMap) const
	Populate kd tree with information about hits in a provided list of clusters. More...

bool	AddClusterToShower (const ClusterEndPointsMap &clusterEndPointsMap, pandora::ClusterList &availableShowerLikeClusters, const pandora::Cluster *const pCluster, pandora::ClusterList &showerCluster) const
	Try to add an available cluster to a given shower cluster, using shower clustering approximation. More...

bool	AddClusterToShower (HitKDTree2D &kdTree, const HitToClusterMap &hitToClusterMap, pandora::ClusterList &availableShowerLikeClusters, const pandora::Cluster *const pCluster, pandora::ClusterList &showerCluster) const
	Try to add an available cluster to a given shower cluster, using cluster hit positions cached in kd tree. More...

EventFeatureInfo	CalculateEventFeatures (const pandora::ClusterList &clusterListU, const pandora::ClusterList &clusterListV, const pandora::ClusterList &clusterListW, const pandora::VertexVector &vertexVector) const
	Calculate the event parameters. More...

void	IncrementShoweryParameters (const pandora::ClusterList &clusterList, unsigned int &nShoweryHits, unsigned int &nHits, float &eventEnergy) const
	Increment the showery hit parameters for a cluster list. More...

bool	IsClusterShowerLike (const pandora::Cluster *const pCluster) const
	Find whether a cluster is shower-like. More...

void	GetLegacyEventShapeFeatures (const pandora::ClusterList &clusterList, float &eventVolume, float &longitudinality) const
	Get the event shape features. More...

void	GetEventShapeFeatures (const ClusterListMap &clusterListMap, float &eventArea, float &longitudinality) const
	Get the event shape features. More...

void	Get2DSpan (const pandora::ClusterList &clusterList, float &xSpan, float &zSpan) const
	Get the coordinate span in one view. More...

void	UpdateSpanCoordinate (const float minPositionCoord, const float maxPositionCoord, pandora::InputFloat &minCoord, pandora::InputFloat &maxCoord) const
	Update the min/max coordinate spans. More...

float	GetCoordinateSpan (const pandora::InputFloat &minCoord, const pandora::InputFloat &maxCoord) const
	Get the coordinate span. More...

void	AddEventFeaturesToVector (const EventFeatureInfo &eventFeatureInfo, LArMvaHelper::MvaFeatureVector &featureVector) const
	Add the event features to a vector in the correct order. More...

void	PopulateVertexFeatureInfoMap (const BeamConstants &beamConstants, const ClusterListMap &clusterListMap, const SlidingFitDataListMap &slidingFitDataListMap, const ShowerClusterListMap &showerClusterListMap, const KDTreeMap &kdTreeMap, const pandora::Vertex *const pVertex, VertexFeatureInfoMap &vertexFeatureInfoMap) const
	Populate the vertex feature info map for a given vertex. More...

void	PopulateInitialScoreList (VertexFeatureInfoMap &vertexFeatureInfoMap, const pandora::Vertex *const pVertex, VertexScoreList &initialScoreList) const
	Populate the initial vertex score list for a given vertex. More...

void	GetBestRegionVertices (VertexScoreList &initialScoreList, pandora::VertexVector &bestRegionVertices) const
	Get the list of top-N separated vertices. More...

void	ProduceTrainingSets (const pandora::VertexVector &vertexVector, const pandora::VertexVector &bestRegionVertices, VertexFeatureInfoMap &vertexFeatureInfoMap, const LArMvaHelper::MvaFeatureVector &eventFeatureList, const KDTreeMap &kdTreeMap) const
	Produce the region and vertex training sets. More...

void	CalculateRPhiScores (pandora::VertexVector &vertexVector, VertexFeatureInfoMap &vertexFeatureInfoMap, const KDTreeMap &kdTreeMap) const
	Calculate the r/phi scores for the vertices in a vector, possibly erasing those that fail the fast score test. More...

std::string	GetInteractionType () const
	Get the interaction type string. More...

const pandora::Vertex *	ProduceTrainingExamples (const pandora::VertexVector &vertexVector, const VertexFeatureInfoMap &vertexFeatureInfoMap, std::bernoulli_distribution &coinFlip, std::mt19937 &generator, const std::string &interactionType, const std::string &trainingOutputFile, const LArMvaHelper::MvaFeatureVector &eventFeatureList, const KDTreeMap &kdTreeMap, const float maxRadius, const bool useRPhi) const
	Produce a set of training examples for a binary classifier. More...

void	GetBestVertex (const pandora::VertexVector &vertexVector, const pandora::Vertex *&pBestVertex, float &bestVertexDr) const
	Use the MC information to get the best vertex from a list. More...

void	GetSharedFeatures (const pandora::Vertex const pVertex1, const pandora::Vertex const pVertex2, const KDTreeMap &kdTreeMap, float &separation, float &axisHits) const
	Calculates the shared features of a pair of vertex candidates. More...

void	IncrementSharedAxisValues (const pandora::CartesianVector pos1, const pandora::CartesianVector pos2, HitKDTree2D &kdTree, float &axisHits) const
	Increments the axis hits information for one view. More...

bool	IsHitInBox (const pandora::CartesianVector &hitPos, const pandora::CartesianVector &point1, const pandora::CartesianVector &point2, const pandora::CartesianVector &point3, const pandora::CartesianVector &point4) const
	Determines whether a hit lies within the box defined by four other positions. More...

void	AddVertexFeaturesToVector (const VertexFeatureInfo &vertexFeatureInfo, LArMvaHelper::MvaFeatureVector &featureVector, const bool useRPhi) const
	Add the vertex features to a vector in the correct order. More...

void	AddSharedFeaturesToVector (const VertexSharedFeatureInfo &vertexSharedFeatureInfo, LArMvaHelper::MvaFeatureVector &featureVector) const
	Add the shared features to a vector in the correct order. More...

void	PopulateFinalVertexScoreList (const VertexFeatureInfoMap &vertexFeatureInfoMap, const pandora::Vertex *const pFavouriteVertex, const pandora::VertexVector &vertexVector, VertexScoreList &finalVertexScoreList) const
	Populate the final vertex score list using the r/phi score to find the best vertex in the vicinity. More...

pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Protected Member Functions inherited from lar_content::VertexSelectionBaseAlgorithm
virtual void	FilterVertexList (const pandora::VertexList *const pInputVertexList, HitKDTree2D &kdTreeU, HitKDTree2D &kdTreeV, HitKDTree2D &kdTreeW, pandora::VertexVector &filteredVertices) const
	Filter the input list of vertices to obtain a reduced number of vertex candidates. More...

virtual void	GetBeamConstants (const pandora::VertexVector &vertexVector, BeamConstants &beamConstants) const
	Get the beam score constants for a provided list of candidate vertices. More...

void	GetClusterLists (const pandora::StringVector &inputClusterListNames, pandora::ClusterList &clusterListU, pandora::ClusterList &clusterListV, pandora::ClusterList &clusterListW) const
	Get the cluster lists. More...

void	CalculateClusterSlidingFits (const pandora::ClusterList &inputClusterList, const unsigned int minClusterCaloHits, const unsigned int slidingFitWindow, SlidingFitDataList &slidingFitDataList) const
	Calculate the cluster sliding fits. More...

float	GetBeamDeweightingScore (const BeamConstants &beamConstants, const pandora::Vertex *const pVertex) const
	Get the beam deweighting score for a vertex. More...

bool	IsBeamModeOn () const
	Whether algorithm is running in beam mode, assuming neutrinos travel in positive z-direction. More...

float	GetVertexEnergy (const pandora::Vertex *const pVertex, const KDTreeMap &kdTreeMap) const
	Calculate the energy of a vertex candidate by summing values from all three planes. More...

float	VertexHitEnergy (const pandora::Vertex *const pVertex, const pandora::HitType hitType, HitKDTree2D &kdTree) const
	Finds the energy of the nearest hit to the vertex candidate in this view. More...

pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Member Functions
void	GetVertexScoreList (const pandora::VertexVector &vertexVector, const BeamConstants &beamConstants, HitKDTree2D &kdTreeU, HitKDTree2D &kdTreeV, HitKDTree2D &kdTreeW, VertexScoreList &vertexScoreList) const
	Get the vertex score list. More...

const pandora::Vertex *	CompareVertices (const pandora::VertexVector &vertexVector, const VertexFeatureInfoMap &vertexFeatureInfoMap, const LArMvaHelper::MvaFeatureVector &eventFeatureList, const KDTreeMap &kdTreeMap, const T &t, const bool useRPhi) const
	Used a binary classifier to compare a set of vertices and pick the best one. More...

Private Attributes
std::string	m_filePathEnvironmentVariable
	The environment variable providing a list of paths to mva files. More...

std::string	m_mvaFileName
	The mva file name. More...

std::string	m_regionMvaName
	The name of the region mva to find. More...

std::string	m_vertexMvaName
	The name of the vertex mva to find. More...

T	m_mvaRegion
	The region mva. More...

T	m_mvaVertex
	The vertex mva. More...

Additional Inherited Members
Public Types inherited from lar_content::TrainedVertexSelectionAlgorithm
typedef std::map< const pandora::Vertex *const, VertexFeatureInfo >	VertexFeatureInfoMap

Public Types inherited from lar_content::VertexSelectionBaseAlgorithm
typedef std::vector< VertexScore >	VertexScoreList

typedef std::vector< SlidingFitData >	SlidingFitDataList

typedef std::vector< ShowerCluster >	ShowerClusterList

typedef KDTreeNodeInfoT< const pandora::CaloHit *, 2 >	HitKDNode2D

typedef std::vector< HitKDNode2D >	HitKDNode2DList

typedef KDTreeLinkerAlgo< const pandora::CaloHit *, 2 >	HitKDTree2D

typedef std::map< pandora::HitType, const pandora::ClusterList & >	ClusterListMap
	Map array of cluster lists for passing to tools. More...

typedef std::map< pandora::HitType, const SlidingFitDataList >	SlidingFitDataListMap
	Map of sliding fit data lists for passing to tools. More...

typedef std::map< pandora::HitType, const ShowerClusterList >	ShowerClusterListMap
	Map of shower cluster lists for passing to tools. More...

typedef std::map< pandora::HitType, const std::reference_wrapper< HitKDTree2D > >	KDTreeMap
	Map array of hit kd trees for passing to tools. More...

typedef MvaFeatureTool< const VertexSelectionBaseAlgorithm const, const pandora::Vertex const, const SlidingFitDataListMap &, const ClusterListMap &, const KDTreeMap &, const ShowerClusterListMap &, const float, float & >	VertexFeatureTool
	The base type for the vertex feature tools. More...

Protected Types inherited from lar_content::TrainedVertexSelectionAlgorithm
typedef std::pair< pandora::CartesianVector, pandora::CartesianVector >	ClusterEndPoints

typedef std::map< const pandora::Cluster *const, ClusterEndPoints >	ClusterEndPointsMap

typedef std::vector< LArMvaHelper::MvaFeatureVector >	FeatureListVector

typedef std::vector< pandora::VertexVector >	VectorOfVertexVectors

typedef KDTreeLinkerAlgo< const pandora::CaloHit *, 2 >	HitKDTree2D

typedef KDTreeNodeInfoT< const pandora::CaloHit *, 2 >	HitKDNode2D

typedef std::vector< HitKDNode2D >	HitKDNode2DList

typedef std::unordered_map< const pandora::CaloHit , const pandora::Cluster >	HitToClusterMap

Protected Attributes inherited from lar_content::TrainedVertexSelectionAlgorithm
VertexFeatureTool::FeatureToolVector	m_featureToolVector
	The feature tool vector. More...

bool	m_trainingSetMode
	Whether to train. More...

bool	m_allowClassifyDuringTraining
	Whether classification is allowed during training. More...

float	m_mcVertexXCorrection
	The correction to the x-coordinate of the MC vertex position. More...

std::string	m_trainingOutputFileRegion
	The training output file for the region mva. More...

std::string	m_trainingOutputFileVertex
	The training output file for the vertex mva. More...

std::string	m_mcParticleListName
	The MC particle list for creating training examples. More...

std::string	m_caloHitListName
	The 2D CaloHit list name. More...

pandora::StringVector	m_inputClusterListNames
	The list of cluster list names. More...

unsigned int	m_minClusterCaloHits
	The min number of hits parameter in the energy score. More...

unsigned int	m_slidingFitWindow
	The layer window for the sliding linear fits. More...

float	m_minShowerSpineLength
	The minimum length at which all are considered to be tracks. More...

float	m_beamDeweightingConstant
	The beam deweighting constant for the initial region score list. More...

float	m_localAsymmetryConstant
	The local asymmetry constant for the initial region score list. More...

float	m_globalAsymmetryConstant
	The global asymmetry constant for the initial region score list. More...

float	m_showerAsymmetryConstant
	The shower asymmetry constant for the initial region score list. More...

float	m_energyKickConstant
	The energy kick constant for the initial region score list. More...

float	m_showerClusteringDistance
	The shower clustering distance. More...

unsigned int	m_minShowerClusterHits
	The minimum number of shower cluster hits. More...

bool	m_useShowerClusteringApproximation
	Whether to use the shower clustering distance approximation. More...

float	m_regionRadius
	The radius for a vertex region. More...

float	m_rPhiFineTuningRadius
	The maximum distance the r/phi tune can move a vertex. More...

float	m_maxTrueVertexRadius
	The maximum distance at which a vertex candidate can be considered the 'true' vertex. More...

bool	m_useRPhiFeatureForRegion
	Whether to use the r/phi feature for the region vertex. More...

bool	m_dropFailedRPhiFastScoreCandidates
	Whether to drop candidates that fail the r/phi fast score test. More...

bool	m_testBeamMode
	Test beam mode. More...

bool	m_legacyEventShapes
	Whether to use the old event shapes calculation. More...

bool	m_legacyVariables
	Whether to only use the old variables. More...

Detailed Description

template<typename T>
class lar_content::MvaVertexSelectionAlgorithm< T >

MvaVertexSelectionAlgorithm class.

Definition at line 37 of file MvaVertexSelectionAlgorithm.h.

Constructor & Destructor Documentation

template<typename T >

lar_content::MvaVertexSelectionAlgorithm< T >::MvaVertexSelectionAlgorithm ( )

Default constructor.

Definition at line 36 of file MvaVertexSelectionAlgorithm.cc.

                                                             :
     TrainedVertexSelectionAlgorithm(),
     m_filePathEnvironmentVariable("FW_SEARCH_PATH")
 {
 }

Member Function Documentation

template<typename T >

const pandora::Vertex * lar_content::MvaVertexSelectionAlgorithm< T >::CompareVertices	(	const pandora::VertexVector &	vertexVector,
		const VertexFeatureInfoMap &	vertexFeatureInfoMap,
		const LArMvaHelper::MvaFeatureVector &	eventFeatureList,
		const KDTreeMap &	kdTreeMap,
		const T &	t,
		const bool	useRPhi
	)		const

private

Used a binary classifier to compare a set of vertices and pick the best one.

Parameters

vertexVector	the vector of vertices
vertexFeatureInfoMap	the vertex feature info map
eventFeatureList	the event feature list
kdTreeMap	the map of 2D hit kd trees
t	the mva
useRPhi	whether to include the r/phi feature

Returns: address of the best vertex

Definition at line 127 of file MvaVertexSelectionAlgorithm.cc.

 {
     const Vertex *pBestVertex(vertexVector.front());
     LArMvaHelper::MvaFeatureVector chosenFeatureList;
 
     VertexFeatureInfo chosenVertexFeatureInfo(vertexFeatureInfoMap.at(pBestVertex));
     this->AddVertexFeaturesToVector(chosenVertexFeatureInfo, chosenFeatureList, useRPhi);
 
     for (const Vertex *const pVertex : vertexVector)
     {
         if (pVertex == pBestVertex)
             continue;
 
         LArMvaHelper::MvaFeatureVector featureList;
         VertexFeatureInfo vertexFeatureInfo(vertexFeatureInfoMap.at(pVertex));
         this->AddVertexFeaturesToVector(vertexFeatureInfo, featureList, useRPhi);
 
         if (!m_legacyVariables)
         {
             LArMvaHelper::MvaFeatureVector sharedFeatureList;
             float separation(0.f), axisHits(0.f);
             this->GetSharedFeatures(pVertex, pBestVertex, kdTreeMap, separation, axisHits);
             VertexSharedFeatureInfo sharedFeatureInfo(separation, axisHits);
             this->AddSharedFeaturesToVector(sharedFeatureInfo, sharedFeatureList);
 
             if (LArMvaHelper::Classify(t, eventFeatureList, featureList, chosenFeatureList, sharedFeatureList))
             {
                 pBestVertex = pVertex;
                 chosenFeatureList = featureList;
             }
         }
         else
         {
             if (LArMvaHelper::Classify(t, eventFeatureList, featureList, chosenFeatureList))
             {
                 pBestVertex = pVertex;
                 chosenFeatureList = featureList;
             }
         }
     }
 
     return pBestVertex;
 }

template<typename T >

void lar_content::MvaVertexSelectionAlgorithm< T >::GetVertexScoreList	(	const pandora::VertexVector &	vertexVector,
		const BeamConstants &	beamConstants,
		HitKDTree2D &	kdTreeU,
		HitKDTree2D &	kdTreeV,
		HitKDTree2D &	kdTreeW,
		VertexScoreList &	vertexScoreList
	)		const

privatevirtual

Get the vertex score list.

Parameters

vertexVector	the vector of vertices
beamConstants	the beam constants
kdTreeU	the hit kd tree for the U view
kdTreeV	the hit kd tree for the V view
kdTreeW	the hit kd tree for the W view
vertexScoreList	the vertex score list to fill

Implements lar_content::TrainedVertexSelectionAlgorithm.

Definition at line 45 of file MvaVertexSelectionAlgorithm.cc.

 {
     ClusterList clustersU, clustersV, clustersW;
     this->GetClusterLists(m_inputClusterListNames, clustersU, clustersV, clustersW);
 
     SlidingFitDataList slidingFitDataListU, slidingFitDataListV, slidingFitDataListW;
     this->CalculateClusterSlidingFits(clustersU, m_minClusterCaloHits, m_slidingFitWindow, slidingFitDataListU);
     this->CalculateClusterSlidingFits(clustersV, m_minClusterCaloHits, m_slidingFitWindow, slidingFitDataListV);
     this->CalculateClusterSlidingFits(clustersW, m_minClusterCaloHits, m_slidingFitWindow, slidingFitDataListW);
 
     ShowerClusterList showerClusterListU, showerClusterListV, showerClusterListW;
     this->CalculateShowerClusterList(clustersU, showerClusterListU);
     this->CalculateShowerClusterList(clustersV, showerClusterListV);
     this->CalculateShowerClusterList(clustersW, showerClusterListW);
 
     // Create maps from hit types to objects for passing to feature tools.
     const ClusterListMap clusterListMap{{TPC_VIEW_U, clustersU}, {TPC_VIEW_V, clustersV}, {TPC_VIEW_W, clustersW}};
 
     const SlidingFitDataListMap slidingFitDataListMap{
         {TPC_VIEW_U, slidingFitDataListU}, {TPC_VIEW_V, slidingFitDataListV}, {TPC_VIEW_W, slidingFitDataListW}};
 
     const ShowerClusterListMap showerClusterListMap{{TPC_VIEW_U, showerClusterListU}, {TPC_VIEW_V, showerClusterListV}, {TPC_VIEW_W, showerClusterListW}};
 
     const KDTreeMap kdTreeMap{{TPC_VIEW_U, kdTreeU}, {TPC_VIEW_V, kdTreeV}, {TPC_VIEW_W, kdTreeW}};
 
     // Calculate the event feature list and the vertex feature map.
     EventFeatureInfo eventFeatureInfo(this->CalculateEventFeatures(clustersU, clustersV, clustersW, vertexVector));
 
     LArMvaHelper::MvaFeatureVector eventFeatureList;
     this->AddEventFeaturesToVector(eventFeatureInfo, eventFeatureList);
 
     VertexFeatureInfoMap vertexFeatureInfoMap;
     for (const Vertex *const pVertex : vertexVector)
     {
         this->PopulateVertexFeatureInfoMap(
             beamConstants, clusterListMap, slidingFitDataListMap, showerClusterListMap, kdTreeMap, pVertex, vertexFeatureInfoMap);
     }
 
     // Use a simple score to get the list of vertices representing good regions.
     VertexScoreList initialScoreList;
     for (const Vertex *const pVertex : vertexVector)
         PopulateInitialScoreList(vertexFeatureInfoMap, pVertex, initialScoreList);
 
     VertexVector bestRegionVertices;
     this->GetBestRegionVertices(initialScoreList, bestRegionVertices);
 
     if (m_trainingSetMode)
         this->ProduceTrainingSets(vertexVector, bestRegionVertices, vertexFeatureInfoMap, eventFeatureList, kdTreeMap);
 
     if ((!m_trainingSetMode || m_allowClassifyDuringTraining) && !bestRegionVertices.empty())
     {
         // Use mva to choose the region.
         const Vertex *const pBestRegionVertex(
             this->CompareVertices(bestRegionVertices, vertexFeatureInfoMap, eventFeatureList, kdTreeMap, m_mvaRegion, m_useRPhiFeatureForRegion));
 
         // Get all the vertices in the best region.
         VertexVector regionalVertices{pBestRegionVertex};
         for (const Vertex *const pVertex : vertexVector)
         {
             if (pVertex == pBestRegionVertex)
                 continue;
 
             if ((pBestRegionVertex->GetPosition() - pVertex->GetPosition()).GetMagnitude() < m_regionRadius)
                 regionalVertices.push_back(pVertex);
         }
 
         this->CalculateRPhiScores(regionalVertices, vertexFeatureInfoMap, kdTreeMap);
 
         if (!regionalVertices.empty())
         {
             // Use mva to choose the vertex and then fine-tune using the RPhi score.
             const Vertex *const pBestVertex(
                 this->CompareVertices(regionalVertices, vertexFeatureInfoMap, eventFeatureList, kdTreeMap, m_mvaVertex, true));
             this->PopulateFinalVertexScoreList(vertexFeatureInfoMap, pBestVertex, vertexVector, vertexScoreList);
         }
     }
 }

template<typename T >

StatusCode lar_content::MvaVertexSelectionAlgorithm< T >::ReadSettings ( const pandora::TiXmlHandle xmlHandle )

protected

Definition at line 175 of file MvaVertexSelectionAlgorithm.cc.

 {
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
         XmlHelper::ReadValue(xmlHandle, "FilePathEnvironmentVariable", m_filePathEnvironmentVariable));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MvaFileName", m_mvaFileName));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "RegionMvaName", m_regionMvaName));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "VertexMvaName", m_vertexMvaName));
 
     // ATTN : Need access to base class member variables at this point, so call read settings prior to end of this function
     PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, TrainedVertexSelectionAlgorithm::ReadSettings(xmlHandle));
 
     if ((!m_trainingSetMode || m_allowClassifyDuringTraining))
     {
         if (m_mvaFileName.empty() || m_regionMvaName.empty() || m_vertexMvaName.empty())
         {
             std::cout << "MvaVertexSelectionAlgorithm: MvaFileName, RegionMvaName and VertexMvaName must be set if training set mode is"
                       << "off or we allow classification during training" << std::endl;
             return STATUS_CODE_INVALID_PARAMETER;
         }
 
         const std::string fullMvaFileName(LArFileHelper::FindFileInPath(m_mvaFileName, m_filePathEnvironmentVariable));
         m_mvaRegion.Initialize(fullMvaFileName, m_regionMvaName);
         m_mvaVertex.Initialize(fullMvaFileName, m_vertexMvaName);
     }
 
     return STATUS_CODE_SUCCESS;
 }