TwoDParticleCreationAlgorithm class. More...

#include <TwoDParticleCreationAlgorithm.h>

Inheritance diagram for lar_content::TwoDParticleCreationAlgorithm:

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

Private Member Functions
pandora::StatusCode	Run ()

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

pandora::StatusCode	CreatePFOs (const pandora::ClusterList *const pClusterList) const
	Create pfos for provided clusters. More...

Private Attributes
std::string	m_inputClusterListNameU
	The input cluster list name for the U view. More...

std::string	m_inputClusterListNameV
	The input cluster list name for the V view. More...

std::string	m_inputClusterListNameW
	The input cluster list name for the W view. More...

std::string	m_outputPfoListName
	The output pfo list name. More...

unsigned int	m_minHitsInCluster
	Min number of hits for clusters to form pfos. More...

float	m_minClusterEnergy
	Min energy for clusters to form pfos. More...

Detailed Description

TwoDParticleCreationAlgorithm class.

Definition at line 19 of file TwoDParticleCreationAlgorithm.h.

Constructor & Destructor Documentation

lar_content::TwoDParticleCreationAlgorithm::TwoDParticleCreationAlgorithm ( )

Default constructor.

Definition at line 18 of file TwoDParticleCreationAlgorithm.cc.

                                                              : m_minHitsInCluster(5), m_minClusterEnergy(0.f)
 {
 }

Member Function Documentation

StatusCode lar_content::TwoDParticleCreationAlgorithm::CreatePFOs ( const pandora::ClusterList *const pClusterList ) const

private

Create pfos for provided clusters.

Parameters

pClusterList address of the cluster list

Definition at line 68 of file TwoDParticleCreationAlgorithm.cc.

 {
     for (ClusterList::const_iterator iter = pClusterList->begin(), iterEnd = pClusterList->end(); iter != iterEnd; ++iter)
     {
         const Cluster *const pCluster = *iter;
 
         if (pCluster->GetNCaloHits() < m_minHitsInCluster)
             continue;
 
         // TODO Finalize particle id and energy measurement here
         const float clusterEnergy(pCluster->GetElectromagneticEnergy());
 
         if (clusterEnergy < m_minClusterEnergy)
             continue;
 
         PandoraContentApi::ParticleFlowObject::Parameters pfoParameters;
 
         // TODO Finalize particle id here
         const ParticleType particleType(PandoraContentApi::GetPlugins(*this)->GetParticleId()->IsMuon(pCluster) ? MU_MINUS : PHOTON);
 
         const ClusterFitResult &fitToAllHitsResult(pCluster->GetFitToAllHitsResult());
 
         if (!fitToAllHitsResult.IsFitSuccessful())
             continue;
 
         // TODO Check remaining parameters
         pfoParameters.m_particleId = particleType;
         pfoParameters.m_charge = 0;
         pfoParameters.m_mass = 0.;
         pfoParameters.m_energy = clusterEnergy;
         pfoParameters.m_momentum = CartesianVector(fitToAllHitsResult.GetDirection() * clusterEnergy);
         pfoParameters.m_clusterList.push_back(pCluster);
 
         const ParticleFlowObject *pPfo(NULL);
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ParticleFlowObject::Create(*this, pfoParameters, pPfo));
     }
 
     return STATUS_CODE_SUCCESS;
 }

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

private

Definition at line 110 of file TwoDParticleCreationAlgorithm.cc.

 {
     PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "OutputPfoListName", m_outputPfoListName));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ClusterListNameU", m_inputClusterListNameU));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ClusterListNameV", m_inputClusterListNameV));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ClusterListNameW", m_inputClusterListNameW));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinHitsInCluster", m_minHitsInCluster));
 
     PANDORA_RETURN_RESULT_IF_AND_IF(
         STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinClusterEnergy", m_minClusterEnergy));
 
     return STATUS_CODE_SUCCESS;
 }

StatusCode lar_content::TwoDParticleCreationAlgorithm::Run ( )

private

Definition at line 24 of file TwoDParticleCreationAlgorithm.cc.

 {
     const PfoList *pPfoList = nullptr;
     std::string pfoListName;
     PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::CreateTemporaryListAndSetCurrent(*this, pPfoList, pfoListName));
 
     if (!m_inputClusterListNameU.empty())
     {
         const ClusterList *pClusterListU = nullptr;
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, m_inputClusterListNameU, pClusterListU));
 
         if (pClusterListU)
             PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, this->CreatePFOs(pClusterListU));
     }
 
     if (!m_inputClusterListNameV.empty())
     {
         const ClusterList *pClusterListV = nullptr;
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, m_inputClusterListNameV, pClusterListV));
 
         if (pClusterListV)
             PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, this->CreatePFOs(pClusterListV));
     }
 
     if (!m_inputClusterListNameW.empty())
     {
         const ClusterList *pClusterListW = nullptr;
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, m_inputClusterListNameW, pClusterListW));
 
         if (pClusterListW)
             PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, this->CreatePFOs(pClusterListW));
     }
 
     if (!pPfoList->empty())
     {
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList<Pfo>(*this, m_outputPfoListName));
         PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<Pfo>(*this, m_outputPfoListName));
     }
 
     return STATUS_CODE_SUCCESS;
 }