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