geo::AuxDetReadoutGeom Class Reference

#include "AuxDetReadoutGeom.h"

Public Member Functions
	AuxDetReadoutGeom (AuxDetReadoutInitializers initializers)
std::size_t	NearestAuxDet (Point_t const &point, std::vector< AuxDetGeo > const &auxDets, double tolerance=0, bool throwIfAbsent=true) const
std::size_t	NearestSensitiveAuxDet (Point_t const &point, std::vector< AuxDetGeo > const &auxDets, double tolerance=0, bool throwIfAbsent=true) const
std::pair< std::size_t, std::size_t >	ChannelToSensitiveAuxDet (std::string const &detName, std::uint32_t channel) const
Point_t	AuxDetChannelToPosition (std::uint32_t channel, std::string const &auxDetName, std::vector< AuxDetGeo > const &auxDets) const

Private Member Functions
std::size_t	DetNameToAuxDet (std::string const &detName) const

Private Attributes
std::map< std::size_t, std::string >	fADGeoToName
	map the AuxDetGeo index to the name More...
std::map< std::string, std::size_t >	fNameToADGeo
	map the names to the AuxDetGeo index More...
std::map< std::size_t, std::vector< chanAndSV > >	fADGeoToChannelAndSV

Detailed Description

Definition at line 42 of file AuxDetReadoutGeom.h.

Constructor & Destructor Documentation

geo::AuxDetReadoutGeom::AuxDetReadoutGeom ( AuxDetReadoutInitializers  initializers )

explicit

Definition at line 29 of file AuxDetReadoutGeom.cxx.

References geo::AuxDetReadoutInitializers::ADGeoToName, fADGeoToChannelAndSV, and fNameToADGeo.

    : fADGeoToName{std::move(initializers.ADGeoToName)}
    , fNameToADGeo{std::move(initializers.NameToADGeo)}
    , fADGeoToChannelAndSV{std::move(initializers.ADGeoToChannelAndSV)}
  {}

Member Function Documentation

geo::Point_t geo::AuxDetReadoutGeom::AuxDetChannelToPosition	(	std::uint32_t	channel,
		std::string const &	auxDetName,
		std::vector< AuxDetGeo > const &	auxDets
	)		const

Definition at line 137 of file AuxDetReadoutGeom.cxx.

References fADGeoToChannelAndSV, and fNameToADGeo.

  {
    // Figure out which detector we are in
    auto ad_it = fNameToADGeo.find(auxDetName);
    if (ad_it == fNameToADGeo.cend()) {
      throw cet::exception("CRTWireReadoutGeom") << "No AuxDetGeo with name " << auxDetName;
    }
    std::size_t const ad = ad_it->second;

    // Get the vector of channel and sensitive volume pairs
    auto it = fADGeoToChannelAndSV.find(ad);
    if (it == fADGeoToChannelAndSV.end()) {
      throw cet::exception("CRTWireReadoutGeom") << "No entry in channel and sensitive volume"
                                                 << " map for AuxDet index " << ad << " bail";
    }

    // Loop over the vector of channel and sensitive volumes to determine the sensitive
    // volume for this channel. Then get the origin of the sensitive volume in the world
    // coordinate system.
    for (auto const& [ch, svindex] : it->second) {
      if (ch == channel) { return auxDets[ad].SensitiveVolume(svindex).GetCenter(); }
    }

    return {};
  }

std::pair< std::size_t, std::size_t > geo::AuxDetReadoutGeom::ChannelToSensitiveAuxDet	(	std::string const &	detName,
		std::uint32_t	channel
	)		const

Definition at line 114 of file AuxDetReadoutGeom.cxx.

References DetNameToAuxDet(), and fADGeoToChannelAndSV.

  {
    std::size_t adGeoIdx = DetNameToAuxDet(detName);

    // look for the index of the sensitive volume for the given channel
    auto it = fADGeoToChannelAndSV.find(adGeoIdx);
    if (it == fADGeoToChannelAndSV.cend()) {
      throw cet::exception("Geometry")
        << "Given AuxDetSensitive channel, " << channel
        << ", cannot be found in vector associated to AuxDetGeo index: " << adGeoIdx
        << ". Vector has size " << it->second.size();
    }

    // get the vector of channels to AuxDetSensitiveGeo index
    if (channel < it->second.size()) return {adGeoIdx, it->second[channel].second};

    throw cet::exception("Geometry") << "Given AuxDetGeo with index " << adGeoIdx
                                     << " does not correspond to any vector of sensitive volumes";
  }

std::size_t geo::AuxDetReadoutGeom::DetNameToAuxDet ( std::string const &  detName ) const

private

Definition at line 167 of file AuxDetReadoutGeom.cxx.

References fADGeoToName.

Referenced by ChannelToSensitiveAuxDet().

  {
    // loop over the map of AuxDet names to Geo object numbers to determine which auxdet
    // we have.  If no name in the map matches the provided string, throw an exception;
    // the list of AuxDetGeo passed as argument is ignored!  Note that fADGeoToName must
    // have been updated by a derived class.
    auto it = std::find_if(fADGeoToName.begin(), fADGeoToName.end(), [&detName](auto const& pr) {
      return pr.second == detName;
    });
    if (it == fADGeoToName.end()) {
      throw cet::exception("Geometry") << "No AuxDetGeo matching name: " << detName;
    }
    return it->first;
  }

std::size_t geo::AuxDetReadoutGeom::NearestAuxDet	(	Point_t const &	point,
		std::vector< AuxDetGeo > const &	auxDets,
		double	tolerance = 0,
		bool	throwIfAbsent = true
	)		const

Definition at line 36 of file AuxDetReadoutGeom.cxx.

References tca::Length().

Referenced by NearestSensitiveAuxDet().

  {
    for (std::size_t a = 0; a < auxDets.size(); ++a) {
      auto const localPoint = auxDets[a].toLocalCoords(point);
      double const HalfCenterWidth = 0.5 * (auxDets[a].HalfWidth1() + auxDets[a].HalfWidth2());

      if (localPoint.Z() >= (-auxDets[a].Length() / 2 - tolerance) &&
          localPoint.Z() <= (auxDets[a].Length() / 2 + tolerance) &&
          localPoint.Y() >= (-auxDets[a].HalfHeight() - tolerance) &&
          localPoint.Y() <= (auxDets[a].HalfHeight() + tolerance) &&
          // if AuxDet a is a box, then HalfSmallWidth = HalfWidth
          localPoint.X() >= (-HalfCenterWidth +
                             localPoint.Z() * (HalfCenterWidth - auxDets[a].HalfWidth2()) /
                               (0.5 * auxDets[a].Length()) -
                             tolerance) &&
          localPoint.X() <= (HalfCenterWidth -
                             localPoint.Z() * (HalfCenterWidth - auxDets[a].HalfWidth2()) /
                               (0.5 * auxDets[a].Length()) +
                             tolerance))
        return a;
    } // for loop over AudDet a

    std::ostringstream msg;
    msg << "Can't find AuxDet for position " << point;

    if (throwIfAbsent) { throw cet::exception("AuxDetReadoutGeom") << msg.str(); }

    mf::LogDebug("AuxDetReadoutGeom") << msg.str();
    return invalid_index;
  }

std::size_t geo::AuxDetReadoutGeom::NearestSensitiveAuxDet	(	Point_t const &	point,
		std::vector< AuxDetGeo > const &	auxDets,
		double	tolerance = 0,
		bool	throwIfAbsent = true
	)		const

Definition at line 71 of file AuxDetReadoutGeom.cxx.

References geo::AuxDetSensitiveGeo::HalfHeight(), geo::AuxDetSensitiveGeo::HalfWidth1(), geo::AuxDetSensitiveGeo::HalfWidth2(), geo::AuxDetSensitiveGeo::Length(), NearestAuxDet(), geo::AuxDetGeo::NSensitiveVolume(), geo::AuxDetGeo::SensitiveVolume(), and geo::AuxDetSensitiveGeo::toLocalCoords().

  {
    std::size_t const auxDetIdx = NearestAuxDet(point, auxDets, tolerance, throwIfAbsent);
    if (auxDetIdx == invalid_index) { return invalid_index; }

    AuxDetGeo const& adg = auxDets[auxDetIdx];
    for (std::size_t a = 0; a < adg.NSensitiveVolume(); ++a) {
      AuxDetSensitiveGeo const& adsg = adg.SensitiveVolume(a);
      auto const localPoint = adsg.toLocalCoords(point);

      double const HalfCenterWidth = 0.5 * (adsg.HalfWidth1() + adsg.HalfWidth2());

      if (localPoint.Z() >= -(adsg.Length() / 2 + tolerance) &&
          localPoint.Z() <= (adsg.Length() / 2 + tolerance) &&
          localPoint.Y() >= -adsg.HalfHeight() - tolerance &&
          localPoint.Y() <= adsg.HalfHeight() + tolerance &&
          // if AuxDet a is a box, then HalfSmallWidth = HalfWidth
          localPoint.X() >=
            -HalfCenterWidth +
              localPoint.Z() * (HalfCenterWidth - adsg.HalfWidth2()) / (0.5 * adsg.Length()) -
              tolerance &&
          localPoint.X() <=
            HalfCenterWidth -
              localPoint.Z() * (HalfCenterWidth - adsg.HalfWidth2()) / (0.5 * adsg.Length()) +
              tolerance)
        return a;
    } // for loop over AuxDetSensitive a

    std::ostringstream msg;
    msg << "Can't find AuxDetSensitive for position " << point;

    if (throwIfAbsent) { throw cet::exception("AuxDetReadoutGeom") << msg.str(); }

    mf::LogDebug("AuxDetReadoutGeom") << msg.str();
    return invalid_index;
  }

Member Data Documentation

std::map<std::size_t, std::vector<chanAndSV> > geo::AuxDetReadoutGeom::fADGeoToChannelAndSV

private

map the AuxDetGeo index to a vector of pairs corresponding to the channel and AuxDetSensitiveGeo index

Definition at line 69 of file AuxDetReadoutGeom.h.

Referenced by AuxDetChannelToPosition(), AuxDetReadoutGeom(), and ChannelToSensitiveAuxDet().

std::map<std::size_t, std::string> geo::AuxDetReadoutGeom::fADGeoToName

private

map the AuxDetGeo index to the name

Definition at line 66 of file AuxDetReadoutGeom.h.

Referenced by DetNameToAuxDet().

std::map<std::string, std::size_t> geo::AuxDetReadoutGeom::fNameToADGeo

private

map the names to the AuxDetGeo index

Definition at line 67 of file AuxDetReadoutGeom.h.

Referenced by AuxDetChannelToPosition(), and AuxDetReadoutGeom().

---

The documentation for this class was generated from the following files:

• larcorealg/v10_00_02/source/larcorealg/Geometry/AuxDetReadoutGeom.h
• larcorealg/v10_00_02/source/larcorealg/Geometry/AuxDetReadoutGeom.cxx