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"

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