lar_content::PcaShowerParticleBuildingAlgorithm Class Reference

PcaShowerParticleBuildingAlgorithm class. More...

#include <PcaShowerParticleBuildingAlgorithm.h>

Inheritance diagram for lar_content::PcaShowerParticleBuildingAlgorithm:

Classes
class	Factory
	Factory class for instantiating algorithm. More...

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

Private Member Functions
void	CreatePfo (const pandora::ParticleFlowObject *const pInputPfo, const pandora::ParticleFlowObject *&pOutputPfo) const
	Create specialised Pfo from an generic input Pfo. More...
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
unsigned int	m_layerFitHalfWindow

Additional Inherited Members
Protected Member Functions inherited from lar_content::CustomParticleCreationAlgorithm
virtual pandora::StatusCode	Run ()

Detailed Description

PcaShowerParticleBuildingAlgorithm class.

Definition at line 21 of file PcaShowerParticleBuildingAlgorithm.h.

Constructor & Destructor Documentation

lar_content::PcaShowerParticleBuildingAlgorithm::PcaShowerParticleBuildingAlgorithm

(

)

Default constructor.

Definition at line 26 of file PcaShowerParticleBuildingAlgorithm.cc.

                                                                       : m_layerFitHalfWindow(20)
{
}

Member Function Documentation

void lar_content::PcaShowerParticleBuildingAlgorithm::CreatePfo ( const pandora::ParticleFlowObject *const  pInputPfo, const pandora::ParticleFlowObject *&  pOutputPfo  ) const

privatevirtual

Create specialised Pfo from an generic input Pfo.

Parameters

pInputPfo	the address of the input Pfo
pOutputPfo	the address of the output Pfo

Implements lar_content::CustomParticleCreationAlgorithm.

Definition at line 32 of file PcaShowerParticleBuildingAlgorithm.cc.

{
    try
    {
        // In cosmic mode, build showers from all daughter pfos, otherwise require that pfo is shower-like
        if (LArPfoHelper::IsNeutrinoFinalState(pInputPfo))
        {
            if (!LArPfoHelper::IsShower(pInputPfo))
                return;
        }
        else
        {
            if (LArPfoHelper::IsFinalState(pInputPfo))
                return;

            if (LArPfoHelper::IsNeutrino(pInputPfo))
                return;
        }

        // Need an input vertex to provide a shower propagation direction
        const Vertex *const pInputVertex = LArPfoHelper::GetVertex(pInputPfo);

        // Run the PCA analysis
        const LArShowerPCA showerPCA(LArPfoHelper::GetPrincipalComponents(pInputPfo, pInputVertex));

        // Build a new pfo
        LArShowerPfoFactory pfoFactory;
        LArShowerPfoParameters pfoParameters;
        pfoParameters.m_particleId = (LArPfoHelper::IsShower(pInputPfo) ? pInputPfo->GetParticleId() : E_MINUS);
        pfoParameters.m_charge = PdgTable::GetParticleCharge(pfoParameters.m_particleId.Get());
        pfoParameters.m_mass = PdgTable::GetParticleMass(pfoParameters.m_particleId.Get());
        pfoParameters.m_energy = 0.f;
        pfoParameters.m_momentum = pInputPfo->GetMomentum();
        pfoParameters.m_propertiesToAdd = pInputPfo->GetPropertiesMap();
        pfoParameters.m_showerVertex = pInputVertex->GetPosition();
        pfoParameters.m_showerCentroid = showerPCA.GetCentroid();
        pfoParameters.m_showerDirection = showerPCA.GetPrimaryAxis();
        pfoParameters.m_showerSecondaryVector = showerPCA.GetSecondaryAxis();
        pfoParameters.m_showerTertiaryVector = showerPCA.GetTertiaryAxis();
        pfoParameters.m_showerEigenValues = showerPCA.GetEigenValues();
        pfoParameters.m_showerLength = showerPCA.GetAxisLengths();
        pfoParameters.m_showerOpeningAngle =
            (showerPCA.GetPrimaryLength() > 0.f ? std::atan(showerPCA.GetSecondaryLength() / showerPCA.GetPrimaryLength()) : 0.f);

        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ParticleFlowObject::Create(*this, pfoParameters, pOutputPfo, pfoFactory));

        const LArShowerPfo *const pLArPfo = dynamic_cast<const LArShowerPfo *>(pOutputPfo);

        if (!pLArPfo)
            throw StatusCodeException(STATUS_CODE_FAILURE);

        // Build a new vertex - TODO: tune vertex position based on PCA results
        const Vertex *pOutputVertex = nullptr;

        PandoraContentApi::Vertex::Parameters vtxParameters;
        vtxParameters.m_position = pInputVertex->GetPosition();
        vtxParameters.m_vertexLabel = pInputVertex->GetVertexLabel();
        vtxParameters.m_vertexType = pInputVertex->GetVertexType();

        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Vertex::Create(*this, vtxParameters, pOutputVertex));
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToPfo<Vertex>(*this, pOutputPfo, pOutputVertex));
    }
    catch (StatusCodeException &statusCodeException)
    {
        if (STATUS_CODE_FAILURE == statusCodeException.GetStatusCode())
            throw statusCodeException;
    }
}

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

privatevirtual

Reimplemented from lar_content::CustomParticleCreationAlgorithm.

Definition at line 103 of file PcaShowerParticleBuildingAlgorithm.cc.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "LayerFitHalfWindow", m_layerFitHalfWindow));

    return CustomParticleCreationAlgorithm::ReadSettings(xmlHandle);
}

Member Data Documentation

unsigned int lar_content::PcaShowerParticleBuildingAlgorithm::m_layerFitHalfWindow

private

Definition at line 43 of file PcaShowerParticleBuildingAlgorithm.h.

---

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

• PcaShowerParticleBuildingAlgorithm.h
• PcaShowerParticleBuildingAlgorithm.cc