lar_content::LArEventTopology::Particle Class Reference

Public Member Functions
	Particle (const pandora::MCParticle *const pRoot)
	Default constructor. More...
	Particle (const pandora::MCParticle *const pRoot, const pandora::CaloHitList &caloHitList)
	Default constructor. More...
virtual	~Particle ()
	Destructor. More...
void	AddChild (Particle *pChild)
	Add a child particle to this particle. More...
void	GetVertices (pandora::CartesianPointVector &vertices) const
	Extract the clear topological vertices from the event. More...
void	Parse (const MCHitMap &mcHitMap)
	Assess this child particle and its descendents to see if it is topologically significant either indepednently or within the context of downstream activity. More...
void	Prune ()
	Prune particles from the hierarchy if they are not deemed topologically significant. More...
const std::string	Print (const MCHitMap &mcHitMap, const std::string &indent) const
	Print the visible hierarchy. More...

Private Attributes
pandora::MCParticleList	m_particles
pandora::CaloHitList	m_caloHits
std::list< Particle * >	m_children
bool	m_fold

Detailed Description

Definition at line 26 of file LArEventTopology.h.

Constructor & Destructor Documentation

lar_content::LArEventTopology::Particle::Particle

(

const pandora::MCParticle *const

pRoot

)

Default constructor.

Referenced by lar_content::LArEventTopology::PruneHierarchy().

lar_content::LArEventTopology::Particle::Particle	(	const pandora::MCParticle *const	pRoot,
		const pandora::CaloHitList &	caloHitList
	)

Default constructor.

Parameters

pRoot	The MC particle
caloHitList	The list of hits associated with this particle

lar_content::LArEventTopology::Particle::~Particle ( )

virtual

Destructor.

Definition at line 149 of file LArEventTopology.cc.

{
    for (Particle *pParticle : m_children)
        delete pParticle;
}

Member Function Documentation

void lar_content::LArEventTopology::Particle::AddChild

(

Particle *

pChild

)

Add a child particle to this particle.

Parameters

pChild

The child particle to add

Definition at line 157 of file LArEventTopology.cc.

Referenced by lar_content::LArEventTopology::Print().

{
    m_children.emplace_back(pChild);
}

void lar_content::LArEventTopology::Particle::GetVertices

(

pandora::CartesianPointVector &

vertices

)

const

Extract the clear topological vertices from the event.

Definition at line 164 of file LArEventTopology.cc.

References util::abs(), and lar_content::LArMCParticleHelper::IsNeutrino().

Referenced by lar_content::LArEventTopology::GetVertices().

{
    const MCParticle *const pMC{m_particles.front()};
    if (LArMCParticleHelper::IsNeutrino(pMC))
    {
        if (std::find(vertices.begin(), vertices.end(), pMC->GetVertex()) == vertices.end())
            vertices.emplace_back(pMC->GetVertex());
    }
    else
    {
        // Consider vertex
        const MCParticle *const pParentMC{pMC->GetParentList().front()};
        if (LArMCParticleHelper::IsNeutrino(pParentMC))
        {
            if (std::find(vertices.begin(), vertices.end(), pMC->GetVertex()) == vertices.end())
                vertices.emplace_back(pMC->GetVertex());
        }
        else
        {
            const CartesianVector &momParent{pParentMC->GetMomentum().GetUnitVector()};
            const CartesianVector &mom{pMC->GetMomentum().GetUnitVector()};
            if ((pMC->GetParticleId() != pParentMC->GetParticleId()) || (mom.GetDotProduct(momParent) < 0.996f))
            {
                // Visible photon vertices occur at the endpoint, not the vertex, so skip photons here
                if (pMC->GetParticleId() != PHOTON && std::find(vertices.begin(), vertices.end(), pMC->GetVertex()) == vertices.end())
                    vertices.emplace_back(pMC->GetVertex());
            }
        }
        // Consider end point
        int nTrackLikeChildren{0};
        for (const Particle *const pChild : m_children)
        {
            const MCParticle *const pChildMC{pChild->m_particles.front()};
            const int pdg{std::abs(pChildMC->GetParticleId())};
            if (!(pdg == PHOTON || pdg == E_MINUS))
                ++nTrackLikeChildren;
        }
        // If only shower children, or if more than 1 track, tag the endpoint as a vertex
        // If  we have exactly 1 track child, then downstream vertex checks will take care of this
        // Also note that the endpoint of a photon actually marks the start of the associated hit deposition
        if (nTrackLikeChildren != 1)
        {
            if (std::find(vertices.begin(), vertices.end(), pMC->GetEndpoint()) == vertices.end())
                vertices.emplace_back(pMC->GetEndpoint());
        }
    }

    for (const Particle *const pChild : m_children)
        pChild->GetVertices(vertices);
}

void lar_content::LArEventTopology::Particle::Parse

(

const MCHitMap &

mcHitMap

)

Assess this child particle and its descendents to see if it is topologically significant either indepednently or within the context of downstream activity.

Parameters

mcHitMap

The map from MC particles to hits

Definition at line 217 of file LArEventTopology.cc.

References util::abs(), and target.

Referenced by lar_content::LArEventTopology::PruneHierarchy().

{
    const MCParticle *const pMC{m_particles.front()};
    const CaloHitList &caloHitList{mcHitMap.find(pMC) != mcHitMap.end() ? mcHitMap.at(pMC) : CaloHitList()};
    int uHits{0}, vHits{0}, wHits{0};
    for (const CaloHit *const pCaloHit : caloHitList)
    {
        switch (pCaloHit->GetHitType())
        {
            case TPC_VIEW_U:
                ++uHits;
                break;
            case TPC_VIEW_V:
                ++vHits;
                break;
            case TPC_VIEW_W:
                ++wHits;
                break;
            default:
                break;
        }
    }

    const int pdg{std::abs(pMC->GetParticleId())};
    const int target{pdg == PHOTON || pdg == E_MINUS ? 8 : 5};
    const int totalHits{uHits + vHits + wHits};
    int goodViews{uHits >= target ? 1 : 0};
    goodViews += vHits >= target ? 1 : 0;
    goodViews += wHits >= target ? 1 : 0;

    if (totalHits < (3 * target) || goodViews < 2)
        m_fold = true;

    for (Particle *const pChild : m_children)
        pChild->Parse(mcHitMap);
}

const std::string lar_content::LArEventTopology::Particle::Print	(	const MCHitMap &	mcHitMap,
		const std::string &	indent
	)		const

Print the visible hierarchy.

Parameters

mcHitMap	The map from MC particles to hits
indent	The level of indent for printing

Returns

The string representation of this particle and its children

Definition at line 273 of file LArEventTopology.cc.

Referenced by lar_content::LArEventTopology::Print().

{
    const MCParticle *const pMC{m_particles.front()};
    const CaloHitList &caloHitList{mcHitMap.find(pMC) != mcHitMap.end() ? mcHitMap.at(pMC) : CaloHitList()};
    int uHits{0}, vHits{0}, wHits{0};
    for (const CaloHit *const pCaloHit : caloHitList)
    {
        switch (pCaloHit->GetHitType())
        {
            case TPC_VIEW_U:
                ++uHits;
                break;
            case TPC_VIEW_V:
                ++vHits;
                break;
            case TPC_VIEW_W:
                ++wHits;
                break;
            default:
                break;
        }
    }

    std::ostringstream os;
    os << indent << pMC->GetParticleId() << ": E: " << pMC->GetEnergy() << " GeV";
    os << " (" << uHits << "," << vHits << "," << wHits << ") " << (m_fold ? "[fold]" : "") << "\n";
    for (const Particle *const pChild : m_children)
        os << pChild->Print(mcHitMap, indent + "  ");

    return os.str();
}

void lar_content::LArEventTopology::Particle::Prune

(

)

Prune particles from the hierarchy if they are not deemed topologically significant.

Definition at line 256 of file LArEventTopology.cc.

Referenced by lar_content::LArEventTopology::PruneHierarchy().

{
    for (Particle *const pChild : m_children)
        pChild->Prune();

    std::list<Particle> toPrune;
    for (auto iter = m_children.begin(); iter != m_children.end();)
    {
        if ((*iter)->m_fold && (*iter)->m_children.empty())
            iter = m_children.erase(iter);
        else
            ++iter;
    }
}

Member Data Documentation

pandora::CaloHitList lar_content::LArEventTopology::Particle::m_caloHits

private

Definition at line 84 of file LArEventTopology.h.

std::list<Particle *> lar_content::LArEventTopology::Particle::m_children

private

Definition at line 85 of file LArEventTopology.h.

bool lar_content::LArEventTopology::Particle::m_fold

private

Definition at line 86 of file LArEventTopology.h.

pandora::MCParticleList lar_content::LArEventTopology::Particle::m_particles

private

Definition at line 83 of file LArEventTopology.h.

---

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

• larpandoracontent/v04_15_01/source/larpandoracontent/LArObjects/LArEventTopology.h
• larpandoracontent/v04_15_01/source/larpandoracontent/LArObjects/LArEventTopology.cc