BeamParticleIdTool.h

Go to the documentation of this file.

#ifndef LAR_BEAM_PARTICLE_ID_TOOL_H
#define LAR_BEAM_PARTICLE_ID_TOOL_H 1

#include "larpandoracontent/LArControlFlow/MasterAlgorithm.h"
#include "larpandoracontent/LArControlFlow/SliceIdBaseTool.h"

namespace lar_content
{

class BeamParticleIdTool : public SliceIdBaseTool
{
public:
    BeamParticleIdTool();

    void SelectOutputPfos(const pandora::Algorithm *const pAlgorithm, const SliceHypotheses &beamSliceHypotheses,
        const SliceHypotheses &crSliceHypotheses, pandora::PfoList &selectedPfos);

private:
    class Plane
    {
    public:
        Plane(const pandora::CartesianVector &normal, const pandora::CartesianVector &point);

        pandora::CartesianVector GetLineIntersection(const pandora::CartesianVector &a0, const pandora::CartesianVector &a) const;

    private:
        pandora::CartesianVector m_unitNormal; 
        pandora::CartesianVector m_point;      
        float m_d;                             
    };

    pandora::StatusCode Initialize();

    void GetSelectedCaloHits(const pandora::CaloHitList &inputCaloHitList, pandora::CaloHitList &outputCaloHitList, float &closestHitToFaceDistance) const;

    void GetTPCIntercepts(const pandora::CartesianVector &a0, const pandora::CartesianVector &majorAxis,
        pandora::CartesianVector &interceptOne, pandora::CartesianVector &interceptTwo) const;

    bool IsContained(const pandora::CartesianVector &spacePoint) const;

    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    typedef std::vector<Plane> PlaneVector;

    bool m_selectAllBeamParticles;            
    bool m_selectOnlyFirstSliceBeamParticles; 
    float m_tpcMinX;                          
    float m_tpcMaxX;                          
    float m_tpcMinY;                          
    float m_tpcMaxY;                          
    float m_tpcMinZ;                          
    float m_tpcMaxZ;                          
    pandora::CartesianVector m_beamTPCIntersection; 
    pandora::CartesianVector m_beamDirection;       
    PlaneVector m_tpcPlanes;                        

    float m_projectionIntersectionCut; 
    float m_closestDistanceCut;        
    float m_angleToBeamCut;            
    float m_selectedFraction;          
    unsigned int m_nSelectedHits;      
};

} // namespace lar_content

#endif // #ifndef LAR_BEAM_PARTICLE_ID_TOOL_H