lar_content::TestBeamParticleCreationAlgorithm Class Reference

TestBeamParticleCreationAlgorithm class. More...

#include "TestBeamParticleCreationAlgorithm.h"

Inheritance diagram for lar_content::TestBeamParticleCreationAlgorithm:

Public Member Functions
	TestBeamParticleCreationAlgorithm ()
	Constructor. More...

Private Member Functions
pandora::StatusCode	Run ()
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
std::string	m_pfoListName
	Input pfo list name. More...
std::string	m_vertexListName
	Input vertex list name. More...
bool	m_keepInteractionVertex
	Retain the vertex for the test beam particle at the low z point. More...
bool	m_keepStartVertex
	Retain the vertex for the test beam particle at the interaction point. More...

Detailed Description

TestBeamParticleCreationAlgorithm class.

Definition at line 19 of file TestBeamParticleCreationAlgorithm.h.

Constructor & Destructor Documentation

lar_content::TestBeamParticleCreationAlgorithm::TestBeamParticleCreationAlgorithm

(

)

Constructor.

Definition at line 22 of file TestBeamParticleCreationAlgorithm.cc.

                                                                     :
    m_pfoListName(""),
    m_vertexListName(""),
    m_keepInteractionVertex(false),
    m_keepStartVertex(true)
{
}

Member Function Documentation

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

private

Definition at line 130 of file TestBeamParticleCreationAlgorithm.cc.

References m_keepInteractionVertex, m_keepStartVertex, m_pfoListName, and m_vertexListName.

{
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "PfoListName", m_pfoListName));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
        "KeepInteractionVertex", m_keepInteractionVertex));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle,
        "KeepStartVertex", m_keepStartVertex));

    if (m_keepStartVertex)
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "VertexListName", m_vertexListName));

    if (m_keepInteractionVertex == m_keepStartVertex)
    {
        std::cout << "TestBeamParticleCreationAlgorithm::ReadSettings - must persist one vertex per pfo" << std::endl;
        return STATUS_CODE_INVALID_PARAMETER;
    }

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TestBeamParticleCreationAlgorithm::Run ( )

private

Definition at line 32 of file TestBeamParticleCreationAlgorithm.cc.

References lar_content::LArPfoHelper::GetCaloHits(), lar_content::LArPfoHelper::GetVertex(), lar_content::LArPfoHelper::IsNeutrino(), m_keepInteractionVertex, m_keepStartVertex, m_pfoListName, m_vertexListName, and max.

{
    const PfoList *pPfoList(nullptr);
    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, m_pfoListName, pPfoList));

    PfoList neutrinoPfos;

    for (const Pfo *const pPfo : *pPfoList)
    {
        if (!LArPfoHelper::IsNeutrino(pPfo))
            continue;

        // ATTN: If the test beam score is not set in the neutrino pfo, set the score to 1 as it has been reconstructed under the neutrino hypothesis
        const PropertiesMap properties(pPfo->GetPropertiesMap());
        PropertiesMap::const_iterator propertiesIter(properties.find("TestBeamScore"));
        const float testBeamScore(propertiesIter != properties.end() ? (*propertiesIter).second : 1.f);

        const PfoList &daughterList(pPfo->GetDaughterPfoList());
        const Pfo *pPrimaryPfo(nullptr);
        CartesianVector positionMinZCaloHit(std::numeric_limits<float>::max(), std::numeric_limits<float>::max(), std::numeric_limits<float>::max());

        for (const Pfo *const pDaughterPfo : daughterList)
        {
            CaloHitList collectedHits;
            LArPfoHelper::GetCaloHits(pDaughterPfo, TPC_3D, collectedHits);

            for (const CaloHit *const pCaloHit : collectedHits)
            {
                if (pCaloHit->GetPositionVector().GetZ() < positionMinZCaloHit.GetZ())
                {
                    positionMinZCaloHit = pCaloHit->GetPositionVector();
                    pPrimaryPfo = pDaughterPfo;
                }
            }
        }

        for (const Pfo *const pPrimaryDaughterPfo : daughterList)
        {
            if (pPrimaryDaughterPfo == pPrimaryPfo)
                continue;

            PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SetPfoParentDaughterRelationship(*this, pPrimaryPfo, pPrimaryDaughterPfo));
        }

        PandoraContentApi::ParticleFlowObject::Metadata pfoMetadata;
        pfoMetadata.m_propertiesToAdd["IsTestBeam"] = 1.f;
        pfoMetadata.m_propertiesToAdd["TestBeamScore"] = testBeamScore;

        // ATTN: If the primary pfo is shower like, the target beam particle is most likely an electron/positron.  If the primary pfo is track like, the target
        // beam particle is most likely a pion as pion interactions are more frequent than proton, kaon and muon interactions in the CERN test beam.
        if (std::abs(pPrimaryPfo->GetParticleId()) != E_MINUS)
            pfoMetadata.m_particleId = PI_PLUS;

        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ParticleFlowObject::AlterMetadata(*this, pPrimaryPfo, pfoMetadata));

        try
        {
            if (m_keepStartVertex)
            {
                if (!m_keepInteractionVertex)
                {
                    const Vertex *const pVertex(LArPfoHelper::GetVertex(pPrimaryPfo));
                    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::RemoveFromPfo(*this, pPrimaryPfo, pVertex));
                    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Delete<Vertex>(*this, pVertex));
                }

                std::string vertexListName;
                const VertexList *pVertexList(nullptr);
                PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::CreateTemporaryListAndSetCurrent(*this, pVertexList, vertexListName));

                PandoraContentApi::Vertex::Parameters parameters;
                parameters.m_position = positionMinZCaloHit;
                parameters.m_vertexLabel = VERTEX_START;
                parameters.m_vertexType = VERTEX_3D;

                const Vertex *pVertex(nullptr);
                PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Vertex::Create(*this, parameters, pVertex));
                PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToPfo(*this, pPrimaryPfo, pVertex));
                PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList<Vertex>(*this, m_vertexListName));
                PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::ReplaceCurrentList<Vertex>(*this, m_vertexListName));
            }
        }
        catch (const StatusCodeException &)
        {
            std::cout << "TestBeamParticleCreationAlgorithm::Run - unable to modify test beam particle vertex" << std::endl;
        }

        neutrinoPfos.push_back(pPfo);
    }

    for (const Pfo *const pNeutrinoPfo : neutrinoPfos)
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Delete<Pfo>(*this, pNeutrinoPfo));

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

bool lar_content::TestBeamParticleCreationAlgorithm::m_keepInteractionVertex

private

Retain the vertex for the test beam particle at the low z point.

Definition at line 34 of file TestBeamParticleCreationAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::TestBeamParticleCreationAlgorithm::m_keepStartVertex

private

Retain the vertex for the test beam particle at the interaction point.

Definition at line 35 of file TestBeamParticleCreationAlgorithm.h.

Referenced by ReadSettings(), and Run().

std::string lar_content::TestBeamParticleCreationAlgorithm::m_pfoListName

private

Input pfo list name.

Definition at line 32 of file TestBeamParticleCreationAlgorithm.h.

Referenced by ReadSettings(), and Run().

std::string lar_content::TestBeamParticleCreationAlgorithm::m_vertexListName

private

Input vertex list name.

Definition at line 33 of file TestBeamParticleCreationAlgorithm.h.

Referenced by ReadSettings(), and Run().

---

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

• larpandoracontent/v03_13_02/source/larpandoracontent/LArThreeDReco/LArEventBuilding/TestBeamParticleCreationAlgorithm.h
• larpandoracontent/v03_13_02/source/larpandoracontent/LArThreeDReco/LArEventBuilding/TestBeamParticleCreationAlgorithm.cc