lar_dl_content::DlHitValidationAlgorithm Class Reference

DlHitValidationlgorithm class. More...

#include <DlHitValidationAlgorithm.h>

Inheritance diagram for lar_dl_content::DlHitValidationAlgorithm:

Public Member Functions
	DlHitValidationAlgorithm ()
	Default constructor. More...
virtual	~DlHitValidationAlgorithm ()

Private Member Functions
pandora::StatusCode	Run ()
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
pandora::StringVector	m_caloHitListNames
	Name of input calo hit list. More...
int	m_confusionU [2][2]
	Confusion matrix for the U view. More...
int	m_confusionV [2][2]
	Confusion matrix for the V view. More...
int	m_confusionW [2][2]
	Confusion matrix for the W view. More...

Detailed Description

DlHitValidationlgorithm class.

Definition at line 19 of file DlHitValidationAlgorithm.h.

Constructor & Destructor Documentation

lar_dl_content::DlHitValidationAlgorithm::DlHitValidationAlgorithm

(

)

Default constructor.

Definition at line 24 of file DlHitValidationAlgorithm.cc.

                                                   : m_confusionU(), m_confusionV(), m_confusionW()
{
}

lar_dl_content::DlHitValidationAlgorithm::~DlHitValidationAlgorithm ( )

virtual

Definition at line 30 of file DlHitValidationAlgorithm.cc.

{
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "u_true_shower", m_confusionU[0][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "u_false_shower", m_confusionU[1][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "u_false_track", m_confusionU[0][1]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "u_true_track", m_confusionU[1][1]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "v_true_shower", m_confusionV[0][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "v_false_shower", m_confusionV[1][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "v_false_track", m_confusionV[0][1]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "v_true_track", m_confusionV[1][1]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "w_true_shower", m_confusionW[0][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "w_false_shower", m_confusionW[1][0]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "w_false_track", m_confusionW[0][1]));
    PANDORA_MONITORING_API(SetTreeVariable(this->GetPandora(), "confusion_tree", "w_true_track", m_confusionW[1][1]));
    PANDORA_MONITORING_API(FillTree(this->GetPandora(), "confusion_tree"));
    try
    {
        PANDORA_MONITORING_API(SaveTree(this->GetPandora(), "confusion_tree", "confusion.root", "UPDATE"));
    }
    catch (const StatusCodeException &)
    {
        std::cout << "DlHitValidationAlgorithm: Unable to write confusion_tree to file" << std::endl;
    }
}

Member Function Documentation

StatusCode lar_dl_content::DlHitValidationAlgorithm::ReadSettings ( const pandora::TiXmlHandle  xmlHandle )

private

Definition at line 111 of file DlHitValidationAlgorithm.cc.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadVectorOfValues(xmlHandle, "CaloHitListNames", m_caloHitListNames));

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_dl_content::DlHitValidationAlgorithm::Run ( )

private

Definition at line 57 of file DlHitValidationAlgorithm.cc.

{
    const int SHOWER_IDX{0}, TRACK_IDX{1};
    for (const std::string listName : m_caloHitListNames)
    {
        const CaloHitList *pCaloHitList(nullptr);
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetList(*this, listName, pCaloHitList));
        const MCParticleList *pMCParticleList(nullptr);
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pMCParticleList));

        const HitType view{pCaloHitList->front()->GetHitType()};

        if (!(view == TPC_VIEW_U || view == TPC_VIEW_V || view == TPC_VIEW_W))
            return STATUS_CODE_NOT_ALLOWED;

        LArMCParticleHelper::PrimaryParameters parameters;
        // Only care about reconstructability with respect to the current view, so skip good view check
        parameters.m_minHitsForGoodView = 0;
        // Turn off max photo propagation for now, only care about killing off daughters of neutrons
        parameters.m_maxPhotonPropagation = std::numeric_limits<float>::max();
        LArMCParticleHelper::MCContributionMap targetMCParticleToHitsMap;
        LArMCParticleHelper::SelectReconstructableMCParticles(
            pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsBeamNeutrinoFinalState, targetMCParticleToHitsMap);

        for (const CaloHit *pCaloHit : *pCaloHitList)
        {
            try
            {
                const MCParticle *const pMCParticle(MCParticleHelper::GetMainMCParticle(pCaloHit));
                const int pdg{std::abs(pMCParticle->GetParticleId())};
                const int truth{(pdg == 11 || pdg == 22) ? SHOWER_IDX : TRACK_IDX};
                const LArCaloHit *pLArCaloHit{dynamic_cast<const LArCaloHit *>(pCaloHit)};
                const float pTrack{pLArCaloHit->GetTrackProbability()};
                const float pShower{pLArCaloHit->GetShowerProbability()};
                const int cls{(pShower > pTrack) ? SHOWER_IDX : TRACK_IDX};
                if (view == TPC_VIEW_U)
                    ++m_confusionU[truth][cls];
                else if (view == TPC_VIEW_V)
                    ++m_confusionV[truth][cls];
                else
                    ++m_confusionW[truth][cls];
            }
            catch (const StatusCodeException &)
            {
                continue;
            }
        }
    }

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

pandora::StringVector lar_dl_content::DlHitValidationAlgorithm::m_caloHitListNames

private

Name of input calo hit list.

Definition at line 33 of file DlHitValidationAlgorithm.h.

int lar_dl_content::DlHitValidationAlgorithm::m_confusionU[2][2]

private

Confusion matrix for the U view.

Definition at line 34 of file DlHitValidationAlgorithm.h.

int lar_dl_content::DlHitValidationAlgorithm::m_confusionV[2][2]

private

Confusion matrix for the V view.

Definition at line 35 of file DlHitValidationAlgorithm.h.

int lar_dl_content::DlHitValidationAlgorithm::m_confusionW[2][2]

private

Confusion matrix for the W view.

Definition at line 36 of file DlHitValidationAlgorithm.h.

---

The documentation for this class was generated from the following files:

• DlHitValidationAlgorithm.h
• DlHitValidationAlgorithm.cc