LArMuonLeadingHelper.h

Go to the documentation of this file.

#ifndef LAR_MUON_LEADING_HELPER_H
#define LAR_MUON_LEADING_HELPER_H 1

#include "Pandora/PandoraInternal.h"

#include "larpandoracontent/LArHelpers/LArMCParticleHelper.h"

namespace lar_content
{
class LArMuonLeadingHelper
{
public:
    class ValidationParameters : public LArMCParticleHelper::PrimaryParameters
    {
    public:
        ValidationParameters();

        float m_maxBremsstrahlungSeparation; 
    };

    typedef std::map<const pandora::MCParticle *, pandora::CaloHitList> LeadingMCParticleToPostBremsstrahlungHitList;

    static MCProcess GetLeadingProcess(const pandora::MCParticle *const pMCParticle);

    static bool IsDeltaRay(const pandora::MCParticle *const pMCParticle);

    static bool IsMichel(const pandora::MCParticle *const pMCParticle);

    static bool IsMuonLeading(const pandora::MCParticle *const pMCParticle);

    static const pandora::MCParticle *GetLeadingParticle(const pandora::MCParticle *const pMCParticle);

    static void SelectReconstructableLeadingParticles(const pandora::MCParticleList *pMCParticleList,
        const pandora::CaloHitList *pCaloHitList, const ValidationParameters &parameters, const pandora::CaloHitList &recoMuonHitList,
        LArMCParticleHelper::MCContributionMap &selectedMCParticlesToHitsMap);

    static void GetPfoMatchContamination(const pandora::MCParticle *const pLeadingParticle, const pandora::CaloHitList &matchedPfoHitList,
        pandora::CaloHitList &parentTrackHits, pandora::CaloHitList &otherTrackHits, pandora::CaloHitList &otherShowerHits);

    static void GetMuonPfoContaminationContribution(const pandora::CaloHitList &cosmicRayPfoHitList,
        const pandora::CaloHitList &leadingMCHitList, pandora::CaloHitList &leadingHitsInParentCosmicRay);

    static float GetClosestDistance(const pandora::Cluster *const pCluster, const pandora::CartesianPointVector &cartesianPointVector);

    static float GetClosestDistance(const pandora::CaloHit *const pCaloHit, const pandora::CartesianPointVector &cartesianPointVector);

    static pandora::StatusCode GetClosestPosition(const pandora::CartesianVector &referencePoint,
        const pandora::CartesianPointVector &cartesianPointVector, const pandora::Cluster *const pCluster, const float maxDistanceToCluster,
        const float maxDistanceToReferencePoint, pandora::CartesianVector &closestPosition);

    static void GetClosestPositions(const pandora::CartesianPointVector &cartesianPointVector1, const pandora::Cluster *const pCluster2,
        pandora::CartesianVector &outputPosition1, pandora::CartesianVector &outputPosition2);

private:
    static void GetMCToLeadingMap(const pandora::MCParticleList *const pMCParticleList, LArMCParticleHelper::MCRelationMap &mcToLeadingMap);

    static void SelectCaloHits(const pandora::CaloHitList *const pCaloHitList, const LArMCParticleHelper::MCRelationMap &mcToTargetMCMap,
        pandora::CaloHitList &selectedCaloHitList, const bool selectInputHits, const float minHitSharingFraction,
        const pandora::CaloHitList &recoMuonHitList, LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList);

    static bool RejectBremsstrahlungHits(
        const pandora::CaloHit *const pCaloHit, LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList);

    static void AddInPostBremsstrahlungHits(const LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList,
        const float maxBremsstrahlungSeparation, LArMCParticleHelper::MCContributionMap &leadingMCToTrueHitListMap);

    static void AddInPostBremsstrahlungHits(const pandora::MCParticle *const pLeadingParticle,
        const LeadingMCParticleToPostBremsstrahlungHitList &leadingMCParticleToPostBremsstrahlungHitList, const float maxBremsstrahlungSeparation,
        LArMCParticleHelper::MCContributionMap &leadingMCToTrueHitListMap, const pandora::HitType tpcView);

    static void SelectLeadingMCParticles(const pandora::MCParticleList *pMCParticleList, pandora::MCParticleVector &selectedParticles);
};

} // namespace lar_content

#endif // #ifndef LAR_MUON_LEADING_HELPER_H