CheatingNeutrinoCreationAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_CHEATING_NEUTRINO_CREATION_ALGORITHM_H
#define LAR_CHEATING_NEUTRINO_CREATION_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

#include "larpandoracontent/LArHelpers/LArMCParticleHelper.h"

#include <unordered_map>

namespace lar_content
{

class CheatingNeutrinoCreationAlgorithm : public pandora::Algorithm
{
public:
    CheatingNeutrinoCreationAlgorithm();

private:
    pandora::StatusCode Run();

    void GetMCNeutrinoVector(pandora::MCParticleVector &mcNeutrinoVector) const;

    void CreateAndSaveNeutrinoPfo(const pandora::MCParticle *const pMCNeutrino, const pandora::ParticleFlowObject *&pNeutrinoPfo) const;

    void AddNeutrinoVertex(const pandora::MCParticle *const pMCNeutrino, const pandora::ParticleFlowObject *const pNeutrinoPfo) const;

    typedef std::unordered_map<const pandora::MCParticle*, const pandora::ParticleFlowObject*> MCParticleToPfoMap;

    void GetMCParticleToDaughterPfoMap(MCParticleToPfoMap &mcParticleToPfoMap) const;

    void CreatePfoHierarchy(const pandora::MCParticle *const pParentMCParticle, const pandora::ParticleFlowObject *const pParentPfo,
        const MCParticleToPfoMap &mcParticleToPfoMap) const;

    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    bool                    m_collapseToPrimaryMCParticles; 

    std::string             m_mcParticleListName;           
    std::string             m_neutrinoPfoListName;          

    std::string             m_vertexListName;               
    pandora::StringVector   m_daughterPfoListNames;         

    float                   m_vertexTolerance;              
};

} // namespace lar_content

#endif // #ifndef LAR_CHEATING_NEUTRINO_CREATION_ALGORITHM_H