lar_content::CheatingNeutrinoDaughterVerticesAlgorithm Class Reference

CheatingNeutrinoDaughterVerticesAlgorithm::Algorithm class. More...

#include "CheatingNeutrinoDaughterVerticesAlgorithm.h"

Inheritance diagram for lar_content::CheatingNeutrinoDaughterVerticesAlgorithm:

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

Private Member Functions
pandora::StatusCode	Run ()
void	GetMCPrimaryMap (LArMCParticleHelper::MCRelationMap &mcPrimaryMap) const
	Get the mapping from mc particle to primary, only required if collapsed mc particle hierarchy specified. More...
void	ProcessRecoNeutrinos (const pandora::PfoList &neutrinoPfos, const LArMCParticleHelper::MCRelationMap &mcPrimaryMap) const
	Process the list of reconstructed neutrinos. More...
void	ProcessDaughterPfo (const pandora::ParticleFlowObject *const pDaughterPfo, const LArMCParticleHelper::MCRelationMap &mcPrimaryMap) const
	Process a daughter pfo. More...
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
bool	m_collapseToPrimaryMCParticles
	Whether to collapse mc particle hierarchies to primary particles. More...
std::string	m_mcParticleListName
	The mc particle list name, required if want to collapse mc particle hierarchy. More...
std::string	m_neutrinoListName
	The input list of pfo list names. More...
std::string	m_vertexListName
	The name of the output cosmic-ray vertex list. More...

Detailed Description

CheatingNeutrinoDaughterVerticesAlgorithm::Algorithm class.

Definition at line 21 of file CheatingNeutrinoDaughterVerticesAlgorithm.h.

Constructor & Destructor Documentation

lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::CheatingNeutrinoDaughterVerticesAlgorithm

(

)

Default constructor.

Definition at line 20 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

                                                                                     :
    m_collapseToPrimaryMCParticles(false)
{
}

Member Function Documentation

void lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::GetMCPrimaryMap ( LArMCParticleHelper::MCRelationMap &  mcPrimaryMap ) const

private

Get the mapping from mc particle to primary, only required if collapsed mc particle hierarchy specified.

Parameters

mcPrimaryMap

to receive the mapping from mc particle to primary

Definition at line 53 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

References lar_content::LArMCParticleHelper::GetMCPrimaryMap(), m_collapseToPrimaryMCParticles, and m_mcParticleListName.

Referenced by Run().

{
    if (m_collapseToPrimaryMCParticles)
    {
        const MCParticleList *pMCParticleList(nullptr);
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, m_mcParticleListName, pMCParticleList));

        LArMCParticleHelper::GetMCPrimaryMap(pMCParticleList, mcPrimaryMap);
    }
}

void lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::ProcessDaughterPfo ( const pandora::ParticleFlowObject *const  pDaughterPfo, const LArMCParticleHelper::MCRelationMap &  mcPrimaryMap  ) const

private

Process a daughter pfo.

Parameters

pDaughterPfo	the address of a daughter pfo
mcPrimaryMap	the mapping from mc particle to primary, only required if collapsed mc particle hierarchy specified

Definition at line 93 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

References lar_content::LArMCParticleHelper::GetMainMCParticle(), m_collapseToPrimaryMCParticles, and m_vertexListName.

Referenced by ProcessRecoNeutrinos().

{
    const MCParticle *pMCParticle(LArMCParticleHelper::GetMainMCParticle(pDaughterPfo));

    if (m_collapseToPrimaryMCParticles)
    {
        LArMCParticleHelper::MCRelationMap::const_iterator primaryIter = mcPrimaryMap.find(pMCParticle);

        if (mcPrimaryMap.end() == primaryIter)
            throw StatusCodeException(STATUS_CODE_NOT_FOUND);

        pMCParticle = primaryIter->second;
    }

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

    PandoraContentApi::Vertex::Parameters parameters;
    parameters.m_position = pMCParticle->GetVertex();
    parameters.m_vertexLabel = VERTEX_INTERACTION;
    parameters.m_vertexType = VERTEX_3D;

    const Vertex *pVertex(nullptr);
    PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::Vertex::Create(*this, parameters, pVertex));

    if (!pVertexList->empty())
    {
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList<Vertex>(*this, m_vertexListName));
        PANDORA_THROW_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::AddToPfo<Vertex>(*this, pDaughterPfo, pVertex));
    }
}

void lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::ProcessRecoNeutrinos ( const pandora::PfoList &  neutrinoPfos, const LArMCParticleHelper::MCRelationMap &  mcPrimaryMap  ) const

private

Process the list of reconstructed neutrinos.

Parameters

neutrinoPfos	the list of neutrino pfos
mcPrimaryMap	the mapping from mc particle to primary, only required if collapsed mc particle hierarchy specified

Definition at line 66 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

References lar_content::LArPfoHelper::GetAllDownstreamPfos(), and ProcessDaughterPfo().

Referenced by Run().

{
    for (const ParticleFlowObject *const pNeutrinoPfo : neutrinoPfos)
    {
        PfoList daughterPfos;
        LArPfoHelper::GetAllDownstreamPfos(pNeutrinoPfo, daughterPfos);

        PfoList::iterator neutrinoIter(std::find(daughterPfos.begin(), daughterPfos.end(), pNeutrinoPfo));

        if (daughterPfos.end() != neutrinoIter)
            daughterPfos.erase(neutrinoIter);

        for (const ParticleFlowObject *const pDaughterPfo : daughterPfos)
        {
            try
            {
                this->ProcessDaughterPfo(pDaughterPfo, mcPrimaryMap);
            }
            catch (const StatusCodeException &)
            {
            }
        }
    }
}

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

private

Definition at line 129 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

References m_collapseToPrimaryMCParticles, m_mcParticleListName, m_neutrinoListName, and m_vertexListName.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "CollapseToPrimaryMCParticles", m_collapseToPrimaryMCParticles));

    if (m_collapseToPrimaryMCParticles)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "MCParticleListName", m_mcParticleListName));
    }

    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NeutrinoPfoListName", m_neutrinoListName));

    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "OutputVertexListName", m_vertexListName));

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::Run ( )

private

Definition at line 27 of file CheatingNeutrinoDaughterVerticesAlgorithm.cc.

References GetMCPrimaryMap(), lar_content::LArPfoHelper::GetRecoNeutrinos(), m_neutrinoListName, and ProcessRecoNeutrinos().

{
    const PfoList *pPfoList(nullptr);
    PANDORA_THROW_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_INITIALIZED, !=, PandoraContentApi::GetList(*this, m_neutrinoListName, pPfoList));

    if (!pPfoList)
    {
        if (PandoraContentApi::GetSettings(*this)->ShouldDisplayAlgorithmInfo())
            std::cout << "CheatingNeutrinoDaughterVerticesAlgorithm: pfo list unavailable." << std::endl;

        return STATUS_CODE_SUCCESS;
    }

    LArMCParticleHelper::MCRelationMap mcPrimaryMap;
    this->GetMCPrimaryMap(mcPrimaryMap);

    PfoList neutrinoPfos;
    LArPfoHelper::GetRecoNeutrinos(pPfoList, neutrinoPfos);

    this->ProcessRecoNeutrinos(neutrinoPfos, mcPrimaryMap);

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

bool lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::m_collapseToPrimaryMCParticles

private

Whether to collapse mc particle hierarchies to primary particles.

Definition at line 57 of file CheatingNeutrinoDaughterVerticesAlgorithm.h.

Referenced by GetMCPrimaryMap(), ProcessDaughterPfo(), and ReadSettings().

std::string lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::m_mcParticleListName

private

The mc particle list name, required if want to collapse mc particle hierarchy.

Definition at line 58 of file CheatingNeutrinoDaughterVerticesAlgorithm.h.

Referenced by GetMCPrimaryMap(), and ReadSettings().

std::string lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::m_neutrinoListName

private

The input list of pfo list names.

Definition at line 60 of file CheatingNeutrinoDaughterVerticesAlgorithm.h.

Referenced by ReadSettings(), and Run().

std::string lar_content::CheatingNeutrinoDaughterVerticesAlgorithm::m_vertexListName

private

The name of the output cosmic-ray vertex list.

Definition at line 61 of file CheatingNeutrinoDaughterVerticesAlgorithm.h.

Referenced by ProcessDaughterPfo(), and ReadSettings().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArCheating/CheatingNeutrinoDaughterVerticesAlgorithm.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArCheating/CheatingNeutrinoDaughterVerticesAlgorithm.cc