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