lar_content::ProtoShowerMatchingTool Class Reference

#include "ProtoShowerMatchingTool.h"

Inheritance diagram for lar_content::ProtoShowerMatchingTool:

Public Member Functions
	ProtoShowerMatchingTool ()
	Default constructor. More...
pandora::StatusCode	Run (const ProtoShowerVector &protoShowerVectorU, const ProtoShowerVector &protoShowerVectorV, const ProtoShowerVector &protoShowerVectorW, ProtoShowerMatchVector &protoShowerMatchVector)

Private Member Functions
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)
bool	ArePathwaysConsistent (const ProtoShower &protoShowerU, const ProtoShower &protoShowerV, const ProtoShower &protoShowerW, Consistency &consistency) const
	Determine whether three 2D connection pathways form a consistent 3D connection pathway. More...
bool	AreShowerStartsConsistent (const ProtoShower &protoShowerU, const ProtoShower &protoShowerV, const ProtoShower &protoShowerW) const
	Determine whether three 2D shower start positions correspond to the same 3D shower start position. More...
bool	AreDirectionsConsistent (const ProtoShower &protoShowerU, const ProtoShower &protoShowerV, const ProtoShower &protoShowerW) const
	Determine whether three 2D initial spine directions correspond to the same 3D initial spine direction. More...
bool	AreDirectionsConsistent (const pandora::CartesianVector &nuVertexU, const pandora::CartesianVector &nuVertexV, const pandora::CartesianVector &nuVertexW, const pandora::CartesianVector &directionU, const pandora::CartesianVector &directionV, const pandora::CartesianVector &directionW) const
	Determine whether three 2D initial spine directions correspond to the same 3D initial spine direction. More...

Private Attributes
unsigned int	m_spineSlidingFitWindow
	The shower spine sliding fit window. More...
float	m_maxXSeparation
	The max. drift-coordinate separation between matched 2D shower start positions. More...
float	m_maxSeparation
	The max. average separation between true and projected 2D shower start positions for a match. More...
float	m_xExtrapolation
	Extrapolation distance in the x-direction. More...
float	m_maxAngularDeviation
	The max. opening angle between true and projected 2D initial directions for a match. More...

Detailed Description

Definition at line 19 of file ProtoShowerMatchingTool.h.

Constructor & Destructor Documentation

lar_content::ProtoShowerMatchingTool::ProtoShowerMatchingTool

(

)

Default constructor.

Definition at line 25 of file ProtoShowerMatchingTool.cc.

                                                 :
    m_spineSlidingFitWindow(20),
    m_maxXSeparation(5.f),
    m_maxSeparation(5.f),
    m_xExtrapolation(5.f),
    m_maxAngularDeviation(5.f)
{
}

Member Function Documentation

bool lar_content::ProtoShowerMatchingTool::AreDirectionsConsistent ( const ProtoShower &  protoShowerU, const ProtoShower &  protoShowerV, const ProtoShower &  protoShowerW  ) const

private

Determine whether three 2D initial spine directions correspond to the same 3D initial spine direction.

Parameters

protoShowerU	the U view ProtoShower
protoShowerV	the V view ProtoShower
protoShowerW	the W view ProtoShower

Returns

whether three 2D initial spine directions correspond to the same 3D initial spine direction

Definition at line 142 of file ProtoShowerMatchingTool.cc.

References f, lar_content::ProtoShower::GetConnectionPathway(), lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerDirection(), lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection(), lar_content::ProtoShower::GetShowerCore(), lar_content::ProtoShower::GetSpineHitList(), lar_content::ConnectionPathway::GetStartDirection(), lar_content::ShowerCore::GetStartPosition(), lar_content::ConnectionPathway::GetStartPosition(), lar_content::LArGeometryHelper::GetWirePitch(), m_maxAngularDeviation, m_spineSlidingFitWindow, m_xExtrapolation, and lar_content::LArGeometryHelper::MergeTwoPositions().

Referenced by ArePathwaysConsistent().

{
    const CartesianVector &directionU1(protoShowerU.GetConnectionPathway().GetStartDirection());
    const CartesianVector &directionV1(protoShowerV.GetConnectionPathway().GetStartDirection());
    const CartesianVector &directionW1(protoShowerW.GetConnectionPathway().GetStartDirection());

    if (this->AreDirectionsConsistent(protoShowerU.GetConnectionPathway().GetStartPosition(), protoShowerV.GetConnectionPathway().GetStartPosition(),
            protoShowerW.GetConnectionPathway().GetStartPosition(), directionU1, directionV1, directionW1))
    {
        return true;
    }
    else
    {
        const bool isDownstream(
            protoShowerW.GetShowerCore().GetStartPosition().GetZ() > protoShowerW.GetConnectionPathway().GetStartPosition().GetZ());

        CartesianPointVector spinePositionsU, spinePositionsV, spinePositionsW;

        for (const CaloHit *const pCaloHit : protoShowerU.GetSpineHitList())
            spinePositionsU.push_back(pCaloHit->GetPositionVector());

        for (const CaloHit *const pCaloHit : protoShowerV.GetSpineHitList())
            spinePositionsV.push_back(pCaloHit->GetPositionVector());

        for (const CaloHit *const pCaloHit : protoShowerW.GetSpineHitList())
            spinePositionsW.push_back(pCaloHit->GetPositionVector());

        const TwoDSlidingFitResult spineFitU(&spinePositionsU, m_spineSlidingFitWindow, LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_U));
        const TwoDSlidingFitResult spineFitV(&spinePositionsV, m_spineSlidingFitWindow, LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_V));
        const TwoDSlidingFitResult spineFitW(&spinePositionsW, m_spineSlidingFitWindow, LArGeometryHelper::GetWirePitch(this->GetPandora(), TPC_VIEW_W));

        const CartesianVector directionU2(isDownstream ? spineFitU.GetGlobalMinLayerDirection() : spineFitU.GetGlobalMaxLayerDirection() * (-1.f));
        const CartesianVector directionV2(isDownstream ? spineFitV.GetGlobalMinLayerDirection() : spineFitV.GetGlobalMaxLayerDirection() * (-1.f));
        const CartesianVector directionW2(isDownstream ? spineFitW.GetGlobalMinLayerDirection() : spineFitW.GetGlobalMaxLayerDirection() * (-1.f));

        return this->AreDirectionsConsistent(protoShowerU.GetConnectionPathway().GetStartPosition(),
            protoShowerV.GetConnectionPathway().GetStartPosition(), protoShowerW.GetConnectionPathway().GetStartPosition(), directionU2,
            directionV2, directionW2);
    }

    return false;
}

bool lar_content::ProtoShowerMatchingTool::AreDirectionsConsistent ( const pandora::CartesianVector &  nuVertexU, const pandora::CartesianVector &  nuVertexV, const pandora::CartesianVector &  nuVertexW, const pandora::CartesianVector &  directionU, const pandora::CartesianVector &  directionV, const pandora::CartesianVector &  directionW  ) const

private

Determine whether three 2D initial spine directions correspond to the same 3D initial spine direction.

Parameters

nuVertexU	the U view projection of the neutrino vertex
nuVertexV	the V view projection of the neutrino vertex
nuVertexW	the W view projection of the neutrino vertex
directionU	the U view initial spine direction
directionU	the V view initial spine direction
directionU	the W view initial spine direction

Returns

whether three 2D initial spine directions correspond to the same 3D initial spine direction

bool lar_content::ProtoShowerMatchingTool::ArePathwaysConsistent ( const ProtoShower &  protoShowerU, const ProtoShower &  protoShowerV, const ProtoShower &  protoShowerW, Consistency &  consistency  ) const

private

Determine whether three 2D connection pathways form a consistent 3D connection pathway.

Parameters

protoShowerU	the U view ProtoShower
protoShowerV	the V view ProtoShower
protoShowerW	the W view ProtoShower
consistency	the basis of the match

Returns

whether three 2D connection pathways form a consistent 3D connection pathway

Definition at line 90 of file ProtoShowerMatchingTool.cc.

References AreDirectionsConsistent(), AreShowerStartsConsistent(), lar_content::DIRECTION, and lar_content::POSITION.

Referenced by Run().

{
    if (this->AreShowerStartsConsistent(protoShowerU, protoShowerV, protoShowerW))
    {
        consistency = Consistency::POSITION;
    }
    else if (this->AreDirectionsConsistent(protoShowerU, protoShowerV, protoShowerW))
    {
        consistency = Consistency::DIRECTION;
    }
    else
    {
        return false;
    }

    return true;
}

bool lar_content::ProtoShowerMatchingTool::AreShowerStartsConsistent ( const ProtoShower &  protoShowerU, const ProtoShower &  protoShowerV, const ProtoShower &  protoShowerW  ) const

private

Determine whether three 2D shower start positions correspond to the same 3D shower start position.

Parameters

protoShowerU	the U view ProtoShower
protoShowerV	the V view ProtoShower
protoShowerW	the W view ProtoShower

Returns

whether three 2D shower start positions correspond to the same 3D shower start position

Definition at line 111 of file ProtoShowerMatchingTool.cc.

References f, lar_content::ProtoShower::GetShowerCore(), lar_content::ShowerCore::GetStartPosition(), m_maxSeparation, m_maxXSeparation, and lar_content::LArGeometryHelper::MergeTwoPositions().

Referenced by ArePathwaysConsistent().

{
    const CartesianVector &showerStartU(protoShowerU.GetShowerCore().GetStartPosition());
    const CartesianVector &showerStartV(protoShowerV.GetShowerCore().GetStartPosition());
    const CartesianVector &showerStartW(protoShowerW.GetShowerCore().GetStartPosition());

    const float xSeparationUV(std::fabs(showerStartU.GetX() - showerStartV.GetX()));
    const float xSeparationUW(std::fabs(showerStartU.GetX() - showerStartW.GetX()));
    const float xSeparationVW(std::fabs(showerStartV.GetX() - showerStartW.GetX()));

    if ((xSeparationUV > m_maxXSeparation) || (xSeparationUW > m_maxXSeparation) || (xSeparationVW > m_maxXSeparation))
        return false;

    float chi2(0.f);
    CartesianVector projectionU(0.f, 0.f, 0.f), projectionV(0.f, 0.f, 0.f), projectionW(0.f, 0.f, 0.f);

    LArGeometryHelper::MergeTwoPositions(this->GetPandora(), TPC_VIEW_V, TPC_VIEW_W, showerStartV, showerStartW, projectionU, chi2);
    LArGeometryHelper::MergeTwoPositions(this->GetPandora(), TPC_VIEW_W, TPC_VIEW_U, showerStartW, showerStartU, projectionV, chi2);
    LArGeometryHelper::MergeTwoPositions(this->GetPandora(), TPC_VIEW_U, TPC_VIEW_V, showerStartU, showerStartV, projectionW, chi2);

    const float separationU((projectionU - showerStartU).GetMagnitude());
    const float separationV((projectionV - showerStartV).GetMagnitude());
    const float separationW((projectionW - showerStartW).GetMagnitude());

    const float metric((separationU + separationV + separationW) / 3.f);

    return (metric < m_maxSeparation);
}

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

private

Definition at line 258 of file ProtoShowerMatchingTool.cc.

References m_maxAngularDeviation, m_maxSeparation, m_maxXSeparation, m_spineSlidingFitWindow, and m_xExtrapolation.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "SpineSlidingFitWindow", m_spineSlidingFitWindow));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxXSeparation", m_maxXSeparation));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxSeparation", m_maxSeparation));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "XExtrapolation", m_xExtrapolation));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxAngularDeviation", m_maxAngularDeviation));

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::ProtoShowerMatchingTool::Run	(	const ProtoShowerVector &	protoShowerVectorU,
		const ProtoShowerVector &	protoShowerVectorV,
		const ProtoShowerVector &	protoShowerVectorW,
		ProtoShowerMatchVector &	protoShowerMatchVector
	)

Definition at line 36 of file ProtoShowerMatchingTool.cc.

References ArePathwaysConsistent(), and lar_content::POSITION.

Referenced by lar_content::ElectronInitialRegionRefinementAlgorithm::RefineShower().

{
    IntVector usedProtoShowersU, usedProtoShowersV, usedProtoShowersW;

    bool added(false);

    for (unsigned int uIndex = 0; uIndex < protoShowerVectorU.size(); ++uIndex)
    {
        added = false;

        if (std::find(usedProtoShowersU.begin(), usedProtoShowersU.end(), uIndex) != usedProtoShowersU.end())
            continue;

        const ProtoShower &protoShowerU(protoShowerVectorU.at(uIndex));

        for (unsigned int vIndex = 0; vIndex < protoShowerVectorV.size(); ++vIndex)
        {
            if (std::find(usedProtoShowersV.begin(), usedProtoShowersV.end(), vIndex) != usedProtoShowersV.end())
                continue;

            const ProtoShower &protoShowerV(protoShowerVectorV.at(vIndex));

            for (unsigned int wIndex = 0; wIndex < protoShowerVectorW.size(); ++wIndex)
            {
                if (std::find(usedProtoShowersW.begin(), usedProtoShowersW.end(), wIndex) != usedProtoShowersW.end())
                    continue;

                const ProtoShower &protoShowerW(protoShowerVectorW.at(wIndex));
                Consistency consistency(Consistency::POSITION);

                if (!this->ArePathwaysConsistent(protoShowerU, protoShowerV, protoShowerW, consistency))
                    continue;

                usedProtoShowersU.push_back(uIndex);
                usedProtoShowersV.push_back(vIndex);
                usedProtoShowersW.push_back(wIndex);

                protoShowerMatchVector.push_back(ProtoShowerMatch(protoShowerU, protoShowerV, protoShowerW, consistency));

                added = true;
                break;
            }

            if (added)
                break;
        }
    }

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

float lar_content::ProtoShowerMatchingTool::m_maxAngularDeviation

private

The max. opening angle between true and projected 2D initial directions for a match.

Definition at line 88 of file ProtoShowerMatchingTool.h.

Referenced by AreDirectionsConsistent(), and ReadSettings().

float lar_content::ProtoShowerMatchingTool::m_maxSeparation

private

The max. average separation between true and projected 2D shower start positions for a match.

Definition at line 86 of file ProtoShowerMatchingTool.h.

Referenced by AreShowerStartsConsistent(), and ReadSettings().

float lar_content::ProtoShowerMatchingTool::m_maxXSeparation

private

The max. drift-coordinate separation between matched 2D shower start positions.

Definition at line 85 of file ProtoShowerMatchingTool.h.

Referenced by AreShowerStartsConsistent(), and ReadSettings().

unsigned int lar_content::ProtoShowerMatchingTool::m_spineSlidingFitWindow

private

The shower spine sliding fit window.

Definition at line 84 of file ProtoShowerMatchingTool.h.

Referenced by AreDirectionsConsistent(), and ReadSettings().

float lar_content::ProtoShowerMatchingTool::m_xExtrapolation

private

Extrapolation distance in the x-direction.

Definition at line 87 of file ProtoShowerMatchingTool.h.

Referenced by AreDirectionsConsistent(), and ReadSettings().

---

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

• larpandoracontent/v04_14_00/source/larpandoracontent/LArShowerRefinement/ProtoShowerMatchingTool.h
• larpandoracontent/v04_14_00/source/larpandoracontent/LArShowerRefinement/ProtoShowerMatchingTool.cc