lar_content::CheatingCosmicRayIdentificationAlg Class Reference

CheatingCosmicRayIdentificationAlg class. More...

#include <CheatingCosmicRayIdentificationAlg.h>

Inheritance diagram for lar_content::CheatingCosmicRayIdentificationAlg:

Public Member Functions
	CheatingCosmicRayIdentificationAlg ()
	Default constructor. More...

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

Private Attributes
std::string	m_inputPfoListName
	The input pfo list name. More...
std::string	m_outputPfoListName
	The output pfo list name. More...
std::string	m_inputDaughterPfoListName
	The input daughter pfo list name (if not specified, will assume same as main input list) More...
std::string	m_outputDaughterPfoListName
	The output daughter pfo list name (if not specified, will assume same as main output list) More...
float	m_maxNeutrinoFraction
	The maximum true neutrino fraction in a particle to be labelled as a cosmic ray. More...

Detailed Description

CheatingCosmicRayIdentificationAlg class.

Definition at line 19 of file CheatingCosmicRayIdentificationAlg.h.

Constructor & Destructor Documentation

lar_content::CheatingCosmicRayIdentificationAlg::CheatingCosmicRayIdentificationAlg

(

)

Default constructor.

Definition at line 22 of file CheatingCosmicRayIdentificationAlg.cc.

                                                                       : m_maxNeutrinoFraction(0.5f)
{
}

Member Function Documentation

StatusCode lar_content::CheatingCosmicRayIdentificationAlg::ReadSettings ( const pandora::TiXmlHandle  xmlHandle )

private

Definition at line 73 of file CheatingCosmicRayIdentificationAlg.cc.

{
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "InputPfoListName", m_inputPfoListName));

    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "OutputPfoListName", m_outputPfoListName));

    m_inputDaughterPfoListName = m_inputPfoListName;
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "InputDaughterPfoListName", m_inputDaughterPfoListName));

    m_outputDaughterPfoListName = m_outputPfoListName;
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "OutputDaughterPfoListName", m_outputDaughterPfoListName));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MaxNeutrinoFraction", m_maxNeutrinoFraction));

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::CheatingCosmicRayIdentificationAlg::Run ( )

private

Definition at line 28 of file CheatingCosmicRayIdentificationAlg.cc.

{
    const PfoList *pPfoList(nullptr);
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_INITIALIZED, !=, PandoraContentApi::GetList(*this, m_inputPfoListName, pPfoList));

    if (!pPfoList)
    {
        if (PandoraContentApi::GetSettings(*this)->ShouldDisplayAlgorithmInfo())
            std::cout << "CheatingCosmicRayIdentificationAlg: pfo list " << m_inputPfoListName << " unavailable." << std::endl;

        return STATUS_CODE_SUCCESS;
    }

    PfoList outputPfoList, outputDaughterPfoList;

    for (const ParticleFlowObject *const pPfo : *pPfoList)
    {
        if (!pPfo->GetParentPfoList().empty())
            continue;

        PfoList downstreamPfos;
        LArPfoHelper::GetAllDownstreamPfos(pPfo, downstreamPfos);

        float thisNeutrinoWeight(0.f), thisTotalWeight(0.f);
        CheatingSliceIdBaseTool::GetTargetParticleWeight(&downstreamPfos, thisNeutrinoWeight, thisTotalWeight, LArMCParticleHelper::IsNeutrino);

        if ((thisTotalWeight < std::numeric_limits<float>::epsilon()) || ((thisNeutrinoWeight / thisTotalWeight) < m_maxNeutrinoFraction))
            outputPfoList.push_back(pPfo);
    }

    if (!outputPfoList.empty())
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::SaveList(*this, m_inputPfoListName, m_outputPfoListName, outputPfoList));

        if (!outputDaughterPfoList.empty())
            PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=,
                PandoraContentApi::SaveList(*this, m_inputDaughterPfoListName, m_outputDaughterPfoListName, outputDaughterPfoList));
    }

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

std::string lar_content::CheatingCosmicRayIdentificationAlg::m_inputDaughterPfoListName

private

The input daughter pfo list name (if not specified, will assume same as main input list)

Definition at line 33 of file CheatingCosmicRayIdentificationAlg.h.

std::string lar_content::CheatingCosmicRayIdentificationAlg::m_inputPfoListName

private

The input pfo list name.

Definition at line 31 of file CheatingCosmicRayIdentificationAlg.h.

float lar_content::CheatingCosmicRayIdentificationAlg::m_maxNeutrinoFraction

private

The maximum true neutrino fraction in a particle to be labelled as a cosmic ray.

Definition at line 35 of file CheatingCosmicRayIdentificationAlg.h.

std::string lar_content::CheatingCosmicRayIdentificationAlg::m_outputDaughterPfoListName

private

The output daughter pfo list name (if not specified, will assume same as main output list)

Definition at line 34 of file CheatingCosmicRayIdentificationAlg.h.

std::string lar_content::CheatingCosmicRayIdentificationAlg::m_outputPfoListName

private

The output pfo list name.

Definition at line 32 of file CheatingCosmicRayIdentificationAlg.h.

---

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

• CheatingCosmicRayIdentificationAlg.h
• CheatingCosmicRayIdentificationAlg.cc