LArSoft: geo::CryostatGeo Class Reference

static bool CoordinateContained(double c, double min, double max, double wiggle=1.)

Returns whether the specified coordinate is in a range.

Definition: BoxBoundedGeo.h:333

bool geo::BoxBoundedGeo::ContainsY	(	double	y,
		double const	wiggle = `1`
	)		const

inlineinherited

Returns whether this TPC contains the specified world y coordinate.

Parameters

y	the absolute ("world") coordinate y
wiggle	expansion factor for the range (see ContainsPosition())

Returns: whether the specified coordinate is in this TPC

See also: ContainsPosition()

Definition at line 171 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::MaxY(), and geo::BoxBoundedGeo::MinY().

Referenced by geo::BoxBoundedGeo::ContainsYZ().

172 { return CoordinateContained(y, MinY(), MaxY(), wiggle); }

Float_t y

Definition: compare.C:6

geo::BoxBoundedGeo::MaxY

double MaxY() const

Returns the world y coordinate of the end of the box.

Definition: BoxBoundedGeo.h:108

static bool CoordinateContained(double c, double min, double max, double wiggle=1.)

Returns whether the specified coordinate is in a range.

Definition: BoxBoundedGeo.h:333

geo::BoxBoundedGeo::MinY

double MinY() const

Returns the world y coordinate of the start of the box.

Definition: BoxBoundedGeo.h:105

bool geo::BoxBoundedGeo::ContainsYZ	(	double	y,
		double	z,
		double const	wiggle = `1`
	)		const

inlineinherited

Returns if TPC contains the specified world y and z coordinates.

Parameters

y	the absolute ("world") coordinate y
z	the absolute ("world") coordinate z
wiggle	expansion factor for the range (see ContainsPosition())

Returns: whether the specified coordinate is in this TPC

See also: ContainsPosition()

Definition at line 192 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::ContainsPosition(), geo::BoxBoundedGeo::ContainsY(), and geo::BoxBoundedGeo::ContainsZ().

Referenced by geo::BoxBoundedGeo::ContainsPosition().

193 { return ContainsY(y, wiggle) && ContainsZ(z, wiggle); }

Float_t y

Definition: compare.C:6

geo::BoxBoundedGeo::ContainsY

Double_t z

Definition: plot.C:279

bool ContainsY(double y, double const wiggle=1) const

Returns whether this TPC contains the specified world y coordinate.

Definition: BoxBoundedGeo.h:171

geo::BoxBoundedGeo::ContainsZ

bool ContainsZ(double z, double const wiggle=1) const

Returns whether this TPC contains the specified world z coordinate.

Definition: BoxBoundedGeo.h:181

bool geo::BoxBoundedGeo::ContainsZ	(	double	z,
		double const	wiggle = `1`
	)		const

inlineinherited

Returns whether this TPC contains the specified world z coordinate.

Parameters

z	the absolute ("world") coordinate z
wiggle	expansion factor for the range (see ContainsPosition())

Returns: whether the specified coordinate is in this TPC

See also: ContainsPosition()

Definition at line 181 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::MaxZ(), and geo::BoxBoundedGeo::MinZ().

Referenced by geo::BoxBoundedGeo::ContainsYZ().

182 { return CoordinateContained(z, MinZ(), MaxZ(), wiggle); }

Double_t z

Definition: plot.C:279

geo::BoxBoundedGeo::MinZ

double MinZ() const

Returns the world z coordinate of the start of the box.

Definition: BoxBoundedGeo.h:120

geo::BoxBoundedGeo::MaxZ

double MaxZ() const

Returns the world z coordinate of the end of the box.

Definition: BoxBoundedGeo.h:123

static bool CoordinateContained(double c, double min, double max, double wiggle=1.)

Returns whether the specified coordinate is in a range.

Definition: BoxBoundedGeo.h:333

static bool geo::BoxBoundedGeo::CoordinateContained	(	double	c,
		double	min,
		double	max,
		double	wiggle = `1.`
	)

inlinestaticinherited

Returns whether the specified coordinate is in a range.

Parameters

c	the coordinate
min	lower boundary of the range
max	upper boundary of the range
wiggle	expansion factor for the range

Returns: whether the specified coordinate is in a range

If the wiggle is larger than 1, the range is expanded by the wiggle factor. If the wiggle is less than 1, the range is shrinked.

Definition at line 333 of file BoxBoundedGeo.h.

Referenced by geo::BoxBoundedGeo::ContainsX(), geo::BoxBoundedGeo::ContainsY(), geo::BoxBoundedGeo::ContainsZ(), geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::InFiducialX(), geo::BoxBoundedGeo::InFiducialY(), and geo::BoxBoundedGeo::InFiducialZ().

       {
         return (c >= (min > 0? min / wiggle: min * wiggle))
           && (c <= (max < 0? max / wiggle: max * wiggle));
       } // CoordinateContained()

static bool geo::BoxBoundedGeo::CoordinateContained	(	double	c,
		double const *	range,
		double	wiggle = `1.`
	)

inlinestaticinherited

Returns whether the specified coordinate is in a range.

Parameters

c	the coordinate
range	pointer to [ lower boundary, upper boundary ] of the range
wiggle	expansion factor for the range

Returns: whether the specified coordinate is in a range

See also: CoordinateContained(double, double, double, double)

If the wiggle is larger than 1, the range is expanded by the wiggle factor. If the wiggle is less than 1, the range is shrinked.

Definition at line 351 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained().

352 { return CoordinateContained(c, range[0], range[1], wiggle); }

static bool CoordinateContained(double c, double min, double max, double wiggle=1.)

Returns whether the specified coordinate is in a range.

Definition: BoxBoundedGeo.h:333

void geo::BoxBoundedGeo::ExtendToInclude	(	Coord_t	x,
		Coord_t	y,
		Coord_t	z
	)

inlineinherited

Extends the current box to also include the specified point.

Parameters

x	x coordinate of the point to include
y	y coordinate of the point to include
z	z coordinate of the point to include

Definition at line 392 of file BoxBoundedGeo.h.

Referenced by geo::PlaneGeo::BoundingBox(), evd_tool::StandardDrawer::DetOutline3D(), and lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry().

393 { ExtendToInclude(geo::Point_t(x, y, z)); }

x

Float_t x

Definition: compare.C:6

Float_t y

Definition: compare.C:6

geo::BoxBoundedGeo::ExtendToInclude

Double_t z

Definition: plot.C:279

void ExtendToInclude(Coord_t x, Coord_t y, Coord_t z)

Extends the current box to also include the specified point.

Definition: BoxBoundedGeo.h:392

geo::Point_t

ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t

Type for representation of position in physical 3D space.

Definition: geo_vectors.h:187

void geo::BoxBoundedGeo::ExtendToInclude ( geo::Point_t const & point )

inlineinherited

Extends the current box to also include the specified point.

Parameters

point coordinates of the point to include

Definition at line 399 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max, geo::BoxBoundedGeo::c_min, geo::BoxBoundedGeo::set_max(), and geo::BoxBoundedGeo::set_min().

       { 
         set_min(c_min, point);
         set_max(c_max, point);
       }

void geo::BoxBoundedGeo::ExtendToInclude ( BoxBoundedGeo const & box )

inlineinherited

Extends the current box to also include the specified one.

Parameters

box	the box to include

It is assumed that the box has its boundaries properly sorted.

Definition at line 411 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max, geo::BoxBoundedGeo::c_min, geo::BoxBoundedGeo::GetIntersections(), geo::BoxBoundedGeo::Max(), geo::BoxBoundedGeo::Min(), geo::BoxBoundedGeo::set_max(), and geo::BoxBoundedGeo::set_min().

       {
         set_min(c_min, box.Min());
         set_max(c_max, box.Max());
       } // ExtendToInclude()

void geo::CryostatGeo::FindOpDet	(	std::vector< const TGeoNode * > &	path,
		unsigned int	depth
	)

private

Definition at line 188 of file CryostatGeo.cxx.

References MakeOpDet(), and OpDetGeoName().

Referenced by CryostatGeo(), and toLocalCoords().

   {
 
     const char* nm = path[depth]->GetName();
     if( (strncmp(nm, OpDetGeoName().c_str(), 6) == 0) ){
       this->MakeOpDet(path,depth);
       return;
     }
 
     //explore the next layer down
     unsigned int deeper = depth+1;
     if(deeper >= path.size()){
       throw cet::exception("BadTGeoNode") << "exceeded maximum TGeoNode depth\n";
     }
 
     const TGeoVolume *v = path[depth]->GetVolume();
     int nd = v->GetNdaughters();
     for(int i = 0; i < nd; ++i){
       path[deeper] = v->GetNode(i);
       this->FindOpDet(path, deeper);
     }
   
   }

void geo::CryostatGeo::FindTPC	(	std::vector< const TGeoNode * > &	path,
		unsigned int	depth
	)

private

Definition at line 105 of file CryostatGeo.cxx.

References MakeTPC().

Referenced by CryostatGeo(), and toLocalCoords().

   {
 
     const char* nm = path[depth]->GetName();
     if( (strncmp(nm, "volTPC", 6) == 0) ){
       this->MakeTPC(path,depth);
       return;
     }
 
     //explore the next layer down
     unsigned int deeper = depth+1;
     if(deeper >= path.size()){
       throw cet::exception("BadTGeoNode") << "exceeded maximum TGeoNode depth\n";
     }
 
     const TGeoVolume *v = path[depth]->GetVolume();
     int nd = v->GetNdaughters();
     for(int i = 0; i < nd; ++i){
       path[deeper] = v->GetNode(i);
       this->FindTPC(path, deeper);
     }
   
   }

geo::TPCID::TPCID_t geo::CryostatGeo::FindTPCAtPosition	(	double const	worldLoc[3],
		double const	wiggle
	)		const

Returns the index of the TPC at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
wiggle	a small factor (like 1+epsilon) to avoid rounding errors

Returns: the TPC index, or geo::TPCID::InvalidID if no TPC is there

Deprecated:: Use PositionToTPCID() instead

Definition at line 236 of file CryostatGeo.cxx.

References geo::TPCID::InvalidID, geo::vect::makePointFromCoords(), PositionToTPCID(), and geo::TPCID::TPC.

Referenced by GetElementPtr(), evd::SimulationDrawer::MCTruth3D(), and OpDet().

   {
     geo::TPCID tpcid
       = PositionToTPCID(geo::vect::makePointFromCoords(worldLoc), wiggle);
     return tpcid? tpcid.TPC: geo::TPCID::InvalidID;
   } // CryostatGeo::FindTPCAtPosition()

geo::Point_t geo::CryostatGeo::GetCenter ( ) const

inline

Returns the geometrical center of the cryostat.

Definition at line 108 of file CryostatGeo.h.

References Boundaries(), and geo::BoxBoundedGeo::Center().

Referenced by PrintCryostatInfo().

109 { return Boundaries().Center(); }

geo::CryostatGeo::Boundaries

geo::BoxBoundedGeo const & Boundaries() const

Returns boundaries of the cryostat (in centimetres).

Definition: CryostatGeo.h:99

geo::BoxBoundedGeo::Center

geo::Point_t Center() const

Returns the center point of the box.

Definition: BoxBoundedGeo.h:141

unsigned int geo::CryostatGeo::GetClosestOpDet ( geo::Point_t const & point ) const

Returns the index of the optical detector in this cryostat closest to point.

Definition at line 339 of file CryostatGeo.cxx.

References geo::OpDetGeo::DistanceToPoint(), max, NOpDet(), and OpDet().

Referenced by GetClosestOpDet(), geo::GeometryCore::GetClosestOpDet(), GetClosestOpDetPtr(), and NOpDet().

                                                                          {
     unsigned int ClosestDet = std::numeric_limits<unsigned int>::max();
     double ClosestDist = std::numeric_limits<double>::max();
     
     for(unsigned int o = 0U; o < NOpDet(); ++o) {
       double const ThisDist = OpDet(o).DistanceToPoint(point); 
       if(ThisDist < ClosestDist) {
         ClosestDist = ThisDist;
         ClosestDet  = o;
       }
     } // for
     return ClosestDet;
   } // CryostatGeo::GetClosestOpDet(geo::Point_t)

unsigned int geo::CryostatGeo::GetClosestOpDet ( double const * point ) const

See also: GetClosestOpDet(geo::Point_t const&) const

Definition at line 354 of file CryostatGeo.cxx.

References GetClosestOpDet(), and geo::vect::makePointFromCoords().

355 { return GetClosestOpDet(geo::vect::makePointFromCoords(point)); }

geo::CryostatGeo::GetClosestOpDet

unsigned int GetClosestOpDet(geo::Point_t const &point) const

Definition: CryostatGeo.cxx:339

geo::vect::makePointFromCoords

GENVECTOR_CONSTEXPR::geo::Point_t makePointFromCoords(Coords &&coords)

Creates a geo::Point_t from its coordinates (see makeFromCoords()).

Definition: geo_vectors_utils.h:1307

geo::OpDetGeo const * geo::CryostatGeo::GetClosestOpDetPtr ( geo::Point_t const & point ) const

Returns the optical detector det in this cryostat nearest to point. If there are no optical detectors, nullptr is returned.

Definition at line 331 of file CryostatGeo.cxx.

References GetClosestOpDet(), max, and OpDet().

Referenced by Boundaries(), and NOpDet().

   {
     unsigned int iOpDet = GetClosestOpDet(point);
     return
       (iOpDet == std::numeric_limits<double>::max())? nullptr: &OpDet(iOpDet);
   }

const TPCGeo& geo::CryostatGeo::GetElement ( TPCID const & tpcid ) const

inline

Alias for TPC().

Definition at line 197 of file CryostatGeo.h.

References TPC().

198 { return TPC(tpcid); }

geo::CryostatGeo::TPC

const TPCGeo & TPC(unsigned int itpc) const

Return the itpc'th TPC in the cryostat.

Definition: CryostatGeo.cxx:174

TPCGeo const* geo::CryostatGeo::GetElementPtr ( TPCID const & tpcid ) const

inline

Alias for TPCPtr().

Definition at line 231 of file CryostatGeo.h.

References FindTPCAtPosition(), PositionToTPC(), PositionToTPCID(), and TPCPtr().

232 { return TPCPtr(tpcid); }

geo::CryostatGeo::TPCPtr

TPCGeo const * TPCPtr(unsigned int itpc) const

Returns the TPC number itpc from this cryostat.

Definition: CryostatGeo.h:217

std::vector< TVector3 > geo::BoxBoundedGeo::GetIntersections	(	TVector3 const &	TrajectoryStart,
		TVector3 const &	TrajectoryDirect
	)		const

inherited

Calculates the entry and exit points of a trajectory on the box surface.

Author: Christoph Rudolf von Rohr (crohr.nosp@m.@fna.nosp@m.l.gov)

Date: July 27th, 2015

Parameters

TrajectoryStart	position of the trajectory source
TrajectoryDirect	direction vector of the trajectory

This member is public since it just gives an output and does not change any member. The algorithm works only for a box shaped active volume with facing walls parallel to axis. If the return std::vector is empty the trajectory does not intersect with the box. Normally the return value should have one (if the trajectory originates in the box) or two (else) entries. If the return value has two entries the first represents the entry point and the second the exit point

Definition at line 118 of file BoxBoundedGeo.cxx.

Referenced by geo::BoxBoundedGeo::ExtendToInclude(), and geo::BoxBoundedGeo::GetIntersections().

   {
     std::vector<TVector3> intersections;
     for (auto const& point: GetIntersections(geo::Point_t(TrajectoryStart), geo::Vector_t(TrajectoryDirect)))
       intersections.emplace_back(point.X(), point.Y(), point.Z());
     return intersections;
   } // GetIntersections(TVector3)

std::vector< geo::Point_t > geo::BoxBoundedGeo::GetIntersections	(	geo::Point_t const &	TrajectoryStart,
		geo::Vector_t const &	TrajectoryDirect
	)		const

inherited

Calculates the entry and exit points of a trajectory on the box surface.

Author: Christoph Rudolf von Rohr (crohr.nosp@m.@fna.nosp@m.l.gov)

Date: July 27th, 2015

Parameters

TrajectoryStart	position of the trajectory source
TrajectoryDirect	direction vector of the trajectory

This member is public since it just gives an output and does not change any member. The algorithm works only for a box shaped active volume with facing walls parallel to axis. If the return std::vector is empty the trajectory does not intersect with the box. Normally the return value should have one (if the trajectory originates in the box) or two (else) entries. If the return value has two entries the first represents the entry point and the second the exit point

Definition at line 27 of file BoxBoundedGeo.cxx.

References geo::vect::bindCoord(), geo::BoxBoundedGeo::c_max, geo::BoxBoundedGeo::c_min, geo::vect::coord(), geo::BoxBoundedGeo::GetIntersections(), geo::BoxBoundedGeo::Max(), and geo::BoxBoundedGeo::Min().

           {
   
     std::vector<geo::Point_t> IntersectionPoints;
     std::vector<double> LineParameters;
     
     // Generate normal vectors and offsets for every plane of the box
     // All normal vectors are headed outwards
     // BUG the double brace syntax is required to work around clang bug 21629
     // (https://bugs.llvm.org/show_bug.cgi?id=21629)
     static std::array<geo::Vector_t, 6U> const NormalVectors = {{
       -geo::Xaxis<geo::Vector_t>(), geo::Xaxis<geo::Vector_t>(), // anode, cathode,
       -geo::Yaxis<geo::Vector_t>(), geo::Yaxis<geo::Vector_t>(), // bottom, top, 
       -geo::Zaxis<geo::Vector_t>(), geo::Zaxis<geo::Vector_t>()  // upstream, downstream
     }};
     // BUG the double brace syntax is required to work around clang bug 21629
     // (https://bugs.llvm.org/show_bug.cgi?id=21629)
     std::array<geo::Point_t, 6U> const NormalVectorOffsets = {{
       geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // Anode offset
       geo::Point_t{ Max().X(), Min().Y(), Min().Z() }, // Cathode
       geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // Bottom
       geo::Point_t{ Min().X(), Max().Y(), Min().Z() }, // Top
       geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // upstream
       geo::Point_t{ Min().X(), Min().Y(), Max().Z() }  // downstream
     }};
     
     // Loop over all surfaces of the box 
     for(unsigned int face_no = 0; face_no < NormalVectors.size(); face_no++)
     {
       // Check if trajectory and surface are not parallel
       if(NormalVectors[face_no].Dot(TrajectoryDirect))
       {
         // Calculate the line parameter for the intersection points
         LineParameters.push_back( NormalVectors[face_no].Dot(NormalVectorOffsets.at(face_no) - TrajectoryStart)
                                 / NormalVectors[face_no].Dot(TrajectoryDirect) );
       }
       else continue;
       
       // Calculate intersection point using the line parameter
       IntersectionPoints.push_back( TrajectoryStart + LineParameters.back()*TrajectoryDirect );
       
       // Coordinate which should be ignored when checking for limits added by Christoph Rudolf von Rohr 05/21/2016
       unsigned int NoCheckCoord;
       
       // Calculate NoCheckCoord out of the face_no
       if(face_no % 2)
       {
           // Convert odd face number to coordinate
           NoCheckCoord = (face_no - 1)/2;
       }
       else
       {
           // Convert even face number to coordinate
           NoCheckCoord = face_no/2;
       }
       
       // Loop over all three space coordinates
       unsigned int coord = 0;
       for(auto extractCoord: geo::vect::coordReaders<geo::Point_t>())
       {
         auto const lastPointCoord = geo::vect::bindCoord(IntersectionPoints.back(), extractCoord);
         auto const minCoord = geo::vect::bindCoord(c_min, extractCoord);
         auto const maxCoord = geo::vect::bindCoord(c_max, extractCoord);
         
         // Changed by Christoph Rudolf von Rohr 05/21/2016
         // Then check if point is not within the surface limits at this coordinate, without looking
         // at the plane normal vector coordinate. We can assume, that our algorithm already found this coordinate correctily.
         // In rare cases, looking for boundaries in this coordinate causes unexpected behavior due to floating point inaccuracies.
         if( coord++ != NoCheckCoord && ((lastPointCoord() > maxCoord()) || (lastPointCoord() < minCoord)) )
         {
           // if off limits, get rid of the useless data and break the coordinate loop
           LineParameters.pop_back();
           IntersectionPoints.pop_back();
           break;
         }
       }// coordinate loop
     }// Surcaces loop
     
     // sort points according to their parameter value (first is entry, second is exit)
     if(LineParameters.size() == 2 && LineParameters.front() > LineParameters.back())
     {
       std::swap(IntersectionPoints.front(),IntersectionPoints.back());
     }
     
     return IntersectionPoints;
   } // GetIntersections()

double geo::CryostatGeo::HalfHeight ( ) const

Half height of the cryostat [cm].

Definition at line 305 of file CryostatGeo.cxx.

References fVolume.

Referenced by geo::GeometryCore::CryostatHalfHeight(), Height(), and InitCryoBoundaries().

   {
     return static_cast<TGeoBBox const*>(fVolume->GetShape())->GetDY();
   }

double geo::CryostatGeo::HalfLength ( ) const

Half height of the cryostat [cm].

Definition at line 311 of file CryostatGeo.cxx.

References fVolume.

Referenced by InitCryoBoundaries(), and Length().

   { 
     return static_cast<TGeoBBox const*>(fVolume->GetShape())->GetDZ();
   }

double geo::BoxBoundedGeo::HalfSizeX ( ) const

inlineinherited

Returns the size from the center to the border on X dimension.

Definition at line 102 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::SizeX().

102 { return SizeX() / 2.0; }

geo::BoxBoundedGeo::SizeX

double SizeX() const

Returns the full size in the X dimension.

Definition: BoxBoundedGeo.h:99

double geo::BoxBoundedGeo::HalfSizeY ( ) const

inlineinherited

Returns the size from the center to the border on Y dimension.

Definition at line 117 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::SizeY().

117 { return SizeY() / 2.0; }

geo::BoxBoundedGeo::SizeY

double SizeY() const

Returns the full size in the Y dimension.

Definition: BoxBoundedGeo.h:114

double geo::BoxBoundedGeo::HalfSizeZ ( ) const

inlineinherited

Returns the size from the center to the border on Z dimension.

Definition at line 132 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::SizeZ().

132 { return SizeZ() / 2.0; }

geo::BoxBoundedGeo::SizeZ

double SizeZ() const

Returns the full size in the Z dimension.

Definition: BoxBoundedGeo.h:129

double geo::CryostatGeo::HalfWidth ( ) const

Half width of the cryostat [cm].

Definition at line 299 of file CryostatGeo.cxx.

References fVolume.

Referenced by geo::GeometryCore::CryostatHalfWidth(), InitCryoBoundaries(), and Width().

   {
     return static_cast<TGeoBBox const*>(fVolume->GetShape())->GetDX();
   }

bool geo::CryostatGeo::HasElement ( unsigned int itpc ) const

inline

Alias for HasTPC().

Definition at line 167 of file CryostatGeo.h.

References HasTPC().

167 { return HasTPC(itpc); }

bool HasTPC(unsigned int itpc) const

Returns whether a TPC with index itpc is present in this cryostat.

Definition: CryostatGeo.h:164

bool geo::CryostatGeo::HasElement ( geo::TPCID const & tpcid ) const

inline

Alias for HasTPC(geo::TPCID const&)

Definition at line 179 of file CryostatGeo.h.

References HasTPC(), and TPC().

179 { return HasTPC(tpcid); }

bool HasTPC(unsigned int itpc) const

Returns whether a TPC with index itpc is present in this cryostat.

Definition: CryostatGeo.h:164

bool geo::CryostatGeo::HasTPC ( unsigned int itpc ) const

inline

Returns whether a TPC with index itpc is present in this cryostat.

Parameters

itpc	index of TPC in this cryostat

Returns: whether the TPC with index itpc is present in this cryostat

Definition at line 164 of file CryostatGeo.h.

References NTPC().

Referenced by HasElement(), geo::GeometryCore::HasTPC(), and TPCPtr().

164 { return itpc < NTPC(); }

geo::CryostatGeo::NTPC

unsigned int NTPC() const

Number of TPCs in this cryostat.

Definition: CryostatGeo.h:155

bool geo::CryostatGeo::HasTPC ( geo::TPCID const & tpcid ) const

inline

Returns whether the TPC in tpcid is present in this cryostat.

Parameters

tpcid full TPC ID

Returns: whether the TPC in tpcid is present in this cryostat

The cryostat number in tpcid is ignored, as it is ignored whether tpcid is invalid.

Definition at line 177 of file CryostatGeo.h.

References HasTPC(), and geo::TPCID::TPC.

Referenced by HasTPC().

177 { return HasTPC(tpcid.TPC); }

geo::CryostatGeo::HalfHeight

bool HasTPC(unsigned int itpc) const

Returns whether a TPC with index itpc is present in this cryostat.

Definition: CryostatGeo.h:164

double geo::CryostatGeo::Height ( ) const

inline

Full height of the cryostat [cm].

Definition at line 89 of file CryostatGeo.h.

References HalfHeight().

Referenced by PrintCryostatInfo().

89 { return 2. * HalfHeight(); }

double HalfHeight() const

Half height of the cryostat [cm].

Definition: CryostatGeo.cxx:305

geo::CryostatID const& geo::CryostatGeo::ID ( ) const

inline

Returns the identifier of this cryostat.

Definition at line 116 of file CryostatGeo.h.

References fID, art::detail::indent(), and PrintCryostatInfo().

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), geo::GeometryCore::FindCryostatAtPosition(), geo::GeometryCore::FindTPCAtPosition(), geo::GeometryCore::GetClosestOpDet(), evd::SimulationDrawer::MCTruth3D(), geo::GeometryCore::PositionToCryostatID(), PositionToTPC(), and PrintCryostatInfo().

116 { return fID; }

geo::CryostatGeo::fID

geo::CryostatID fID

ID of this cryostat.

Definition: CryostatGeo.h:413

bool geo::BoxBoundedGeo::InFiducialX	(	double	x,
		double	neg_margin,
		double	pos_margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world x coordinate.

Parameters

x	the absolute ("world") coordinate x
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns: whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialX(x, 5., 8.);

on a TPC with x from 0 to 100 (cm, presumably) will return true for all x between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume. Note that x is by definition the drift direction, and a reconstructed x typically depends on an assumption respect to the event time.

Definition at line 243 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::MaxX(), and geo::BoxBoundedGeo::MinX().

Referenced by geo::BoxBoundedGeo::InFiducialX().

       {
         return CoordinateContained(x, MinX() + neg_margin, MaxX() - pos_margin);
       }

bool geo::BoxBoundedGeo::InFiducialX	(	double	x,
		double	margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world x coordinate.

See also: InFiducialX(double, double, double) const

Margins are symmetric.

Definition at line 253 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::InFiducialX().

254 { return InFiducialX(x, margin, margin); }

x

Float_t x

Definition: compare.C:6

geo::BoxBoundedGeo::InFiducialX

bool InFiducialX(double x, double neg_margin, double pos_margin) const

Returns whether TPC fiducial volume contains world x coordinate.

Definition: BoxBoundedGeo.h:243

bool geo::BoxBoundedGeo::InFiducialY	(	double	y,
		double	neg_margin,
		double	pos_margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world y coordinate.

Parameters

y	the absolute ("world") coordinate y
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns: whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialY(y, 5., 8.);

on a TPC with y from 0 to 100 (cm, presumably) will return true for all y between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume.

Definition at line 274 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::MaxY(), and geo::BoxBoundedGeo::MinY().

Referenced by geo::BoxBoundedGeo::InFiducialY().

       {
         return CoordinateContained(y, MinY() + neg_margin, MaxY() - pos_margin);
       }

bool geo::BoxBoundedGeo::InFiducialY	(	double	y,
		double	margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world y coordinate.

See also: InFiducialY(double, double, double) const

Margins are symmetric.

Definition at line 284 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::InFiducialY().

285 { return InFiducialY(y, margin, margin); }

geo::BoxBoundedGeo::InFiducialY

Float_t y

Definition: compare.C:6

bool InFiducialY(double y, double neg_margin, double pos_margin) const

Returns whether TPC fiducial volume contains world y coordinate.

Definition: BoxBoundedGeo.h:274

bool geo::BoxBoundedGeo::InFiducialZ	(	double	z,
		double	neg_margin,
		double	pos_margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world z coordinate.

Parameters

z	the absolute ("world") coordinate y
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns: whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialZ(z, 5., 8.);

on a TPC with z from 0 to 100 (cm, presumably) will return true for all z between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume.

Definition at line 305 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), geo::BoxBoundedGeo::MaxZ(), and geo::BoxBoundedGeo::MinZ().

Referenced by geo::BoxBoundedGeo::InFiducialZ().

       {
         return CoordinateContained(z, MinZ() + neg_margin, MaxZ() - pos_margin);
       }

bool geo::BoxBoundedGeo::InFiducialZ	(	double	z,
		double	margin
	)		const

inlineinherited

Returns whether TPC fiducial volume contains world z coordinate.

See also: InFiducialZ(double, double, double) const

Margins are symmetric.

Definition at line 315 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::CoordinateContained(), and geo::BoxBoundedGeo::InFiducialZ().

316 { return InFiducialZ(z, margin, margin); }

geo::BoxBoundedGeo::InFiducialZ

Double_t z

Definition: plot.C:279

bool InFiducialZ(double z, double neg_margin, double pos_margin) const

Returns whether TPC fiducial volume contains world z coordinate.

Definition: BoxBoundedGeo.h:305

void geo::CryostatGeo::InitCryoBoundaries ( )

private

Fill the boundary information of the cryostat.

Definition at line 358 of file CryostatGeo.cxx.

References HalfHeight(), HalfLength(), HalfWidth(), geo::BoxBoundedGeo::SetBoundaries(), toWorldCoords(), and Volume().

Referenced by CryostatGeo(), and toLocalCoords().

                                        {
     
     // check that this is indeed a box
     if (!dynamic_cast<TGeoBBox*>(Volume()->GetShape())) {
       // at initialisation time we don't know yet our real ID
       throw cet::exception("CryostatGeo") << "Cryostat is not a box! (it is a "
         << Volume()->GetShape()->IsA()->GetName() << ")\n";
     }
     
     // get the half width, height, etc of the cryostat
     const double halflength = HalfLength();
     const double halfwidth  = HalfWidth();
     const double halfheight = HalfHeight();
     
     SetBoundaries(
       toWorldCoords(LocalPoint_t{ -halfwidth, -halfheight, -halflength }),
       toWorldCoords(LocalPoint_t{ +halfwidth, +halfheight, +halflength })
       );
     
   } // CryostatGeo::InitCryoBoundaries()

double geo::CryostatGeo::Length ( ) const

inline

Length of the cryostat [cm].

Definition at line 91 of file CryostatGeo.h.

References HalfLength().

Referenced by geo::GeometryCore::CryostatLength(), and PrintCryostatInfo().

91 { return 2. * HalfLength(); }

geo::CryostatGeo::HalfLength

double HalfLength() const

Half height of the cryostat [cm].

Definition: CryostatGeo.cxx:311

void geo::CryostatGeo::LocalToWorld	(	const double *	cryo,
		double *	world
	)		const

inline

Transform point from local cryostat frame to world frame.

Definition at line 321 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformation< StoredMatrix >::LocalToWorld().

Referenced by geo::sortCryoStandard().

322 { fTrans.LocalToWorld(cryo, world); }

geo::LocalTransformation::LocalToWorld

void LocalToWorld(double const *local, double *world) const

Transforms a point from local frame to world frame.

geo::LocalTransformation::LocalToWorld

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

template<typename Point >

Point geo::CryostatGeo::LocalToWorld ( Point const & local ) const

inline

Transform point from local cryostat frame to world frame.

Deprecated:: This method breaks the distinction between local and global vectors, since input and output vectors share the same type; use the "official" vector types geo::Point_t, geo::CryostatGeo::LocalPoint_t, geo::Vector_t and geo::CryostatGeo::LocalVector_t, and then use the method geo::CryostatGeo::toWorldCoords() instead.

Definition at line 333 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformation< StoredMatrix >::LocalToWorld().

334 { return fTrans.LocalToWorld(local); }

void LocalToWorld(double const *local, double *world) const

Transforms a point from local frame to world frame.

geo::LocalTransformation::LocalToWorldVect

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

void geo::CryostatGeo::LocalToWorldVect	(	const double *	cryo,
		double *	world
	)		const

inline

Transform direction vector from local to world.

Definition at line 341 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformation< StoredMatrix >::LocalToWorldVect().

342 { fTrans.LocalToWorldVect(cryo, world); }

void LocalToWorldVect(double const *local, double *world) const

Transforms a vector from local frame to world frame.

geo::CryostatGeo::fVolume

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

void geo::CryostatGeo::MakeOpDet	(	std::vector< const TGeoNode * > &	path,
		int	depth
	)

private

Definition at line 214 of file CryostatGeo.cxx.

References fOpDets.

Referenced by FindOpDet(), and toLocalCoords().

   {
     fOpDets.emplace_back(path, depth);
   }

void geo::CryostatGeo::MakeTPC	(	std::vector< const TGeoNode * > &	path,
		int	depth
	)

private

Definition at line 131 of file CryostatGeo.cxx.

References fTPCs.

Referenced by FindTPC(), and toLocalCoords().

   {
     fTPCs.emplace_back(path, depth);
   }

double geo::CryostatGeo::Mass ( ) const

inline

Mass of the cryostat.

Definition at line 93 of file CryostatGeo.h.

References fVolume.

93 { return fVolume->Weight(); }

TGeoVolume * fVolume

Total volume of cryostat, called volCryostat in GDML file.

Definition: CryostatGeo.h:411

geo::Point_t geo::BoxBoundedGeo::Max ( ) const

inlineinherited

Returns the corner point with the largest coordinates.

Definition at line 138 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max.

Referenced by geo::BoxBoundedGeo::ExtendToInclude(), geo::BoxBoundedGeo::GetIntersections(), geo::GeometryCore::MaterialName(), PrintCryostatInfo(), and geo::TPCGeo::PrintTPCInfo().

138 { return c_max; }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterX(), geo::BoxBoundedGeo::ContainsX(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialX(), evd::SimulationDrawer::MCTruthOrtho(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), geo::BoxBoundedGeo::SizeX(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_max

maximum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:469

unsigned int geo::CryostatGeo::MaxPlanes ( ) const

Returns the largest number of planes among the TPCs in this cryostat.

Definition at line 279 of file CryostatGeo.cxx.

References fTPCs, geo::TPCGeo::Nplanes(), and TPC().

Referenced by ems::MultiEMShowers::convCluster(), ems::EMShower3D::LinkCandidates(), PositionToTPC(), and PrintCryostatInfo().

                                             {
     unsigned int maxPlanes = 0;
     for (geo::TPCGeo const& TPC: fTPCs) {
       unsigned int maxPlanesInTPC = TPC.Nplanes();
       if (maxPlanesInTPC > maxPlanes) maxPlanes = maxPlanesInTPC;
     } // for
     return maxPlanes;
   } // CryostatGeo::MaxPlanes()

unsigned int geo::CryostatGeo::MaxWires ( ) const

Returns the largest number of wires among the TPCs in this cryostat.

Definition at line 289 of file CryostatGeo.cxx.

References fTPCs, geo::TPCGeo::MaxWires(), and TPC().

Referenced by PositionToTPC(), and PrintCryostatInfo().

                                            {
     unsigned int maxWires = 0;
     for (geo::TPCGeo const& TPC: fTPCs) {
       unsigned int maxWiresInTPC = TPC.MaxWires();
       if (maxWiresInTPC > maxWires) maxWires = maxWiresInTPC;
     } // for
     return maxWires;
   } // CryostatGeo::MaxWires()

double geo::BoxBoundedGeo::MaxX ( ) const

inlineinherited

Returns the world x coordinate of the end of the box.

Definition at line 93 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max.

93 { return c_max.X(); }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterY(), geo::BoxBoundedGeo::ContainsY(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialY(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), tca::SaveCRInfo(), geo::BoxBoundedGeo::SizeY(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_max

maximum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:469

double geo::BoxBoundedGeo::MaxY ( ) const

inlineinherited

Returns the world y coordinate of the end of the box.

Definition at line 108 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max.

108 { return c_max.Y(); }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterZ(), geo::BoxBoundedGeo::ContainsZ(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialZ(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), tca::SaveCRInfo(), geo::BoxBoundedGeo::SizeZ(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_max

maximum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:469

double geo::BoxBoundedGeo::MaxZ ( ) const

inlineinherited

Returns the world z coordinate of the end of the box.

Definition at line 123 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max.

123 { return c_max.Z(); }

Coords_t c_max

maximum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:469

geo::Point_t geo::BoxBoundedGeo::Min ( ) const

inlineinherited

Returns the corner point with the smallest coordinates.

Definition at line 135 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_min.

Referenced by geo::BoxBoundedGeo::ExtendToInclude(), geo::BoxBoundedGeo::GetIntersections(), geo::GeometryCore::MaterialName(), PrintCryostatInfo(), and geo::TPCGeo::PrintTPCInfo().

135 { return c_min; }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterX(), geo::BoxBoundedGeo::ContainsX(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialX(), pma::Track3D::InitFromMiddle(), evd::SimulationDrawer::MCTruthOrtho(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), evgen::ActiveVolumeVertexSampler::sample_vertex_pos(), geo::BoxBoundedGeo::SizeX(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_min

minimum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:468

double geo::BoxBoundedGeo::MinX ( ) const

inlineinherited

Returns the world x coordinate of the start of the box.

Definition at line 90 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_min.

90 { return c_min.X(); }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterY(), geo::BoxBoundedGeo::ContainsY(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialY(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), tca::SaveCRInfo(), geo::BoxBoundedGeo::SizeY(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_min

minimum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:468

double geo::BoxBoundedGeo::MinY ( ) const

inlineinherited

Returns the world y coordinate of the start of the box.

Definition at line 105 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_min.

105 { return c_min.Y(); }

Referenced by Boundaries(), geo::BoxBoundedGeo::CenterZ(), geo::BoxBoundedGeo::ContainsZ(), evd_tool::ICARUSDrawer::DetOutline3D(), evd_tool::StandardDrawer::DrawBoxBoundedGeoOutline(), lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), mvapid::MVAAlg::GetDetectorEdges(), geo::BoxBoundedGeo::InFiducialZ(), pma::Node3D::Node3D(), evgen::NucleonDecay::produce(), evgen::NeutronOsc::produce(), tca::SaveCRInfo(), geo::BoxBoundedGeo::SizeZ(), TPCarea(), and geo::GeometryCore::WorldBox().

Coords_t c_min

minimum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:468

double geo::BoxBoundedGeo::MinZ ( ) const

inlineinherited

Returns the world z coordinate of the start of the box.

Definition at line 120 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_min.

120 { return c_min.Z(); }

Coords_t c_min

minimum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:468

unsigned int geo::CryostatGeo::NElements ( ) const

inline

Alias for NTPC().

Definition at line 157 of file CryostatGeo.h.

References fTPCs.

Referenced by geo::GeometryCore::NTPC().

157 { return fTPCs.size(); }

geo::CryostatGeo::fOpDets

TPCList_t fTPCs

List of tpcs in this cryostat.

Definition: CryostatGeo.h:409

unsigned int geo::CryostatGeo::NOpDet ( ) const

inline

Number of optical detectors in this TPC.

Definition at line 295 of file CryostatGeo.h.

References fOpDets, GetClosestOpDet(), GetClosestOpDetPtr(), and OpDet().

Referenced by GetClosestOpDet(), geo::GeometryCore::NOpDets(), geo::GeometryCore::OpDetFromCryo(), geo::GeometryCore::OpDetGeoFromOpDet(), and PrintCryostatInfo().

295 { return fOpDets.size(); }

std::vector< OpDetGeo > fOpDets

List of opdets in this cryostat.

Definition: CryostatGeo.h:410

unsigned int geo::CryostatGeo::NTPC ( ) const

inline

Number of TPCs in this cryostat.

Definition at line 155 of file CryostatGeo.h.

References fTPCs.

Referenced by util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadoutGeometry::Construct(), evd_tool::ICARUSDrawer::DetOutline3D(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), vertex::FeatureVertexFinder::Find3dVtxFrom2dClusterVtxCand(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), corner::CornerFinderAlg::get_feature_points_fast(), groupTPCsByDriftDir(), HasTPC(), geo::ChannelMapStandardAlg::Initialize(), trkf::SpacePointAlg::makeSpacePoints(), PrintCryostatInfo(), vertex::FeatureVertexFinder::produce(), larg4::LArG4::produce(), evd::SimulationDrawer::SimulationDrawer(), geo::GeometryCore::TotalNTPC(), trkf::SpacePointAlg::update(), and UpdateAfterSorting().

155 { return fTPCs.size(); }

geo::CryostatGeo::fOpDetGeoName

TPCList_t fTPCs

List of tpcs in this cryostat.

Definition: CryostatGeo.h:409

const OpDetGeo & geo::CryostatGeo::OpDet ( unsigned int iopdet ) const

Return the iopdet'th optical detector in the cryostat.

Definition at line 220 of file CryostatGeo.cxx.

References FindTPCAtPosition(), and fOpDets.

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::FillFlashProperties(), opdet::FlashHypothesisAnaAlg::FillOpDetPositions(), GetClosestOpDet(), GetClosestOpDetPtr(), NOpDet(), and geo::GeometryCore::OpDetGeoFromOpDet().

   {
     if(iopdet >= fOpDets.size()){
       throw cet::exception("OpDetOutOfRange") << "Request for non-existant OpDet " 
                                               << iopdet;
     }
 
     return fOpDets[iopdet];
   }

std::string geo::CryostatGeo::OpDetGeoName ( ) const

inline

Get name of opdet geometry element.

Definition at line 311 of file CryostatGeo.h.

References fOpDetGeoName.

Referenced by FindOpDet(), and geo::GeometryCore::OpDetGeoName().

311 { return fOpDetGeoName; }

std::string fOpDetGeoName

Name of opdet geometry elements in gdml.

Definition: CryostatGeo.h:412

TPCGeo const & geo::CryostatGeo::PositionToTPC	(	geo::Point_t const &	point,
		double	wiggle
	)		const

Returns the ID of the TPC at specified location.

Parameters

point	the location (world reference frame)
wiggle	a small factor (like 1+epsilon) to avoid rounding errors

Returns: the ID of the TPC at the specified point (invalid ID if none)

Definition at line 257 of file CryostatGeo.cxx.

References ID(), and PositionToTPCptr().

Referenced by GetElementPtr(), PositionToTPC(), and PositionToTPCID().

   {
     geo::TPCGeo const* tpc = PositionToTPCptr(point, wiggle);
     if (!tpc) {
       throw cet::exception("CryostatGeo")
         << "Can't find any TPC for position " << point << " within " << ID()
         << "\n";
     }
     return *tpc;
   }

TPCGeo const& geo::CryostatGeo::PositionToTPC	(	double const	worldLoc[3],
		double	wiggle
	)		const

inline

Returns the ID of the TPC at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
wiggle	a small factor (like 1+epsilon) to avoid rounding errors

Returns: the ID of the TPC at the specified point (invalid ID if none)

Definition at line 268 of file CryostatGeo.h.

References geo::vect::makePointFromCoords(), MaxPlanes(), MaxWires(), PositionToTPC(), and PositionToTPCptr().

269 { return PositionToTPC(geo::vect::makePointFromCoords(worldLoc), wiggle); }

geo::CryostatGeo::PositionToTPC

TPCGeo const & PositionToTPC(geo::Point_t const &point, double wiggle) const

Returns the ID of the TPC at specified location.

Definition: CryostatGeo.cxx:257

geo::vect::makePointFromCoords

GENVECTOR_CONSTEXPR::geo::Point_t makePointFromCoords(Coords &&coords)

Creates a geo::Point_t from its coordinates (see makeFromCoords()).

Definition: geo_vectors_utils.h:1307

geo::TPCID geo::CryostatGeo::PositionToTPCID	(	geo::Point_t const &	point,
		double	wiggle
	)		const

Returns the ID of the TPC at specified location.

Parameters

point	3D coordinates of the point (world reference frame)
wiggle	a small factor (like 1+epsilon) to avoid rounding errors

Returns: the ID of the TPC at the specified point (invalid ID if none)

Definition at line 247 of file CryostatGeo.cxx.

References geo::TPCGeo::ID(), PositionToTPC(), and PositionToTPCptr().

Referenced by FindTPCAtPosition(), geo::GeometryCore::FindTPCAtPosition(), and GetElementPtr().

   { 
     geo::TPCGeo const* tpc = PositionToTPCptr(point, wiggle);
     return tpc? tpc->ID(): geo::TPCID{};
   }

geo::TPCGeo const * geo::CryostatGeo::PositionToTPCptr	(	geo::Point_t const &	point,
		double	wiggle
	)		const

Returns a pointer to the TPC at specified location.

Parameters

point	position in space [cm]
wiggle	a small factor (like 1+ε) to avoid rounding errors

Returns: a pointer to the geo::TPCGeo at point (nullptr if none)

Definition at line 270 of file CryostatGeo.cxx.

References TPCs().

Referenced by PositionToTPC(), PositionToTPCID(), and geo::GeometryCore::PositionToTPCptr().

   {
     for (auto const& tpc: TPCs())
       if (tpc.ContainsPosition(point, wiggle)) return &tpc;
     return nullptr;
   } // CryostatGeo::PositionToTPCptr()

template<typename Stream >

void geo::CryostatGeo::PrintCryostatInfo	(	Stream &&	out,
		std::string	indent = `""`,
		unsigned int	verbosity = `1`
	)		const

Prints information about this cryostat.

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 cryostat ID
1 _(default)_: also center and size
2: also number of TPCs, optical detectors, and maximum wires per plane
and of planes for TPC
3: also information on bounding box

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

Definition at line 423 of file CryostatGeo.h.

References BoundingBox(), GetCenter(), Height(), ID(), Length(), geo::BoxBoundedGeo::Max(), MaxPlanes(), MaxWires(), geo::BoxBoundedGeo::Min(), NOpDet(), NTPC(), and Width().

Referenced by ID().

         {
   
   //----------------------------------------------------------------------------
   out << "Cryostat " << std::string(ID());
   
   if (verbosity-- <= 0) return; // 0
   
   //----------------------------------------------------------------------------
   out
     << " (" << Width() << " x " << Height() << " x " << Length() << ") cm^3 at "
       << GetCenter();
   
   if (verbosity-- <= 0) return; // 1
   
   //----------------------------------------------------------------------------
   
   out << "\n" << indent
     << "hosts " << NTPC() << " TPCs (largest number of planes: " << MaxPlanes()
       << ", of wires: " << MaxWires() << ") and "
       << NOpDet() << " optical detectors"
     ;
   
   if (verbosity-- <= 0) return; // 2
   
   //----------------------------------------------------------------------------
   // print also the containing box
   geo::BoxBoundedGeo const& box = BoundingBox();
   out << "\n" << indent
     << "bounding box: " << box.Min() << " -- " << box.Max();
   
 //  if (verbosity-- <= 0) return; // 3
   
   //----------------------------------------------------------------------------
 } // geo::CryostatGeo::PrintCryostatInfo()

static void geo::BoxBoundedGeo::set_max	(	Coord_t &	var,
		Coord_t	value
	)

inlinestaticinherited

Sets var to value if value is larger than the current var value.

Definition at line 446 of file BoxBoundedGeo.h.

References fhicl::detail::atom::value().

Referenced by geo::BoxBoundedGeo::ExtendToInclude().

447 { if (value > var) var = value; }

fhicl::detail::atom::value

std::string value(boost::any const &)

Definition: printing_helpers.cc:97

static void geo::BoxBoundedGeo::set_max	(	Coords_t &	var,
		geo::Point_t const &	value
	)

inlinestaticinherited

Sets each coordinate of var to the one in value if the latter is larger.

Definition at line 458 of file BoxBoundedGeo.h.

       {
         if (value.X() > var.X()) var.SetX(value.X());
         if (value.Y() > var.Y()) var.SetY(value.Y());
         if (value.Z() > var.Z()) var.SetZ(value.Z());
       }

static void geo::BoxBoundedGeo::set_min	(	Coord_t &	var,
		Coord_t	value
	)

inlinestaticinherited

Sets var to value if value is smaller than the current var value.

Definition at line 442 of file BoxBoundedGeo.h.

References fhicl::detail::atom::value().

Referenced by geo::BoxBoundedGeo::ExtendToInclude().

443 { if (value < var) var = value; }

fhicl::detail::atom::value

std::string value(boost::any const &)

Definition: printing_helpers.cc:97

static void geo::BoxBoundedGeo::set_min	(	Coords_t &	var,
		geo::Point_t const &	value
	)

inlinestaticinherited

Sets each coordinate of var to the one in value if the latter is smaller.

Definition at line 450 of file BoxBoundedGeo.h.

       {
         if (value.X() < var.X()) var.SetX(value.X());
         if (value.Y() < var.Y()) var.SetY(value.Y());
         if (value.Z() < var.Z()) var.SetZ(value.Z());
       }

void geo::BoxBoundedGeo::SetBoundaries	(	Coord_t	x_min,
		Coord_t	x_max,
		Coord_t	y_min,
		Coord_t	y_max,
		Coord_t	z_min,
		Coord_t	z_max
	)

inlineinherited

Sets the boundaries in world coordinates as specified.

Parameters

x_min	lower x coordinate
x_max	upper x coordinate
y_min	lower y coordinate
y_max	upper y coordinate
z_min	lower z coordinate
z_max	upper z coordinate

Definition at line 367 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max, geo::BoxBoundedGeo::c_min, and geo::BoxBoundedGeo::SortCoordinates().

Referenced by geo::PlaneGeo::BoundingBox(), InitCryoBoundaries(), and geo::TPCGeo::InitTPCBoundaries().

       {
         c_min.SetXYZ(x_min, y_min, z_min);
         c_max.SetXYZ(x_max, y_max, z_max);
         SortCoordinates();
       }

void geo::BoxBoundedGeo::SetBoundaries	(	Coords_t	lower,
		Coords_t	upper
	)

inlineinherited

Sets the boundaries in world coordinates as specified.

Parameters

lower	lower coordinates (x, y, z)
upper	upper coordinates (x, y, z)

Definition at line 383 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::c_max, geo::BoxBoundedGeo::c_min, and geo::BoxBoundedGeo::SortCoordinates().

384 { c_min = lower; c_max = upper; SortCoordinates(); }

Coords_t c_max

maximum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:469

geo::BoxBoundedGeo::SortCoordinates

Coords_t c_min

minimum coordinates (x, y, z)

Definition: BoxBoundedGeo.h:468

void SortCoordinates()

Makes sure each coordinate of the minimum point is smaller than maximum.

Definition: BoxBoundedGeo.cxx:127

double geo::BoxBoundedGeo::SizeX ( ) const

inlineinherited

Returns the full size in the X dimension.

Definition at line 99 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::MaxX(), and geo::BoxBoundedGeo::MinX().

Referenced by geo::BoxBoundedGeo::HalfSizeX().

99 { return MaxX() - MinX(); }

geo::BoxBoundedGeo::MinX

double MinX() const

Returns the world x coordinate of the start of the box.

Definition: BoxBoundedGeo.h:90

geo::BoxBoundedGeo::MaxX

double MaxX() const

Returns the world x coordinate of the end of the box.

Definition: BoxBoundedGeo.h:93

double geo::BoxBoundedGeo::SizeY ( ) const

inlineinherited

Returns the full size in the Y dimension.

Definition at line 114 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::MaxY(), and geo::BoxBoundedGeo::MinY().

Referenced by geo::BoxBoundedGeo::HalfSizeY().

114 { return MaxY() - MinY(); }

geo::BoxBoundedGeo::MaxY

double MaxY() const

Returns the world y coordinate of the end of the box.

Definition: BoxBoundedGeo.h:108

geo::BoxBoundedGeo::MinY

double MinY() const

Returns the world y coordinate of the start of the box.

Definition: BoxBoundedGeo.h:105

double geo::BoxBoundedGeo::SizeZ ( ) const

inlineinherited

Returns the full size in the Z dimension.

Definition at line 129 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::MaxZ(), and geo::BoxBoundedGeo::MinZ().

Referenced by geo::BoxBoundedGeo::HalfSizeZ().

129 { return MaxZ() - MinZ(); }

geo::BoxBoundedGeo::MinZ

double MinZ() const

Returns the world z coordinate of the start of the box.

Definition: BoxBoundedGeo.h:120

geo::BoxBoundedGeo::MaxZ

double MaxZ() const

Returns the world z coordinate of the end of the box.

Definition: BoxBoundedGeo.h:123

void geo::CryostatGeo::SortSubVolumes ( geo::GeoObjectSorter const & sorter )

Method to sort TPCGeo objects.

Definition at line 139 of file CryostatGeo.cxx.

References fTPCs, util::SortByPointers(), geo::TPCGeo::SortSubVolumes(), geo::GeoObjectSorter::SortTPCs(), and TPC().

Referenced by toLocalCoords().

   {
     //
     // TPCs
     //
     util::SortByPointers
       (fTPCs, [&sorter](auto& coll){ sorter.SortTPCs(coll); });
     
     for (geo::TPCGeo& TPC: fTPCs) { 
       TPC.SortSubVolumes(sorter);
     } // for TPCs
     
     //
     // optical detectors
     //
     
     // sorting of optical detectors happens elsewhere
     
   } // CryostatGeo::SortSubVolumes()

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

inline

Transform point from world frame to local cryostat frame.

Definition at line 365 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformationGeo< StoredMatrix, LocalPoint, LocalVector >::toLocalCoords().

366 { return fTrans.toLocalCoords(world); }

geo::LocalTransformationGeo::toLocalCoords

LocalPoint_t toLocalCoords(GlobalPoint_t const &world) const

Transforms a point from world frame to local frame.

Definition: LocalTransformationGeo.h:133

geo::LocalTransformationGeo::toLocalCoords

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

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

inline

Transform direction vector from world to local.

Definition at line 373 of file CryostatGeo.h.

References FindOpDet(), FindTPC(), fTrans, InitCryoBoundaries(), MakeOpDet(), MakeTPC(), SortSubVolumes(), geo::LocalTransformationGeo< StoredMatrix, LocalPoint, LocalVector >::toLocalCoords(), and UpdateAfterSorting().

374 { return fTrans.toLocalCoords(world); }

LocalPoint_t toLocalCoords(GlobalPoint_t const &world) const

Transforms a point from world frame to local frame.

Definition: LocalTransformationGeo.h:133

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

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

inline

Transform point from local cryostat frame to world frame.

Definition at line 337 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformationGeo< StoredMatrix, LocalPoint, LocalVector >::toWorldCoords().

Referenced by InitCryoBoundaries().

338 { return fTrans.toWorldCoords(local); }

geo::LocalTransformationGeo::toWorldCoords

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

GlobalPoint_t toWorldCoords(LocalPoint_t const &local) const

Transforms a point from local frame to world frame.

Definition: LocalTransformationGeo.h:103

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

inline

Transform direction vector from local to world.

Definition at line 345 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformationGeo< StoredMatrix, LocalPoint, LocalVector >::toWorldCoords().

346 { return fTrans.toWorldCoords(local); }

geo::LocalTransformationGeo::toWorldCoords

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

GlobalPoint_t toWorldCoords(LocalPoint_t const &local) const

Transforms a point from local frame to world frame.

Definition: LocalTransformationGeo.h:103

const TPCGeo & geo::CryostatGeo::TPC ( unsigned int itpc ) const

Return the itpc'th TPC in the cryostat.

Exceptions

cet::exception (category "TPCOutOfRange") if no such TPC

Definition at line 174 of file CryostatGeo.cxx.

References TPCPtr().

Referenced by geo::buildDriftVolumes(), pma::ProjectionMatchingAlg::buildShowerSeg(), lar_pandora::LArPandoraOutput::CalculateT0(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), lar_pandora::LArPandoraInput::CreatePandoraHits2D(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), evd_tool::ICARUSDrawer::DetOutline3D(), hit::DisambigCheater::DisambigCheater(), evd::RecoBaseDrawer::DrawTrack2D(), cluster::EndPointAlg::EndPoint(), cluster::HoughBaseAlg::FastTransform(), trkf::SpacePointAlg::fillComplexSpacePoint(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), corner::CornerFinderAlg::get_feature_points_fast(), GetElement(), trkf::BezierTrack::GetProjectedPointUVWT(), trkf::BezierTrack::GetProjectedPointUVWX(), trkf::BezierTrack::GetTrackPitch(), lar_pandora::LArPandoraInput::GetTrueX0(), lar_pandora::PFParticleHitDumper::GetUVW(), groupTPCsByDriftDir(), HasElement(), cluster::fuzzyClusterAlg::InitFuzzy(), geo::ChannelMapStandardAlg::Initialize(), trkf::SpacePointAlg::makeSpacePoints(), MaxPlanes(), MaxWires(), cluster::ClusterCheater::produce(), vertex::HarrisVertexFinder::produce(), vertex::FeatureVertexFinder::produce(), evd::SimulationDrawer::SimulationDrawer(), SortSubVolumes(), apa::APAGeometryAlg::ThreeChanPos(), TPC(), geo::GeometryCore::TPC(), cluster::HoughBaseAlg::Transform(), apa::DisambigAlg::TrivialDisambig(), DUNE::NeutrinoTrackingEff::truthLength(), trkf::SpacePointAlg::update(), apa::DisambigAlg::UseEndPts(), and trkf::VertexFitAlg::VertexFit().

   {
     TPCGeo const* pTPC = TPCPtr(itpc);
     if(!pTPC){
       throw cet::exception("TPCOutOfRange") << "Request for non-existant TPC "
                                             << itpc << "\n";
     }
 
     return *pTPC;
   }

const TPCGeo& geo::CryostatGeo::TPC ( TPCID const & tpcid ) const

inline

Returns the TPC in tpcid from this cryostat.

Parameters

tpcid full TPC ID

Returns: a constant reference to the TPC in tpcid

Exceptions

cet::exception (category "TPCOutOfRange") if no such TPC

The cryostat number in tpcid is ignored, as it is ignored whether tpcid is invalid.

Definition at line 194 of file CryostatGeo.h.

References TPC(), and geo::TPCID::TPC.

195 { return TPC(tpcid.TPC); }

geo::CryostatGeo::TPC

const TPCGeo & TPC(unsigned int itpc) const

Return the itpc'th TPC in the cryostat.

Definition: CryostatGeo.cxx:174

TPCGeo const* geo::CryostatGeo::TPCPtr ( unsigned int itpc ) const

inline

Returns the TPC number itpc from this cryostat.

Parameters

itpc	the number of local TPC

Returns: a constant pointer to the TPC, or nullptr if it does not exist

Definition at line 217 of file CryostatGeo.h.

References fTPCs, and HasTPC().

Referenced by GetElementPtr(), TPC(), TPCPtr(), and geo::GeometryCore::TPCPtr().

218 { return HasTPC(itpc)? &(fTPCs[itpc]): nullptr; }

TPCList_t fTPCs

List of tpcs in this cryostat.

Definition: CryostatGeo.h:409

bool HasTPC(unsigned int itpc) const

Returns whether a TPC with index itpc is present in this cryostat.

Definition: CryostatGeo.h:164

TPCGeo const* geo::CryostatGeo::TPCPtr ( TPCID const & tpcid ) const

inline

Returns the TPC in tpcid from this cryostat.

Parameters

tpcid full TPC ID

Returns: a constant pointer to the TPC, or nullptr if it does not exist

The cryostat number in tpcid is ignored, as it is ignored whether tpcid is invalid.

Definition at line 228 of file CryostatGeo.h.

References geo::TPCID::TPC, and TPCPtr().

229 { return TPCPtr(tpcid.TPC); }

geo::CryostatGeo::TPCPtr

TPCGeo const * TPCPtr(unsigned int itpc) const

Returns the TPC number itpc from this cryostat.

Definition: CryostatGeo.h:217

auto const& geo::CryostatGeo::TPCs ( ) const

inline

Returns a container with references to all TPCs.

Returns: a container with references to all TPCs

The returned value can be used in a range-for loop like:

for (geo::TPCGeo const& tpc: cryo.TPCs()) { ... }

Definition at line 210 of file CryostatGeo.h.

References fTPCs.

Referenced by PositionToTPCptr().

210 { return fTPCs; }

geo::CryostatGeo::fVolume

TPCList_t fTPCs

List of tpcs in this cryostat.

Definition: CryostatGeo.h:409

void geo::CryostatGeo::UpdateAfterSorting ( geo::CryostatID cryoid )

Performs all needed updates after geometry has sorted the cryostats.

Definition at line 161 of file CryostatGeo.cxx.

References fID, fTPCs, and NTPC().

Referenced by toLocalCoords().

                                                            {
     
     // update the cryostat ID
     fID = cryoid;
     
     // trigger all the TPCs to update as well
     for (unsigned int tpc = 0; tpc < NTPC(); ++tpc)
       fTPCs[tpc].UpdateAfterSorting(geo::TPCID(fID, tpc));
     
   } // CryostatGeo::UpdateAfterSorting()

const TGeoVolume* geo::CryostatGeo::Volume ( ) const

inline

Pointer to ROOT's volume descriptor.

Definition at line 95 of file CryostatGeo.h.

References fVolume.

Referenced by InitCryoBoundaries().

95 { return fVolume; }

TGeoVolume * fVolume

Total volume of cryostat, called volCryostat in GDML file.

Definition: CryostatGeo.h:411

double geo::CryostatGeo::Width ( ) const

inline

Full width of the cryostat [cm].

Definition at line 87 of file CryostatGeo.h.

References HalfWidth().

Referenced by PrintCryostatInfo().

87 { return 2. * HalfWidth(); }

geo::CryostatGeo::HalfWidth

double HalfWidth() const

Half width of the cryostat [cm].

Definition: CryostatGeo.cxx:299

void geo::CryostatGeo::WorldToLocal	(	const double *	world,
		double *	cryo
	)		const

inline

Transform point from world frame to local cryostat frame.

Definition at line 349 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformation< StoredMatrix >::WorldToLocal().

350 { fTrans.WorldToLocal(world, cryo); }

geo::LocalTransformation::WorldToLocal

void WorldToLocal(double const *world, double *local) const

Transforms a point from world frame to local frame.

geo::LocalTransformation::WorldToLocal

LocalTransformation_t fTrans

Cryostat-to-world transformation.

Definition: CryostatGeo.h:408

template<typename Point >

Point geo::CryostatGeo::WorldToLocal ( Point const & world ) const

inline

Transform point from world frame to local cryostat frame.

Deprecated:: This method breaks the distinction between local and global vectors, since input and output vectors share the same type; use the "official" vector types geo::Point_t, geo::CryostatGeo::LocalPoint_t, geo::Vector_t and geo::CryostatGeo::LocalVector_t, and then use the method geo::CryostatGeo::toLocalCoords() instead.

Definition at line 361 of file CryostatGeo.h.

References fTrans, and geo::LocalTransformation< StoredMatrix >::WorldToLocal().

362 { return fTrans.WorldToLocal(world); }

void WorldToLocal(double const *world, double *local) const

Transforms a point from world frame to local frame.