lar_content::LArParticleIdPlugins::LArMuonId Class Reference

LArMuonId class. More...

#include "LArParticleIdPlugins.h"

Inheritance diagram for lar_content::LArParticleIdPlugins::LArMuonId:

Public Member Functions
	LArMuonId ()
	Default constructor. More...
bool	IsMatch (const pandora::Cluster *const pCluster) const
bool	IsMatch (const pandora::ParticleFlowObject *const pPfo) const

Private Member Functions
float	GetMuonTrackWidth (const TwoDSlidingFitResult &twoDSlidingFitResult) const
	Get the muon track width estimator for a provided sliding fit result. More...
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
unsigned int	m_layerFitHalfWindow
	Layer fit half window, used for calculating sliding muon track width. More...
float	m_minLayerOccupancy
	Min layer occupancy for for muon identification. More...
float	m_maxTrackWidth
	Max muon track width estimator for muon identification. More...
float	m_trackResidualQuantile
	Track residual quantile, used for calculating muon track width. More...
unsigned int	m_minClustersPassingId
	Match pfo if at sufficient clusters in pfo pass the cluster particle id logic. More...

Detailed Description

LArMuonId class.

Definition at line 29 of file LArParticleIdPlugins.h.

Constructor & Destructor Documentation

lar_content::LArParticleIdPlugins::LArMuonId::LArMuonId

(

)

Default constructor.

Definition at line 29 of file LArParticleIdPlugins.cc.

References lar_content::LArClusterHelper::GetLayerOccupancy(), GetMuonTrackWidth(), lar_content::LArGeometryHelper::GetWireZPitch(), IsMatch(), m_layerFitHalfWindow, m_maxTrackWidth, and m_minLayerOccupancy.

                                         :
    m_layerFitHalfWindow(20),
    m_minLayerOccupancy(0.75f),
    m_maxTrackWidth(0.5f),
    m_trackResidualQuantile(0.8f),
    m_minClustersPassingId(2)
{
}

Member Function Documentation

float lar_content::LArParticleIdPlugins::LArMuonId::GetMuonTrackWidth ( const TwoDSlidingFitResult &  twoDSlidingFitResult ) const

private

Get the muon track width estimator for a provided sliding fit result.

Parameters

twoDSlidingFitResult

the sliding fit result

Returns

the muon track width estimator

Definition at line 80 of file LArParticleIdPlugins.cc.

References f, lar_content::TwoDSlidingFitResult::GetCluster(), lar_content::TwoDSlidingFitResult::GetLayer(), lar_content::TwoDSlidingFitResult::GetLayerFitResultMap(), lar_content::TwoDSlidingFitResult::GetLocalPosition(), and m_trackResidualQuantile.

Referenced by LArMuonId().

{
    FloatVector residuals;
    const OrderedCaloHitList &orderedCaloHitList(twoDSlidingFitResult.GetCluster()->GetOrderedCaloHitList());
    const LayerFitResultMap &layerFitResultMap(twoDSlidingFitResult.GetLayerFitResultMap());

    for (OrderedCaloHitList::const_iterator iter = orderedCaloHitList.begin(); iter != orderedCaloHitList.end(); ++iter)
    {
        for (CaloHitList::const_iterator hitIter = iter->second->begin(), hitIterEnd = iter->second->end(); hitIter != hitIterEnd; ++hitIter)
        {
            float rL(0.f), rT(0.f);
            twoDSlidingFitResult.GetLocalPosition((*hitIter)->GetPositionVector(), rL, rT);
            const int layer(twoDSlidingFitResult.GetLayer(rL));

            LayerFitResultMap::const_iterator fitResultIter = layerFitResultMap.find(layer);

            if (layerFitResultMap.end() == fitResultIter)
                continue;

            const double fitT(fitResultIter->second.GetFitT());
            const double gradient(fitResultIter->second.GetGradient());
            const double residualSquared((fitT - rT) * (fitT - rT) / (1. + gradient * gradient)); // angular correction (note: this is cheating!)
            residuals.push_back(residualSquared);
        }
    }

    if (residuals.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    std::sort(residuals.begin(), residuals.end());
    const float theQuantile(residuals[m_trackResidualQuantile * residuals.size()]);

    return std::sqrt(theQuantile);
}

bool lar_content::LArParticleIdPlugins::LArMuonId::IsMatch

(

const pandora::Cluster *const

pCluster

)

const

Referenced by IsMatch(), and LArMuonId().

bool lar_content::LArParticleIdPlugins::LArMuonId::IsMatch

(

const pandora::ParticleFlowObject *const

pPfo

)

const

Definition at line 56 of file LArParticleIdPlugins.cc.

References lar_content::LArPfoHelper::GetTwoDClusterList(), IsMatch(), and m_minClustersPassingId.

{
    ClusterList clusterList;
    LArPfoHelper::GetTwoDClusterList(pPfo, clusterList);

    if (clusterList.empty())
        return false;

    unsigned int nClustersPassing(0);

    for (const Cluster *const pCluster : clusterList)
    {
        if (this->IsMatch(pCluster))
            ++nClustersPassing;
    }

    if (nClustersPassing < m_minClustersPassingId)
        return false;

    return true;
}

StatusCode lar_content::LArParticleIdPlugins::LArMuonId::ReadSettings ( const pandora::TiXmlHandle  xmlHandle )

private

Definition at line 117 of file LArParticleIdPlugins.cc.

References m_layerFitHalfWindow, m_maxTrackWidth, m_minClustersPassingId, m_minLayerOccupancy, and m_trackResidualQuantile.

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

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinLayerOccupancy", m_minLayerOccupancy));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxTrackWidth", m_maxTrackWidth));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "TrackResidualQuantile", m_trackResidualQuantile));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinClustersPassingId", m_minClustersPassingId));

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

unsigned int lar_content::LArParticleIdPlugins::LArMuonId::m_layerFitHalfWindow

private

Layer fit half window, used for calculating sliding muon track width.

Definition at line 52 of file LArParticleIdPlugins.h.

Referenced by LArMuonId(), and ReadSettings().

float lar_content::LArParticleIdPlugins::LArMuonId::m_maxTrackWidth

private

Max muon track width estimator for muon identification.

Definition at line 54 of file LArParticleIdPlugins.h.

Referenced by LArMuonId(), and ReadSettings().

unsigned int lar_content::LArParticleIdPlugins::LArMuonId::m_minClustersPassingId

private

Match pfo if at sufficient clusters in pfo pass the cluster particle id logic.

Definition at line 56 of file LArParticleIdPlugins.h.

Referenced by IsMatch(), and ReadSettings().

float lar_content::LArParticleIdPlugins::LArMuonId::m_minLayerOccupancy

private

Min layer occupancy for for muon identification.

Definition at line 53 of file LArParticleIdPlugins.h.

Referenced by LArMuonId(), and ReadSettings().

float lar_content::LArParticleIdPlugins::LArMuonId::m_trackResidualQuantile

private

Track residual quantile, used for calculating muon track width.

Definition at line 55 of file LArParticleIdPlugins.h.

Referenced by GetMuonTrackWidth(), and ReadSettings().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArPlugins/LArParticleIdPlugins.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArPlugins/LArParticleIdPlugins.cc