AuxDetGeometryCore.h

Go to the documentation of this file.

#ifndef GEO_AUXDETGEOMETRYCORE_H
#define GEO_AUXDETGEOMETRYCORE_H

// LArSoft libraries
#include "larcorealg/Geometry/AuxDetChannelMapAlg.h"
#include "larcorealg/Geometry/AuxDetGeo.h"
#include "larcoreobj/SimpleTypesAndConstants/geo_vectors.h"

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

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

namespace geo {

  // Forward declarations within namespace.
  class AuxDetSensitiveGeo;

  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);

    // 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 const& ROOTFile() const { return fROOTfile; }

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

    std::string const& 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(Point_t const& worldLoc, double tolerance = 0) const;

    void FindAuxDetSensitiveAtPosition(Point_t const& worldLoc,
                                       size_t& adg,
                                       size_t& sv,
                                       double tolerance = 0) const;

    AuxDetGeo const& PositionToAuxDet(Point_t const& worldLoc,
                                      unsigned int& ad,
                                      double tolerance = 0) const;

    AuxDetSensitiveGeo const& PositionToAuxDetSensitive(Point_t const& worldLoc,
                                                        size_t& ad,
                                                        size_t& sv,
                                                        double tolerance = 0) const;

    uint32_t PositionToAuxDetChannel(Point_t const& worldLoc, size_t& ad, size_t& sv) const;

    Point_t AuxDetChannelToPosition(std::string const& auxDetName, uint32_t channel) const;

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

    AuxDetSensitiveGeo const& ChannelToAuxDetSensitive(
      std::string const& auxDetName,
      uint32_t 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::unique_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 ClearGeometry();

    AuxDetGeometryData_t fGeoData; 

    std::string fDetectorName;              
    std::string fGDMLfile;                  
    std::string fROOTfile;                  
    fhicl::ParameterSet fBuilderParameters; 
    std::unique_ptr<AuxDetChannelMapAlg const>
      fChannelMapAlg; 
  };                  // class GeometryCore

} // namespace geo

#endif // GEO_AUXDETGEOMETRYCORE_H