lar_content::CosmicRayTaggingTool::CRCandidate Class Reference

Class to encapsulate the logic required determine if a Pfo should or shouldn't be tagged as a cosmic ray. More...

Public Member Functions
	CRCandidate (const pandora::Pandora &pandora, const pandora::ParticleFlowObject *const pPfo, const unsigned int sliceId)
	Constructor. More...

Public Attributes
const pandora::ParticleFlowObject *const	m_pPfo
	Address of the candidate Pfo. More...
unsigned int	m_sliceId
	Slice ID. More...
bool	m_canFit
	If there are a sufficient number of 3D hits to perform a fitting. More...
pandora::CartesianVector	m_endPoint1
	First fitted end point in 3D. More...
pandora::CartesianVector	m_endPoint2
	Second fitted end point in 3D. More...
double	m_length
	Straight line length of the linear fit. More...
double	m_curvature
	Measure of the curvature of the track. More...
double	m_theta
	Direction made with vertical. More...

Private Member Functions
void	CalculateFitVariables (const ThreeDSlidingFitResult &slidingFitResult)
	Calculate all variables which require a fit. More...

Detailed Description

Class to encapsulate the logic required determine if a Pfo should or shouldn't be tagged as a cosmic ray.

Definition at line 38 of file CosmicRayTaggingTool.h.

Constructor & Destructor Documentation

lar_content::CosmicRayTaggingTool::CRCandidate::CRCandidate	(	const pandora::Pandora &	pandora,
		const pandora::ParticleFlowObject *const	pPfo,
		const unsigned int	sliceId
	)

Constructor.

Parameters

pandora	the relevant pandora instance
pPfo	the address of the candidate pfo
slice	the slice id

Definition at line 521 of file CosmicRayTaggingTool.cc.

                                                                                                                                     :
    m_pPfo(pPfo),
    m_sliceId(sliceId),
    m_canFit(false),
    m_endPoint1(std::numeric_limits<float>::max(), std::numeric_limits<float>::max(), std::numeric_limits<float>::max()),
    m_endPoint2(std::numeric_limits<float>::max(), std::numeric_limits<float>::max(), std::numeric_limits<float>::max()),
    m_length(std::numeric_limits<float>::max()),
    m_curvature(std::numeric_limits<float>::max()),
    m_theta(std::numeric_limits<float>::max())
{
    ClusterList clusters3D;
    LArPfoHelper::GetThreeDClusterList(pPfo, clusters3D);

    if (!clusters3D.empty() && (clusters3D.front()->GetNCaloHits() > 15)) // TODO Configurable
    {
        m_canFit = true;
        const LArTPC *const pFirstLArTPC(pandora.GetGeometry()->GetLArTPCMap().begin()->second);
        const ThreeDSlidingFitResult slidingFitResult(clusters3D.front(), 5, pFirstLArTPC->GetWirePitchW()); // TODO Configurable
        this->CalculateFitVariables(slidingFitResult);
    }
}

Member Function Documentation

void lar_content::CosmicRayTaggingTool::CRCandidate::CalculateFitVariables ( const ThreeDSlidingFitResult &  slidingFitResult )

private

Calculate all variables which require a fit.

Parameters

slidingFitResult

the three dimensional sliding fit result

Definition at line 545 of file CosmicRayTaggingTool.cc.

{
    m_endPoint1 = slidingFitResult.GetGlobalMinLayerPosition();
    m_endPoint2 = slidingFitResult.GetGlobalMaxLayerPosition();
    m_length = (m_endPoint2 - m_endPoint1).GetMagnitude();

    if (std::fabs(m_length) > std::numeric_limits<float>::epsilon())
        m_theta = std::fabs(m_endPoint2.GetY() - m_endPoint1.GetY()) / m_length;

    const float layerPitch(slidingFitResult.GetFirstFitResult().GetLayerPitch());

    CartesianPointVector directionList;
    for (int i = slidingFitResult.GetMinLayer(); i < slidingFitResult.GetMaxLayer(); ++i)
    {
        CartesianVector direction(0.f, 0.f, 0.f);
        if (STATUS_CODE_SUCCESS == slidingFitResult.GetGlobalFitDirection(static_cast<float>(i) * layerPitch, direction))
            directionList.push_back(direction);
    }

    CartesianVector meanDirection(0.f, 0.f, 0.f);
    for (const CartesianVector &direction : directionList)
        meanDirection += direction;

    if (!directionList.empty() > 0)
        meanDirection *= 1.f / static_cast<float>(directionList.size());

    m_curvature = 0.f;
    for (const CartesianVector &direction : directionList)
        m_curvature += (direction - meanDirection).GetMagnitude();

    if (!directionList.empty() > 0)
        m_curvature *= 1.f / static_cast<float>(directionList.size());
}

Member Data Documentation

bool lar_content::CosmicRayTaggingTool::CRCandidate::m_canFit

If there are a sufficient number of 3D hits to perform a fitting.

Definition at line 52 of file CosmicRayTaggingTool.h.

double lar_content::CosmicRayTaggingTool::CRCandidate::m_curvature

Measure of the curvature of the track.

Definition at line 56 of file CosmicRayTaggingTool.h.

pandora::CartesianVector lar_content::CosmicRayTaggingTool::CRCandidate::m_endPoint1

First fitted end point in 3D.

Definition at line 53 of file CosmicRayTaggingTool.h.

pandora::CartesianVector lar_content::CosmicRayTaggingTool::CRCandidate::m_endPoint2

Second fitted end point in 3D.

Definition at line 54 of file CosmicRayTaggingTool.h.

double lar_content::CosmicRayTaggingTool::CRCandidate::m_length

Straight line length of the linear fit.

Definition at line 55 of file CosmicRayTaggingTool.h.

const pandora::ParticleFlowObject* const lar_content::CosmicRayTaggingTool::CRCandidate::m_pPfo

Address of the candidate Pfo.

Definition at line 50 of file CosmicRayTaggingTool.h.

unsigned int lar_content::CosmicRayTaggingTool::CRCandidate::m_sliceId

Slice ID.

Definition at line 51 of file CosmicRayTaggingTool.h.

double lar_content::CosmicRayTaggingTool::CRCandidate::m_theta

Direction made with vertical.

Definition at line 57 of file CosmicRayTaggingTool.h.

---

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

• larpandoracontent/larpandoracontent/LArControlFlow/CosmicRayTaggingTool.h
• larpandoracontent/larpandoracontent/LArControlFlow/CosmicRayTaggingTool.cc