LArPandoraInput.h

Go to the documentation of this file.

#ifndef LAR_PANDORA_INPUT_H
#define LAR_PANDORA_INPUT_H 1

namespace detinfo {
  class DetectorPropertiesData;
}

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

#include "larpandoracontent/LArObjects/LArMCParticle.h"

namespace pandora {
  class Pandora;
}

namespace lar_pandora {

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

      const pandora::Pandora* m_pPrimaryPandora; 
      bool m_useHitWidths;                       
      bool m_useBirksCorrection;                 
      bool m_useActiveBoundingBox;               
      int m_uidOffset;                           
      int m_hitCounterOffset;                    
      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 art::Event& evt,
                                    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:
    typedef std::map<std::string, lar_content::MCProcess> MCProcessMap;

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

    static void GetTrueStartAndEndPoints(const geo::TPCID& ref_tpcid,
                                         const art::Ptr<simb::MCParticle>& particle,
                                         int& startT,
                                         int& endT);

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

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

    static void FillMCProcessMap(MCProcessMap& processMap);
  };

} // namespace lar_pandora

#endif // #ifndef LAR_PANDORA_INPUT_H