lar_content::SimplePCAThreeDClusteringTool Class Reference

#include "SimplePCAThreeDClusteringTool.h"

Inheritance diagram for lar_content::SimplePCAThreeDClusteringTool:

Public Member Functions
	SimplePCAThreeDClusteringTool ()

Private Member Functions
bool	Run (const pandora::CaloHitList &inputCaloHitList, std::vector< pandora::CaloHitList > &outputCaloHitListsVector)
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle)

Detailed Description

Definition at line 18 of file SimplePCAThreeDClusteringTool.h.

Constructor & Destructor Documentation

lar_content::SimplePCAThreeDClusteringTool::SimplePCAThreeDClusteringTool

(

)

Definition at line 20 of file SimplePCAThreeDClusteringTool.cc.

{
}

Member Function Documentation

StatusCode lar_content::SimplePCAThreeDClusteringTool::ReadSettings ( const pandora::TiXmlHandle  )

private

Definition at line 111 of file SimplePCAThreeDClusteringTool.cc.

{
    return STATUS_CODE_SUCCESS;
}

bool lar_content::SimplePCAThreeDClusteringTool::Run ( const pandora::CaloHitList &  inputCaloHitList, std::vector< pandora::CaloHitList > &  outputCaloHitListsVector  )

privatevirtual

Implements lar_content::ClusteringTool.

Definition at line 26 of file SimplePCAThreeDClusteringTool.cc.

References f, and GetIntersection().

{
    // Begin with a PCA
    CartesianVector centroid(0.f, 0.f, 0.f);
    LArPcaHelper::EigenVectors eigenVecs;
    LArPcaHelper::EigenValues eigenValues(0.f, 0.f, 0.f);
    LArPcaHelper::RunPca(inputCaloHitList, centroid, eigenValues, eigenVecs);

    // By convention, the primary axis has a positive z-component - define ortho directions
    const CartesianVector orthoDirection1(eigenVecs.at(1).GetZ() > 0.f ? eigenVecs.at(1) : eigenVecs.at(1) * -1.f);

    //Lists for hits that have a positive and negative projection on the secondary axis
    CaloHitList posCaloHitList, negCaloHitList;

    for (const CaloHit *const pCaloHit3D : inputCaloHitList)
    {
        const CartesianVector pCaloHit3DPosition{pCaloHit3D->GetPositionVector()};

        if ((pCaloHit3DPosition - centroid).GetDotProduct(orthoDirection1) < 0)
        {
            negCaloHitList.push_back(pCaloHit3D);
        }
        else
        {
            posCaloHitList.push_back(pCaloHit3D);
        }
    }

    outputCaloHitListsVector.push_back(posCaloHitList);
    outputCaloHitListsVector.push_back(negCaloHitList);

    //Now, run PCA independently on pos and neg hit lists, to refine split

    //Get pos list PCA
    CartesianVector centroidPos(0.f, 0.f, 0.f);
    LArPcaHelper::EigenVectors eigenVecsPos;
    LArPcaHelper::EigenValues eigenValuesPos(0.f, 0.f, 0.f);
    LArPcaHelper::RunPca(posCaloHitList, centroidPos, eigenValuesPos, eigenVecsPos);

    // By convention, the primary axis has a positive z-component.
    const CartesianVector axisDirectionPos(eigenVecsPos.at(0).GetZ() > 0.f ? eigenVecsPos.at(0) : eigenVecsPos.at(0) * -1.f);

    //Get neg list PCA
    CartesianVector centroidNeg(0.f, 0.f, 0.f);
    LArPcaHelper::EigenVectors eigenVecsNeg;
    LArPcaHelper::EigenValues eigenValuesNeg(0.f, 0.f, 0.f);
    LArPcaHelper::RunPca(negCaloHitList, centroidNeg, eigenValuesNeg, eigenVecsNeg);

    // By convention, the primary axis has a positive z-component.
    const CartesianVector axisDirectionNeg(eigenVecsNeg.at(0).GetZ() > 0.f ? eigenVecsNeg.at(0) : eigenVecsNeg.at(0) * -1.f);

    //Get intersection point of the two new principal axes
    CartesianVector intersectionPoint(0.f, 0.f, 0.f);
    float displacementPos(0.f), displacementNeg(0.f);

    try
    {
        LArPointingClusterHelper::GetIntersection(
            centroidPos, axisDirectionPos, centroidNeg, axisDirectionNeg, intersectionPoint, displacementPos, displacementNeg);
    }
    catch (const StatusCodeException &)
    {
        std::cout << "Exception caught! Cannot get intersection between positive and negative hit lists primary PCA axes!" << std::endl;
    }

    //Clear pos and neg lists
    posCaloHitList.clear();
    negCaloHitList.clear();

    //Loop over original hit list, check whether it is within smaller cone of pos or neg axis, attach to relevant list
    for (const CaloHit *const pCaloHit3D : inputCaloHitList)
    {
        const float cosConeAxisPos = axisDirectionPos.GetCosOpeningAngle(pCaloHit3D->GetPositionVector() - intersectionPoint);
        const float cosConeAxisNeg = axisDirectionNeg.GetCosOpeningAngle(pCaloHit3D->GetPositionVector() - intersectionPoint);
        if (cosConeAxisPos > cosConeAxisNeg)
            posCaloHitList.push_back(pCaloHit3D);
        else
            negCaloHitList.push_back(pCaloHit3D);
    }

    return true;
}

---

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

• larpandoracontent/v04_14_00/source/larpandoracontent/LArReclustering/SimplePCAThreeDClusteringTool.h
• larpandoracontent/v04_14_00/source/larpandoracontent/LArReclustering/SimplePCAThreeDClusteringTool.cc