LArSoft: geo::TPCGeo Class Reference

geo::PlaneGeo::ComposePoint

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

template<typename Point = DefaultPoint_t>

Point geo::TPCGeo::ComposePoint	(	double	distance,
		Projection_t const &	proj
	)		const

inline

Returns the 3D point from composition of projection and distance.

Template Parameters

Point type of point to be returned

Parameters

distance	distance of the target point from the reference plane
proj	projection of the target point on the reference plane

Returns: the 3D point from composition of projection and distance

See also: DecomposePoint()

The returned point is the reference point, translated by two 3D vectors:

a vector parallel to the plane normal, with norm the input distance
a vector lying on the plane, whose projection via PointProjection() gives the input projection

The choice of the projection reference point embodies the same convention used in PointProjection() and DecomposePoint().

Definition at line 480 of file TPCGeo.h.

References geo::PlaneGeo::ComposePoint(), proj, and ReferencePlane().

481 { return ReferencePlane().ComposePoint<Point>(distance, proj); }

Point ComposePoint(WireDecomposedVector_t const &decomp) const

Returns the 3D vector from composition of projection and distance.

Definition: PlaneGeo.h:818

dcel2d::Point

std::tuple< double, double, const reco::ClusterHit3D * > Point

Definitions used by the VoronoiDiagram algorithm.

Definition: DCEL.h:34

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

proj

Float_t proj

Definition: plot.C:34

double geo::TPCGeo::ComputeDriftDistance ( ) const

private

Computes the distance between the cathode and the last wire plane (last respect to the sorting order).

Definition at line 480 of file TPCGeo.cxx.

References geo::PlaneGeo::DistanceFromPlane(), fPlanes, and GetCathodeCenter().

Referenced by DriftDistance().

                                             {
     
     //
     // 1. find the center of the face of the TPC opposite to the anode
     // 2. compute the distance of it from the last wire plane
     //
     
     geo::PlaneGeo const& plane = fPlanes.back();
     return std::abs(plane.DistanceFromPlane(GetCathodeCenter()));
     
   } // TPCGeo::ComputeDriftDistance()

bool geo::BoxBoundedGeo::ContainsPosition	(	geo::Point_t const &	point,
		double	wiggle = `1.0`
	)		const

inlineinherited

Returns whether this volume contains the specified point.

Parameters

point	the point [cm]
wiggle	expansion factor for the range

Returns: whether the specified coordinate is in this volume

If the wiggle is larger than 1, each size of the volume is expanded by the wiggle factor. If the wiggle is less than 1, each size is shrinked.

Definition at line 207 of file BoxBoundedGeo.h.

References geo::BoxBoundedGeo::ContainsX(), and geo::BoxBoundedGeo::ContainsYZ().

Referenced by geo::BoxBoundedGeo::ContainsPosition(), geo::BoxBoundedGeo::ContainsYZ(), and geo::GeometryCore::MaterialName().

       {
         return ContainsX(point.X(), wiggle)
           && ContainsYZ(point.Y(), point.Z(), wiggle);
       } // ContainsPosition()

bool geo::BoxBoundedGeo::ContainsPosition	(	TVector3 const &	point,
		double	wiggle = `1.0`
	)		const

inherited

See also: ContainsPosition(geo::Point_t const&, double) const.

Definition at line 18 of file BoxBoundedGeo.cxx.

References geo::BoxBoundedGeo::ContainsPosition(), and geo::vect::toPoint().

19 { return ContainsPosition(geo::vect::toPoint(point), wiggle); }

geo::BoxBoundedGeo::ContainsPosition

::geo::Point_t toPoint(Point const &p)

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

bool ContainsPosition(geo::Point_t const &point, double wiggle=1.0) const

Returns whether this volume contains the specified point.

Definition: BoxBoundedGeo.h:207

bool geo::BoxBoundedGeo::ContainsPosition	(	double const *	point,
		double	wiggle = `1.0`
	)		const

inherited

See also: ContainsPosition(geo::Point_t const&, double) const.

Definition at line 23 of file BoxBoundedGeo.cxx.

References geo::vect::makePointFromCoords().

24 { return ContainsPosition(geo::vect::makePointFromCoords(point), wiggle); }

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::BoxBoundedGeo::ContainsPosition

bool ContainsPosition(geo::Point_t const &point, double wiggle=1.0) const

Returns whether this volume contains the specified point.

Definition: BoxBoundedGeo.h:207

bool geo::BoxBoundedGeo::ContainsX	(	double	x,
		double const	wiggle = `1`
	)		const

inlineinherited

Returns whether this TPC contains the specified world x coordinate.

Parameters

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

Returns: whether the specified coordinate is in this TPC

See also: ContainsPosition()

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 161 of file BoxBoundedGeo.h.

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

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

162 { return CoordinateContained(x, MinX(), MaxX(), wiggle); }

x

Float_t x

Definition: compare.C:6

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

geo::BoxBoundedGeo::CoordinateContained

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

geo::BoxBoundedGeo::CoordinateContained

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

geo::BoxBoundedGeo::CoordinateContained

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::TPCGeo::CoordinateContained	(	double	c,
		double	min,
		double	max,
		double	wiggle = `1.`
	)

inlinestatic

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 610 of file TPCGeo.h.

Referenced by CoordinateContained().

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

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

inlinestatic

Definition at line 617 of file TPCGeo.h.

References CoordinateContained(), FindPlane(), and MakePlane().

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

geo::TPCGeo::CoordinateContained

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

Returns whether the specified coordinate is in a range.

Definition: TPCGeo.h:610

DecomposedVector_t geo::TPCGeo::DecomposePoint ( geo::Point_t const & point ) const

inline

Decomposes a 3D point in two components.

Parameters

point the point to be decomposed

Returns: the two components of point, on the plane and orthogonal to it

The point is decomposed in:

a component orthogonal to the reference plane, expressed as a signed real number
a component lying on the reference plane, expressed as a 2D vector

The distance is from the reference plane (DistanceFromReferencePlane()). The projection on the plane is obtained following the same convention as PointProjection().

Definition at line 384 of file TPCGeo.h.

References geo::PlaneGeo::DecomposePointWidthDepth(), and ReferencePlane().

Referenced by DecomposePoint().

385 { return ReferencePlane().DecomposePointWidthDepth(point); }

geo::PlaneGeo::DecomposePointWidthDepth

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

WDDecomposedVector_t DecomposePointWidthDepth(geo::Point_t const &point) const

Decomposes a 3D point in two components.

Definition: PlaneGeo.h:878

DecomposedVector_t geo::TPCGeo::DecomposePoint ( TVector3 const & point ) const

inline

Decomposes a 3D point in two components.

Parameters

point the point to be decomposed

Returns: the two components of point, on the plane and orthogonal to it

The point is decomposed in:

a component orthogonal to the reference plane, expressed as a signed real number
a component lying on the reference plane, expressed as a 2D vector

The distance is from the reference plane (DistanceFromReferencePlane()). The projection on the plane is obtained following the same convention as PointProjection().

Definition at line 386 of file TPCGeo.h.

References DecomposePoint(), and geo::vect::toPoint().

387 { return DecomposePoint(geo::vect::toPoint(point)); }

geo::TPCGeo::DecomposePoint

::geo::Point_t toPoint(Point const &p)

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

DecomposedVector_t DecomposePoint(geo::Point_t const &point) const

Decomposes a 3D point in two components.

Definition: TPCGeo.h:384

short int geo::TPCGeo::DetectDriftDirection ( ) const

Returns the expected drift direction based on geometry.

The return value is coded as follow:

+1: positive x
+2: positive y
+3: positive z
-1: negative x
-2: negative y
-3: negative z
0: other (or algorithm failed)

The current implementation is based on the assumption that electrons in the middle of TPC will drift toward the wire planes, and it "never fails".

Definition at line 181 of file TPCGeo.cxx.

References geo::PlaneGeo::GetBoxCenter(), GetCenter(), ID(), Nplanes(), and Plane().

Referenced by GetCathodeCenterImpl(), nnet::TrainingDataAlg::getProjection(), pma::PMAlgTracker::mergeCoLinear(), ResetDriftDirection(), and toLocalCoords().

                                                {
     
     //
     // 1. determine the drift axis
     // 2. determine the drift direction on it
     // 
     // We assume that all the planes cover most of the TPC face; therefore,
     // the centre of the plane and the one of the TPC should be very close
     // to each other, when projected on the same drift plane.
     // Here we find which is the largest coordinate difference.
     
     if (Nplanes() == 0) {
       // chances are that we get this because stuff is not initialised yet,
       // and then even the ID might be wrong
       throw cet::exception("TPCGeo")
         << "DetectDriftDirection(): no planes in TPC " << std::string(ID())
         << "\n";
     }
     
     auto const TPCcenter = GetCenter();
     auto const PlaneCenter = Plane(0).GetBoxCenter(); // any will do
     
     auto const driftVector = PlaneCenter - TPCcenter; // approximation!
     
     if ((std::abs(driftVector.X()) > std::abs(driftVector.Y()))
       && (std::abs(driftVector.X()) > std::abs(driftVector.Z())))
     {
       // x is the solution
       return (driftVector.X() > 0)? +1: -1;
     }
     else if (std::abs(driftVector.Y()) > std::abs(driftVector.Z()))
     {
       // y is the man
       return (driftVector.Y() > 0)? +2: -2;
     }
     else {
       // z is the winner
       return (driftVector.Z() > 0)? +3: -3;
     }
     
   } // TPCGeo::DetectDriftDirection()

double geo::TPCGeo::DistanceFromReferencePlane ( geo::Point_t const & point ) const

inline

Returns the distance of the specified point from reference plane.

Parameters

point a point in world coordinates [cm]

Returns: the signed distance from the plane

The distance is defined positive if the point lies in the inner side with respect to the reference plane.

It should match the drift distance from that plane.

Definition at line 343 of file TPCGeo.h.

References geo::PlaneGeo::DistanceFromPlane(), and ReferencePlane().

Referenced by DistanceFromReferencePlane().

344 { return ReferencePlane().DistanceFromPlane(point); }

geo::PlaneGeo::DistanceFromPlane

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

double DistanceFromPlane(geo::Point_t const &point) const

Returns the distance of the specified point from the wire plane.

Definition: PlaneGeo.h:557

double geo::TPCGeo::DistanceFromReferencePlane ( TVector3 const & point ) const

inline

Returns the plane used for reference by projection methods.

Definition at line 345 of file TPCGeo.h.

References DistanceFromReferencePlane(), and geo::vect::toPoint().

346 { return DistanceFromReferencePlane(geo::vect::toPoint(point)); }

geo::TPCGeo::DistanceFromReferencePlane

double DistanceFromReferencePlane(geo::Point_t const &point) const

Returns the distance of the specified point from reference plane.

Definition: TPCGeo.h:343

geo::TPCGeo::fDriftDirection

::geo::Point_t toPoint(Point const &p)

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

template<typename Vector >

Vector geo::TPCGeo::DriftDir ( ) const

Returns the direction of the drift (vector pointing toward the planes).

Definition at line 687 of file TPCGeo.h.

References fDriftDir.

Referenced by checkTPCcoords(), DriftDirection(), groupTPCsByDriftDir(), and UpdateAfterSorting().

688 { return geo::vect::convertTo<Vector>(fDriftDir); }

geo::TPCGeo::fDriftDir

geo::Vector_t fDriftDir

Direction electrons drift along.

Definition: TPCGeo.h:650

DriftDirection_t geo::TPCGeo::DriftDirection ( ) const

inline

Returns an enumerator value describing the drift direction.

Definition at line 127 of file TPCGeo.h.

References DriftDir(), and fDriftDirection.

Referenced by util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadout::DriftIonizationElectrons(), and detsim::SimDriftElectrons::produce().

127 { return fDriftDirection; }

DriftDirection_t fDriftDirection

Direction of the electron drift in the TPC.

Definition: TPCGeo.h:635

double geo::TPCGeo::DriftDistance ( ) const

inline

Drift distance is defined as the distance between the last anode plane and the opposite face of the TPC, in centimeters.

Definition at line 135 of file TPCGeo.h.

References ComputeDriftDistance(), fPlanes, and Nplanes().

Referenced by pma::PMAlgCosmicTagger::GetDimensions(), and PrintTPCInfo().

135 { return ComputeDriftDistance(); }

geo::TPCGeo::ComputeDriftDistance

double ComputeDriftDistance() const

Definition: TPCGeo.cxx:480

void geo::TPCGeo::DriftPoint	(	geo::Point_t &	position,
		double	distance
	)		const

inline

Shifts the position of an electron drifted by a distance.

Parameters

position	_(modified)_ the position of the electron
distance	drift distance to shift the electron by [cm]

See also: geo::Plane::DriftPoint(), ReferencePlane()

This operation is delegated to the reference plane (see geo::Plane::DriftPoint()).

Definition at line 360 of file TPCGeo.h.

References geo::PlaneGeo::DriftPoint(), and ReferencePlane().

361 { ReferencePlane().DriftPoint(position, distance); }

geo::PlaneGeo::DriftPoint

void DriftPoint(geo::Point_t &position, double distance) const

Shifts the position of an electron drifted by a distance.

Definition: PlaneGeo.h:579

geo::PlaneGeo::DriftPoint

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

void geo::TPCGeo::DriftPoint	(	TVector3 &	position,
		double	distance
	)		const

inline

Shifts the position of an electron drifted by a distance.

Parameters

position	_(modified)_ the position of the electron
distance	drift distance to shift the electron by [cm]

See also: geo::Plane::DriftPoint(), ReferencePlane()

This operation is delegated to the reference plane (see geo::Plane::DriftPoint()).

Definition at line 362 of file TPCGeo.h.

References geo::PlaneGeo::DriftPoint(), and ReferencePlane().

363 { ReferencePlane().DriftPoint(position, distance); }

void DriftPoint(geo::Point_t &position, double distance) const

Shifts the position of an electron drifted by a distance.

Definition: PlaneGeo.h:579

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

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::TPCGeo::FindPlane	(	GeoNodePath_t &	path,
		size_t	depth
	)

private

Definition at line 148 of file TPCGeo.cxx.

References MakePlane().

Referenced by CoordinateContained(), and TPCGeo().

   {
 
     const char* nm = path[depth]->GetName();
     if( (strncmp(nm, "volTPCPlane", 11) == 0) ){
       this->MakePlane(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->FindPlane(path, deeper);
     }
   
   }

geo::PlaneGeo const& geo::TPCGeo::FirstPlane ( ) const

inline

Returns the first wire plane (the closest to TPC center).

Definition at line 225 of file TPCGeo.h.

References fPlanes.

Referenced by ReferencePlane().

225 { return fPlanes[0]; }

geo::TPCGeo::fActiveCenter

std::vector< PlaneGeo > fPlanes

List of planes in this plane.

Definition: TPCGeo.h:632

template<typename Point = DefaultPoint_t>

Point geo::TPCGeo::GetActiveVolumeCenter ( ) const

inline

Returns the center of the TPC active volume in world coordinates [cm].

Definition at line 247 of file TPCGeo.h.

References fActiveCenter.

Referenced by evd_tool::StandardDrawer::DrawActiveTPCoutline(), and evd::SimulationDrawer::SimulationDrawer().

248 { return geo::vect::convertTo<Point>(fActiveCenter); }

geo::Point_t fActiveCenter

Center of the active volume, in world coordinates [cm].

Definition: TPCGeo.h:638

template<typename Point = DefaultPoint_t>

Point geo::TPCGeo::GetCathodeCenter ( ) const

inline

Returns the center of the active volume face opposite to the wire planes [cm]

Definition at line 253 of file TPCGeo.h.

References GetCathodeCenterImpl().

Referenced by geo::DriftPartitions::computeCoverage(), and ComputeDriftDistance().

254 { return geo::vect::convertTo<Point>(GetCathodeCenterImpl()); }

geo::TPCGeo::GetCathodeCenterImpl

geo::Point_t GetCathodeCenterImpl() const

Definition: TPCGeo.cxx:360

geo::Point_t geo::TPCGeo::GetCathodeCenterImpl ( ) const

private

Definition at line 360 of file TPCGeo.cxx.

References ActiveHalfHeight(), ActiveHalfWidth(), ActiveLength(), DetectDriftDirection(), and geo::vect::Xcoord().

Referenced by GetCathodeCenter().

                                                 {
     
     //
     // 1. find the center of the face of the TPC opposite to the anode
     // 2. compute the distance of it from the last wire plane
     //
     
     //
     // find the cathode center
     //
     geo::Point_t cathodeCenter = GetActiveVolumeCenter<geo::Point_t>();
     switch (DetectDriftDirection()) {
       case -1:
         geo::vect::Xcoord(cathodeCenter) += ActiveHalfWidth();
       //  cathodeCenter.SetX(cathodeCenter.X() + ActiveHalfWidth());
         break;
       case +1:
         cathodeCenter.SetX(cathodeCenter.X() - ActiveHalfWidth());
         break;
       case -2:
         cathodeCenter.SetY(cathodeCenter.Y() + ActiveHalfHeight());
         break;
       case +2:
         cathodeCenter.SetY(cathodeCenter.Y() - ActiveHalfHeight());
         break;
       case -3:
         cathodeCenter.SetZ(cathodeCenter.Z() + ActiveLength() / 2.0);
         break;
       case +3:
         cathodeCenter.SetZ(cathodeCenter.Z() - ActiveLength() / 2.0);
         break;
       case 0:
       default:
         // in this case, a better algorithm is probably needed
         throw cet::exception("TPCGeo")
           << "CathodeCenter(): Can't determine the cathode plane (code="
           << DetectDriftDirection() << ")\n";
     } // switch
     return cathodeCenter;
     
   } // TPCGeo::GetCathodeCenterImpl()

template<typename Point >

Point geo::TPCGeo::GetCenter ( ) const

Returns the center of the TPC volume in world coordinates [cm].

Definition at line 693 of file TPCGeo.h.

References toWorldCoords().

Referenced by DetectDriftDirection(), evd_tool::ICARUSDrawer::DetOutline3D(), LastPlane(), ResetDriftDirection(), evd::SimulationDrawer::SimulationDrawer(), and SortPlanes().

                                  {
   
   // convert the origin (default constructed TVector)
   return geo::vect::convertTo<Point>(toWorldCoords(LocalPoint_t{ 0.0, 0.0, 0.0 }));
   
 } // geo::TPCGeo::GetCenter()

const PlaneGeo& geo::TPCGeo::GetElement ( PlaneID const & planeid ) const

inline

Returns the plane in planeid from this TPC.

Parameters

planeid full plane ID

Returns: a constant reference to the plane in planeid

Exceptions

cet::exception (category "PlaneOutOfRange") if no such plane

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 194 of file TPCGeo.h.

References Plane().

195 { return Plane(planeid); }

geo::TPCGeo::Plane

PlaneGeo const & Plane(geo::View_t view) const

Return the plane in the tpc with View_t view.

Definition: TPCGeo.cxx:298

PlaneGeo const* geo::TPCGeo::GetElementPtr ( PlaneID const & planeid ) const

inline

Returns the plane in planeid from this TPC.

Parameters

planeid full plane ID

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

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 217 of file TPCGeo.h.

References PlanePtr(), and SmallestPlane().

218 { return PlanePtr(planeid); }

geo::TPCGeo::PlanePtr

PlaneGeo const * PlanePtr(unsigned int iplane) const

Returns the plane number iplane from this TPC.

Definition: TPCGeo.h:203

template<typename Point = geo::Point_t>

Point geo::TPCGeo::GetFrontFaceCenter ( ) const

inline

Returns the center of the active TPC volume side facing negative z.

Definition at line 258 of file TPCGeo.h.

References GetFrontFaceCenterImpl().

259 { return geo::vect::convertTo<Point>(GetFrontFaceCenterImpl()); }

geo::TPCGeo::GetFrontFaceCenterImpl

geo::Point_t GetFrontFaceCenterImpl() const

Definition: TPCGeo.cxx:404

geo::Point_t geo::TPCGeo::GetFrontFaceCenterImpl ( ) const

private

Definition at line 404 of file TPCGeo.cxx.

References ActiveBoundingBox().

Referenced by GetFrontFaceCenter().

                                                   {
     auto const& activeBox = ActiveBoundingBox();
     return { activeBox.CenterX(), activeBox.CenterY(), activeBox.MinZ() };
   } // TPCGeo::GetFrontFaceCenterImpl()

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::TPCGeo::HalfHeight ( ) const

inline

Height is associated with y coordinate [cm].

Definition at line 103 of file TPCGeo.h.

References fHalfHeight.

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), larg::LArSimChannelAna::ensureHists(), hit::MagDriftAna::ensureHists(), pma::PMAlgCosmicTagger::GetDimensions(), pma::PMAlgStitching::GetTPCXOffsets(), Height(), InitTPCBoundaries(), and evd::SimulationDrawer::SimulationDrawer().

103 { return fHalfHeight; }

geo::TPCGeo::fHalfHeight

double fHalfHeight

Half height of total volume.

Definition: TPCGeo.h:644

double geo::TPCGeo::HalfLength ( ) const

inline

Length is associated with z coordinate [cm].

Definition at line 109 of file TPCGeo.h.

References fLength.

Referenced by InitTPCBoundaries().

109 { return fLength / 2.0; }

geo::TPCGeo::fLength

double fLength

Length of total volume.

Definition: TPCGeo.h:645

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::TPCGeo::HalfWidth ( ) const

inline

Width is associated with x coordinate [cm].

Definition at line 99 of file TPCGeo.h.

References fHalfWidth.

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), larg::LArSimChannelAna::ensureHists(), hit::MagDriftAna::ensureHists(), pma::PMAlgCosmicTagger::GetDimensions(), pma::PMAlgStitching::GetTPCXOffsets(), InitTPCBoundaries(), evd::SimulationDrawer::SimulationDrawer(), and Width().

99 { return fHalfWidth; }

geo::TPCGeo::fHalfWidth

double fHalfWidth

Half width of total volume.

Definition: TPCGeo.h:643

bool geo::TPCGeo::HasElement ( unsigned int iplane ) const

inline

Returns whether a plane with index iplane is present in this TPC.

Parameters

iplane index of plane in this TPC

Returns: whether the plane with index iplane is present in this TPC

Definition at line 156 of file TPCGeo.h.

References HasPlane().

156 { return HasPlane(iplane); }

bool HasPlane(unsigned int iplane) const

Returns whether a plane with index iplane is present in this TPC.

Definition: TPCGeo.h:155

bool geo::TPCGeo::HasElement ( geo::PlaneID const & planeid ) const

inline

Returns whether the plane in planeid is present in this TPC.

Parameters

planeid full plane ID

Returns: whether the plane in planeid is present in this TPC

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 170 of file TPCGeo.h.

References HasPlane(), and Plane().

171 { return HasPlane(planeid); }

bool HasPlane(unsigned int iplane) const

Returns whether a plane with index iplane is present in this TPC.

Definition: TPCGeo.h:155

bool geo::TPCGeo::HasPlane ( unsigned int iplane ) const

inline

Returns whether a plane with index iplane is present in this TPC.

Parameters

iplane index of plane in this TPC

Returns: whether the plane with index iplane is present in this TPC

Definition at line 155 of file TPCGeo.h.

References Nplanes().

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), pma::VtxCandidate::ComputeMse2D(), pma::PMAlgStitching::GetTPCXOffsets(), pma::ProjectionMatchingAlg::guideEndpoints(), HasElement(), HasPlane(), geo::GeometryCore::HasPlane(), pma::PMAlgTracker::matchCluster(), pma::Node3D::Node3D(), PlanePtr(), trkf::PMAlgTrajFitter::produce(), trkf::PMAlgTrackMaker::produce(), and pma::ProjectionMatchingAlg::TestTrk().

155 { return iplane < Nplanes(); }

geo::TPCGeo::Nplanes

unsigned int Nplanes() const

Number of planes in this tpc.

Definition: TPCGeo.h:145

bool geo::TPCGeo::HasPlane ( geo::PlaneID const & planeid ) const

inline

Returns whether the plane in planeid is present in this TPC.

Parameters

planeid full plane ID

Returns: whether the plane in planeid is present in this TPC

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 168 of file TPCGeo.h.

References HasPlane(), and geo::PlaneID::Plane.

169 { return HasPlane(planeid.Plane); }

geo::BoxBoundedGeo::InFiducialX

bool HasPlane(unsigned int iplane) const

Returns whether a plane with index iplane is present in this TPC.

Definition: TPCGeo.h:155

double geo::TPCGeo::Height ( ) const

inline

Height is associated with y coordinate [cm].

Definition at line 105 of file TPCGeo.h.

References HalfHeight().

Referenced by PrintTPCInfo().

105 { return 2.0 * HalfHeight(); }

geo::TPCGeo::HalfHeight

double HalfHeight() const

Height is associated with y coordinate [cm].

Definition: TPCGeo.h:103

template<typename Vector = DefaultVector_t>

decltype(auto) geo::TPCGeo::HeightDir ( ) const

inline

Returns the direction Height() is measured on.

Definition at line 120 of file TPCGeo.h.

References fHeightDir.

120 { return geo::vect::convertTo<Vector>(fHeightDir); }

geo::TPCGeo::fHeightDir

geo::Vector_t fHeightDir

Direction height refers to.

Definition: TPCGeo.h:648

geo::TPCID const& geo::TPCGeo::ID ( ) const

inline

Returns the identifier of this TPC.

Definition at line 278 of file TPCGeo.h.

References fID.

Referenced by checkTPCcoords(), DetectDriftDirection(), evd_tool::ICARUSDrawer::DetOutline3D(), geo::GeometryCore::PositionToTPC(), geo::CryostatGeo::PositionToTPCID(), geo::GeometryCore::PositionToTPCID(), and PrintTPCInfo().

278 { return fID; }

geo::TPCGeo::fID

geo::TPCID fID

ID of this TPC.

Definition: TPCGeo.h:654

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

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::TPCGeo::InitTPCBoundaries ( )

private

Recomputes the TPC boundary.

Definition at line 494 of file TPCGeo.cxx.

References ActiveHalfHeight(), ActiveHalfLength(), ActiveHalfWidth(), fActiveBox, HalfHeight(), HalfLength(), HalfWidth(), geo::BoxBoundedGeo::SetBoundaries(), and toWorldCoords().

Referenced by TPCGeo().

                                  {
     // note that this assumes no rotations of the TPC
     // (except for rotations of a flat angle around one of the three main axes);
     // to avoid this, we should transform the six vertices
     // rather than just the centre
     
     // we rely on the asumption that the center of TPC is at the local origin
     SetBoundaries(
       toWorldCoords(LocalPoint_t(-HalfWidth(), -HalfHeight(), -HalfLength())),
       toWorldCoords(LocalPoint_t(+HalfWidth(), +HalfHeight(), +HalfLength()))
       );
     
     // the center of the active volume may be elsewhere than the local origin:
     auto const& activeCenter = GetActiveVolumeCenter<geo::Point_t>();
     fActiveBox.SetBoundaries(
       activeCenter.X() - ActiveHalfWidth(),  activeCenter.X() + ActiveHalfWidth(),
       activeCenter.Y() - ActiveHalfHeight(), activeCenter.Y() + ActiveHalfHeight(),
       activeCenter.Z() - ActiveHalfLength(), activeCenter.Z() + ActiveHalfLength()
       );
 
     
   } // CryostatGeo::InitTPCBoundaries()

geo::PlaneGeo const& geo::TPCGeo::LastPlane ( ) const

inline

Returns the last wire plane (the farther from TPC center).

Definition at line 228 of file TPCGeo.h.

References fPlanes, GetCenter(), MaxWires(), Nplanes(), and Views().

Referenced by geo::DriftPartitions::computeCoverage().

228 { return fPlanes[Nplanes() - 1]; }

geo::TPCGeo::Nplanes

unsigned int Nplanes() const

Number of planes in this tpc.

Definition: TPCGeo.h:145

Referenced by DUNE::NeutrinoShowerEff::beginJob(), DUNE::NeutrinoTrackingEff::beginJob(), DUNE::MuonTrackingEff::beginJob(), tca::FillWireHitRange(), cluster::ClusterCrawlerAlg::FindHammerClusters(), pma::PMAlgCosmicTagger::GetDimensions(), pma::PMAlgStitching::GetTPCXOffsets(), tca::InsideTPC(), evd::RecoBaseDrawer::OpFlashOrtho(), evd::Ortho3DPad::Ortho3DPad(), evgen::NDKGen::produce(), geo::sortTPCStandard(), trkf::CCTrackMaker::TagCosmics(), apa::DisambigAlg::TrivialDisambig(), apa::DisambigAlg::UseEndPts(), and cluster::ClusterCrawlerAlg::VtxMatch().

std::vector< PlaneGeo > fPlanes

List of planes in this plane.

Definition: TPCGeo.h:632

double geo::TPCGeo::Length ( ) const

inline

Length is associated with z coordinate [cm].

Definition at line 107 of file TPCGeo.h.

References fLength.

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), larg::LArSimChannelAna::ensureHists(), hit::MagDriftAna::ensureHists(), pma::PMAlgCosmicTagger::GetDimensions(), pma::PMAlgStitching::GetTPCXOffsets(), PrintTPCInfo(), evd::SimulationDrawer::SimulationDrawer(), and apa::APAGeometryAlg::ThreeChanPos().

107 { return fLength; }

geo::TPCGeo::fLength

double fLength

Length of total volume.

Definition: TPCGeo.h:645

template<typename Vector = DefaultVector_t>

decltype(auto) geo::TPCGeo::LengthDir ( ) const

inline

Returns the direction Length() is measured on.

Definition at line 124 of file TPCGeo.h.

References fLengthDir.

124 { return geo::vect::convertTo<Vector>(fLengthDir); }

geo::TPCGeo::fLengthDir

geo::Vector_t fLengthDir

Direction length refers to.

Definition: TPCGeo.h:649

void geo::TPCGeo::LocalToWorld	(	const double *	tpc,
		double *	world
	)		const

inline

Transform point from local TPC frame to world frame.

Definition at line 490 of file TPCGeo.h.

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

491 { fTrans.LocalToWorld(tpc, 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

TPC-to-world transformation.

Definition: TPCGeo.h:630

TVector3 geo::TPCGeo::LocalToWorld ( const TVector3 & local ) const

inline

Transform point from local TPC frame to world frame.

Definition at line 494 of file TPCGeo.h.

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

495 { return fTrans.LocalToWorld<TVector3>(local); }

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

Transforms a point from local frame to world frame.

geo::LocalTransformation::LocalToWorldVect

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

void geo::TPCGeo::LocalToWorldVect	(	const double *	tpc,
		double *	world
	)		const

inline

Transform direction vector from local to world.

Definition at line 502 of file TPCGeo.h.

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

503 { fTrans.LocalToWorldVect(tpc, world); }

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

Transforms a vector from local frame to world frame.

geo::LocalTransformation::LocalToWorldVect

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

TVector3 geo::TPCGeo::LocalToWorldVect ( const TVector3 & local ) const

inline

Transform direction vector from local to world.

Definition at line 506 of file TPCGeo.h.

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

507 { return fTrans.LocalToWorldVect<TVector3>(local); }

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

Transforms a vector from local frame to world frame.

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

void geo::TPCGeo::MakePlane	(	GeoNodePath_t &	path,
		size_t	depth
	)

private

Definition at line 174 of file TPCGeo.cxx.

References fPlanes.

Referenced by CoordinateContained(), and FindPlane().

   {
     fPlanes.emplace_back(path, depth);
   }

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(), geo::CryostatGeo::PrintCryostatInfo(), and PrintTPCInfo().

138 { return c_max; }

Referenced by geo::CryostatGeo::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::TPCGeo::MaxWires ( ) const

Returns the largest number of wires among the planes in this TPC.

Definition at line 332 of file TPCGeo.cxx.

References fPlanes.

Referenced by LastPlane(), geo::CryostatGeo::MaxWires(), and PrintTPCInfo().

                                       {
     unsigned int maxWires = 0;
     for (geo::PlaneGeo const& plane: fPlanes) {
       unsigned int maxWiresInPlane = plane.Nwires();
       if (maxWiresInPlane > maxWires) maxWires = maxWiresInPlane;
     } // for
     return maxWires;
   } // TPCGeo::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 geo::CryostatGeo::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 geo::CryostatGeo::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(), geo::CryostatGeo::PrintCryostatInfo(), and PrintTPCInfo().

135 { return c_min; }

Referenced by geo::CryostatGeo::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 geo::CryostatGeo::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 geo::CryostatGeo::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::TPCGeo::NElements ( ) const

inline

Number of planes in this tpc.

Definition at line 146 of file TPCGeo.h.

References fPlanes.

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

146 { return fPlanes.size(); }

Referenced by calo::LinearEnergyAlg::CalculateEnergy(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), DetectDriftDirection(), hit::DisambigCheater::DisambigCheater(), DriftDistance(), larg4::LArVoxelReadout::DriftIonizationElectrons(), trkf::SpacePointAlg::fillComplexSpacePoint(), tca::FillWireHitRange(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), tca::FindShowers3D(), trkf::BezierTrack::GetProjectedPointUVWT(), trkf::BezierTrack::GetProjectedPointUVWX(), trkf::BezierTrack::GetTrackPitch(), groupByDriftCoordinate(), HasPlane(), geo::ChannelMapStandardAlg::Initialize(), LastPlane(), trkf::SpacePointAlg::makeSpacePoints(), tca::TruthMatcher::MatchTruth(), geo::CryostatGeo::MaxPlanes(), tca::MergePFPTjs(), geo::GeometryCore::Print(), PrintTPCInfo(), detsim::SimDriftElectrons::produce(), vertex::FeatureVertexFinder::produce(), cluster::ClusterCrawlerAlg::RunCrawler(), tca::TrajClusterAlg::RunTrajClusterAlg(), SortSubVolumes(), trkf::CosmicTrackerAlg::TrackTrajectory(), trkf::SpacePointAlg::update(), UpdateAfterSorting(), UpdatePlaneViewCache(), pma::Node3D::UpdateProj2D(), trkf::VertexFitAlg::VertexFit(), cluster::ClusterCrawlerAlg::VtxMatch(), and geo::GeometryCore::WireAngleToVertical().

std::vector< PlaneGeo > fPlanes

List of planes in this plane.

Definition: TPCGeo.h:632

unsigned int geo::TPCGeo::Nplanes ( ) const

inline

Number of planes in this tpc.

Definition at line 145 of file TPCGeo.h.

145 { return fPlanes.size(); }

Referenced by hit::MagDriftAna::analyze(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), evd::RecoBaseDrawer::Cluster2D(), pma::CmToWireDrift(), lar_pandora::LArPandoraInput::CreatePandoraHits2D(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), DetectDriftDirection(), hit::DisambigCheater::DisambigCheater(), evd::RecoBaseDrawer::DrawTrack2D(), evd::RecoBaseDrawer::DrawTrackVertexAssns2D(), larg4::LArVoxelReadout::DriftIonizationElectrons(), cluster::EndPointAlg::EndPoint(), evd::RecoBaseDrawer::EndPoint2D(), evd::RecoBaseDrawer::Event2D(), cluster::HoughBaseAlg::FastTransform(), corner::CornerFinderAlg::get_feature_points_fast(), pma::ProjectionMatchingAlg::GetCloseHits(), GetElement(), pma::GetProjectionToPlane(), trkf::BezierTrack::GetTrackPitch(), lar_pandora::PFParticleHitDumper::GetUVW(), HasElement(), geo::ChannelMapStandardAlg::Initialize(), lar_pandora::LArPandoraGeometry::LoadGeometry(), trkf::SpacePointAlg::makeSpacePoints(), Plane(), geo::GeometryCore::Print(), cluster::ClusterCheater::produce(), vertex::HarrisVertexFinder::produce(), evd::RecoBaseDrawer::Prong2D(), evd::RecoBaseDrawer::RecoBaseDrawer(), ResetDriftDirection(), geo::GeometryCore::ThirdPlane_dTdW(), geo::GeometryCore::ThirdPlaneSlope(), cluster::HoughBaseAlg::Transform(), trkf::SpacePointAlg::update(), pma::Node3D::UpdateProj2D(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc_test(), geo::GeometryCore::WireAngleToVertical(), pma::WireDriftToCm(), and WirePitch().

std::vector< PlaneGeo > fPlanes

List of planes in this plane.

Definition: TPCGeo.h:632

const PlaneGeo & geo::TPCGeo::Plane ( geo::View_t view ) const

Return the plane in the tpc with View_t view.

Definition at line 298 of file TPCGeo.cxx.

References fPlanes, fViewToPlaneNumber, and geo::PlaneID::InvalidID.

   {
     geo::PlaneID::PlaneID_t const p = fViewToPlaneNumber[size_t(view)];
     if (p == geo::PlaneID::InvalidID) {
       throw cet::exception("TPCGeo")
         << "TPCGeo[" << ((void*) this) << "]::Plane(): no plane for view #"
         << (size_t) view << "\n";
     }
     return fPlanes[p];
   } // TPCGeo::Plane(geo::View_t)

const PlaneGeo & geo::TPCGeo::Plane ( unsigned int iplane ) const

Return the iplane'th plane in the TPC.

Exceptions

cet::exception (category "PlaneOutOfRange") if no such plane

Definition at line 287 of file TPCGeo.cxx.

References PlanePtr().

   {
     geo::PlaneGeo const* pPlane = PlanePtr(iplane);
     if (!pPlane){
       throw cet::exception("PlaneOutOfRange") << "Request for non-existant plane " << iplane << "\n";
     }
 
     return *pPlane;
   }

const PlaneGeo& geo::TPCGeo::Plane ( PlaneID const & planeid ) const

inline

Returns the plane in planeid from this TPC.

Parameters

planeid full plane ID

Returns: a constant reference to the plane in planeid

Exceptions

cet::exception (category "PlaneOutOfRange") if no such plane

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 192 of file TPCGeo.h.

References Plane(), and geo::PlaneID::Plane.

193 { return Plane(planeid.Plane); }

geo::TPCGeo::Plane

PlaneGeo const & Plane(geo::View_t view) const

Return the plane in the tpc with View_t view.

Definition: TPCGeo.cxx:298

double geo::TPCGeo::Plane0Pitch ( unsigned int p ) const

Returns the center of the TPC volume in world coordinates [cm].

Definition at line 354 of file TPCGeo.cxx.

References fPlane0Pitch.

Referenced by ActiveBoundingBox(), larg4::LArVoxelReadout::DriftIonizationElectrons(), groupByDriftCoordinate(), detsim::SimDriftElectrons::produce(), and trkf::SpacePointAlg::update().

   {
     return fPlane0Pitch[p];
   }

const double * geo::TPCGeo::PlaneLocation ( unsigned int p ) const

Returns the coordinates of the center of the specified plane [cm].

Definition at line 412 of file TPCGeo.cxx.

References fPlaneLocation.

Referenced by ActiveBoundingBox(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadout::DriftIonizationElectrons(), pma::PMAlgStitching::GetTPCXOffsets(), detsim::SimDriftElectrons::produce(), DUNE::NeutrinoTrackingEff::truthLength(), and trkf::SpacePointAlg::update().

   {
     return &fPlaneLocation[p][0];
   }

double geo::TPCGeo::PlanePitch	(	unsigned int	p1 = `0`,
		unsigned int	p2 = `1`
	)		const

Returns the center of the TPC volume in world coordinates [cm].

Definition at line 418 of file TPCGeo.cxx.

References fPlane0Pitch.

Referenced by ActiveBoundingBox(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadout::DriftIonizationElectrons(), geo::GeometryCore::PlanePitch(), and detsim::SimDriftElectrons::produce().

   {
     return std::abs(fPlane0Pitch[p2] - fPlane0Pitch[p1]);
   }

PlaneGeo const* geo::TPCGeo::PlanePtr ( unsigned int iplane ) const

inline

Returns the plane number iplane from this TPC.

Parameters

iplane the number of local plane

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

Definition at line 203 of file TPCGeo.h.

References fPlanes, and HasPlane().

Referenced by GetElementPtr(), Plane(), PlanePtr(), and geo::GeometryCore::PlanePtr().

204 { return HasPlane(iplane)? &(fPlanes[iplane]): nullptr; }

bool HasPlane(unsigned int iplane) const

Returns whether a plane with index iplane is present in this TPC.

Definition: TPCGeo.h:155

std::vector< PlaneGeo > fPlanes

List of planes in this plane.

Definition: TPCGeo.h:632

PlaneGeo const* geo::TPCGeo::PlanePtr ( PlaneID const & planeid ) const

inline

Returns the plane in planeid from this TPC.

Parameters

planeid full plane ID

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

The cryostat and TPC numbers in planeid are ignored, as it is ignored whether planeid is invalid.

Definition at line 215 of file TPCGeo.h.

References geo::PlaneID::Plane, and PlanePtr().

216 { return PlanePtr(planeid.Plane); }

geo::TPCGeo::PlanePtr

PlaneGeo const * PlanePtr(unsigned int iplane) const

Returns the plane number iplane from this TPC.

Definition: TPCGeo.h:203

Projection_t geo::TPCGeo::PointProjection ( geo::Point_t const & point ) const

inline

Returns the projection of the specified point on the plane.

Parameters

point the 3D point to be projected, in world coordinates

Returns: a 2D vector representing the projection of point on the plane

The returned vector is a 2D vector expressing the projection of the point (from world coordinates) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. The origin point is returned by ProjectionReferencePoint(). All coordinates are in centimeters.

Definition at line 417 of file TPCGeo.h.

References Projection().

418 { return Projection(point); }

Projection_t Projection(geo::Point_t const &point) const

Returns the projection of the specified point on the plane.

Definition: TPCGeo.h:415

Projection_t geo::TPCGeo::PointProjection ( TVector3 const & point ) const

inline

Returns the projection of the specified point on the plane.

Parameters

point the 3D point to be projected, in world coordinates

Returns: a 2D vector representing the projection of point on the plane

The returned vector is a 2D vector expressing the projection of the point (from world coordinates) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. The origin point is returned by ProjectionReferencePoint(). All coordinates are in centimeters.

Definition at line 419 of file TPCGeo.h.

References Projection(), and geo::vect::toPoint().

420 { return Projection(geo::vect::toPoint(point)); }

::geo::Point_t toPoint(Point const &p)

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

geo::PlaneGeo::PointWidthDepthProjection

Projection_t Projection(geo::Point_t const &point) const

Returns the projection of the specified point on the plane.

Definition: TPCGeo.h:415

template<typename Stream >

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

Prints information about this TPC.

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 TPC ID
1 _(default)_: also center and size
2: also drift direction, cathode position and number of planes
3: also maximum number of wires per plane
4: also information on main direction
5: also information on bounding box

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

Definition at line 703 of file TPCGeo.h.

References ActiveBoundingBox(), ActiveHeight(), ActiveLength(), ActiveWidth(), BoundingBox(), DriftDistance(), Height(), ID(), Length(), geo::BoxBoundedGeo::Max(), MaxWires(), geo::BoxBoundedGeo::Min(), Nplanes(), and Width().

Referenced by geo::part::PartitionDataDescriber< geo::TPCGeo >::PartitionDataDescriber(), geo::GeometryCore::Print(), and toLocalCoords().

         {
   
   //----------------------------------------------------------------------------
   out << "TPC " << std::string(ID());
   
   if (verbosity-- <= 0) return; // 0
   
   //----------------------------------------------------------------------------
   out
     << " (" << Width() << " x " << Height() << " x " << Length() << ") cm^3 at "
       << GetCenter<geo::Point_t>();
   
   if (verbosity-- <= 0) return; // 1
   
   //----------------------------------------------------------------------------
   
   out << "\n" << indent
     << "drift direction " << DriftDir<geo::Vector_t>()
       << " from cathode around " << GetCathodeCenter<geo::Point_t>()
       << " through " << DriftDistance() << " cm toward "
       << Nplanes() << " wire planes"
     ;
   
   if (verbosity-- <= 0) return; // 2
   
   //----------------------------------------------------------------------------
   out << "\n" << indent
     << "maximum wires on any plane: " << MaxWires();
     
   if (verbosity-- <= 0) return; // 3
   
   //----------------------------------------------------------------------------
   out << "\n" << indent
     << "active volume ("
       << ActiveWidth() << " x " << ActiveHeight() << " x " << ActiveLength()
       << ") cm^3, front face at " << GetFrontFaceCenter<geo::Point_t>()
       << " cm;"
     << "\n" << indent << "main directions:"
       << " width " << WidthDir<geo::Vector_t>()
       << " height " << HeightDir<geo::Vector_t>()
       << " length " << LengthDir<geo::Vector_t>()
     ;
     
   if (verbosity-- <= 0) return; // 4
   
   //----------------------------------------------------------------------------
   // print also the containing box
   geo::BoxBoundedGeo const& box = BoundingBox();
   out << "\n" << indent
     << "bounding box: " << box.Min() << " -- " << box.Max();
   
 //  if (verbosity-- <= 0) return; // 5
   
   //----------------------------------------------------------------------------
   // print also the active box
   geo::BoxBoundedGeo const& activeBox = ActiveBoundingBox();
   out << "\n" << indent
     << "active volume box: " << activeBox.Min() << " -- "  << activeBox.Max();
   
 //  if (verbosity-- <= 0) return; // 6
   
   //----------------------------------------------------------------------------
 } // geo::TPCGeo::PrintTPCInfo()

Projection_t geo::TPCGeo::Projection ( geo::Point_t const & point ) const

inline

Returns the projection of the specified point on the plane.

Parameters

point the 3D point to be projected, in world coordinates

Returns: a 2D vector representing the projection of point on the plane

The returned vector is a 2D vector expressing the projection of the point (from world coordinates) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. The origin point is returned by ProjectionReferencePoint(). All coordinates are in centimeters.

Definition at line 415 of file TPCGeo.h.

References geo::PlaneGeo::PointWidthDepthProjection(), and ReferencePlane().

Referenced by PointProjection(), and VectorProjection().

416 { return ReferencePlane().PointWidthDepthProjection(point); }

WidthDepthProjection_t PointWidthDepthProjection(geo::Point_t const &point) const

Returns the projection of the specified point on the plane.

Definition: PlaneGeo.h:897

geo::PlaneGeo::VectorWidthDepthProjection

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

Projection_t geo::TPCGeo::Projection ( geo::Vector_t const & v ) const

inline

Returns the projection of the specified vector on the plane.

Parameters

v	the 3D vector to be projected, in world units [cm]

Returns: a 2D vector representing the projection of v on the plane

The returned vector is a 2D vector expressing the projection of the vector (from world units) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. All coordinates are in centimeters.

Definition at line 436 of file TPCGeo.h.

References ReferencePlane(), and geo::PlaneGeo::VectorWidthDepthProjection().

437 { return ReferencePlane().VectorWidthDepthProjection(v); }

WidthDepthProjection_t VectorWidthDepthProjection(geo::Vector_t const &v) const

Returns the projection of the specified vector on the plane.

Definition: PlaneGeo.h:917

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

template<typename Point = DefaultPoint_t>

Point geo::TPCGeo::ProjectionReferencePoint ( ) const

inline

Returns the reference point used by PointProjection().

Template Parameters

Point type of point to be returned

The returned point is such that its decomposition results in a null projection and a 0 distance from the plane.

Definition at line 398 of file TPCGeo.h.

References geo::PlaneGeo::GetCenter(), and ReferencePlane().

399 { return ReferencePlane().GetCenter<Point>(); }

dcel2d::Point

std::tuple< double, double, const reco::ClusterHit3D * > Point

Definitions used by the VoronoiDiagram algorithm.

Definition: DCEL.h:34

geo::PlaneGeo::GetCenter

Point GetCenter() const

Returns the centre of the wire plane in world coordinates [cm].

Definition: PlaneGeo.h:426

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

template<typename Vector = DefaultPoint_t>

Vector geo::TPCGeo::RefDepthDir ( ) const

inline

Return the direction of reference plane depth.

Template Parameters

Vector type of vector to return (current default: TVector3)

See also: RefWidthDir(), DriftDir()

The precise definition of the vector is arbitrary, but it is defined to lie on the wire plane and so that RefWidthDir(), RefDepthDir() and a vector opposite to DriftDir() make a orthonormal base. That base (width, depth, normal) is guaranteed to be positive defined.

Definition at line 329 of file TPCGeo.h.

References geo::PlaneGeo::DepthDir(), and ReferencePlane().

Referenced by geo::buildDriftVolumes().

329 { return ReferencePlane().DepthDir<Vector>(); }

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

geo::PlaneGeo::DepthDir

Vector DepthDir() const

Return the direction of plane depth.

Definition: PlaneGeo.h:219

geo::PlaneGeo const& geo::TPCGeo::ReferencePlane ( ) const

inline

Returns the plane used for reference by projection methods.

Definition at line 299 of file TPCGeo.h.

References FirstPlane().

Referenced by ComposePoint(), DecomposePoint(), DistanceFromReferencePlane(), DriftPoint(), Projection(), ProjectionReferencePoint(), RefDepthDir(), ReferencePlaneID(), and RefWidthDir().

299 { return FirstPlane(); }

geo::TPCGeo::FirstPlane

geo::PlaneGeo const & FirstPlane() const

Returns the first wire plane (the closest to TPC center).

Definition: TPCGeo.h:225

geo::PlaneID const& geo::TPCGeo::ReferencePlaneID ( ) const

inline

Returns the ID of the plane used for reference by projection methods.

Definition at line 302 of file TPCGeo.h.

References geo::PlaneGeo::ID(), and ReferencePlane().

303 { return ReferencePlane().ID(); }

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

geo::PlaneGeo::ID

geo::PlaneID const & ID() const

Returns the identifier of this plane.

Definition: PlaneGeo.h:190

template<typename Vector = DefaultPoint_t>

Vector geo::TPCGeo::RefWidthDir ( ) const

inline

Return the direction of reference plane width.

Template Parameters

Vector type of vector to return (current default: TVector3)

See also: RefDepthDir(), DriftDir()

The precise definition of the vector is arbitrary, but it is defined to lie on the wire plane and so that RefWidthDir(), RefDepthDir() and a vector opposite to DriftDir() make a orthonormal base. That base (width, depth, normal) is guaranteed to be positive defined.

Definition at line 316 of file TPCGeo.h.

References ReferencePlane(), and geo::PlaneGeo::WidthDir().

Referenced by geo::buildDriftVolumes().

316 { return ReferencePlane().WidthDir<Vector>(); }

fhicl::detail::atom::value

geo::PlaneGeo const & ReferencePlane() const

Returns the plane used for reference by projection methods.

Definition: TPCGeo.h:299

geo::PlaneGeo::WidthDir

Vector WidthDir() const

Return the direction of plane width.

Definition: PlaneGeo.h:207

void geo::TPCGeo::ResetDriftDirection ( )

private

Recomputes the drift direction; needs planes to have been initialised.

Definition at line 432 of file TPCGeo.cxx.

References DetectDriftDirection(), e, fDriftDir, fDriftDirection, geo::PlaneGeo::GetBoxCenter(), GetCenter(), geo::kNeg, geo::kNegX, geo::kPos, geo::kPosX, geo::kUnknownDrift, Plane(), geo::vect::round01(), geo::Xaxis(), geo::Yaxis(), and geo::Zaxis().

Referenced by TPCGeo().

                                    {
     
     auto const driftDirCode = DetectDriftDirection();
     switch (driftDirCode) {
       case +1:
         fDriftDirection = geo::kPosX; // this is the same as kPos!
         fDriftDir = geo::Xaxis();
         break;
       case -1:
         fDriftDirection = geo::kNegX; // this is the same as kNeg!
         fDriftDir = -geo::Xaxis();
         break;
       case +2:
         fDriftDir = geo::Yaxis();
         fDriftDirection = geo::kPos;
         break;
       case -2:
         fDriftDir = -geo::Yaxis();
         fDriftDirection = geo::kNeg;
         break;
       case +3:
         fDriftDir = geo::Zaxis();
         fDriftDirection = geo::kPos;
         break;
       case -3:
         fDriftDir = -geo::Zaxis();
         fDriftDirection = geo::kNeg;
         break;
       default:
         // TPC ID is likely not yet set
         fDriftDirection = kUnknownDrift;
         
         // we estimate the drift direction roughly from the geometry
         fDriftDir = Plane(0).GetBoxCenter() - GetCenter();
         
         mf::LogError("TPCGeo")
           << "Unable to detect drift direction (result: " << driftDirCode
           << ", drift: ( " << fDriftDir.X() << " ; " << fDriftDir.Y() << " ; "
           << fDriftDir.Z() << " )";
         break;
     } // switch
     
     geo::vect::round01(fDriftDir, 1e-4);
     
   } // TPCGeo::ResetDriftDirection()

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

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(), geo::CryostatGeo::InitCryoBoundaries(), and 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

geo::PlaneGeo const & geo::TPCGeo::SmallestPlane ( ) const

Returns the wire plane with the smallest surface.

Definition at line 311 of file TPCGeo.cxx.

References fPlanes.

Referenced by GetElementPtr().

                                                  {
     
     //
     // Returns the plane with the smallest width x depth. No nonsense here.
     //
     
     auto iPlane = fPlanes.begin(), pend = fPlanes.end();
     auto smallestPlane = iPlane;
     double smallestSurface = smallestPlane->Width() * smallestPlane->Depth();
     while (++iPlane != pend) {
       double const surface = iPlane->Width() * iPlane->Depth();
       if (surface > smallestSurface) continue;
       smallestSurface = surface;
       smallestPlane = iPlane;
     } // while
     return *smallestPlane;
     
   } // TPCGeo::SmallestPlane()

void geo::TPCGeo::SortPlanes ( std::vector< geo::PlaneGeo > & planes ) const

private

Sorts (in place) the specified PlaneGeo objects by drift distance.

Definition at line 535 of file TPCGeo.cxx.

References geo::vect::dot(), GetCenter(), and geo::vect::normalize().

Referenced by SortSubVolumes().

                                                               {
     //
     // Sort planes by increasing drift distance.
     // 
     // This function should work in bootstrap mode, relying on least things as
     // possible. Therefore we compute here a proxy of the drift axis.
     //
     
     //
     // determine the drift axis (or close to): from TPC center to plane center
     // 
     
     // Instead of using the plane center, which might be not available yet,
     // we use the plane box center, which only needs the geometry description
     // to be available.
     // We use the first plane -- it does not make any difference.
     decltype(auto) TPCcenter = GetCenter();
     auto driftAxis
       = geo::vect::normalize(planes[0].GetBoxCenter() - TPCcenter);
     
     //
     // associate each plane with its distance from the center of TPC
     //
     decltype(auto) center = GetCenter();
     // pair: <plane pointer,drift distance>
     std::vector<std::pair<geo::PlaneGeo*, double>> planesWithDistance;
     planesWithDistance.reserve(planes.size());
     for (geo::PlaneGeo& plane: planes) {
       double const driftDistance
         = geo::vect::dot(plane.GetBoxCenter() - TPCcenter, driftAxis);
       planesWithDistance.emplace_back(&plane, driftDistance);
     } // for
     
     //
     // sort by distance
     //
     std::sort(planesWithDistance.begin(), planesWithDistance.end(), 
       [](auto const& a, auto const& b){ return a.second < b.second; }
       );
     
     //
     // extract the result
     //
     std::vector<geo::PlaneGeo> sortedPlanes;
     sortedPlanes.reserve(planesWithDistance.size());
     for (auto const& pair: planesWithDistance)
       sortedPlanes.push_back(std::move(*(pair.first)));
     
     planes = std::move(sortedPlanes);
     
   } // TPCGeo::SortPlanes()

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

Apply sorting to the PlaneGeo objects.

Definition at line 225 of file TPCGeo.cxx.

References fPlane0Pitch, fPlaneLocation, fPlanes, fViewToPlaneNumber, geo::PlaneID::InvalidID, geo::kUnknown, Nplanes(), geo::origin(), and SortPlanes().

Referenced by geo::CryostatGeo::SortSubVolumes(), and toLocalCoords().

   {
     SortPlanes(fPlanes);
     
     double origin[3] = {0.};
  
     // set the plane pitch for this TPC
     double xyz[3]  = {0.};
     fPlanes[0].LocalToWorld(origin,xyz);
     double xyz1[3] = {0.};
     fPlaneLocation.clear();
     fPlaneLocation.resize(fPlanes.size());
     for(unsigned int i = 0; i < fPlaneLocation.size(); ++i) fPlaneLocation[i].resize(3);
     fPlane0Pitch.clear();
     fPlane0Pitch.resize(this->Nplanes(), 0.);
     for(size_t p = 0; p < this->Nplanes(); ++p){
       fPlanes[p].LocalToWorld(origin,xyz1);
       if(p > 0) fPlane0Pitch[p] = fPlane0Pitch[p-1] + std::abs(xyz1[0]-xyz[0]);
       else      fPlane0Pitch[p] = 0.;
       xyz[0] = xyz1[0];
       fPlaneLocation[p][0] = xyz1[0];
       fPlaneLocation[p][1] = xyz1[1];
       fPlaneLocation[p][2] = xyz1[2];
     }
 
     // the PlaneID_t cast convert InvalidID into a rvalue (non-reference);
     // leaving it a reference would cause C++ to treat it as such,
     // that can't be because InvalidID is a static member constant without an address
     // (it is not defined in any translation unit, just declared in header)
     fViewToPlaneNumber.resize
       (1U + (size_t) geo::kUnknown, (geo::PlaneID::PlaneID_t) geo::PlaneID::InvalidID);
     for(size_t p = 0; p < this->Nplanes(); ++p)
       fViewToPlaneNumber[(size_t) fPlanes[p].View()] = p;
 
     for(size_t p = 0; p < fPlanes.size(); ++p) fPlanes[p].SortWires(sorter);
     
   }

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

inline

Transform point from world frame to local TPC frame.

Definition at line 522 of file TPCGeo.h.

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

523 { 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

TPC-to-world transformation.

Definition: TPCGeo.h:630

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

inline

Transform direction vector from world to local.

Definition at line 534 of file TPCGeo.h.

References DetectDriftDirection(), fTrans, art::detail::indent(), PrintTPCInfo(), SortSubVolumes(), geo::LocalTransformationGeo< StoredMatrix, LocalPoint, LocalVector >::toLocalCoords(), and UpdateAfterSorting().

535 { 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

geo::TPCGeo::fTotalVolume

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

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

inline

Half width (associated with x coordinate) of active TPC volume [cm].

Definition at line 112 of file TPCGeo.h.

References fTotalVolume.

112 { return fTotalVolume; }

TGeoVolume * fTotalVolume

Total volume of TPC, called volTPC in GDML file.

Definition: TPCGeo.h:634

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

inline

Transform point from local TPC frame to world frame.

Definition at line 498 of file TPCGeo.h.

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

Referenced by GetCenter(), and InitTPCBoundaries().

499 { return fTrans.toWorldCoords(local); }

geo::LocalTransformationGeo::toWorldCoords

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

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::TPCGeo::toWorldCoords ( LocalVector_t const & local ) const

inline

Transform direction vector from local to world.

Definition at line 510 of file TPCGeo.h.

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

511 { return fTrans.toWorldCoords(local); }

geo::LocalTransformationGeo::toWorldCoords

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

GlobalPoint_t toWorldCoords(LocalPoint_t const &local) const

Transforms a point from local frame to world frame.

Definition: LocalTransformationGeo.h:103

void geo::TPCGeo::UpdateAfterSorting ( geo::TPCID tpcid )

Performs all updates after cryostat has sorted TPCs.

Definition at line 265 of file TPCGeo.cxx.

References DriftDir(), e, fID, fPlanes, lar::util::makeVector3DComparison(), Nplanes(), and UpdatePlaneViewCache().

Referenced by toLocalCoords().

                                                 {
     
     // reset the ID
     fID = tpcid;
     
     // ask the planes to update; also check
     
     for (unsigned int plane = 0; plane < Nplanes(); ++plane) {
       fPlanes[plane].UpdateAfterSorting(geo::PlaneID(fID, plane), *this);
       
       // check that the plane normal is opposite to the TPC drift direction
       assert(lar::util::makeVector3DComparison(1e-5)
         .equal(-(fPlanes[plane].GetNormalDirection()), DriftDir()));
       
     } // for
     
     UpdatePlaneViewCache();
     
   } // TPCGeo::UpdateAfterSorting()

void geo::TPCGeo::UpdatePlaneViewCache ( )

private

Refills the plane vs. view cache of the TPC.

Definition at line 519 of file TPCGeo.cxx.

References fPlanes, fViewToPlaneNumber, geo::PlaneID::InvalidID, geo::kUnknown, and Nplanes().

Referenced by UpdateAfterSorting().

                                     {
     
     // the PlaneID_t cast convert InvalidID into a rvalue (non-reference);
     // leaving it a reference would cause C++ to treat it as such,
     // that can't be because InvalidID is a static member constant without an address
     // (it is not defined in any translation unit, just declared in header)
     fViewToPlaneNumber.clear();
     fViewToPlaneNumber.resize
       (1U + (size_t) geo::kUnknown, (geo::PlaneID::PlaneID_t) geo::PlaneID::InvalidID);
     for(size_t p = 0; p < Nplanes(); ++p)
       fViewToPlaneNumber[(size_t) fPlanes[p].View()] = p;
     
   } // TPCGeo::UpdatePlaneViewCache()

Projection_t geo::TPCGeo::VectorProjection ( geo::Vector_t const & v ) const

inline

Returns the projection of the specified vector on the plane.

Parameters

v	the 3D vector to be projected, in world units [cm]

Returns: a 2D vector representing the projection of v on the plane

The returned vector is a 2D vector expressing the projection of the vector (from world units) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. All coordinates are in centimeters.

Definition at line 438 of file TPCGeo.h.

References Projection().

439 { return Projection(v); }

Projection_t Projection(geo::Point_t const &point) const

Returns the projection of the specified point on the plane.

Definition: TPCGeo.h:415

Projection_t geo::TPCGeo::VectorProjection ( TVector3 const & v ) const

inline

Returns the projection of the specified vector on the plane.

Parameters

v	the 3D vector to be projected, in world units [cm]

Returns: a 2D vector representing the projection of v on the plane

The returned vector is a 2D vector expressing the projection of the vector (from world units) on the reference plane. The vector is expressed as $ ( w, d ) $ , components following the width direction (RefWidthDir()) and the depth direction (RefDepthDir()) respectively. All coordinates are in centimeters.

Definition at line 440 of file TPCGeo.h.

References Projection(), and geo::vect::toVector().

441 { return Projection(geo::vect::toVector(v)); }

geo::vect::toVector

::geo::Vector_t toVector(Vector const &v)

Convert the specified vector into a geo::Vector_t.

Definition: geo_vectors_utils.h:1286

geo::LocalTransformation::WorldToLocal

Projection_t Projection(geo::Point_t const &point) const

Returns the projection of the specified point on the plane.

Definition: TPCGeo.h:415

std::set< geo::View_t > geo::TPCGeo::Views ( ) const

Returns a set of all views covered in this TPC.

Definition at line 343 of file TPCGeo.cxx.

References fPlanes, and geo::PlaneGeo::View().

Referenced by LastPlane().

                                         {
     std::set<geo::View_t> views;
     std::transform(fPlanes.cbegin(), fPlanes.cend(),
       std::inserter(views, views.begin()), std::mem_fn(&geo::PlaneGeo::View));
     return views;
   } // TPCGeo::Views()

double geo::TPCGeo::Width ( ) const

inline

Width is associated with x coordinate [cm].

Definition at line 101 of file TPCGeo.h.

References HalfWidth().

Referenced by PrintTPCInfo().

101 { return 2.0 * HalfWidth(); }

geo::TPCGeo::HalfWidth

double HalfWidth() const

Width is associated with x coordinate [cm].

Definition: TPCGeo.h:99

template<typename Vector = DefaultVector_t>

decltype(auto) geo::TPCGeo::WidthDir ( ) const

inline

Returns the direction Width() is measured on.

Definition at line 116 of file TPCGeo.h.

References fWidthDir.

116 { return geo::vect::convertTo<Vector>(fWidthDir); }

geo::TPCGeo::fWidthDir

geo::Vector_t fWidthDir

Direction width refers to.

Definition: TPCGeo.h:647

double geo::TPCGeo::WirePitch ( unsigned plane = 0 ) const

Returns the center of the TPC volume in world coordinates [cm].

Definition at line 426 of file TPCGeo.cxx.

References Plane(), and geo::PlaneGeo::WirePitch().

Referenced by ActiveBoundingBox(), hit::MagDriftAna::analyze(), geo::ChannelMapStandardAlg::Initialize(), and trkf::SpacePointAlg::update().

   { 
     return this->Plane(plane).WirePitch();
   }

void geo::TPCGeo::WorldToLocal	(	const double *	world,
		double *	tpc
	)		const

inline

Transform point from world frame to local TPC frame.

Definition at line 514 of file TPCGeo.h.

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

515 { fTrans.WorldToLocal(world, tpc); }

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

Transforms a point from world frame to local frame.

geo::LocalTransformation::WorldToLocal

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

TVector3 geo::TPCGeo::WorldToLocal ( TVector3 const & world ) const

inline

Transform point from world frame to local TPC frame.

Definition at line 518 of file TPCGeo.h.

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

519 { return fTrans.WorldToLocal<TVector3>(world); }

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

Transforms a point from world frame to local frame.

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

void geo::TPCGeo::WorldToLocalVect	(	const double *	world,
		double *	tpc
	)		const

inline

Transform direction vector from world to local.

Definition at line 526 of file TPCGeo.h.

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

527 { fTrans.WorldToLocalVect(world, tpc); }

geo::LocalTransformation::WorldToLocalVect

LocalTransformation_t fTrans

TPC-to-world transformation.

Definition: TPCGeo.h:630

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

Transforms a vector from world frame to local frame.

TVector3 geo::TPCGeo::WorldToLocalVect ( TVector3 const & world ) const

inline

Transform direction vector from world to local.

Definition at line 530 of file TPCGeo.h.

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

531 { return fTrans.WorldToLocalVect<TVector3>(world); }