lar_content::TransverseAssociationAlgorithm::LArTransverseCluster Class Reference

LArTransverseCluster class. More...

Public Member Functions
	LArTransverseCluster (const pandora::Cluster *const pSeedCluster, const pandora::ClusterVector &associatedClusters)
	Constructor. More...
const pandora::Cluster *	GetSeedCluster () const
	Constructor. More...
const pandora::ClusterVector &	GetAssociatedClusters () const
	Get the associated cluster vector. More...
const pandora::CartesianVector &	GetInnerVertex () const
	Get the inner vertex position. More...
const pandora::CartesianVector &	GetOuterVertex () const
	Get the outer vertex position. More...
const pandora::CartesianVector &	GetDirection () const
	Get the direction. More...

Private Attributes
const pandora::Cluster *	m_pSeedCluster
pandora::ClusterVector	m_associatedClusters
pandora::CartesianVector	m_innerVertex
pandora::CartesianVector	m_outerVertex
pandora::CartesianVector	m_direction

Detailed Description

LArTransverseCluster class.

Definition at line 40 of file TransverseAssociationAlgorithm.h.

Constructor & Destructor Documentation

lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::LArTransverseCluster	(	const pandora::Cluster *const	pSeedCluster,
		const pandora::ClusterVector &	associatedClusters
	)

Constructor.

Parameters

pSeedCluster
associatedClusters

Definition at line 780 of file TransverseAssociationAlgorithm.cc.

                                                                                                                                                   :
    m_pSeedCluster(pSeedCluster),
    m_associatedClusters(associatedClusters),
    m_innerVertex(0.f, 0.f, 0.f),
    m_outerVertex(0.f, 0.f, 0.f),
    m_direction(0.f, 0.f, 0.f)
{
    double Swzz(0.), Swxx(0.), Swzx(0.), Swz(0.), Swx(0.), Sw(0.);
    double minX(std::numeric_limits<double>::max());
    double maxX(-std::numeric_limits<double>::max());

    ClusterList clusterList(1, pSeedCluster);
    clusterList.insert(clusterList.end(), associatedClusters.begin(), associatedClusters.end());

    for (ClusterList::const_iterator iterI = clusterList.begin(), iterEndI = clusterList.end(); iterI != iterEndI; ++iterI)
    {
        for (OrderedCaloHitList::const_iterator iterJ = (*iterI)->GetOrderedCaloHitList().begin(),
                                                iterEndJ = (*iterI)->GetOrderedCaloHitList().end();
             iterJ != iterEndJ; ++iterJ)
        {
            for (CaloHitList::const_iterator iterK = iterJ->second->begin(), iterEndK = iterJ->second->end(); iterK != iterEndK; ++iterK)
            {
                const CaloHit *const pCaloHit = *iterK;

                if (pCaloHit->GetPositionVector().GetX() < minX)
                    minX = pCaloHit->GetPositionVector().GetX();

                if (pCaloHit->GetPositionVector().GetX() > maxX)
                    maxX = pCaloHit->GetPositionVector().GetX();

                Swzz += pCaloHit->GetPositionVector().GetZ() * pCaloHit->GetPositionVector().GetZ();
                Swxx += pCaloHit->GetPositionVector().GetX() * pCaloHit->GetPositionVector().GetX();
                Swzx += pCaloHit->GetPositionVector().GetZ() * pCaloHit->GetPositionVector().GetX();
                Swz += pCaloHit->GetPositionVector().GetZ();
                Swx += pCaloHit->GetPositionVector().GetX();
                Sw += 1.;
            }
        }
    }

    if (Sw > 0.f)
    {
        const double averageX(Swx / Sw);
        const double averageZ(Swz / Sw);

        if (Sw * Swxx - Swx * Swx > 0.)
        {
            double m((Sw * Swzx - Swx * Swz) / (Sw * Swxx - Swx * Swx));
            double px(1. / std::sqrt(1. + m * m));
            double pz(m / std::sqrt(1. + m * m));

            m_innerVertex.SetValues(static_cast<float>(minX), 0.f, static_cast<float>(averageZ + m * (minX - averageX)));
            m_outerVertex.SetValues(static_cast<float>(maxX), 0.f, static_cast<float>(averageZ + m * (maxX - averageX)));
            m_direction.SetValues(static_cast<float>(px), 0.f, static_cast<float>(pz));
        }
        else
        {
            m_innerVertex.SetValues(static_cast<float>(averageX), 0.f, static_cast<float>(averageZ));
            m_outerVertex.SetValues(static_cast<float>(averageX), 0.f, static_cast<float>(averageZ));
            m_direction.SetValues(1.f, 0.f, 0.f);
        }
    }
    else
    {
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);
    }
}

Member Function Documentation

const pandora::ClusterVector & lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::GetAssociatedClusters ( ) const

inline

Get the associated cluster vector.

Returns

the associated cluster vector

Definition at line 365 of file TransverseAssociationAlgorithm.h.

{
    return m_associatedClusters;
}

const pandora::CartesianVector & lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::GetDirection ( ) const

inline

Get the direction.

Returns

the direction

Definition at line 386 of file TransverseAssociationAlgorithm.h.

{
    return m_direction;
}

const pandora::CartesianVector & lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::GetInnerVertex ( ) const

inline

Get the inner vertex position.

Returns

the inner vertex position

Definition at line 372 of file TransverseAssociationAlgorithm.h.

{
    return m_innerVertex;
}

const pandora::CartesianVector & lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::GetOuterVertex ( ) const

inline

Get the outer vertex position.

Returns

the outer vertex position

Definition at line 379 of file TransverseAssociationAlgorithm.h.

{
    return m_outerVertex;
}

const pandora::Cluster * lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::GetSeedCluster ( ) const

inline

Constructor.

Returns

the address of the seed cluster

Definition at line 358 of file TransverseAssociationAlgorithm.h.

{
    return m_pSeedCluster;
}

Member Data Documentation

pandora::ClusterVector lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::m_associatedClusters

private

Definition at line 88 of file TransverseAssociationAlgorithm.h.

pandora::CartesianVector lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::m_direction

private

Definition at line 91 of file TransverseAssociationAlgorithm.h.

pandora::CartesianVector lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::m_innerVertex

private

Definition at line 89 of file TransverseAssociationAlgorithm.h.

pandora::CartesianVector lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::m_outerVertex

private

Definition at line 90 of file TransverseAssociationAlgorithm.h.

const pandora::Cluster* lar_content::TransverseAssociationAlgorithm::LArTransverseCluster::m_pSeedCluster

private

Definition at line 87 of file TransverseAssociationAlgorithm.h.

---

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

• TransverseAssociationAlgorithm.h
• TransverseAssociationAlgorithm.cc