MuonLeadingEventValidationAlgorithm.h

Go to the documentation of this file.

 
#ifndef LAR_MUON_LEADING_EVENT_VALIDATION_ALGORITHM_H
#define LAR_MUON_LEADING_EVENT_VALIDATION_ALGORITHM_H 1
 
#include "Pandora/Algorithm.h"
 
#include "larpandoracontent/LArHelpers/LArMuonLeadingHelper.h"
 
#include "larpandoracontent/LArMonitoring/EventValidationBaseAlgorithm.h"
 
#ifdef MONITORING
#include "PandoraMonitoringApi.h"
#endif
 
#include <map>
 
namespace lar_content
{
class MuonLeadingEventValidationAlgorithm : public EventValidationBaseAlgorithm
{
public:
    MuonLeadingEventValidationAlgorithm();
 
    virtual ~MuonLeadingEventValidationAlgorithm();
 
private:
    void FillValidationInfo(const pandora::MCParticleList *const pMCParticleList, const pandora::CaloHitList *const pCaloHitList,
        const pandora::PfoList *const pPfoList, ValidationInfo &validationInfo) const;
 
    void GetRecoCosmicRayHits(const pandora::MCParticleList *const pMCParticleList, const pandora::CaloHitList *const pCaloHitList,
        const pandora::PfoList *const pPfoList, pandora::CaloHitList &recoCosmicRayHitList) const;
 
    void PerformUnfoldedMatching(const pandora::MCParticleList *const pMCParticleList, const pandora::CaloHitList *const pCaloHitList,
        const pandora::PfoList *const pPfoList, const pandora::CaloHitList &recoCosmicRayHitList, const float minHitSharingFraction,
        ValidationInfo &validationInfo) const;
 
    void RemoveIncorrectlyReconstructedCosmicRays(const pandora::MCParticleList *const pMCParticleList,
        const pandora::CaloHitList *const pCaloHitList, const pandora::PfoList *const pPfoList, ValidationInfo &validationInfo) const;
 
    void DetermineIncorrectlyReconstructedCosmicRays(const pandora::MCParticleList *const pMCParticleList, const pandora::CaloHitList *const pCaloHitList,
        const pandora::PfoList *const pPfoList, pandora::MCParticleList &incorrectlyReconstructedCosmicRays) const;
 
    void ProcessOutput(const ValidationInfo &validationInfo, const bool useInterpretedMatching, const bool printToScreen, const bool fillTree) const;
 
#ifdef MONITORING
    void PrintHits(const pandora::CaloHitList caloHitList, const bool isCR) const;
#endif
 
#ifdef MONITORING
    void PrintHits(const pandora::CaloHitList totalCaloHitList, const pandora::CaloHitList otherShowerCaloHitList,
        const pandora::CaloHitList otherTrackCaloHitList, const pandora::CaloHitList parentTrackCaloHitList, const std::string &stringTag) const;
#endif
 
    void PrintHits(const pandora::CaloHitList totalCaloHitList, const pandora::CaloHitList leadingCaloHitList, const std::string &stringTag) const;
 
    void FillContaminationHitsDistance(const pandora::CaloHitList &contaminationHits, const pandora::CaloHitList &leadingMCHits,
        pandora::FloatVector &bestMatchContaminationHitsDistance) const;
 
    void GetHitsOfType(const pandora::CaloHitList &inputList, const pandora::HitType hitType, pandora::CaloHitList &outputList) const;
 
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);
 
    LArMuonLeadingHelper::ValidationParameters m_validationParameters; 
    bool m_removeRecoCosmicRayHits;   
    bool m_deltaRayMode;              
    bool m_michelMode;                
    int m_cosmicRaysToSkip;           
    bool m_visualize;                 
    bool m_ignoreIncorrectCosmicRays; 
    bool m_writeRawMatchesToTree;     
    std::vector<int> m_deltaRayIDs;   
};
 
} // namespace lar_content
 
#endif // LAR_MUON_LEADING_EVENT_VALIDATION_ALGORITHM_H