geo::AuxDetGeo Class Reference

#include <AuxDetGeo.h>

Classes
struct	AuxDetGeoCoordinatesTag

Public Types
using	AuxDetSensitiveList_t = std::vector< geo::AuxDetSensitiveGeo >
	Type of list of sensitive volumes. More...
Types for geometry-local reference vectors.
These types represents points and displacement vectors in the reference frame defined in the auxiliary detector geometry box from the GDML geometry description. No alias is explicitly defined for the LArSoft global vector types, geo::Point_t and geo::Vector_t. Remember the LocalPoint_t and LocalVector_t vectors from different instances of geo::AuxDetGeo have the same type but are not compatible.
using	LocalPoint_t = geo::Point3DBase_t< AuxDetGeoCoordinatesTag >
	Type of points in the local GDML auxiliary detector frame. More...
using	LocalVector_t = geo::Vector3DBase_t< AuxDetGeoCoordinatesTag >
	Type of displacement vectors in the local GDML auxiliary detector frame. More...

Public Member Functions
	AuxDetGeo (TGeoNode const &node, geo::TransformationMatrix &&trans, AuxDetSensitiveList_t &&sensitive)
void	GetCenter (double *xyz, double localz=0.0) const
	Return the center position of an AuxDet. More...
geo::Point_t	GetCenter (double localz=0.0) const
	Returns the geometric center of the sensitive volume. More...
geo::Vector_t	GetNormalVector () const
	Returns the unit normal vector to the detector. More...
void	GetNormalVector (double *xyzDir) const
	Fills the unit normal vector to the detector. More...
std::string	Name () const
AuxDetSensitiveGeo const &	SensitiveVolume (size_t sv) const
size_t	NSensitiveVolume () const
void	SortSubVolumes (geo::GeoObjectSorter const &sorter)
template<typename Stream >
void	PrintAuxDetInfo (Stream &&out, std::string indent="", unsigned int verbosity=1) const
	Prints information about this auxiliary detector. More...
std::string	AuxDetInfo (std::string indent="", unsigned int verbosity=1) const
	Returns a string with auxiliary detector information. More...
Box geometry
double	Length () const
double	HalfWidth1 () const
double	HalfWidth2 () const
double	HalfHeight () const
const TGeoVolume *	TotalVolume () const
geo::Length_t	DistanceToPoint (geo::Point_t const &point) const
	Returns the distance of point from the center of the detector. More...
geo::Length_t	DistanceToPoint (double const *point) const
Coordinate transformation
void	LocalToWorld (const double *auxdet, double *world) const
	Transform point from local auxiliary detector frame to world frame. More...
geo::Point_t	toWorldCoords (LocalPoint_t const &local) const
	Transform point from local auxiliary detector frame to world frame. More...
void	LocalToWorldVect (const double *auxdet, double *world) const
	Transform direction vector from local to world. More...
geo::Vector_t	toWorldCoords (LocalVector_t const &local) const
	Transform direction vector from local to world. More...
void	WorldToLocal (const double *world, double *auxdet) const
	Transform point from world frame to local auxiliary detector frame. More...
LocalPoint_t	toLocalCoords (geo::Point_t const &world) const
	Transform point from world frame to local auxiliary detector frame. More...
void	WorldToLocalVect (const double *world, double *auxdet) const
	Transform direction vector from world to local. More...
LocalVector_t	toLocalCoords (geo::Vector_t const &world) const
	Transform direction vector from world to local. More...
Access to the sensitive volumes in the detector
std::size_t	FindSensitiveVolume (geo::Point_t const &point) const
std::size_t	FindSensitiveVolume (double const worldLoc[3]) const
AuxDetSensitiveGeo const &	PositionToSensitiveVolume (geo::Point_t const &point, size_t &sv) const
AuxDetSensitiveGeo const &	PositionToSensitiveVolume (double const worldLoc[3], size_t &sv) const

Static Public Attributes
static constexpr unsigned int	MaxVerbosity = 4
	Maximum verbosity supported by PrintAuxDetInfo(). More...

Private Types
using	LocalTransformation_t = geo::LocalTransformationGeo< geo::TransformationMatrix, LocalPoint_t, LocalVector_t >

Private Member Functions
void	InitShapeSize ()
	Extracts the size of the detector from the geometry information. More...

Private Attributes
const TGeoVolume *	fTotalVolume
	Total volume of AuxDet, called vol*. More...
LocalTransformation_t	fTrans
	Auxiliary detector-to-world transformation. More...
double	fLength
	length of volume, along z direction in local More...
double	fHalfWidth1
	1st half width of volume, at -z/2 in local coordinates More...
double	fHalfWidth2
	2nd half width (width1==width2 for boxes), at +z/2 More...
double	fHalfHeight
	half height of volume More...
std::vector< AuxDetSensitiveGeo >	fSensitive
	sensitive volumes in the detector More...

Detailed Description

Definition at line 38 of file AuxDetGeo.h.

Member Typedef Documentation

using geo::AuxDetGeo::AuxDetSensitiveList_t = std::vector<geo::AuxDetSensitiveGeo>

Type of list of sensitive volumes.

Definition at line 42 of file AuxDetGeo.h.

using geo::AuxDetGeo::LocalPoint_t = geo::Point3DBase_t<AuxDetGeoCoordinatesTag>

Type of points in the local GDML auxiliary detector frame.

Definition at line 64 of file AuxDetGeo.h.

using geo::AuxDetGeo::LocalTransformation_t = geo::LocalTransformationGeo <geo::TransformationMatrix, LocalPoint_t, LocalVector_t>

private

Definition at line 221 of file AuxDetGeo.h.

using geo::AuxDetGeo::LocalVector_t = geo::Vector3DBase_t<AuxDetGeoCoordinatesTag>

Type of displacement vectors in the local GDML auxiliary detector frame.

Definition at line 67 of file AuxDetGeo.h.

Constructor & Destructor Documentation

geo::AuxDetGeo::AuxDetGeo	(	TGeoNode const &	node,
		geo::TransformationMatrix &&	trans,
		AuxDetSensitiveList_t &&	sensitive
	)

Definition at line 32 of file AuxDetGeo.cxx.

    : fTotalVolume(node.GetVolume())
    , fTrans(std::move(trans))
    , fSensitive(std::move(sensitive))
  {

    if(!fTotalVolume)
      throw cet::exception("AuxDetGeo") << "cannot find AuxDet volume\n";

    MF_LOG_DEBUG("Geometry") << "detector total  volume is " << fTotalVolume->GetName();

    // if there are no sensitive volumes then this aux det
    // could be from an older gdml file than the introduction of AuxDetSensitiveGeo
    // in that case assume the full AuxDetGeo is sensitive and copy its information
    // into a single AuxDetSensitive
    if (fSensitive.empty())
      fSensitive.emplace_back(node, geo::TransformationMatrix(fTrans.Matrix()));

    InitShapeSize();

  }

Member Function Documentation

std::string geo::AuxDetGeo::AuxDetInfo	(	std::string	indent = "",
		unsigned int	verbosity = 1
	)		const

Returns a string with auxiliary detector information.

See Also

PrintAuxDetInfo()

The arguments and provided information are the same as in PrintAuxDetInfo().

Definition at line 150 of file AuxDetGeo.cxx.

  {
    std::ostringstream sstr;
    PrintAuxDetInfo(sstr, indent, verbosity);
    return sstr.str();
  } // AuxDetGeo::AuxDetInfo()

geo::Length_t geo::AuxDetGeo::DistanceToPoint ( geo::Point_t const &  point ) const

inline

Returns the distance of point from the center of the detector.

Definition at line 111 of file AuxDetGeo.h.

      { return (point - GetCenter()).R(); }

geo::Length_t geo::AuxDetGeo::DistanceToPoint

(

double const *

point

)

const

Definition at line 87 of file AuxDetGeo.cxx.

    { return DistanceToPoint(geo::vect::makePointFromCoords(point)); }

std::size_t geo::AuxDetGeo::FindSensitiveVolume

(

geo::Point_t const &

point

)

const

Definition at line 91 of file AuxDetGeo.cxx.

  {
    for(std::size_t a = 0; a < fSensitive.size(); ++a) {
      auto const& sensVol = SensitiveVolume(a);

      auto const local = sensVol.toLocalCoords(point);

      double const HalfCenterWidth = sensVol.HalfCenterWidth();

      double const deltaWidth
        = local.Z()*(HalfCenterWidth-sensVol.HalfWidth2())/sensVol.HalfLength();

      if (local.Z() >= - sensVol.HalfLength()       &&
          local.Z() <=   sensVol.HalfLength()       &&
          local.Y() >= - sensVol.HalfHeight()     &&
          local.Y() <=   sensVol.HalfHeight()     &&
          // if SensitiveVolume a is a box, then HalfSmallWidth = HalfWidth
          local.X() >= - HalfCenterWidth + deltaWidth &&
          local.X() <=   HalfCenterWidth - deltaWidth
        )  return a;

    }// for loop over AuxDetSensitive a

    throw cet::exception("AuxDetGeo")
      << "Can't find AuxDetSensitive for position " << point << "\n";
    // the following is not very useful right now...
    return std::numeric_limits<std::size_t>::max();
  } // AuxDetGeo::FindSensitiveVolume(geo::Point_t)

std::size_t geo::AuxDetGeo::FindSensitiveVolume

(

double const

worldLoc[3]

)

const

Deprecated:

Use the version with geo::Point_t argument instead

Definition at line 121 of file AuxDetGeo.cxx.

    { return FindSensitiveVolume(geo::vect::makePointFromCoords(worldPos)); }

void geo::AuxDetGeo::GetCenter	(	double *	xyz,
		double	localz = 0.0
	)		const

Return the center position of an AuxDet.

Parameters

xyz	_(output)_ the returned location: { x, y, z } [cm]
localz	(default: 0) distance along the length of the volume (z) [cm]

Definition at line 62 of file AuxDetGeo.cxx.

                                                            {
    auto const& center = GetCenter(localz);
    xyz[0] = center.X();
    xyz[1] = center.Y();
    xyz[2] = center.Z();
  } // AuxDetGeo::GetCenter(double*)

geo::Point_t geo::AuxDetGeo::GetCenter

(

double

localz = 0.0

)

const

Returns the geometric center of the sensitive volume.

Parameters

localz

(default: 0) distance from the center along the length of the volume (z) [cm]

Returns

the geometric center of the sensitive volume [cm]

Definition at line 58 of file AuxDetGeo.cxx.

    { return toWorldCoords(LocalPoint_t{ 0.0, 0.0, localz }); }

geo::Vector_t geo::AuxDetGeo::GetNormalVector

(

)

const

Returns the unit normal vector to the detector.

Definition at line 72 of file AuxDetGeo.cxx.

    { return toWorldCoords(geo::Zaxis<LocalVector_t>()); }

void geo::AuxDetGeo::GetNormalVector

(

double *

xyzDir

)

const

Fills the unit normal vector to the detector.

Definition at line 78 of file AuxDetGeo.cxx.

                                                      {
    auto const& norm = GetNormalVector();
    xyzDir[0] = norm.X();
    xyzDir[1] = norm.Y();
    xyzDir[2] = norm.Z();
  } // AuxDetGeo::GetNormalVector(double*)

double geo::AuxDetGeo::HalfHeight ( ) const

inline

Definition at line 105 of file AuxDetGeo.h.

{ return fHalfHeight;  }

double geo::AuxDetGeo::HalfWidth1 ( ) const

inline

Definition at line 103 of file AuxDetGeo.h.

{ return fHalfWidth1;  }

double geo::AuxDetGeo::HalfWidth2 ( ) const

inline

Definition at line 104 of file AuxDetGeo.h.

{ return fHalfWidth2;  }

void geo::AuxDetGeo::InitShapeSize ( )

private

Extracts the size of the detector from the geometry information.

Definition at line 159 of file AuxDetGeo.cxx.

                                {
    // set the ends depending on whether the shape is a box or trapezoid
    std::string volName(fTotalVolume->GetName());
    if( volName.find("Trap") != std::string::npos ) {

      //       Small Width
      //          ____          Height is the thickness
      //         /    \     T     of the trapezoid
      //        /      \    |
      //       /        \   | Length
      //      /__________\  _
      //         Width
      fHalfHeight      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDy1(); // same as Dy2()
      fLength          = 2.0*((TGeoTrd2*)fTotalVolume->GetShape())->GetDz();
      fHalfWidth1      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDx1(); // at -Dz
      fHalfWidth2      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDx2(); // at +Dz
    }
    else {
      fHalfWidth1      =     ((TGeoBBox*)fTotalVolume->GetShape())->GetDX();
      fHalfHeight      =     ((TGeoBBox*)fTotalVolume->GetShape())->GetDY();
      fLength          = 2.0*((TGeoBBox*)fTotalVolume->GetShape())->GetDZ();
      fHalfWidth2      = fHalfWidth1;
    }
  } // AuxDetGeo::InitShapeSize()

double geo::AuxDetGeo::Length ( ) const

inline

Definition at line 102 of file AuxDetGeo.h.

{ return fLength;      }

void geo::AuxDetGeo::LocalToWorld ( const double *  auxdet, double *  world  ) const

inline

Transform point from local auxiliary detector frame to world frame.

Definition at line 120 of file AuxDetGeo.h.

      { fTrans.LocalToWorld(auxdet, world); }

void geo::AuxDetGeo::LocalToWorldVect ( const double *  auxdet, double *  world  ) const

inline

Transform direction vector from local to world.

Definition at line 128 of file AuxDetGeo.h.

      { fTrans.LocalToWorldVect(auxdet, world); }

std::string geo::AuxDetGeo::Name ( ) const

inline

Definition at line 154 of file AuxDetGeo.h.

{ return fTotalVolume->GetName(); }

size_t geo::AuxDetGeo::NSensitiveVolume ( ) const

inline

Definition at line 172 of file AuxDetGeo.h.

{ return fSensitive.size(); }

AuxDetSensitiveGeo const & geo::AuxDetGeo::PositionToSensitiveVolume	(	geo::Point_t const &	point,
		size_t &	sv
	)		const

Definition at line 127 of file AuxDetGeo.cxx.

  {
    sv = FindSensitiveVolume(point);
    if (sv == std::numeric_limits<std::size_t>::max()) {
      throw cet::exception("AuxDetGeo")
        << "Can't find AuxDetSensitiveGeo for position " << point << "\n";
    }
    return SensitiveVolume(sv);
  }

AuxDetSensitiveGeo const & geo::AuxDetGeo::PositionToSensitiveVolume	(	double const	worldLoc[3],
		size_t &	sv
	)		const

Deprecated:

Use the version with geo::Point_t argument instead

Definition at line 139 of file AuxDetGeo.cxx.

    { return PositionToSensitiveVolume(geo::vect::makePointFromCoords(worldLoc), sv); }

template<typename Stream >

void geo::AuxDetGeo::PrintAuxDetInfo	(	Stream &&	out,
		std::string	indent = "",
		unsigned int	verbosity = 1
	)		const

Prints information about this auxiliary detector.

Template Parameters

Stream

type of output stream to use

Parameters

out	stream to send the information to
indent	prepend each line with this string
verbosity	amount of information printed

Note that the first line out the output is not indented.

Verbosity levels

• 0: only detector name
• 1 _(default)_: also center
• 2: also size
• 3: also number of sensitive detectors
• 4: also normal direction

The constant MaxVerbosity is set to the highest supported verbosity level.

Definition at line 245 of file AuxDetGeo.h.

        {

  //----------------------------------------------------------------------------
  out << "\"" << Name() << "\"";

  if (verbosity-- <= 0) return; // 0

  //----------------------------------------------------------------------------
  out << " centered at " << GetCenter() << " cm";

  if (verbosity-- <= 0) return; // 1

  //----------------------------------------------------------------------------
  lar::util::RealComparisons<double> coordIs(1e-4);
  out << ", size ( " << (2.0 * HalfWidth1());
  if (coordIs.nonEqual(HalfWidth1(), HalfWidth2()))
    out << "/" << (2.0 * HalfWidth2());
  out << " x " << (2.0 * HalfHeight()) << " x " << Length() << " ) cm";

  if (verbosity-- <= 0) return; // 2

  //----------------------------------------------------------------------------
  out << "\n" << indent
    << "with ";
  switch (NSensitiveVolume()) {
    case 0:  out << "no sensitive volume"; break;
    case 1:  out << "1 sensitive volume"; break;
    default: out << NSensitiveVolume() << " sensitive volumes"; break;
  } // switch

  if (verbosity-- <= 0) return; // 3

  //----------------------------------------------------------------------------
  out << ", normal facing " << GetNormalVector();

//  if (verbosity-- <= 0) return; // 4

  //----------------------------------------------------------------------------

} // geo::AuxDetGeo::PrintAuxDetInfo()

AuxDetSensitiveGeo const& geo::AuxDetGeo::SensitiveVolume ( size_t  sv ) const

inline

Definition at line 171 of file AuxDetGeo.h.

{ return fSensitive[sv];   }

void geo::AuxDetGeo::SortSubVolumes

(

geo::GeoObjectSorter const &

sorter

)

Definition at line 143 of file AuxDetGeo.cxx.

  {
    sorter.SortAuxDetSensitive(fSensitive);
  }

LocalPoint_t geo::AuxDetGeo::toLocalCoords ( geo::Point_t const &  world ) const

inline

Transform point from world frame to local auxiliary detector frame.

Definition at line 140 of file AuxDetGeo.h.

      { return fTrans.toLocalCoords(world); }

LocalVector_t geo::AuxDetGeo::toLocalCoords ( geo::Vector_t const &  world ) const

inline

Transform direction vector from world to local.

Definition at line 148 of file AuxDetGeo.h.

      { return fTrans.toLocalCoords(world); }

const TGeoVolume* geo::AuxDetGeo::TotalVolume ( ) const

inline

Definition at line 106 of file AuxDetGeo.h.

{ return fTotalVolume; }

geo::Point_t geo::AuxDetGeo::toWorldCoords ( LocalPoint_t const &  local ) const

inline

Transform point from local auxiliary detector frame to world frame.

Definition at line 124 of file AuxDetGeo.h.

      { return fTrans.toWorldCoords(local); }

geo::Vector_t geo::AuxDetGeo::toWorldCoords ( LocalVector_t const &  local ) const

inline

Transform direction vector from local to world.

Definition at line 132 of file AuxDetGeo.h.

      { return fTrans.toWorldCoords(local); }

void geo::AuxDetGeo::WorldToLocal ( const double *  world, double *  auxdet  ) const

inline

Transform point from world frame to local auxiliary detector frame.

Definition at line 136 of file AuxDetGeo.h.

      { fTrans.WorldToLocal(world, auxdet); }

void geo::AuxDetGeo::WorldToLocalVect ( const double *  world, double *  auxdet  ) const

inline

Transform direction vector from world to local.

Definition at line 144 of file AuxDetGeo.h.

      { fTrans.WorldToLocalVect(world, auxdet); }

Member Data Documentation

double geo::AuxDetGeo::fHalfHeight

private

half height of volume

Definition at line 228 of file AuxDetGeo.h.

double geo::AuxDetGeo::fHalfWidth1

private

1st half width of volume, at -z/2 in local coordinates

Definition at line 226 of file AuxDetGeo.h.

double geo::AuxDetGeo::fHalfWidth2

private

2nd half width (width1==width2 for boxes), at +z/2

Definition at line 227 of file AuxDetGeo.h.

double geo::AuxDetGeo::fLength

private

length of volume, along z direction in local

Definition at line 225 of file AuxDetGeo.h.

std::vector<AuxDetSensitiveGeo> geo::AuxDetGeo::fSensitive

private

sensitive volumes in the detector

Definition at line 229 of file AuxDetGeo.h.

const TGeoVolume* geo::AuxDetGeo::fTotalVolume

private

Total volume of AuxDet, called vol*.

Definition at line 223 of file AuxDetGeo.h.

LocalTransformation_t geo::AuxDetGeo::fTrans

private

Auxiliary detector-to-world transformation.

Definition at line 224 of file AuxDetGeo.h.

constexpr unsigned int geo::AuxDetGeo::MaxVerbosity = 4

static

Maximum verbosity supported by PrintAuxDetInfo().

Definition at line 215 of file AuxDetGeo.h.

---

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

• AuxDetGeo.h
• AuxDetGeo.cxx