lar_content::EventWritingAlgorithm Class Reference

EventWritingAlgorithm class. More...

#include "EventWritingAlgorithm.h"

Inheritance diagram for lar_content::EventWritingAlgorithm:

Public Member Functions
	EventWritingAlgorithm ()
	Default constructor. More...
	~EventWritingAlgorithm ()
	Destructor. More...

Private Member Functions
pandora::StatusCode	Initialize ()
pandora::StatusCode	Run ()
bool	PassNuanceCodeFilter () const
	Whether current event passes nuance code filter. More...
bool	PassMCParticleFilter () const
	Whether current event passes mc particle constituent filter. More...
bool	PassNeutrinoVertexFilter () const
	Whether current event passes neutrino vertex position filter (e.g. fiducial volume cut) More...
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)

Private Attributes
pandora::FileType	m_geometryFileType
	The geometry file type. More...
pandora::FileType	m_eventFileType
	The event file type. More...
pandora::FileWriter *	m_pEventFileWriter
	Address of the event file writer. More...
pandora::FileWriter *	m_pGeometryFileWriter
	Address of the geometry file writer. More...
bool	m_shouldWriteGeometry
	Whether to write geometry to a specified file. More...
bool	m_writtenGeometry
	Whether geometry has been written. More...
std::string	m_geometryFileName
	Name of the output geometry file. More...
bool	m_shouldWriteEvents
	Whether to write events to a specified file. More...
std::string	m_eventFileName
	Name of the output event file. More...
bool	m_shouldWriteMCRelationships
	Whether to write mc relationship information to the events file. More...
bool	m_shouldWriteTrackRelationships
	Whether to write track relationship information to the events file. More...
bool	m_shouldOverwriteEventFile
	Whether to overwrite existing event file with specified name, or append. More...
bool	m_shouldOverwriteGeometryFile
	Whether to overwrite existing geometry file with specified name, or append. More...
bool	m_useLArCaloHits
	Whether to write lar calo hits, or standard pandora calo hits. More...
unsigned int	m_larCaloHitVersion
	LArCaloHit version for LArCaloHitFactory. More...
bool	m_useLArMCParticles
	Whether to write lar mc particles, or standard pandora mc particles. More...
bool	m_shouldFilterByNuanceCode
	Whether to filter output by nuance code. More...
int	m_filterNuanceCode
	The filter nuance code (required if specify filter by nuance code) More...
bool	m_shouldFilterByMCParticles
	Whether to filter output by mc particle constituents. More...
bool	m_neutrinoInducedOnly
	Whether to consider only mc particles that were neutrino induced. More...
unsigned int	m_matchingMinPrimaryHits
	The minimum number of mc primary hits used in matching scheme. More...
unsigned int	m_nNonNeutrons
	The requested number of mc primaries that are not neutrons. More...
unsigned int	m_nMuons
	The requested number of mc primaries that are muons. More...
unsigned int	m_nElectrons
	The requested number of mc primaries that are electrons. More...
unsigned int	m_nProtons
	The requested number of mc primaries that are protons. More...
unsigned int	m_nPhotons
	The requested number of mc primaries that are photons. More...
unsigned int	m_nChargedPions
	The requested number of mc primaries that are charged pions. More...
bool	m_shouldFilterByNeutrinoVertex
	Whether to filter output by neutrino vertex position (e.g. fiducial volume cut) More...
float	m_detectorHalfLengthX
	Half length of detector in x dimension. More...
float	m_detectorHalfLengthY
	Half length of detector in y dimension. More...
float	m_detectorHalfLengthZ
	Half length of detector in z dimension. More...
float	m_coordinateOffsetX
	Origin offset (from detector corner) in x dimension. More...
float	m_coordinateOffsetY
	Origin offset (from detector corner) in y dimension. More...
float	m_coordinateOffsetZ
	Origin offset (from detector corner) in z dimension. More...
float	m_selectedBorderX
	Required distance from detector edge in x dimension. More...
float	m_selectedBorderY
	Required distance from detector edge in y dimension. More...
float	m_selectedBorderZ
	Required distance from detector edge in z dimension. More...

Detailed Description

EventWritingAlgorithm class.

Definition at line 28 of file EventWritingAlgorithm.h.

Constructor & Destructor Documentation

lar_content::EventWritingAlgorithm::EventWritingAlgorithm

(

)

Default constructor.

Definition at line 27 of file EventWritingAlgorithm.cc.

                                             :
    m_geometryFileType(UNKNOWN_FILE_TYPE),
    m_eventFileType(UNKNOWN_FILE_TYPE),
    m_pEventFileWriter(nullptr),
    m_pGeometryFileWriter(nullptr),
    m_shouldWriteGeometry(false),
    m_writtenGeometry(false),
    m_shouldWriteEvents(true),
    m_shouldWriteMCRelationships(true),
    m_shouldWriteTrackRelationships(true),
    m_shouldOverwriteEventFile(false),
    m_shouldOverwriteGeometryFile(false),
    m_useLArCaloHits(true),
    m_larCaloHitVersion(1),
    m_useLArMCParticles(true),
    m_shouldFilterByNuanceCode(false),
    m_filterNuanceCode(0),
    m_shouldFilterByMCParticles(false),
    m_neutrinoInducedOnly(false),
    m_matchingMinPrimaryHits(15),
    m_nNonNeutrons(0),
    m_nMuons(0),
    m_nElectrons(0),
    m_nProtons(0),
    m_nPhotons(0),
    m_nChargedPions(0),
    m_shouldFilterByNeutrinoVertex(false),
    m_detectorHalfLengthX(-1.f),
    m_detectorHalfLengthY(-1.f),
    m_detectorHalfLengthZ(-1.f),
    m_coordinateOffsetX(0.f),
    m_coordinateOffsetY(0.f),
    m_coordinateOffsetZ(0.f),
    m_selectedBorderX(-1.f),
    m_selectedBorderY(-1.f),
    m_selectedBorderZ(-1.f)
{
}

lar_content::EventWritingAlgorithm::~EventWritingAlgorithm

(

)

Destructor.

Definition at line 68 of file EventWritingAlgorithm.cc.

References m_pEventFileWriter, and m_pGeometryFileWriter.

{
    delete m_pEventFileWriter;
    delete m_pGeometryFileWriter;
}

Member Function Documentation

StatusCode lar_content::EventWritingAlgorithm::Initialize ( )

private

Definition at line 76 of file EventWritingAlgorithm.cc.

References m_eventFileName, m_eventFileType, m_geometryFileName, m_geometryFileType, m_larCaloHitVersion, m_pEventFileWriter, m_pGeometryFileWriter, m_shouldOverwriteEventFile, m_shouldOverwriteGeometryFile, m_shouldWriteEvents, m_shouldWriteGeometry, m_useLArCaloHits, and m_useLArMCParticles.

{
    if (m_shouldWriteGeometry)
    {
        const FileMode fileMode(m_shouldOverwriteGeometryFile ? OVERWRITE : APPEND);

        if (BINARY == m_geometryFileType)
        {
            m_pGeometryFileWriter = new BinaryFileWriter(this->GetPandora(), m_geometryFileName, fileMode);
        }
        else if (XML == m_geometryFileType)
        {
            m_pGeometryFileWriter = new XmlFileWriter(this->GetPandora(), m_geometryFileName, fileMode);
        }
        else
        {
            return STATUS_CODE_FAILURE;
        }
    }

    if (m_shouldWriteEvents)
    {
        const FileMode fileMode(m_shouldOverwriteEventFile ? OVERWRITE : APPEND);

        if (BINARY == m_eventFileType)
        {
            m_pEventFileWriter = new BinaryFileWriter(this->GetPandora(), m_eventFileName, fileMode);
        }
        else if (XML == m_eventFileType)
        {
            m_pEventFileWriter = new XmlFileWriter(this->GetPandora(), m_eventFileName, fileMode);
        }
        else
        {
            return STATUS_CODE_FAILURE;
        }

        if (m_useLArCaloHits)
            m_pEventFileWriter->SetFactory(new LArCaloHitFactory(m_larCaloHitVersion));

        if (m_useLArMCParticles)
            m_pEventFileWriter->SetFactory(new LArMCParticleFactory);
    }

    return STATUS_CODE_SUCCESS;
}

bool lar_content::EventWritingAlgorithm::PassMCParticleFilter ( ) const

private

Whether current event passes mc particle constituent filter.

Returns

boolean

Definition at line 179 of file EventWritingAlgorithm.cc.

References util::abs(), lar_content::LArMCParticleHelper::IsBeamNeutrinoFinalState(), lar_content::LArMCParticleHelper::IsBeamParticle(), lar_content::LArMCParticleHelper::IsCosmicRay(), m_nChargedPions, m_nElectrons, m_neutrinoInducedOnly, m_nMuons, m_nNonNeutrons, m_nPhotons, m_nProtons, lar_content::LArMCParticleHelper::SelectReconstructableMCParticles(), and lar_content::LArMCParticleHelper::SortByMomentum().

Referenced by Run().

{
    const MCParticleList *pMCParticleList = nullptr;
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pMCParticleList));

    const CaloHitList *pCaloHitList = nullptr;
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pCaloHitList));

    LArMCParticleHelper::PrimaryParameters parameters;
    LArMCParticleHelper::MCContributionMap mcParticlesToGoodHitsMap;
    LArMCParticleHelper::SelectReconstructableMCParticles(
        pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsBeamNeutrinoFinalState, mcParticlesToGoodHitsMap);

    if (!m_neutrinoInducedOnly)
    {
        LArMCParticleHelper::SelectReconstructableMCParticles(
            pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsBeamParticle, mcParticlesToGoodHitsMap);
        LArMCParticleHelper::SelectReconstructableMCParticles(
            pMCParticleList, pCaloHitList, parameters, LArMCParticleHelper::IsCosmicRay, mcParticlesToGoodHitsMap);
    }

    unsigned int nNonNeutrons(0), nMuons(0), nElectrons(0), nProtons(0), nPhotons(0), nChargedPions(0);

    MCParticleVector mcPrimaryVector;
    for (const auto &mapEntry : mcParticlesToGoodHitsMap)
        mcPrimaryVector.push_back(mapEntry.first);
    std::sort(mcPrimaryVector.begin(), mcPrimaryVector.end(), LArMCParticleHelper::SortByMomentum);

    for (const MCParticle *const pMCPrimary : mcPrimaryVector)
    {
        const unsigned int particleId(std::abs(pMCPrimary->GetParticleId()));
        if (NEUTRON != particleId)
            ++nNonNeutrons;

        if (MU_MINUS == particleId)
            ++nMuons;
        else if (E_MINUS == particleId)
            ++nElectrons;
        else if (PROTON == particleId)
            ++nProtons;
        else if (PHOTON == particleId)
            ++nPhotons;
        else if (PI_PLUS == particleId)
            ++nChargedPions;
    }

    if ((nNonNeutrons == m_nNonNeutrons) && (nMuons == m_nMuons) && (nElectrons == m_nElectrons) && (nProtons == m_nProtons) &&
        (nPhotons == m_nPhotons) && (nChargedPions == m_nChargedPions))
    {
        return true;
    }

    return false;
}

bool lar_content::EventWritingAlgorithm::PassNeutrinoVertexFilter ( ) const

private

Whether current event passes neutrino vertex position filter (e.g. fiducial volume cut)

Returns

boolean

Definition at line 236 of file EventWritingAlgorithm.cc.

References lar_content::LArMCParticleHelper::GetTrueNeutrinos(), m_coordinateOffsetX, m_coordinateOffsetY, m_coordinateOffsetZ, m_detectorHalfLengthX, m_detectorHalfLengthY, m_detectorHalfLengthZ, m_selectedBorderX, m_selectedBorderY, and m_selectedBorderZ.

Referenced by Run().

{
    const MCParticleList *pMCParticleList = nullptr;
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pMCParticleList));

    MCParticleVector mcNeutrinoVector;
    LArMCParticleHelper::GetTrueNeutrinos(pMCParticleList, mcNeutrinoVector);

    for (const MCParticle *const pMCNeutrino : mcNeutrinoVector)
    {
        const CartesianVector &neutrinoInteractionVertex(pMCNeutrino->GetEndpoint());

        if ((neutrinoInteractionVertex.GetX() < (m_detectorHalfLengthX - m_coordinateOffsetX - m_selectedBorderX)) &&
            (neutrinoInteractionVertex.GetX() > (-m_coordinateOffsetX + m_selectedBorderX)) &&
            (neutrinoInteractionVertex.GetY() < (m_detectorHalfLengthY - m_coordinateOffsetY - m_selectedBorderY)) &&
            (neutrinoInteractionVertex.GetY() > (-m_coordinateOffsetY + m_selectedBorderY)) &&
            (neutrinoInteractionVertex.GetZ() < (m_detectorHalfLengthZ - m_coordinateOffsetZ - m_selectedBorderZ)) &&
            (neutrinoInteractionVertex.GetZ() > (-m_coordinateOffsetZ + m_selectedBorderZ)))
        {
            return true;
        }
    }

    return false;
}

bool lar_content::EventWritingAlgorithm::PassNuanceCodeFilter ( ) const

private

Whether current event passes nuance code filter.

Returns

boolean

Definition at line 158 of file EventWritingAlgorithm.cc.

References lar_content::LArMCParticle::GetNuanceCode(), lar_content::LArMCParticleHelper::GetTrueNeutrinos(), and m_filterNuanceCode.

Referenced by Run().

{
    const MCParticleList *pMCParticleList = nullptr;
    PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pMCParticleList));

    MCParticleVector mcNeutrinoList;
    LArMCParticleHelper::GetTrueNeutrinos(pMCParticleList, mcNeutrinoList);

    for (const MCParticle *const pMCNeutrino : mcNeutrinoList)
    {
        const LArMCParticle *const pLArMCNeutrino = dynamic_cast<const LArMCParticle *>(pMCNeutrino);

        if (pLArMCNeutrino && (pLArMCNeutrino->GetNuanceCode() == m_filterNuanceCode))
            return true;
    }

    return false;
}

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

private

Definition at line 264 of file EventWritingAlgorithm.cc.

References m_coordinateOffsetX, m_coordinateOffsetY, m_coordinateOffsetZ, m_detectorHalfLengthX, m_detectorHalfLengthY, m_detectorHalfLengthZ, m_eventFileName, m_eventFileType, m_filterNuanceCode, m_geometryFileName, m_geometryFileType, m_larCaloHitVersion, m_matchingMinPrimaryHits, m_nChargedPions, m_nElectrons, m_neutrinoInducedOnly, m_nMuons, m_nNonNeutrons, m_nPhotons, m_nProtons, m_selectedBorderX, m_selectedBorderY, m_selectedBorderZ, m_shouldFilterByMCParticles, m_shouldFilterByNeutrinoVertex, m_shouldFilterByNuanceCode, m_shouldOverwriteEventFile, m_shouldOverwriteGeometryFile, m_shouldWriteEvents, m_shouldWriteGeometry, m_shouldWriteMCRelationships, m_shouldWriteTrackRelationships, m_useLArCaloHits, and m_useLArMCParticles.

{
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ShouldWriteGeometry", m_shouldWriteGeometry));

    if (m_shouldWriteGeometry)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "GeometryFileName", m_geometryFileName));

        std::string fileExtension(m_geometryFileName.substr(m_geometryFileName.find_last_of(".")));
        std::transform(fileExtension.begin(), fileExtension.end(), fileExtension.begin(), ::tolower);

        if (std::string(".xml") == fileExtension)
        {
            m_geometryFileType = XML;
        }
        else if (std::string(".pndr") == fileExtension)
        {
            m_geometryFileType = BINARY;
        }
        else
        {
            std::cout << "EventReadingAlgorithm: Unknown geometry file type specified " << std::endl;
            return STATUS_CODE_INVALID_PARAMETER;
        }
    }

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "ShouldWriteEvents", m_shouldWriteEvents));

    if (m_shouldWriteEvents)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "EventFileName", m_eventFileName));

        std::string fileExtension(m_eventFileName.substr(m_eventFileName.find_last_of(".")));
        std::transform(fileExtension.begin(), fileExtension.end(), fileExtension.begin(), ::tolower);

        if (std::string(".xml") == fileExtension)
        {
            m_eventFileType = XML;
        }
        else if (std::string(".pndr") == fileExtension)
        {
            m_eventFileType = BINARY;
        }
        else
        {
            std::cout << "EventReadingAlgorithm: Unknown event file type specified " << std::endl;
            return STATUS_CODE_INVALID_PARAMETER;
        }
    }

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldWriteMCRelationships", m_shouldWriteMCRelationships));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldWriteTrackRelationships", m_shouldWriteTrackRelationships));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldOverwriteEventFile", m_shouldOverwriteEventFile));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldOverwriteGeometryFile", m_shouldOverwriteGeometryFile));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "UseLArCaloHits", m_useLArCaloHits));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "LArCaloHitVersion", m_larCaloHitVersion));

    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "UseLArMCParticles", m_useLArMCParticles));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldFilterByNuanceCode", m_shouldFilterByNuanceCode));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldFilterByMCParticles", m_shouldFilterByMCParticles));

    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=,
        XmlHelper::ReadValue(xmlHandle, "ShouldFilterByNeutrinoVertex", m_shouldFilterByNeutrinoVertex));

    if (m_shouldFilterByNuanceCode)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "FilterNuanceCode", m_filterNuanceCode));
    }

    if (m_shouldFilterByMCParticles)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NeutrinoInducedOnly", m_neutrinoInducedOnly));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "MatchingMinPrimaryHits", m_matchingMinPrimaryHits));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NNonNeutrons", m_nNonNeutrons));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NMuons", m_nMuons));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NElectrons", m_nElectrons));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NProtons", m_nProtons));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NPhotons", m_nPhotons));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "NChargedPions", m_nChargedPions));
    }

    if (m_shouldFilterByNeutrinoVertex)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "DetectorHalfLengthX", m_detectorHalfLengthX));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "DetectorHalfLengthY", m_detectorHalfLengthY));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "DetectorHalfLengthZ", m_detectorHalfLengthZ));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "CoordinateOffsetX", m_coordinateOffsetX));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "CoordinateOffsetY", m_coordinateOffsetY));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "CoordinateOffsetZ", m_coordinateOffsetZ));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "SelectedBorderX", m_selectedBorderX));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "SelectedBorderY", m_selectedBorderY));
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, XmlHelper::ReadValue(xmlHandle, "SelectedBorderZ", m_selectedBorderZ));
    }

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::EventWritingAlgorithm::Run ( )

private

Definition at line 125 of file EventWritingAlgorithm.cc.

References m_pEventFileWriter, m_pGeometryFileWriter, m_shouldFilterByMCParticles, m_shouldFilterByNeutrinoVertex, m_shouldFilterByNuanceCode, m_shouldWriteEvents, m_shouldWriteGeometry, m_shouldWriteMCRelationships, m_shouldWriteTrackRelationships, m_writtenGeometry, PassMCParticleFilter(), PassNeutrinoVertexFilter(), and PassNuanceCodeFilter().

{
    // ATTN Should complete geometry creation in LArSoft begin job, but some channel status service functionality unavailable at that point
    if (!m_writtenGeometry && m_pGeometryFileWriter && m_shouldWriteGeometry)
    {
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, m_pGeometryFileWriter->WriteGeometry());
        m_writtenGeometry = true;
    }

    bool matchNuanceCode(!m_shouldFilterByNuanceCode || this->PassNuanceCodeFilter());
    bool matchParticles(!m_shouldFilterByMCParticles || this->PassMCParticleFilter());
    bool matchNeutrinoVertexPosition(!m_shouldFilterByNeutrinoVertex || this->PassNeutrinoVertexFilter());

    if (matchNuanceCode && matchParticles && matchNeutrinoVertexPosition && m_pEventFileWriter && m_shouldWriteEvents)
    {
        const CaloHitList *pCaloHitList = nullptr;
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pCaloHitList));

        const TrackList *pTrackList = nullptr;
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pTrackList));

        const MCParticleList *pMCParticleList = nullptr;
        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=, PandoraContentApi::GetCurrentList(*this, pMCParticleList));

        PANDORA_RETURN_RESULT_IF(STATUS_CODE_SUCCESS, !=,
            m_pEventFileWriter->WriteEvent(*pCaloHitList, *pTrackList, *pMCParticleList, m_shouldWriteMCRelationships, m_shouldWriteTrackRelationships));
    }

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

float lar_content::EventWritingAlgorithm::m_coordinateOffsetX

private

Origin offset (from detector corner) in x dimension.

Definition at line 108 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_coordinateOffsetY

private

Origin offset (from detector corner) in y dimension.

Definition at line 109 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_coordinateOffsetZ

private

Origin offset (from detector corner) in z dimension.

Definition at line 110 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_detectorHalfLengthX

private

Half length of detector in x dimension.

Definition at line 105 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_detectorHalfLengthY

private

Half length of detector in y dimension.

Definition at line 106 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_detectorHalfLengthZ

private

Half length of detector in z dimension.

Definition at line 107 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

std::string lar_content::EventWritingAlgorithm::m_eventFileName

private

Name of the output event file.

Definition at line 79 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

pandora::FileType lar_content::EventWritingAlgorithm::m_eventFileType

private

The event file type.

Definition at line 69 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

int lar_content::EventWritingAlgorithm::m_filterNuanceCode

private

The filter nuance code (required if specify filter by nuance code)

Definition at line 92 of file EventWritingAlgorithm.h.

Referenced by PassNuanceCodeFilter(), and ReadSettings().

std::string lar_content::EventWritingAlgorithm::m_geometryFileName

private

Name of the output geometry file.

Definition at line 76 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

pandora::FileType lar_content::EventWritingAlgorithm::m_geometryFileType

private

The geometry file type.

Definition at line 68 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_larCaloHitVersion

private

LArCaloHit version for LArCaloHitFactory.

Definition at line 88 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_matchingMinPrimaryHits

private

The minimum number of mc primary hits used in matching scheme.

Definition at line 96 of file EventWritingAlgorithm.h.

Referenced by ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nChargedPions

private

The requested number of mc primaries that are charged pions.

Definition at line 102 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nElectrons

private

The requested number of mc primaries that are electrons.

Definition at line 99 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_neutrinoInducedOnly

private

Whether to consider only mc particles that were neutrino induced.

Definition at line 95 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nMuons

private

The requested number of mc primaries that are muons.

Definition at line 98 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nNonNeutrons

private

The requested number of mc primaries that are not neutrons.

Definition at line 97 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nPhotons

private

The requested number of mc primaries that are photons.

Definition at line 101 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

unsigned int lar_content::EventWritingAlgorithm::m_nProtons

private

The requested number of mc primaries that are protons.

Definition at line 100 of file EventWritingAlgorithm.h.

Referenced by PassMCParticleFilter(), and ReadSettings().

pandora::FileWriter* lar_content::EventWritingAlgorithm::m_pEventFileWriter

private

Address of the event file writer.

Definition at line 71 of file EventWritingAlgorithm.h.

Referenced by Initialize(), Run(), and ~EventWritingAlgorithm().

pandora::FileWriter* lar_content::EventWritingAlgorithm::m_pGeometryFileWriter

private

Address of the geometry file writer.

Definition at line 72 of file EventWritingAlgorithm.h.

Referenced by Initialize(), Run(), and ~EventWritingAlgorithm().

float lar_content::EventWritingAlgorithm::m_selectedBorderX

private

Required distance from detector edge in x dimension.

Definition at line 111 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_selectedBorderY

private

Required distance from detector edge in y dimension.

Definition at line 112 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

float lar_content::EventWritingAlgorithm::m_selectedBorderZ

private

Required distance from detector edge in z dimension.

Definition at line 113 of file EventWritingAlgorithm.h.

Referenced by PassNeutrinoVertexFilter(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_shouldFilterByMCParticles

private

Whether to filter output by mc particle constituents.

Definition at line 94 of file EventWritingAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldFilterByNeutrinoVertex

private

Whether to filter output by neutrino vertex position (e.g. fiducial volume cut)

Definition at line 104 of file EventWritingAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldFilterByNuanceCode

private

Whether to filter output by nuance code.

Definition at line 91 of file EventWritingAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldOverwriteEventFile

private

Whether to overwrite existing event file with specified name, or append.

Definition at line 84 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_shouldOverwriteGeometryFile

private

Whether to overwrite existing geometry file with specified name, or append.

Definition at line 85 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_shouldWriteEvents

private

Whether to write events to a specified file.

Definition at line 78 of file EventWritingAlgorithm.h.

Referenced by Initialize(), ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldWriteGeometry

private

Whether to write geometry to a specified file.

Definition at line 74 of file EventWritingAlgorithm.h.

Referenced by Initialize(), ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldWriteMCRelationships

private

Whether to write mc relationship information to the events file.

Definition at line 81 of file EventWritingAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_shouldWriteTrackRelationships

private

Whether to write track relationship information to the events file.

Definition at line 82 of file EventWritingAlgorithm.h.

Referenced by ReadSettings(), and Run().

bool lar_content::EventWritingAlgorithm::m_useLArCaloHits

private

Whether to write lar calo hits, or standard pandora calo hits.

Definition at line 87 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_useLArMCParticles

private

Whether to write lar mc particles, or standard pandora mc particles.

Definition at line 89 of file EventWritingAlgorithm.h.

Referenced by Initialize(), and ReadSettings().

bool lar_content::EventWritingAlgorithm::m_writtenGeometry

private

Whether geometry has been written.

Definition at line 75 of file EventWritingAlgorithm.h.

Referenced by Run().

---

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

• larpandoracontent/v04_08_01/source/larpandoracontent/LArPersistency/EventWritingAlgorithm.h
• larpandoracontent/v04_08_01/source/larpandoracontent/LArPersistency/EventWritingAlgorithm.cc