AuxDetGeometryCore.h

Go to the documentation of this file.

#ifndef GEO_AUXDETGEOMETRYCORE_H
#define GEO_AUXDETGEOMETRYCORE_H

// LArSoft libraries

// Framework and infrastructure libraries
#include "fhiclcpp/ParameterSet.h"

// ROOT libraries
#include <TVector3.h>

// C/C++ standard libraries
#include <cstddef> // size_t
#include <string>
#include <vector>
#include <memory> // std::shared_ptr<>

#include "larcorealg/Geometry/AuxDetChannelMapAlg.h"

// ROOT class prototypes
class TGeoManager;
class TGeoNode;
class TGeoMaterial;


namespace geo {
  
  
  // Forward declarations within namespace.
  class AuxDetGeo;
  class AuxDetSensitiveGeo;
  class AuxDetGeometryCore;
  
  
  struct AuxDetGeometryData_t {
    
    using AuxDetList_t = std::vector<AuxDetGeo*>;
    
    AuxDetList_t   auxDets;   
    
  }; // AuxDetGeometryData_t
  
  
  class AuxDetGeometryCore {
  public:
    
    using AuxDetList_t = AuxDetGeometryData_t::AuxDetList_t;    
    
    AuxDetGeometryCore(fhicl::ParameterSet const& pset);
    
    ~AuxDetGeometryCore();
   
    // You shall not copy or move or assign me!
    AuxDetGeometryCore(AuxDetGeometryCore const&) = delete;
    AuxDetGeometryCore(AuxDetGeometryCore&&) = delete;
    AuxDetGeometryCore& operator= (AuxDetGeometryCore const&) = delete;
    AuxDetGeometryCore& operator= (AuxDetGeometryCore&&) = delete;
    
    
    std::string ROOTFile() const { return fROOTfile; }

    std::string GDMLFile() const { return fGDMLfile; }
            

    std::string DetectorName() const { return fDetectorName; }    
    
    //
    // object description and information
    //
        
    
    //
    // group features
    //
    
    unsigned int NAuxDets() const { return AuxDets().size(); }
    
    unsigned int NAuxDetSensitive(size_t const& aid) const;
    
    //
    // access
    //
    
    std::vector<AuxDetGeo*> const& AuxDetGeoVec() const { return AuxDets(); }
    
    AuxDetGeo const& AuxDet(unsigned int const ad = 0) const;
    
    unsigned int FindAuxDetAtPosition(double const worldLoc[3]) const;
    
    void  FindAuxDetSensitiveAtPosition(double const worldLoc[3],
                                        size_t     & adg,
                                        size_t     & sv) const;
    
    AuxDetGeo const& PositionToAuxDet(double const worldLoc[3], 
                                      unsigned int &ad) const;
    
    const AuxDetSensitiveGeo& PositionToAuxDetSensitive(double const worldLoc[3],
                                                        size_t     & ad,
                                                        size_t     & sv) const;
    
    uint32_t                 PositionToAuxDetChannel(double const worldLoc[3],
                                                     size_t     & ad,
                                                     size_t     & sv) const;
    TVector3                 AuxDetChannelToPosition(uint32_t    const& channel,
                                                     std::string const& auxDetName) const;


    const AuxDetGeo&         ChannelToAuxDet(std::string const& auxDetName,
                                             uint32_t    const& channel) const; // return the AuxDetGeo for the given detector 
                                                                                // name and channel

    const AuxDetSensitiveGeo& ChannelToAuxDetSensitive(std::string const& auxDetName,
                                                       uint32_t    const& channel) const; // return the AuxDetSensitiveGeo for the given
    
    
    void LoadGeometryFile(std::string gdmlfile, std::string rootfile);
    
    
    bool hasAuxDetChannelMap() const { return bool(fChannelMapAlg); }
    
    
    void ApplyChannelMap(std::shared_ptr<geo::AuxDetChannelMapAlg> pChannelMap);
    
    
  protected:

    geo::AuxDetChannelMapAlg const* AuxDetChannelMap() const { return fChannelMapAlg.get(); }
    
    AuxDetList_t&       AuxDets()       { return fGeoData.auxDets; }
    AuxDetList_t const& AuxDets() const { return fGeoData.auxDets; }
    
  private:
    
    void FindAuxDet(std::vector<const TGeoNode*>& path, unsigned int depth);
    
    void MakeAuxDet(std::vector<const TGeoNode*>& path, int depth);
    
    void ClearGeometry();
    
    AuxDetGeometryData_t fGeoData;  
    
    std::string    fDetectorName;   
    std::string    fGDMLfile;       
    std::string    fROOTfile;       
    std::shared_ptr<const geo::AuxDetChannelMapAlg> fChannelMapAlg;  
  }; // class GeometryCore
  
} // namespace geo




#endif // GEO_AUXDETGEOMETRYCORE_H