LArPandoraInput.h

Go to the documentation of this file.

#ifndef LAR_PANDORA_INPUT_H
#define LAR_PANDORA_INPUT_H 1

#include "lardata/ArtDataHelper/MVAReader.h"

#include "larpandora/LArPandoraInterface/ILArPandora.h"
#include "larpandora/LArPandoraInterface/LArPandoraHelper.h"
#include "larpandora/LArPandoraInterface/LArPandoraGeometry.h"

namespace lar_pandora
{

class LArPandoraInput
{
public:
    class Settings
    {
    public:
        Settings();

        const pandora::Pandora *m_pPrimaryPandora;          
        bool                    m_useHitWidths;             
        bool                    m_useBirksCorrection;       
        int                     m_uidOffset;                
        double                  m_dx_cm;                    
        double                  m_int_cm;                   
        double                  m_rad_cm;                   
        double                  m_dEdX_mip;                 
        double                  m_mips_max;                 
        double                  m_mips_if_negative;         
        double                  m_mips_to_gev;              
        double                  m_recombination_factor;     
    };

    static void CreatePandoraHits2D(const Settings &settings, const LArDriftVolumeMap &driftVolumeMap, const HitVector &hitVector, IdToHitMap &idToHitMap);

    static void CreatePandoraLArTPCs(const Settings &settings, const LArDriftVolumeList &driftVolumeList);

    static void CreatePandoraDetectorGaps(const Settings &settings, const LArDriftVolumeList &driftVolumeList,
        const LArDetectorGapList &listOfGaps);

    static void CreatePandoraReadoutGaps(const Settings &settings, const LArDriftVolumeMap &driftVolumeMap);

    static void CreatePandoraMCParticles(const Settings &settings, const MCTruthToMCParticles &truthToParticles,
        const MCParticlesToMCTruth &particlesToTruth, const RawMCParticleVector &generatorMCParticleVector);

    static void FindPrimaryParticles(const RawMCParticleVector &mcParticleVector, std::map<const simb::MCParticle, bool> &primaryMCParticleMap);

    static bool IsPrimaryMCParticle(const art::Ptr<simb::MCParticle> &mcParticle, std::map<const simb::MCParticle, bool> &primaryMCParticleMap);

    static void CreatePandoraMCLinks2D(const Settings &settings, const HitMap &hitMap, const HitsToTrackIDEs &hitToParticleMap);

private:
    static void GetTrueStartAndEndPoints(const Settings &settings, const art::Ptr<simb::MCParticle> &particle, int &startT, int &endT);

    static void GetTrueStartAndEndPoints(const unsigned int cstat, const unsigned int tpc, const art::Ptr<simb::MCParticle> &particle,
        int &startT, int &endT);

    static float GetTrueX0(const art::Ptr<simb::MCParticle> &particle, const int nT);

    static double GetMips(const Settings &settings, const double hit_Charge, const geo::View_t hit_View);
};

} // namespace lar_pandora

#endif // #ifndef LAR_PANDORA_INPUT_H