#include "WireReadoutStandardGeom.h"

Inheritance diagram for geo::WireReadoutStandardGeom:

Public Types
using	Segment = std::pair< Point_t, Point_t >

Public Member Functions
	WireReadoutStandardGeom (fhicl::ParameterSet const &pset, GeometryCore const *geom, std::unique_ptr< WireReadoutSorter > sorter)

std::vector< WireID >	ChannelToWire (raw::ChannelID_t channel) const override

unsigned int	Nchannels () const override
	Returns the total number of channels present (not necessarily contiguous) More...

unsigned int	Nchannels (readout::ROPID const &ropid) const override
	Returns the number of channels in the specified ROP. More...

std::set< PlaneID > const &	PlaneIDs () const override
	Returns a list of the plane IDs in the whole detector. More...

unsigned int	MaxPlanes () const
	Returns the largest number of planes among all TPCs in this detector. More...

unsigned int	Nviews () const
	Returns the number of views (different wire orientations) More...

std::set< View_t > const &	Views () const
	Returns a list of possible views in the detector. More...

std::set< View_t >	Views (TPCID const &tpcid) const

void	WireEndPoints (WireID const &wireid, double xyzStart, double xyzEnd) const
	Fills two arrays with the coordinates of the wire end points. More...

double	ThirdPlane_dTdW (PlaneID const &pid1, double slope1, PlaneID const &pid2, double slope2, PlaneID const &output_plane) const
	Returns dT/dW on the third plane, given it in the other two. More...

double	ThirdPlane_dTdW (PlaneID const &pid1, double slope1, PlaneID const &pid2, double slope2) const
	Returns dT/dW on the third plane, given it in the other two. More...

double	Plane0Pitch (TPCID const &tpcid, unsigned int p1) const

double	PlanePitch (TPCID const &tpcid, unsigned int p1=0, unsigned int p2=1) const


double	WireCoordinate (double YPos, double ZPos, PlaneID const &planeID) const override
	Returns the index of the wire nearest to the specified position. More...


WireID	NearestWireID (Point_t const &worldPos, PlaneID const &planeID) const override
	Returns the ID of the wire nearest to the specified position. More...


raw::ChannelID_t	PlaneWireToChannel (WireID const &wireID) const override
	Returns the channel ID a wire is connected to. More...

TPC set mapping
unsigned int	NTPCsets (readout::CryostatID const &cryoid) const override
	Returns the total number of TPC sets in the specified cryostat. More...

unsigned int	MaxTPCsets () const override
	Returns the largest number of TPC sets any cryostat in the detector has. More...

bool	HasTPCset (readout::TPCsetID const &tpcsetid) const override

readout::TPCsetID	TPCtoTPCset (TPCID const &tpcid) const override
	Returns the ID of the TPC set the specified TPC belongs to. More...

std::vector< TPCID >	TPCsetToTPCs (readout::TPCsetID const &tpcsetid) const override
	Returns a list of ID of TPCs belonging to the specified TPC set. More...

TPCID	FirstTPCinTPCset (readout::TPCsetID const &tpcsetid) const override
	Returns the ID of the first TPC belonging to the specified TPC set. More...

Readout plane mapping
unsigned int	NROPs (readout::TPCsetID const &tpcsetid) const override
	Returns the total number of ROPs in the specified TPC set. More...

unsigned int	MaxROPs () const override
	Returns the largest number of ROPs a TPC set in the detector has. More...

bool	HasROP (readout::ROPID const &ropid) const override

readout::ROPID	WirePlaneToROP (PlaneID const &planeid) const override
	Returns the ID of the ROP planeid belongs to, or invalid if none. More...

std::vector< PlaneID >	ROPtoWirePlanes (readout::ROPID const &ropid) const override
	Returns a list of ID of wire planes belonging to the specified ROP. More...

std::vector< TPCID >	ROPtoTPCs (readout::ROPID const &ropid) const override
	Returns a list of ID of TPCs the specified ROP spans. More...

readout::ROPID	ChannelToROP (raw::ChannelID_t channel) const override
	Returns the ID of the ROP the channel belongs to (invalid if none) More...

raw::ChannelID_t	FirstChannelInROP (readout::ROPID const &ropid) const override
	Returns the ID of the first channel in the specified readout plane. More...

PlaneID	FirstWirePlaneInROP (readout::ROPID const &ropid) const override
	Returns the ID of the first plane belonging to the specified ROP. More...


bool	HasPlane (PlaneID const &planeid) const
	Returns whether we have the specified plane. More...

bool	HasElement (PlaneID const &planeid) const
	Returns whether we have the specified plane. More...


bool	HasElement (WireID const &wireid) const
	Returns whether we have the specified wire. More...

bool	HasWire (WireID const &wireid) const
	Returns whether we have the specified wire. More...


unsigned int	Nplanes (TPCID const &tpcid=details::tpc_zero) const
	Returns the total number of planes in the specified TPC. More...

unsigned int	NElements (TPCID const &tpcid) const
	Returns the total number of planes in the specified TPC. More...

unsigned int	NSiblingElements (PlaneID const &planeid) const
	Returns the total number of planes in the specified TPC. More...

Wire details
unsigned int	NElements (PlaneID const &planeid) const
	Returns the total number of wires in the specified plane. More...

unsigned int	NSiblingElements (WireID const &wireid) const
	Returns the total number of wires in the specified plane. More...

unsigned int	MaxWires () const
	Returns the total number of wires in the specified plane. More...

unsigned int	Nwires (PlaneID const &planeid) const
	Returns the total number of wires in the specified plane. More...


PlaneGeo const &	FirstPlane (TPCID const &tpcid) const
	Returns the first plane of the specified TPC. More...


PlaneGeo const &	Plane (TPCID const &tpcid, View_t view) const
	Returns the specified wire. More...

PlaneGeo const &	Plane (PlaneID const &planeid) const
	Returns the specified wire. More...

PlaneGeo const &	GetElement (PlaneID const &planeid) const
	Returns the specified wire. More...


WireGeo const &	GetElement (WireID const &wireid) const
	Returns the specified wire. More...

WireGeo const &	Wire (WireID const &wireid) const
	Returns the specified wire. More...


PlaneGeo const *	PlanePtr (PlaneID const &planeid) const
	Returns the specified plane. More...

PlaneGeo const *	GetElementPtr (PlaneID const &planeid) const
	Returns the specified plane. More...


WireGeo const *	GetElementPtr (WireID const &wireid) const
	Returns the specified wire. More...

WireGeo const *	WirePtr (WireID const &wireid) const
	Returns the specified wire. More...


Segment	WireEndPoints (WireID const &wireid) const
	Returns a segment whose ends are the wire end points. More...

Wire access and information
double	WireAngleToVertical (View_t view, TPCID const &tpcid) const
	Returns the angle of the wires in the specified view from vertical. More...


bool	WireIDsIntersect (WireID const &wid1, WireID const &wid2, Point_t &intersection) const
	Computes the intersection between two wires. More...


std::optional< WireIDIntersection >	WireIDsIntersect (WireID const &wid1, WireID const &wid2) const
	Computes the intersection between two wires. More...

Wire geometry queries
Please note the differences between functions: ChannelsIntersect(), WireIDsIntersect() and IntersectionPoint() all calculate wires intersection using the same equation. ChannelsIntersect() and WireIdsIntersect() will return true if the two wires cross, return false if they don't. IntersectionPoint() does not check if the two wires cross.
PlaneID	ThirdPlane (PlaneID const &pid1, PlaneID const &pid2) const
	Returns the plane that is not in the specified arguments. More...

double	ThirdPlaneSlope (PlaneID const &pid1, double slope1, PlaneID const &pid2, double slope2, PlaneID const &output_plane) const
	Returns the slope on the third plane, given it in the other two. More...

double	ThirdPlaneSlope (PlaneID const &pid1, double slope1, PlaneID const &pid2, double slope2) const
	Returns the slope on the third plane, given it in the other two. More...

double	ThirdPlaneSlope (PlaneID::PlaneID_t plane1, double slope1, PlaneID::PlaneID_t plane2, double slope2, TPCID const &tpcid) const
	Returns the slope on the third plane, given it in the other two. More...

TPC channel mapping
std::vector< raw::ChannelID_t >	ChannelsInTPCs () const
	Returns an std::vector<ChannelID_t> in all TPCs in a TPCSet. More...

virtual bool	HasChannel (raw::ChannelID_t channel) const
	Returns whether the specified channel is valid This default implementation assumes all channels up to Nchannels() valid. More...

View_t	View (raw::ChannelID_t const channel) const
	Returns the view (wire orientation) on the specified TPC channel. More...

View_t	View (readout::ROPID const &ropid) const
	Returns the view of the channels in the specified readout plane. More...

SigType_t	SignalType (PlaneID const &pid) const
	Returns the type of signal on the channels of specified TPC plane. More...

SigType_t	SignalType (raw::ChannelID_t const channel) const
	Return the signal type of the specified channel. More...

SigType_t	SignalType (readout::ROPID const &ropid) const
	Return the signal type on the specified readout plane. More...

raw::ChannelID_t	NearestChannel (Point_t const &worldLoc, PlaneID const &planeid) const
	Returns the ID of the channel nearest to the specified position. More...

std::optional< WireIDIntersection >	ChannelsIntersect (raw::ChannelID_t c1, raw::ChannelID_t c2) const
	Returns an intersection point of two channels. More...

Optical detector channel mapping
virtual unsigned int	NOpChannels (unsigned int NOpDets) const
	Returns the number of optical channels contained in some detectors. More...

unsigned int	NOpChannels () const
	Number of electronics channels for all the optical detectors. More...

unsigned int	MaxOpChannel () const
	Largest optical channel number. More...

virtual unsigned int	MaxOpChannel (unsigned int NOpDets) const
	Returns the number of optical channels contained in some detectors. More...

bool	IsValidOpChannel (int opChannel) const
	Is this a valid OpChannel number? More...

virtual bool	IsValidOpChannel (unsigned int opChannel, unsigned int NOpDets) const
	Returns whether the ID identifies a valid optical detector channel. More...

virtual unsigned int	NOpHardwareChannels (unsigned int opDet) const
	Returns the number of channels in the specified optical detectors. More...

virtual unsigned int	OpChannel (unsigned int detNum, unsigned int hwchannel=0) const
	Returns the channel ID of the specified hardware channel. More...

virtual unsigned int	OpDetFromOpChannel (unsigned int opChannel) const
	Returns the optical detector the specified optical channel belongs. More...

OpDetGeo const &	OpDetGeoFromOpChannel (unsigned int opChannel) const
	Returns the optical detector the specified optical channel belongs. More...

virtual unsigned int	HardwareChannelFromOpChannel (unsigned int opChannel) const
	Returns the hardware channel number of specified optical channel. More...

TPC set mapping
readout::TPCsetID	FindTPCsetAtPosition (Point_t const &worldLoc) const
	Returns the total number of TPC sets in the specified cryostat. More...

Testing (not in the interface)
std::vector< std::vector< std::vector< raw::ChannelID_t > > > const &	FirstChannelInNextPlane () const
	Retrieve the private fFirstChannelInNextPlane vector for testing. More...

std::vector< std::vector< std::vector< raw::ChannelID_t > > > const &	FirstChannelInThisPlane () const
	Retrieve the private fFirstChannelInThisPlane vector for testing. More...

Static Public Member Functions
static auto const &	start (Segment const &s)

static auto const &	finish (Segment const &s)

static double	ComputeThirdPlaneSlope (double angle1, double slope1, double angle2, double slope2, double angle_target)
	Returns the slope on the third plane, given it in the other two. More...

static double	ComputeThirdPlane_dTdW (double angle1, double pitch1, double dTdW1, double angle2, double pitch2, double dTdW2, double angle_target, double pitch_target)
	Returns the slope on the third plane, given it in the other two. More...

Protected Types
template<typename T >
using	TPCInfoMap_t = std::vector< std::vector< T >>
	Data type for per-TPC information. More...

template<typename T >
using	PlaneInfoMap_t = TPCInfoMap_t< std::vector< T >>
	Data type for per-plane information. More...

Protected Member Functions
template<typename T >
T const *	GetElementPtr (PlaneInfoMap_t< T > const &map, PlaneID const &id) const
	Returns a pointer to the specified element, or nullptr if invalid. More...

virtual SigType_t	SignalTypeForROPIDImpl (readout::ROPID const &ropid) const
	Return the signal type on the specified readout plane. More...

template<typename T >
T const &	AccessElement (PlaneInfoMap_t< T > const &map, PlaneID const &id) const
	Returns the specified element of the plane map. More...

template<typename T >
size_t	AccessElementSize (PlaneInfoMap_t< T > const &map, TPCID const &id) const
	Returns the number of elements in the specified TPC of the plane map. More...

Internal structure data access
These functions allow access to the XxxInfoMap_t types based on geometry element IDs. They are strictly internal.
template<typename T >
T const &	AccessElement (TPCInfoMap_t< T > const &map, TPCID const &id) const
	Returns the specified element of the TPC map. More...

template<typename T >
size_t	AccessElementSize (TPCInfoMap_t< T > const &map, CryostatID const &id) const
	Returns the number of elements in the specified cryostat of the TPC map. More...

template<typename T >
bool	isValidElement (TPCInfoMap_t< T > const &map, CryostatID const &id) const
	Returns whether the ID specifies a valid entry. More...

template<typename T >
bool	isValidElement (TPCInfoMap_t< T > const &map, TPCID const &id) const
	Returns the specified element of the TPC map. More...


template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, CryostatID const &id) const
	Returns whether the ID specifies a valid entry. More...

template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, TPCID const &id) const
	Returns whether the ID specifies a valid entry. More...

template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, PlaneID const &id) const
	Returns whether the ID specifies a valid entry. More...

Protected Attributes
PlaneInfoMap_t< raw::ChannelID_t >	fFirstChannelInThisPlane

PlaneInfoMap_t< raw::ChannelID_t >	fFirstChannelInNextPlane

Private Member Functions
SigType_t	SignalTypeForChannelImpl (raw::ChannelID_t const channel) const override
	Return the signal type of the specified channel. More...

unsigned int	WireCount (PlaneID const &id) const
	Retrieved the wire cound for the specified plane ID. More...

unsigned int	MaxTPCs () const
	Returns the largest number of TPCs in a single cryostat. More...

Private Attributes
unsigned int	fNcryostat
	number of cryostats in the detector More...

unsigned int	fNchannels
	number of channels in the detector More...

raw::ChannelID_t	fTopChannel
	book keeping highest channel # More...

std::vector< unsigned int >	fNTPC
	number of TPCs in each cryostat More...

std::set< View_t >	fViews
	vector of the views present in the detector More...

std::set< PlaneID >	fPlaneIDs
	vector of the PlaneIDs present in the detector More...

PlaneInfoMap_t< float >	fFirstWireProj

PlaneInfoMap_t< float >	fOrthVectorsY
	Unit vectors orthogonal to wires in. More...

PlaneInfoMap_t< float >	fOrthVectorsZ

PlaneInfoMap_t< float >	fWireCounts

TPCInfoMap_t< unsigned int >	fNPlanes

PlaneInfoMap_t< unsigned int >	fPlaneBaselines

PlaneInfoMap_t< unsigned int >	fWiresPerPlane

Detailed Description

Definition at line 21 of file WireReadoutStandardGeom.h.

Member Typedef Documentation

template<typename T >

using geo::WireReadoutGeom::PlaneInfoMap_t = TPCInfoMap_t<std::vector<T>>

protectedinherited

Data type for per-plane information.

Definition at line 982 of file WireReadoutGeom.h.

using geo::WireReadoutGeom::Segment = std::pair<Point_t, Point_t>

inherited

Definition at line 62 of file WireReadoutGeom.h.

template<typename T >

using geo::WireReadoutGeom::TPCInfoMap_t = std::vector<std::vector<T>>

protectedinherited

Data type for per-TPC information.

Definition at line 978 of file WireReadoutGeom.h.

Constructor & Destructor Documentation

geo::WireReadoutStandardGeom::WireReadoutStandardGeom	(	fhicl::ParameterSet const &	pset,
		GeometryCore const *	geom,
		std::unique_ptr< WireReadoutSorter >	sorter
	)

Definition at line 57 of file WireReadoutStandardGeom.cxx.

References geo::WireGeo::CosThetaZ(), geo::WireReadoutGeom::fFirstChannelInNextPlane, geo::WireReadoutGeom::fFirstChannelInThisPlane, fFirstWireProj, fNchannels, fNcryostat, fNPlanes, fNTPC, fOrthVectorsY, fOrthVectorsZ, fPlaneBaselines, fPlaneIDs, fTopChannel, fWireCounts, fWiresPerPlane, fhicl::ParameterSet::get(), geo::WireGeo::GetCenter(), MF_LOG_DEBUG, geo::WireReadoutGeom::Nplanes(), and geo::WireGeo::SinThetaZ().

     : WireReadoutGeom{geom,
                       std::make_unique<WireReadoutGeomBuilderStandard>(
                         pset.get<fhicl::ParameterSet>("Builder", {})),
                       std::move(sorter)}
   {
     fNcryostat = geom->Ncryostats();
 
     mf::LogInfo("WireReadoutStandardGeom") << "Initializing standard wire-readout geometry...";
 
     fNTPC.resize(fNcryostat);
     fWireCounts.resize(fNcryostat);
     fNPlanes.resize(fNcryostat);
     fFirstWireProj.resize(fNcryostat);
     fOrthVectorsY.resize(fNcryostat);
     fOrthVectorsZ.resize(fNcryostat);
     fPlaneBaselines.resize(fNcryostat);
     fWiresPerPlane.resize(fNcryostat);
     fFirstChannelInNextPlane.resize(fNcryostat);
     fFirstChannelInThisPlane.resize(fNcryostat);
     fPlaneIDs.clear();
     fTopChannel = 0;
 
     int RunningTotal = 0;
 
     for (auto const& cryo : geom->Iterate<CryostatGeo>()) {
       auto const cs = cryo.ID().Cryostat;
       fNTPC[cs] = cryo.NTPC();
 
       // Size up all the vectors
       fWireCounts[cs].resize(fNTPC[cs]);
       fFirstWireProj[cs].resize(fNTPC[cs]);
       fOrthVectorsY[cs].resize(fNTPC[cs]);
       fOrthVectorsZ[cs].resize(fNTPC[cs]);
       fPlaneBaselines[cs].resize(fNTPC[cs]);
       fWiresPerPlane[cs].resize(fNTPC[cs]);
       fNPlanes[cs].resize(fNTPC[cs]);
       fFirstChannelInThisPlane[cs].resize(fNTPC[cs]);
       fFirstChannelInNextPlane[cs].resize(fNTPC[cs]);
 
       for (auto const& TPC : geom->Iterate<TPCGeo>(cryo.ID())) {
         auto const TPCCount = TPC.ID().TPC;
         unsigned int PlanesThisTPC = Nplanes(TPC.ID());
         fWireCounts[cs][TPCCount].resize(PlanesThisTPC);
         fFirstWireProj[cs][TPCCount].resize(PlanesThisTPC);
         fOrthVectorsY[cs][TPCCount].resize(PlanesThisTPC);
         fOrthVectorsZ[cs][TPCCount].resize(PlanesThisTPC);
         fNPlanes[cs][TPCCount] = PlanesThisTPC;
         for (auto const& plane : Iterate<PlaneGeo>(TPC.ID())) {
           fPlaneIDs.insert(plane.ID());
           auto const PlaneCount = plane.ID().Plane;
 
           double ThisWirePitch = plane.WirePitch();
           fWireCounts[cs][TPCCount][PlaneCount] = plane.Nwires();
 
           WireGeo const& firstWire = plane.Wire(0);
           double const sth = firstWire.SinThetaZ(), cth = firstWire.CosThetaZ();
 
           auto WireCenter1 = firstWire.GetCenter();
           auto WireCenter2 = plane.Wire(1).GetCenter();
 
           // figure out if we need to flip the orthogonal vector
           // (should point from wire n -> n+1)
           double OrthY = cth, OrthZ = -sth;
           if (((WireCenter2.Y() - WireCenter1.Y()) * OrthY +
                (WireCenter2.Z() - WireCenter1.Z()) * OrthZ) < 0) {
             OrthZ *= -1;
             OrthY *= -1;
           }
 
           // Overall we are trying to build an expression that looks like
           //  int NearestWireNumber = round((worldPos.OrthVector - FirstWire.OrthVector)/WirePitch);
           // That runs as fast as humanly possible.
           // We predivide everything by the wire pitch so we don't do this in the loop.
           //
           // Putting this together into the useful constants we will use later per plane and tpc:
           fOrthVectorsY[cs][TPCCount][PlaneCount] = OrthY / ThisWirePitch;
           fOrthVectorsZ[cs][TPCCount][PlaneCount] = OrthZ / ThisWirePitch;
 
           fFirstWireProj[cs][TPCCount][PlaneCount] =
             WireCenter1.Y() * OrthY + WireCenter1.Z() * OrthZ;
           fFirstWireProj[cs][TPCCount][PlaneCount] /= ThisWirePitch;
 
           // now to count up wires in each plane and get first channel in each plane
           int WiresThisPlane = plane.Nwires();
           fWiresPerPlane[cs].at(TPCCount).push_back(WiresThisPlane);
           fPlaneBaselines[cs].at(TPCCount).push_back(RunningTotal);
 
           RunningTotal += WiresThisPlane;
 
           fFirstChannelInThisPlane[cs].at(TPCCount).push_back(fTopChannel);
           fTopChannel += WiresThisPlane;
           fFirstChannelInNextPlane[cs].at(TPCCount).push_back(fTopChannel);
 
         } // end loop over planes
       }   // end loop over TPCs
     }     // end loop over cryostats
 
     // calculate the total number of channels in the detector
     fNchannels = fTopChannel;
 
     MF_LOG_DEBUG("ChannelMapStandard") << "# of channels is " << fNchannels;
   }

Member Function Documentation

template<typename T >

T const& geo::WireReadoutGeom::AccessElement	(	TPCInfoMap_t< T > const &	map,
		TPCID const &	id
	)		const

inlineprotectedinherited

Returns the specified element of the TPC map.

Definition at line 994 of file WireReadoutGeom.h.

Referenced by FirstChannelInROP(), NROPs(), WireCoordinate(), and WireCount().

     {
       return map[id.Cryostat][id.TPC];
     }

template<typename T >

T const& geo::WireReadoutGeom::AccessElement	(	PlaneInfoMap_t< T > const &	map,
		PlaneID const &	id
	)		const

inlineprotectedinherited

Returns the specified element of the plane map.

Definition at line 1022 of file WireReadoutGeom.h.

     {
       return map[id.Cryostat][id.TPC][id.Plane];
     }

template<typename T >

size_t geo::WireReadoutGeom::AccessElementSize	(	TPCInfoMap_t< T > const &	map,
		CryostatID const &	id
	)		const

inlineprotectedinherited

Returns the number of elements in the specified cryostat of the TPC map.

Definition at line 1001 of file WireReadoutGeom.h.

     {
       return map[id.Cryostat].size();
     }

template<typename T >

size_t geo::WireReadoutGeom::AccessElementSize	(	PlaneInfoMap_t< T > const &	map,
		TPCID const &	id
	)		const

inlineprotectedinherited

Returns the number of elements in the specified TPC of the plane map.

Definition at line 1029 of file WireReadoutGeom.h.

     {
       return map[id.Cryostat][id.TPC].size();
     }

std::optional< WireIDIntersection > geo::WireReadoutGeom::ChannelsIntersect	(	raw::ChannelID_t	c1,
		raw::ChannelID_t	c2
	)		const

inherited

Returns an intersection point of two channels.

Parameters

c1	one channel ID
c2	the other channel ID

Returns: WireIDIntersection object if intersection point was found; otherwise std::nullopt

Todo:

what happens for channels from different TPCs?

what happens for channels with multiple intersection points?

Deprecated:: This is clearly not APA-aware

Definition at line 543 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ChannelToWire(), and geo::WireReadoutGeom::WireIDsIntersect().

Referenced by util::GeometryUtilities::GetProjectedPoint(), util::GeometryUtilities::GetYZ(), and geo::WireReadoutGeom::HasChannel().

   {
     // [GP] these errors should be exceptions, and this function is deprecated because it
     // violates interoperability
     std::vector<WireID> chan1wires = ChannelToWire(c1);
     if (chan1wires.empty()) {
       mf::LogError("ChannelsIntersect")
         << "1st channel " << c1 << " maps to no wire (is it a real one?)";
       return std::nullopt;
     }
     std::vector<WireID> chan2wires = ChannelToWire(c2);
     if (chan2wires.empty()) {
       mf::LogError("ChannelsIntersect")
         << "2nd channel " << c2 << " maps to no wire (is it a real one?)";
       return std::nullopt;
     }
 
     if (chan1wires.size() > 1) {
       mf::LogWarning("ChannelsIntersect")
         << "1st channel " << c1 << " maps to " << chan2wires.size() << " wires; using the first!";
       return std::nullopt;
     }
     if (chan2wires.size() > 1) {
       mf::LogError("ChannelsIntersect")
         << "2nd channel " << c2 << " maps to " << chan2wires.size() << " wires; using the first!";
       return std::nullopt;
     }
 
     return WireIDsIntersect(chan1wires[0], chan2wires[0]);
   }

std::vector< raw::ChannelID_t > geo::WireReadoutGeom::ChannelsInTPCs ( ) const

inherited

Returns an std::vector<ChannelID_t> in all TPCs in a TPCSet.

Definition at line 462 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::Nchannels(), geo::WireReadoutGeom::PlaneWireToChannel(), and geo::WireReadoutGeom::TPCsetToTPCs().

Referenced by geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     std::vector<raw::ChannelID_t> channels;
     channels.reserve(Nchannels());
 
     for (auto const& ts : Iterate<readout::TPCsetID>()) {
       for (auto const t : TPCsetToTPCs(ts)) {
         for (auto const& wire : Iterate<WireID>(t)) {
           channels.push_back(PlaneWireToChannel(wire));
         }
       }
     }
     std::sort(channels.begin(), channels.end());
     auto last = std::unique(channels.begin(), channels.end());
     channels.erase(last, channels.end());
     return channels;
   }

readout::ROPID geo::WireReadoutStandardGeom::ChannelToROP ( raw::ChannelID_t channel ) const

overridevirtual

Returns the ID of the ROP the channel belongs to (invalid if none)

Implements geo::WireReadoutGeom.

Definition at line 423 of file WireReadoutStandardGeom.cxx.

References ChannelToWire(), raw::isValidChannelID(), and WirePlaneToROP().

   {
     if (!raw::isValidChannelID(channel)) return {}; // invalid ROP returned
 
     // which wires does the channel cover?
     std::vector<WireID> wires = ChannelToWire(channel);
 
     // - none:
     if (wires.empty()) return {}; // invalid ROP returned
 
     // - one: maps its plane ID into a ROP ID
     return WirePlaneToROP(wires[0]);
   }

std::vector< WireID > geo::WireReadoutStandardGeom::ChannelToWire ( raw::ChannelID_t channel ) const

overridevirtual

Returns a list of TPC wires connected to the specified readout channel ID

Exceptions

cet::exception (category: "Geometry") if non-existent channel

Implements geo::WireReadoutGeom.

Definition at line 164 of file WireReadoutStandardGeom.cxx.

References geo::WireReadoutGeom::fFirstChannelInNextPlane, geo::WireReadoutGeom::fFirstChannelInThisPlane, fNcryostat, fNTPC, and fTopChannel.

Referenced by ChannelToROP().

   {
     std::vector<WireID> AllSegments;
     unsigned int cstat = 0;
     unsigned int tpc = 0;
     unsigned int plane = 0;
     unsigned int wire = 0;
 
     // first check if this channel ID is legal
     if (channel > fTopChannel)
       throw cet::exception("Geometry") << "ILLEGAL CHANNEL ID for channel " << channel << "\n";
 
     // then go find which plane, tpc and cryostat it is in from the information we stored earlier
     bool foundWid(false);
     for (unsigned int csloop = 0; csloop != fNcryostat; ++csloop) {
       for (unsigned int tpcloop = 0; tpcloop != fNTPC[csloop]; ++tpcloop) {
         for (unsigned int planeloop = 0;
              planeloop != fFirstChannelInNextPlane[csloop][tpcloop].size();
              ++planeloop) {
           if (channel < fFirstChannelInNextPlane[csloop][tpcloop][planeloop]) {
             cstat = csloop;
             tpc = tpcloop;
             plane = planeloop;
             wire = channel - fFirstChannelInThisPlane[cstat][tpcloop][planeloop];
             foundWid = true;
             break;
           }
           if (foundWid) break;
         } // end plane loop
         if (foundWid) break;
       } // end tpc loop
       if (foundWid) break;
     } // end cryostat loop
 
     WireID CodeWire(cstat, tpc, plane, wire);
 
     AllSegments.push_back(CodeWire);
 
     return AllSegments;
   }

double geo::WireReadoutGeom::ComputeThirdPlane_dTdW	(	double	angle1,
		double	pitch1,
		double	dTdW1,
		double	angle2,
		double	pitch2,
		double	dTdW2,
		double	angle_target,
		double	pitch_target
	)

staticinherited

Returns the slope on the third plane, given it in the other two.

Parameters

angle1	angle or the wires on the first plane
pitch1	wire pitch on the first plane
dTdW1	slope in dt/dw units as observed on the first plane
angle2	angle or the wires on the second plane
pitch2	wire pitch on the second plane
dTdW2	slope in dt/dw units as observed on the second plane
angle_target	angle or the wires on the target plane
pitch_target	wire pitch on the target plane

Returns: dt/dw slope as measured on the third plane, or 999 if infinity

The input slope must be specified in dt/dw non-homogeneous coordinates.

This function will return a small slope if both input slopes are small.

Definition at line 730 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ComputeThirdPlaneSlope().

Referenced by geo::WireReadoutGeom::ThirdPlane_dTdW(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     // we need to convert dt/dw into homogeneous coordinates, and then back;
     // slope = [dT * (TDCperiod / driftVelocity)] / [dW * wirePitch]
     // The coefficient of dT is assumed to be the same for all the planes, and it finally
     // cancels out. Pitches cancel out only if they are all the same.
     return pitch_target *
            ComputeThirdPlaneSlope(angle1, dTdW1 / pitch1, angle2, dTdW2 / pitch2, angle_target);
   }

double geo::WireReadoutGeom::ComputeThirdPlaneSlope	(	double	angle1,
		double	slope1,
		double	angle2,
		double	slope2,
		double	angle_target
	)

staticinherited

Returns the slope on the third plane, given it in the other two.

Parameters

angle1	angle or the wires on the first plane
slope1	slope as observed on the first plane
angle2	angle or the wires on the second plane
slope2	slope as observed on the second plane
angle_target	angle or the wires on the target plane

Returns: the slope as measure on the third plane, or 999 if infinity

This function will return a small slope if both input slopes are small.

Definition at line 702 of file WireReadoutGeom.cxx.

References util::abs().

Referenced by geo::WireReadoutGeom::ComputeThirdPlane_dTdW(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     // note that, if needed, the trigonometric functions can be pre-calculated.
 
     // Can't resolve very small slopes
     if ((std::abs(slope1) < 0.001) && (std::abs(slope2)) < 0.001) return 0.001;
 
     // We need the "wire coordinate direction" for each plane.  This is perpendicular to
     // the wire orientation.
     double slope3 = 0.001;
     if (std::abs(slope1) > 0.001 && std::abs(slope2) > 0.001) {
       slope3 =
         (+(1. / slope1) * std::sin(angle3 - angle2) - (1. / slope2) * std::sin(angle3 - angle1)) /
         std::sin(angle1 - angle2);
     }
     if (slope3 != 0.)
       slope3 = 1. / slope3;
     else
       slope3 = 999.;
 
     return slope3;
   }

readout::TPCsetID geo::WireReadoutGeom::FindTPCsetAtPosition ( Point_t const & worldLoc ) const

inherited

Returns the total number of TPC sets in the specified cryostat.

Parameters

cryoid cryostat ID

Returns: number of TPC sets in the cryostat, or 0 if no cryostat found

Definition at line 576 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fGeom, geo::GeometryCore::FindTPCAtPosition(), and geo::WireReadoutGeom::TPCtoTPCset().

Referenced by geo::WireReadoutGeom::HasChannel().

   {
     return TPCtoTPCset(fGeom->FindTPCAtPosition(worldLoc));
   }

static auto const& geo::WireReadoutGeom::finish ( Segment const & s )

inlinestaticinherited

Definition at line 64 of file WireReadoutGeom.h.

References geo::Iterable< details::ReadoutIterationPolicy, details::ToReadoutGeometryElement >::begin(), geo::Iterable< details::ReadoutIterationPolicy, details::ToReadoutGeometryElement >::end(), geo::WireReadoutGeom::HasPlane(), and geo::Iterable< details::ReadoutIterationPolicy, details::ToReadoutGeometryElement >::Iterate().

Referenced by geo::WireReadoutGeom::WireIDsIntersect().

64 { return s.second; }

std::vector<std::vector<std::vector<raw::ChannelID_t> > > const& geo::WireReadoutGeom::FirstChannelInNextPlane ( ) const

inlineinherited

Retrieve the private fFirstChannelInNextPlane vector for testing.

Definition at line 936 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::fFirstChannelInThisPlane.

     {
       return fFirstChannelInThisPlane;
     }

raw::ChannelID_t geo::WireReadoutStandardGeom::FirstChannelInROP ( readout::ROPID const & ropid ) const

overridevirtual

Returns the ID of the first channel in the specified readout plane.

Parameters

ropid ID of the readout plane

Returns: ID of first channel, or raw::InvalidChannelID if ID is invalid

Note that this check is performed on the validity of the readout plane ID, that does not necessarily imply that the readout plane specified by the ID actually exists. Check if the ROP exists with HasROP(). The behaviour for non-existing readout planes is undefined.

Implements geo::WireReadoutGeom.

Definition at line 438 of file WireReadoutStandardGeom.cxx.

References geo::WireReadoutGeom::AccessElement(), fPlaneBaselines, raw::InvalidChannelID, and geo::CryostatID::isValid.

   {
     if (!ropid.isValid) return raw::InvalidChannelID;
     return (raw::ChannelID_t)AccessElement(fPlaneBaselines, ConvertROPtoWirePlane(ropid));
   }

std::vector<std::vector<std::vector<raw::ChannelID_t> > > const& geo::WireReadoutGeom::FirstChannelInThisPlane ( ) const

inlineinherited

Retrieve the private fFirstChannelInThisPlane vector for testing.

Definition at line 942 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::fFirstChannelInNextPlane, geo::WireReadoutGeom::SignalTypeForChannelImpl(), and geo::WireReadoutGeom::SignalTypeForROPIDImpl().

     {
       return fFirstChannelInNextPlane;
     }

PlaneGeo const & geo::WireReadoutGeom::FirstPlane ( TPCID const & tpcid ) const

inherited

Returns the first plane of the specified TPC.

Parameters

tpc	ID of the TPC

Returns: a constant reference to the specified plane

Exceptions

cet::exception (WireReadoutGeom category) if plane not present

Definition at line 202 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::PlanePtr().

Referenced by larg4::LArVoxelReadout::DriftIonizationElectrons(), trkf::TrackLineFitAlg::TrkLineFit(), DUNE::NeutrinoTrackingEff::truthLength(), and geo::WireReadoutGeom::Views().

   {
     if (auto const* plane_ptr = PlanePtr({tpcid, 0})) { return *plane_ptr; }
     throw cet::exception("WireReadoutGeom")
       << "TPC with ID " << tpcid << " does not have a first plane.";
   }

TPCID geo::WireReadoutStandardGeom::FirstTPCinTPCset ( readout::TPCsetID const & tpcsetid ) const

overridevirtual

Returns the ID of the first TPC belonging to the specified TPC set.

Implements geo::WireReadoutGeom.

Definition at line 362 of file WireReadoutStandardGeom.cxx.

Referenced by NROPs().

   {
     return ConvertTPCsetToTPC(tpcsetid);
   }

PlaneID geo::WireReadoutStandardGeom::FirstWirePlaneInROP ( readout::ROPID const & ropid ) const

overridevirtual

Returns the ID of the first plane belonging to the specified ROP.

Implements geo::WireReadoutGeom.

Definition at line 445 of file WireReadoutStandardGeom.cxx.

Referenced by Nchannels(), and ROPtoWirePlanes().

   {
     return ConvertROPtoWirePlane(ropid);
   }

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

inlineinherited

Returns the specified wire.

Parameters

tpcid	ID of the TPC
view	number of the requested view
planeid	ID of the plane

Returns: a constant reference to the specified plane

Exceptions

cet::exception	(`GeometryCore` category) if cryostat not present
cet::exception	(`TPCOutOfRange` category) if no such TPC
cet::exception	(`PlaneOutOfRange` category) if no such plane

The GetElement() method is overloaded and its return depends on the type of ID.

Definition at line 152 of file WireReadoutGeom.h.

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

152 { return Plane(planeid); }

geo::WireReadoutGeom::Plane

PlaneGeo const & Plane(TPCID const &tpcid, View_t view) const

Returns the specified wire.

Definition: WireReadoutGeom.cxx:210

WireGeo const& geo::WireReadoutGeom::GetElement ( WireID const & wireid ) const

inlineinherited

Returns the specified wire.

Parameters

wireid ID of the wire

Returns: a constant reference to the specified wire

Exceptions

cet::exception if not found

The GetElement() method is overloaded and its return depends on the type of ID.

Definition at line 226 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::Wire(), and geo::WireReadoutGeom::WireEndPoints().

226 { return Wire(wireid); }

geo::WireReadoutGeom::Wire

WireGeo const & Wire(WireID const &wireid) const

Returns the specified wire.

Definition: WireReadoutGeom.cxx:270

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

inlineinherited

Returns the specified plane.

Parameters

planeid plane ID

Returns: a constant pointer to the specified plane, or nullptr if none

The GetElementPtr() method is overloaded and its return depends on the type of ID.

Definition at line 164 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::MaxWires(), geo::WireReadoutGeom::NElements(), geo::WireReadoutGeom::Nwires(), and geo::WireReadoutGeom::PlanePtr().

Referenced by PlaneWireToChannel().

164 { return PlanePtr(planeid); }

geo::WireReadoutGeom::PlanePtr

PlaneGeo const * PlanePtr(PlaneID const &planeid) const

Returns the specified plane.

Definition: WireReadoutGeom.cxx:231

WireGeo const* geo::WireReadoutGeom::GetElementPtr ( WireID const & wireid ) const

inlineinherited

Returns the specified wire.

Parameters

wireid wire ID

Returns: a constant pointer to the specified wire, or nullptr if none

The GetElementPtr() method is overloaded and its return depends on the type of ID.

Definition at line 213 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::Wire(), and geo::WireReadoutGeom::WirePtr().

213 { return WirePtr(wireid); }

geo::WireReadoutGeom::WirePtr

WireGeo const * WirePtr(WireID const &wireid) const

Returns the specified wire.

Definition: WireReadoutGeom.cxx:263

template<typename T >

T const* geo::WireReadoutGeom::GetElementPtr	(	PlaneInfoMap_t< T > const &	map,
		PlaneID const &	id
	)		const

inlineprotectedinherited

Returns a pointer to the specified element, or nullptr if invalid.

Definition at line 1055 of file WireReadoutGeom.h.

     {
       if (id.Cryostat >= map.size()) return nullptr;
       auto const& cryo_map = map[id.Cryostat];
       if (id.TPC >= cryo_map.size()) return nullptr;
       auto const& TPC_map = cryo_map[id.TPC];
       if (id.Plane >= TPC_map.size()) return nullptr;
       auto const& plane_map = TPC_map[id.Plane];
       return &plane_map;
     } // GetElementPtr()

unsigned int geo::WireReadoutGeom::HardwareChannelFromOpChannel ( unsigned int opChannel ) const

virtualinherited

Returns the hardware channel number of specified optical channel.

Parameters

opChannel the optical detector channel being queried

Returns: the optical detector the specified optical channel belongs to

If the specified optical channel is invalid, behaviour is undefined.

Definition at line 426 of file WireReadoutGeom.cxx.

Referenced by geo::WireReadoutGeom::HasChannel(), and geo::WireReadoutGeom::IsValidOpChannel().

   {
     return 0;
   }

virtual bool geo::WireReadoutGeom::HasChannel ( raw::ChannelID_t channel ) const

inlinevirtualinherited

Returns whether the specified channel is valid This default implementation assumes all channels up to Nchannels() valid.

Definition at line 567 of file WireReadoutGeom.h.

     {
       return raw::isValidChannelID(channel) ? channel < Nchannels() : false;
     }

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

inlineinherited

Returns whether we have the specified plane.

The HasElement() method is overloaded and its meaning depends on the type of ID.

Definition at line 79 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::HasPlane(), geo::WireReadoutGeom::MaxPlanes(), geo::WireReadoutGeom::Nplanes(), and geo::details::tpc_zero.

79 { return HasPlane(planeid); }

geo::WireReadoutGeom::HasPlane

bool HasPlane(PlaneID const &planeid) const

Returns whether we have the specified plane.

Definition: WireReadoutGeom.cxx:188

bool geo::WireReadoutGeom::HasElement ( WireID const & wireid ) const

inlineinherited

Returns whether we have the specified wire.

The HasElement() method is overloaded and its meaning depends on the type of ID.

Definition at line 201 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::HasWire(), and geo::WireReadoutGeom::WirePtr().

201 { return HasWire(wireid); }

geo::WireReadoutGeom::HasWire

bool HasWire(WireID const &wireid) const

Returns whether we have the specified wire.

Definition: WireReadoutGeom.cxx:256

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

inherited

Returns whether we have the specified plane.

The HasElement() method is overloaded and its meaning depends on the type of ID.

Definition at line 188 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::PlanePtr().

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), geo::WireReadoutGeom::finish(), pma::ProjectionMatchingAlg::guideEndpoints(), geo::WireReadoutGeom::HasElement(), pma::PMAlgTracker::matchCluster(), pma::Node3D::Node3D(), trkf::PMAlgTrajFitter::produce(), and pma::ProjectionMatchingAlg::TestTrk_().

   {
     return PlanePtr(planeid) != nullptr;
   }

bool geo::WireReadoutStandardGeom::HasROP ( readout::ROPID const & ropid ) const

overridevirtual

Returns whether we have the specified ROP

Returns: whether the readout plane is valid and exists

Implements geo::WireReadoutGeom.

Definition at line 394 of file WireReadoutStandardGeom.cxx.

References NROPs(), and readout::ROPID::ROP.

Referenced by Nchannels().

   {
     return ropid.ROP < NROPs(ropid);
   }

bool geo::WireReadoutStandardGeom::HasTPCset ( readout::TPCsetID const & tpcsetid ) const

overridevirtual

Returns whether we have the specified TPC set

Returns: whether the TPC set is valid and exists

Implements geo::WireReadoutGeom.

Definition at line 342 of file WireReadoutStandardGeom.cxx.

References NTPCsets(), and readout::TPCsetID::TPCset.

Referenced by NROPs().

   {
     return tpcsetid.TPCset < NTPCsets(tpcsetid);
   }

bool geo::WireReadoutGeom::HasWire ( WireID const & wireid ) const

inherited

Returns whether we have the specified wire.

The HasElement() method is overloaded and its meaning depends on the type of ID.

Definition at line 256 of file WireReadoutGeom.cxx.

References geo::PlaneGeo::HasWire(), and geo::WireReadoutGeom::PlanePtr().

Referenced by tca::FillmAllTraj(), geo::WireReadoutGeom::HasElement(), geo::WireReadoutGeom::NSiblingElements(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), pma::ProjectionMatchingAlg::validate_on_adc_test(), and geo::WireReadoutGeom::WireIDIntersectionCheck().

   {
     PlaneGeo const* pPlane = PlanePtr(wireid);
     return pPlane ? pPlane->HasWire(wireid) : false;
   }

template<typename T >

bool geo::WireReadoutGeom::isValidElement	(	TPCInfoMap_t< T > const &	map,
		CryostatID const &	id
	)		const

inlineprotectedinherited

Returns whether the ID specifies a valid entry.

Definition at line 1009 of file WireReadoutGeom.h.

Referenced by geo::WireReadoutGeom::isValidElement().

     {
       return id.Cryostat < map.size();
     }

template<typename T >

bool geo::WireReadoutGeom::isValidElement	(	TPCInfoMap_t< T > const &	map,
		TPCID const &	id
	)		const

inlineprotectedinherited

Returns the specified element of the TPC map.

Definition at line 1014 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::isValidElement().

     {
       return isValidElement(map, id.asCryostatID()) && (id.TPC < map[id.Cryostat].size());
     }

template<typename T >

bool geo::WireReadoutGeom::isValidElement	(	PlaneInfoMap_t< T > const &	map,
		CryostatID const &	id
	)		const

inlineprotectedinherited

Returns whether the ID specifies a valid entry.

Definition at line 1037 of file WireReadoutGeom.h.

     {
       return id.Cryostat < map.size();
     }

template<typename T >

bool geo::WireReadoutGeom::isValidElement	(	PlaneInfoMap_t< T > const &	map,
		TPCID const &	id
	)		const

inlineprotectedinherited

Returns whether the ID specifies a valid entry.

Definition at line 1042 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::isValidElement().

     {
       return isValidElement(map, id.asCryostatID()) && (id.TPC < map[id.Cryostat].size());
     }

template<typename T >

bool geo::WireReadoutGeom::isValidElement	(	PlaneInfoMap_t< T > const &	map,
		PlaneID const &	id
	)		const

inlineprotectedinherited

Returns whether the ID specifies a valid entry.

Definition at line 1047 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::isValidElement().

     {
       return isValidElement(map, id.asTPCID()) && (id.Plane < AccessSize(map, id.asTPCID()));
     }

bool geo::WireReadoutGeom::IsValidOpChannel ( int opChannel ) const

inherited

Is this a valid OpChannel number?

Definition at line 493 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fGeom, and geo::GeometryCore::NOpDets().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), geo::WireReadoutGeom::HasChannel(), cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison(), and opdet::RunHitFinder().

   {
     return IsValidOpChannel(opChannel, fGeom->NOpDets());
   }

bool geo::WireReadoutGeom::IsValidOpChannel	(	unsigned int	opChannel,
		unsigned int	NOpDets
	)		const

virtualinherited

Returns whether the ID identifies a valid optical detector channel.

Parameters

opChannel	channel number
NOpDets	number of optical detectors in the detector

Returns: whether opChannel would be a valid channel

The specification of the number of optical channels reflects the logic described in NOpChannel().

Definition at line 432 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::HardwareChannelFromOpChannel(), geo::WireReadoutGeom::NOpChannels(), geo::WireReadoutGeom::NOpHardwareChannels(), and geo::WireReadoutGeom::OpDetFromOpChannel().

   {
     // Check channel number
     if (opChannel >= NOpChannels(NOpDets)) return false;
 
     // Check opdet number
     unsigned int opdet = OpDetFromOpChannel(opChannel);
     if (opdet >= NOpDets) return false;
 
     // Check hardware channel number
     unsigned int hChan = HardwareChannelFromOpChannel(opChannel);
     return hChan < NOpHardwareChannels(opdet);
   }

unsigned int geo::WireReadoutGeom::MaxOpChannel ( ) const

inherited

Largest optical channel number.

Definition at line 487 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fGeom, and geo::GeometryCore::NOpDets().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), opdet::ConstructFlash(), geo::WireReadoutGeom::HasChannel(), and cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison().

   {
     return MaxOpChannel(fGeom->NOpDets());
   }

unsigned int geo::WireReadoutGeom::MaxOpChannel ( unsigned int NOpDets ) const

virtualinherited

Returns the number of optical channels contained in some detectors.

Parameters

NOpDets number of optical detectors

Returns: optical channels contained in NOpDets detectors

This function returns the first optical channel ID larger than the last channel ID in a detector with NOpDets optical detectors (with the same logic as NOpChannels()). For example, in a detector with 32 channels with contiguous IDs starting at 0, this function would return 32. If the channels started with ID 1, this function would instead return 33 and if there were a 16 channel gap, so that valid channels are from 0 to 15 and from 32 to 47, this function would return 48.

Definition at line 394 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::NOpChannels().

   {
     // By default just return the number of optical detectos
     return NOpChannels(NOpDets);
   }

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

inherited

Returns the largest number of planes among all TPCs in this detector.

Definition at line 166 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fPlanes.

Referenced by pma::PMAlgTracker::build(), ShowerRecoTools::ShowerTrajPointdEdx::CalculateElement(), shower::EMShowerAlg::CheckIsolatedHits_(), img::DataProviderAlg::DataProviderAlg(), trkf::TrackKalmanFitter::doFitWork(), geo::details::ReadoutIterationPolicy::EndPlaneID(), geo::WireReadoutGeom::HasElement(), shower::EMShowerAlg::MakeShower(), shower::EMShowerAlg::MakeSpacePoints(), geo::WireReadoutGeom::Nviews(), shower::EMShowerAlg::OrderShowerHits(), pma::PMAlgTracker::PMAlgTracker(), shower::EMShowerAlg::RelativeWireWidth_(), trkf::TrackKalmanFitter::sortOutput(), and shower::EMShowerAlg::WorstPlane_().

   {
     unsigned int maxPlanes = 0;
     for (auto const& [tpc_id, planes] : fPlanes) {
       unsigned int maxPlanesInTPC = planes.size();
       if (maxPlanesInTPC > maxPlanes) maxPlanes = maxPlanesInTPC;
     }
     return maxPlanes;
   }

unsigned int geo::WireReadoutStandardGeom::MaxROPs ( ) const

overridevirtual

Returns the largest number of ROPs a TPC set in the detector has.

Implements geo::WireReadoutGeom.

Definition at line 384 of file WireReadoutStandardGeom.cxx.

References fNPlanes.

   {
     unsigned int max = 0;
     for (auto const& cryo_tpc : fNPlanes)
       for (unsigned int nPlanes : cryo_tpc)
         if (nPlanes > max) max = nPlanes;
     return max;
   }

unsigned int geo::WireReadoutStandardGeom::MaxTPCs ( ) const

private

Returns the largest number of TPCs in a single cryostat.

Definition at line 368 of file WireReadoutStandardGeom.cxx.

References fNTPC.

Referenced by MaxTPCsets(), and WireCount().

   {
     unsigned int max = 0;
     for (unsigned int nTPCs : fNTPC)
       if (nTPCs > max) max = nTPCs;
     return max;
   }

unsigned int geo::WireReadoutStandardGeom::MaxTPCsets ( ) const

overridevirtual

Returns the largest number of TPC sets any cryostat in the detector has.

Implements geo::WireReadoutGeom.

Definition at line 336 of file WireReadoutStandardGeom.cxx.

References MaxTPCs().

   {
     return MaxTPCs();
   }

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

inherited

Returns the total number of wires in the specified plane.

Parameters

planeid plane ID

Returns: number of wires in specified plane, or 0 if no plane found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 177 of file WireReadoutGeom.cxx.

Referenced by cluster::BlurredClusteringAlg::ConvertRecobHitsToVector(), geo::details::ReadoutIterationPolicy::EndWireID(), and geo::WireReadoutGeom::GetElementPtr().

   {
     unsigned int maxWires = 0;
     for (PlaneGeo const& plane : Iterate<PlaneGeo>()) {
       unsigned int maxWiresInPlane = plane.Nwires();
       if (maxWiresInPlane > maxWires) maxWires = maxWiresInPlane;
     }
     return maxWires;
   }

raw::ChannelID_t geo::WireReadoutStandardGeom::Nchannels ( ) const

overridevirtual

Returns the total number of channels present (not necessarily contiguous)

Implements geo::WireReadoutGeom.

Definition at line 206 of file WireReadoutStandardGeom.cxx.

References fNchannels.

   {
     return fNchannels;
   }

unsigned int geo::WireReadoutStandardGeom::Nchannels ( readout::ROPID const & ropid ) const

overridevirtual

Returns the number of channels in the specified ROP.

Returns: number of channels in the specified ROP, 0 if non-existent

Implements geo::WireReadoutGeom.

Definition at line 212 of file WireReadoutStandardGeom.cxx.

References FirstWirePlaneInROP(), HasROP(), and WireCount().

   {
     if (!HasROP(ropid)) return 0;
     // The number of channels matches the number of wires. Life is easy.
     return WireCount(FirstWirePlaneInROP(ropid));
   }

raw::ChannelID_t geo::WireReadoutGeom::NearestChannel	(	Point_t const &	worldLoc,
		PlaneID const &	planeid
	)		const

inherited

Returns the ID of the channel nearest to the specified position.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the wire plane the channel must belong to

Returns: the ID of the channel, or raw::InvalidChannelID if invalid wire

Bug:: on invalid wire, a geo::InvalidWireError exception is thrown

Definition at line 524 of file WireReadoutGeom.cxx.

References raw::InvalidChannelID, geo::PlaneGeo::NearestWireID(), geo::WireReadoutGeom::Plane(), and geo::WireReadoutGeom::PlaneWireToChannel().

Referenced by larg4::LArVoxelReadout::DriftIonizationElectrons(), geo::WireReadoutGeom::HasChannel(), and apa::DisambigAlg::TrivialDisambig().

   {
     // This method is supposed to return a channel number rather than a wire number.
     // Perform the conversion here (although, maybe faster if we deal in wire numbers
     // rather than channel numbers?)
 
     // NOTE on failure both NearestChannel() and upstream:
     // * according to documentation, should return invalid channel
     // * in the actual code throw an exception because of a BUG
     //
     // The following implementation automatically becomes in fact compliant to the
     // documentation if upstreams becomes compliant to.  When that happens, just delete
     // this comment.
     WireID const wireID = Plane(planeid).NearestWireID(worldPos);
     return wireID ? PlaneWireToChannel(wireID) : raw::InvalidChannelID;
   }

WireID geo::WireReadoutStandardGeom::NearestWireID	(	Point_t const &	worldPos,
		PlaneID const &	planeID
	)		const

overridevirtual

Returns the ID of the wire nearest to the specified position.

Parameters

worldPos	position to be tested
planeID	plane containing the wire

Returns: the ID of the wire closest to worldPos in the specified plane

Exceptions

InvalidWireIDError the ID found is not present in the detector

See also: WireCoordinate(double, double, geo::PlaneID const&)

The plane is required to be valid and exist in the detector. Otherwise, the behaviour is undefined. An exception is thrown if the wire that would be the closest is actually not present; but no check is performed whether the specified position is outside the wire plane: wires are extrapolated to be infinitely long. In other words, the result can be trusted only as long as the position is within the specified wire plane.

Implements geo::WireReadoutGeom.

Definition at line 229 of file WireReadoutStandardGeom.cxx.

References WireCoordinate(), and WireCount().

   {
 
     // This part is the actual calculation of the nearest wire number, where we assume
     //  uniform wire pitch and angle within a wireplane
 
     // add 0.5 to have the correct rounding
     int NearestWireNumber = int(0.5 + WireCoordinate(worldPos.Y(), worldPos.Z(), planeID));
 
     // If we are outside of the wireplane range, throw an exception
     // (this response maintains consistency with the previous
     // implementation based on geometry lookup)
     if (NearestWireNumber < 0 || (unsigned int)NearestWireNumber >= WireCount(planeID)) {
       int wireNumber = NearestWireNumber; // save for the output
 
       if (NearestWireNumber < 0)
         NearestWireNumber = 0;
       else
         NearestWireNumber = WireCount(planeID) - 1;
 
       throw InvalidWireError("Geometry", planeID, wireNumber, NearestWireNumber)
         << "Can't Find Nearest Wire for position (" << worldPos.X() << "," << worldPos.Y() << ","
         << worldPos.Z() << ")"
         << " in plane " << std::string(planeID) << " approx wire number # " << wireNumber
         << " (capped from " << NearestWireNumber << ")\n";
     }
 
     return WireID(planeID, (WireID::WireID_t)NearestWireNumber);
   }

unsigned int geo::WireReadoutGeom::NElements ( TPCID const & tpcid ) const

inlineinherited

Returns the total number of planes in the specified TPC.

Parameters

tpcid TPC ID

Returns: number of planes in specified TPC, or 0 if no TPC found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 103 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::Nplanes().

Referenced by geo::WireReadoutGeom::GetElementPtr().

103 { return Nplanes(tpcid); }

geo::WireReadoutGeom::Nplanes

unsigned int Nplanes(TPCID const &tpcid=details::tpc_zero) const

Returns the total number of planes in the specified TPC.

Definition: WireReadoutGeom.cxx:194

unsigned int geo::WireReadoutGeom::NElements ( PlaneID const & planeid ) const

inherited

Returns the total number of wires in the specified plane.

Parameters

planeid plane ID

Returns: number of wires in specified plane, or 0 if no plane found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

unsigned int geo::WireReadoutGeom::NOpChannels ( unsigned int NOpDets ) const

virtualinherited

Returns the number of optical channels contained in some detectors.

Parameters

NOpDets number of optical detectors

Returns: optical channels contained in NOpDets detectors

This function returns how many channels can be expected to be present in a detector with NOpDets optical detectors. This is an upper limit, as not all channels have necessarily to be present.

For example: if a detector has four channels per optical detector, the returned value will be four times the argument NOpDets. If there is a single channel on each optical detector, the return value will be the value NOpDets (this is in fact the fallback implementation). If each optical detector can have anywhere between 2 and 12 channels, the returned value is 12 times NOpDets, and it will be an overestimation of the number of channels.

Definition at line 387 of file WireReadoutGeom.cxx.

Referenced by opdet::OptDetDigitizer::produce().

   {
     // By default just return the number of optical detectos
     return NOpDets;
   }

unsigned int geo::WireReadoutGeom::NOpChannels ( ) const

inherited

Number of electronics channels for all the optical detectors.

Definition at line 481 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fGeom, and geo::GeometryCore::NOpDets().

Referenced by geo::WireReadoutGeom::HasChannel(), geo::WireReadoutGeom::IsValidOpChannel(), and geo::WireReadoutGeom::MaxOpChannel().

   {
     return NOpChannels(fGeom->NOpDets());
   }

unsigned int geo::WireReadoutGeom::NOpHardwareChannels ( unsigned int opDet ) const

virtualinherited

Returns the number of channels in the specified optical detectors.

Parameters

opDet ID of the chosen optical detector

Returns: optical channels contained in optical detector with ID opDet

This function returns how many channels are actually present in the optical detector with the specified ID.

For example: if a detector has four channels per optical detector, the returned value will be four, regardless opDet, and . If there is a single channel on each optical detector, the return value will be 1, again ignoring opDet (this is in fact the fallback implementation). If each optical detector can have anywhere between 2 and 12 channels, the returned value will be 2, 12, etc., that is the exact number of channels in opDet.

Although implementations are encouraged to return 0 on invalid optical detectors, the actual return value in that case is undefined.

Definition at line 401 of file WireReadoutGeom.cxx.

Referenced by geo::WireReadoutGeom::HasChannel(), and geo::WireReadoutGeom::IsValidOpChannel().

   {
     // By default, 1 channel per optical detector
     return 1;
   }

unsigned int geo::WireReadoutGeom::Nplanes ( TPCID const & tpcid = details::tpc_zero ) const

inherited

Returns the total number of planes in the specified TPC.

Parameters

tpcid TPC ID

Returns: number of planes in specified TPC, or 0 if no TPC found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 194 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fPlanes.

   {
     auto it = fPlanes.find(tpcid);
     if (it == fPlanes.cend()) { return 0; }
     return it->second.size();
   }

unsigned int geo::WireReadoutStandardGeom::NROPs ( readout::TPCsetID const & tpcsetid ) const

overridevirtual

Returns the total number of ROPs in the specified TPC set.

Parameters

tpcsetid TPC set ID

Returns: number of readout planes in the TPC sets, or 0 if ID is invalid

Note that this methods explicitly check the existence of the TPC set.

In this mapping, planes have independent readout and there is one wire plane in each readout plane and one readout plane for each wire plane.

Implements geo::WireReadoutGeom.

Definition at line 377 of file WireReadoutStandardGeom.cxx.

References geo::WireReadoutGeom::AccessElement(), FirstTPCinTPCset(), fNPlanes, and HasTPCset().

Referenced by HasROP().

   {
     if (!HasTPCset(tpcsetid)) return 0;
     return AccessElement(fNPlanes, FirstTPCinTPCset(tpcsetid));
   }

unsigned int geo::WireReadoutGeom::NSiblingElements ( PlaneID const & planeid ) const

inlineinherited

Returns the total number of planes in the specified TPC.

Parameters

tpcid TPC ID

Returns: number of planes in specified TPC, or 0 if no TPC found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 104 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::Nplanes(), and geo::WireReadoutGeom::Nviews().

Referenced by geo::details::ReadoutIterationPolicy::NSiblings().

104 { return Nplanes(planeid); }

geo::WireReadoutGeom::Nplanes

unsigned int Nplanes(TPCID const &tpcid=details::tpc_zero) const

Returns the total number of planes in the specified TPC.

Definition: WireReadoutGeom.cxx:194

unsigned int geo::WireReadoutGeom::NSiblingElements ( WireID const & wireid ) const

inlineinherited

Returns the total number of wires in the specified plane.

Parameters

planeid plane ID

Returns: number of wires in specified plane, or 0 if no plane found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 186 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::HasWire(), and geo::WireReadoutGeom::Nwires().

186 { return Nwires(wireid); }

geo::WireReadoutGeom::Nwires

unsigned int Nwires(PlaneID const &planeid) const

Returns the total number of wires in the specified plane.

Definition: WireReadoutGeom.cxx:249

unsigned int geo::WireReadoutStandardGeom::NTPCsets ( readout::CryostatID const & cryoid ) const

overridevirtual

Returns the total number of TPC sets in the specified cryostat.

Parameters

cryoid cryostat ID

Returns: number of TPC sets in the cryostat, or 0 if no cryostat found

In this mapping, TPCs have independent readout and there is one TPC in each TPC set and one TPC set for each TPC.

Implements geo::WireReadoutGeom.

Definition at line 329 of file WireReadoutStandardGeom.cxx.

References geo::CryostatID::Cryostat, fNTPC, and geo::CryostatID::isValid.

Referenced by HasTPCset().

   {
     // return the same number as the number of TPCs
     return (cryoid.isValid && cryoid.Cryostat < fNTPC.size()) ? fNTPC[cryoid.Cryostat] : 0;
   }

unsigned int geo::WireReadoutGeom::Nviews ( ) const

inherited

Returns the number of views (different wire orientations)

Returns the number of different views, or wire orientations, in the detector.

The function assumes that all TPCs in all cryostats of a detector have the same number of planes, which should be a safe assumption.

Todo:: Change return type to size_t

Definition at line 160 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::MaxPlanes().

Referenced by geo::WireReadoutGeom::NSiblingElements().

   {
     return MaxPlanes();
   }

unsigned int geo::WireReadoutGeom::Nwires ( PlaneID const & planeid ) const

inherited

Returns the total number of wires in the specified plane.

Parameters

planeid plane ID

Returns: number of wires in specified plane, or 0 if no plane found

The NElements() and NSiblingElements() methods are overloaded and their return depends on the type of ID.

Todo:: Change return type to size_t

Definition at line 249 of file WireReadoutGeom.cxx.

References geo::PlaneGeo::NElements(), and geo::WireReadoutGeom::PlanePtr().

Referenced by lar_cluster3d::StandardHit3DBuilder::BuildChannelStatusVec(), lar_cluster3d::SnippetHit3DBuilder::BuildChannelStatusVec(), trkf::CCTrackMaker::FillWireHitRange(), tca::FillWireHitRange(), geo::WireReadoutGeom::GetElementPtr(), util::GeometryUtilities::GetYZ(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire_(), corner::CornerFinderAlg::GrabWires(), corner::CornerFinderAlg::InitializeGeometry(), lar_cluster3d::StandardHit3DBuilder::NearestWireID(), lar_cluster3d::SnippetHit3DBuilder::NearestWireID(), geo::WireReadoutGeom::NSiblingElements(), util::GeometryUtilities::SelectPolygonHitList(), img::DataProviderAlg::setWireDriftData(), and geo::WireReadoutGeom::WireIDIntersectionCheck().

   {
     PlaneGeo const* pPlane = PlanePtr(planeid);
     return pPlane ? pPlane->NElements() : 0;
   }

unsigned int geo::WireReadoutGeom::OpChannel	(	unsigned int	detNum,
		unsigned int	hwchannel = `0`
	)		const

virtualinherited

Returns the channel ID of the specified hardware channel.

Parameters

detNum	optical detector ID
hwchannel	hardware channel within the specified optical detector

Returns: ID of the channel identified by detector and hardware channel IDs

If the input IDs identify a non-existing channel, the result is undefined.

Definition at line 408 of file WireReadoutGeom.cxx.

Referenced by geo::WireReadoutGeom::HasChannel().

   {
     return detNum;
   }

unsigned int geo::WireReadoutGeom::OpDetFromOpChannel ( unsigned int opChannel ) const

virtualinherited

Returns the optical detector the specified optical channel belongs.

Parameters

opChannel the optical detector channel being queried

Returns: the optical detector the specified optical channel belongs to

If the specified optical channel is invalid, behaviour is undefined.

Definition at line 414 of file WireReadoutGeom.cxx.

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), geo::WireReadoutGeom::HasChannel(), geo::WireReadoutGeom::IsValidOpChannel(), geo::WireReadoutGeom::OpDetGeoFromOpChannel(), and cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison().

   {
     return opChannel;
   }

OpDetGeo const & geo::WireReadoutGeom::OpDetGeoFromOpChannel ( unsigned int opChannel ) const

inherited

Returns the optical detector the specified optical channel belongs.

Parameters

opChannel the optical detector channel being queried

Returns: the optical detector geometry object the specified optical channel belongs to

If the specified optical channel is invalid, behaviour is undefined.

Definition at line 420 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fGeom, geo::WireReadoutGeom::OpDetFromOpChannel(), and geo::GeometryCore::OpDetGeoFromOpDet().

Referenced by opdet::GetHitGeometryInfo(), and geo::WireReadoutGeom::HasChannel().

   {
     return fGeom->OpDetGeoFromOpDet(OpDetFromOpChannel(opChannel));
   }

PlaneGeo const & geo::WireReadoutGeom::Plane	(	TPCID const &	tpcid,
		View_t	view
	)		const

inherited

Returns the specified wire.

Parameters

tpcid	ID of the TPC
view	number of the requested view
planeid	ID of the plane

Returns: a constant reference to the specified plane

Exceptions

cet::exception	(`GeometryCore` category) if cryostat not present
cet::exception	(`TPCOutOfRange` category) if no such TPC
cet::exception	(`PlaneOutOfRange` category) if no such plane

The GetElement() method is overloaded and its return depends on the type of ID.

Definition at line 210 of file WireReadoutGeom.cxx.

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

   {
     auto it = fPlanes.find(tpcid);
     if (it == fPlanes.cend()) {
       throw cet::exception("WireReadoutGeom") << "No TPC with ID " << tpcid << " supported.";
     }
     for (PlaneGeo const& plane : it->second) {
       if (plane.View() == view) { return plane; }
     }
     throw cet::exception("WireReadoutGeom")
       << "TPCGeo " << tpcid << " has no plane for view #" << (size_t)view << "\n";
   }

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

inherited

Returns the specified wire.

Parameters

tpcid	ID of the TPC
view	number of the requested view
planeid	ID of the plane

Returns: a constant reference to the specified plane

Exceptions

cet::exception	(`GeometryCore` category) if cryostat not present
cet::exception	(`TPCOutOfRange` category) if no such TPC
cet::exception	(`PlaneOutOfRange` category) if no such plane

The GetElement() method is overloaded and its return depends on the type of ID.

Definition at line 224 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::PlanePtr().

   {
     if (auto const* plane_ptr = PlanePtr(planeid)) { return *plane_ptr; }
     throw cet::exception("WireReadoutGeom") << "Plane with ID " << planeid << " is not supported.";
   }

double geo::WireReadoutGeom::Plane0Pitch	(	TPCID const &	tpcid,
		unsigned int	p1
	)		const

inherited

Definition at line 124 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fPlane0Pitch.

Referenced by larg4::LArVoxelReadout::DriftIonizationElectrons(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     return fPlane0Pitch.at(tpcid)[p];
   }

std::set< PlaneID > const & geo::WireReadoutStandardGeom::PlaneIDs ( ) const

overridevirtual

Returns a list of the plane IDs in the whole detector.

Implements geo::WireReadoutGeom.

Definition at line 323 of file WireReadoutStandardGeom.cxx.

References fPlaneIDs.

   {
     return fPlaneIDs;
   }

double geo::WireReadoutGeom::PlanePitch	(	TPCID const &	tpcid,
		unsigned int	p1 = `0`,
		unsigned int	p2 = `1`
	)		const

inherited

Definition at line 130 of file WireReadoutGeom.cxx.

References util::abs(), and geo::WireReadoutGeom::fPlane0Pitch.

Referenced by detinfo::DetectorPropertiesStandard::DataFor(), larg4::LArVoxelReadout::DriftIonizationElectrons(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     auto const& plane0_pitch = fPlane0Pitch.at(tpcid);
     return std::abs(plane0_pitch[p2] - plane0_pitch[p1]);
   }

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

inherited

Returns the specified plane.

Parameters

planeid plane ID

Returns: a constant pointer to the specified plane, or nullptr if none

The GetElementPtr() method is overloaded and its return depends on the type of ID.

Definition at line 231 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fPlanes, n, geo::PlaneID::parentID(), and geo::PlaneID::Plane.

Referenced by geo::WireReadoutGeom::FirstPlane(), geo::WireReadoutGeom::GetElement(), geo::details::getElementPtr(), geo::WireReadoutGeom::GetElementPtr(), geo::WireReadoutGeom::HasPlane(), geo::WireReadoutGeom::HasWire(), geo::WireReadoutGeom::Nwires(), geo::WireReadoutGeom::Plane(), and geo::WireReadoutGeom::WirePtr().

   {
     auto const [tpc_id, plane] = std::make_pair(planeid.parentID(), planeid.Plane);
     auto it = fPlanes.find(tpc_id);
     if (it == fPlanes.cend()) { return nullptr; }
     if (std::size_t const n = it->second.size(); n <= plane) { return nullptr; }
     return &it->second[plane];
   }

raw::ChannelID_t geo::WireReadoutStandardGeom::PlaneWireToChannel ( WireID const & wireID ) const

overridevirtual

Returns the channel ID a wire is connected to.

Parameters

wireID ID of the wire

Returns: the ID of the channel

See also: PlaneWireToChannel(geo::WireID const&)

Behaviour on an invalid or not present wires is undefined.

Implements geo::WireReadoutGeom.

Definition at line 273 of file WireReadoutStandardGeom.cxx.

References fPlaneBaselines, geo::WireReadoutGeom::GetElementPtr(), raw::InvalidChannelID, and geo::WireID::Wire.

   {
     unsigned int const* pBaseLine = GetElementPtr(fPlaneBaselines, wireID);
     // This is the actual lookup part - first make sure coordinates are legal
     if (pBaseLine) {
       // if the channel has legal coordinates, its ID is given by the wire
       // number above the number of wires in lower planes, tpcs and cryostats
       return *pBaseLine + wireID.Wire;
     }
     else {
       // if the coordinates were bad, throw an exception
       throw cet::exception("WireReadoutStandardGeom")
         << "NO CHANNEL FOUND for " << std::string(wireID);
     }
 
     // made it here, that shouldn't happen, return raw::InvalidChannelID
     mf::LogWarning("WireReadoutStandardGeom")
       << "should not be at the point in the function, returning "
       << "invalid channel";
     return raw::InvalidChannelID;
   }

std::vector< TPCID > geo::WireReadoutStandardGeom::ROPtoTPCs ( readout::ROPID const & ropid ) const

overridevirtual

Returns a list of ID of TPCs the specified ROP spans.

Parameters

ropid ID of the readout plane

Returns: the list of TPC IDs, empty if readout plane ID is invalid

In this mapping, readout planes and wire planes are mapped one-to-one. The returned list contains always one entry, unless the specified readout plane ID is invalid, in which case the list is empty. Note that this check is performed on the validity of the readout plane ID, that does not necessarily imply that the readout plane specified by the ID actually exists. Check if the ROP exists with HasROP(). The behaviour on non-existing readout planes is undefined.

Implements geo::WireReadoutGeom.

Definition at line 414 of file WireReadoutStandardGeom.cxx.

References readout::ROPID::asTPCsetID(), and geo::CryostatID::isValid.

   {
     std::vector<TPCID> IDs;
     // we take the TPC set of the ROP and convert it straight into a TPC ID
     if (ropid.isValid) IDs.emplace_back(ConvertTPCsetToTPC(ropid.asTPCsetID()));
     return IDs;
   }

std::vector< PlaneID > geo::WireReadoutStandardGeom::ROPtoWirePlanes ( readout::ROPID const & ropid ) const

overridevirtual

Returns a list of ID of wire planes belonging to the specified ROP.

Parameters

ropid ID of the readout plane to convert into wire planes

Returns: the list of wire plane IDs, empty if readout plane ID is invalid

In this mapping, readout planes and wire planes are mapped one-to-one. The returned list contains always one entry, unless the specified readout plane ID is invalid, in which case the list is empty. Note that this check is performed on the validity of the readout plane ID, that does not necessarily imply that the readout plane specified by the ID actually exists.

Implements geo::WireReadoutGeom.

Definition at line 406 of file WireReadoutStandardGeom.cxx.

References FirstWirePlaneInROP(), and geo::CryostatID::isValid.

   {
     std::vector<PlaneID> IDs;
     if (ropid.isValid) IDs.emplace_back(FirstWirePlaneInROP(ropid));
     return IDs;
   }

SigType_t geo::WireReadoutGeom::SignalType ( PlaneID const & pid ) const

inherited

Returns the type of signal on the channels of specified TPC plane.

Parameters

plane TPC plane ID

Returns: the type of signal on the specified plane, or geo::kMysteryType

Assumes that all the channels on the plane have the same signal type.

Todo:: verify that kMysteryType is returned on invalid plane

Definition at line 499 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::WirePlaneToROP().

Referenced by geo::WireReadoutDumper::dumpTPCplane(), reco3d::TripletFinder::FillHitMap(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire_(), geo::WireReadoutGeom::HasChannel(), shwf::ShowerReco::LongTransEnergy(), caldata::CalWire::produce(), hit::RFFHitFinderAlg::Run(), hit::CCHitFinderAlg::RunCCHitFinder(), and geo::WireReadoutGeom::SignalTypeForROPIDImpl().

   {
     // map wire plane -> readout plane -> first channel, then use SignalType(channel)
 
     auto const ropid = WirePlaneToROP(pid);
     if (!ropid.isValid) {
       throw cet::exception("WireReadoutGeom") << "SignalType(): Mapping of wire plane "
                                               << std::string(pid) << " to readout plane failed!\n";
     }
     return SignalType(ropid);
   }

SigType_t geo::WireReadoutGeom::SignalType ( raw::ChannelID_t const channel ) const

inherited

Return the signal type of the specified channel.

Parameters

channel ID of the channel

Returns: signal type of the channel, or geo::kMysteryType if not known

On any type of error (e.g., invalid or unknown channel ID), geo::kMysteryType is returned.

Definition at line 446 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::SignalTypeForChannelImpl().

   {
     return SignalTypeForChannelImpl(channel);
   }

SigType_t geo::WireReadoutGeom::SignalType ( readout::ROPID const & ropid ) const

inherited

Return the signal type on the specified readout plane.

Parameters

ropid ID of the readout plane

Returns: signal type on the plane, or geo::kMysteryType if not known

If the readout plane ID is marked invalid, geo::kMysteryType is returned. If the readout plane is not marked invalid, but it does not match an existing readout plane, the result is undefined.

The default implementation uses readout plane to channel mapping. Other implementation may decide to do the opposite.

Definition at line 451 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::SignalTypeForROPIDImpl().

   {
     return SignalTypeForROPIDImpl(ropid);
   }

SigType_t geo::WireReadoutStandardGeom::SignalTypeForChannelImpl ( raw::ChannelID_t const channel ) const

overrideprivatevirtual

Return the signal type of the specified channel.

Parameters

channel ID of the channel

Returns: signal type of the channel, or geo::kMysteryType if not known

On any type of error (e.g., invalid or unknown channel ID), geo::kMysteryType is returned.

Implements geo::WireReadoutGeom.

Definition at line 296 of file WireReadoutStandardGeom.cxx.

References geo::WireReadoutGeom::fFirstChannelInNextPlane, geo::WireReadoutGeom::fFirstChannelInThisPlane, fNchannels, fNPlanes, fNTPC, geo::kCollection, geo::kInduction, and geo::kMysteryType.

   {
 
     // still assume one cryostat for now -- faster
     unsigned int nChanPerTPC = fNchannels / fNTPC[0];
     // casting wil trunc towards 0 -- faster than floor
     unsigned int tpc = channel / nChanPerTPC;
     //need number of planes to know Collection
     unsigned int PlanesThisTPC = fNPlanes[0][tpc];
 
     SigType_t sigt = kMysteryType;
     if ((channel >= fFirstChannelInThisPlane[0][tpc][0]) &&
         (channel < fFirstChannelInNextPlane[0][tpc][PlanesThisTPC - 2])) {
       sigt = kInduction;
     }
     else if ((channel >= fFirstChannelInThisPlane[0][tpc][PlanesThisTPC - 1]) &&
              (channel < fFirstChannelInNextPlane[0][tpc][PlanesThisTPC - 1])) {
       sigt = kCollection;
     }
     else
       mf::LogWarning("BadChannelSignalType")
         << "Channel " << channel << " not given signal type." << std::endl;
 
     return sigt;
   }

SigType_t geo::WireReadoutGeom::SignalTypeForROPIDImpl ( readout::ROPID const & ropid ) const

protectedvirtualinherited

Return the signal type on the specified readout plane.

Parameters

ropid ID of the readout plane

Returns: signal type on the plane, or geo::kMysteryType if not known

If the readout plane ID is marked invalid, geo::kMysteryType is returned. If the readout plane is not marked invalid, but it does not match an existing readout plane, the result is undefined.

The default implementation uses readout plane to channel mapping. Other implementation may decide to do the opposite.

Definition at line 456 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::FirstChannelInROP(), and geo::WireReadoutGeom::SignalType().

Referenced by geo::WireReadoutGeom::FirstChannelInThisPlane(), and geo::WireReadoutGeom::SignalType().

   {
     return SignalType(FirstChannelInROP(ropid));
   }

static auto const& geo::WireReadoutGeom::start ( Segment const & s )

inlinestaticinherited

Definition at line 63 of file WireReadoutGeom.h.

Referenced by geo::WireReadoutGeom::WireEndPoints(), and geo::WireReadoutGeom::WireIDsIntersect().

63 { return s.first; }

PlaneID geo::WireReadoutGeom::ThirdPlane	(	PlaneID const &	pid1,
		PlaneID const &	pid2
	)		const

inherited

Returns the plane that is not in the specified arguments.

Parameters

pid1	a plane
pid2	another plane

Returns: the ID to the third plane

Exceptions

cet::exception	(category: "GeometryCore") if other than 3 planes
cet::exception	(category: "GeometryCore") if pid1 and pid2 match

This function requires a geometry with exactly three planes. If the two input planes are not on the same TPC, the result is undefined.

Definition at line 610 of file WireReadoutGeom.cxx.

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

Referenced by geo::WireReadoutGeom::ThirdPlane_dTdW(), geo::WireReadoutGeom::ThirdPlaneSlope(), and geo::WireReadoutGeom::WireEndPoints().

   {
     // how many planes in the TPC pid1 belongs to:
     unsigned const int nPlanes = Nplanes(pid1);
     if (nPlanes != 3) {
       throw cet::exception("GeometryCore")
         << "ThirdPlane() supports only TPCs with 3 planes, and I see " << nPlanes << " instead\n";
     }
 
     PlaneID::PlaneID_t target_plane = nPlanes;
     for (PlaneID::PlaneID_t iPlane = 0; iPlane < nPlanes; ++iPlane) {
       if ((iPlane == pid1.Plane) || (iPlane == pid2.Plane)) continue;
       if (target_plane != nPlanes) {
         throw cet::exception("GeometryCore")
           << "ThirdPlane() found too many planes that are not " << std::string(pid1) << " nor "
           << std::string(pid2) << "! (first " << target_plane << ", then " << iPlane << ")\n";
       } // if we had a target already
       target_plane = iPlane;
     } // for
     if (target_plane == nPlanes) {
       throw cet::exception("GeometryCore")
         << "ThirdPlane() can't find a plane that is not " << std::string(pid1) << " nor "
         << std::string(pid2) << "!\n";
     }
 
     return PlaneID(pid1, target_plane);
   }

double geo::WireReadoutGeom::ThirdPlane_dTdW	(	PlaneID const &	pid1,
		double	slope1,
		PlaneID const &	pid2,
		double	slope2,
		PlaneID const &	output_plane
	)		const

inherited

Returns dT/dW on the third plane, given it in the other two.

Parameters

pid1	ID of the plane of the first dT/dW
dTdW1	dT/dW as seen on the first plane
pid2	ID of the plane of the second dT/dW
dTdW2	dT/dW as seen on the second plane
output_plane	ID of the plane on which to calculate the slope

Returns: dT/dW on the third plane, or -999. if dT/dW would be infinity

Exceptions

cet::exception	(category: "GeometryCore") if different TPC
cet::exception	(category: "GeometryCore") if same plane
cet::exception	(category: "GeometryCore") if other than 3 planes

Given a dT/dW as projected in two planes, returns the dT/dW as projected in the third plane. The dT/dW are defined in time ticks/wide number units.

Definition at line 663 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ComputeThirdPlane_dTdW(), geo::PlaneGeo::PhiZ(), geo::WireReadoutGeom::Plane(), and geo::PlaneGeo::WirePitch().

Referenced by geo::WireReadoutGeom::ThirdPlane_dTdW(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     CheckIndependentPlanesOnSameTPC(pid1, pid2, "ThirdPlane_dTdW()");
 
     double angle[3], pitch[3];
     PlaneGeo const* const planes[3] = {&Plane(pid1), &Plane(pid2), &Plane(output_plane)};
 
     // We need wire pitch and "wire coordinate direction" for each plane.  The latter is
     // perpendicular to the wire orientation.  PlaneGeo::PhiZ() defines the right
     // orientation too.
     for (size_t i = 0; i < 3; ++i) {
       angle[i] = planes[i]->PhiZ();
       pitch[i] = planes[i]->WirePitch();
     }
 
     return ComputeThirdPlane_dTdW(
       angle[0], pitch[0], slope1, angle[1], pitch[1], slope2, angle[2], pitch[2]);
   }

double geo::WireReadoutGeom::ThirdPlane_dTdW	(	PlaneID const &	pid1,
		double	slope1,
		PlaneID const &	pid2,
		double	slope2
	)		const

inherited

Returns dT/dW on the third plane, given it in the other two.

Parameters

pid1	ID of the plane of the first dT/dW
dTdW1	dT/dW as seen on the first plane
pid2	ID of the plane of the second dT/dW
dTdW2	dT/dW as seen on the second plane

Returns: dT/dW on the third plane, or -999. if dT/dW would be infinity

Exceptions

cet::exception	(category: "GeometryCore") if different TPC
cet::exception	(category: "GeometryCore") if same plane
cet::exception	(category: "GeometryCore") if other than 3 planes

Given a dT/dW as projected in two planes, returns the dT/dW as projected in the third plane. This function is a shortcut assuming exactly three wire planes in the TPC, in which case the output plane is chosen as the one that is neither of the input planes.

Definition at line 687 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ThirdPlane(), and geo::WireReadoutGeom::ThirdPlane_dTdW().

   {
     return ThirdPlane_dTdW(pid1, slope1, pid2, slope2, ThirdPlane(pid1, pid2));
   }

double geo::WireReadoutGeom::ThirdPlaneSlope	(	PlaneID const &	pid1,
		double	slope1,
		PlaneID const &	pid2,
		double	slope2,
		PlaneID const &	output_plane
	)		const

inherited

Returns the slope on the third plane, given it in the other two.

Parameters

pid1	ID of the plane of the first slope
slope1	slope as seen on the first plane
pid2	ID of the plane of the second slope
slope2	slope as seen on the second plane
output_plane	ID of the plane on which to calculate the slope

Returns: the slope on the third plane, or -999. if slope would be infinity

Exceptions

cet::exception	(category: "GeometryCore") if different TPC
cet::exception	(category: "GeometryCore") if input planes match

Given a slope as projected in two planes, returns the slope as projected in the specified output plane. The slopes are defined in uniform units; they should be computed as distance ratios (or tangent of a geometrical angle; the formula is still valid using dt/dw directly in case of equal wire pitch in all planes and uniform drift velocity.

Definition at line 639 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ComputeThirdPlaneSlope(), and geo::WireReadoutGeom::Plane().

Referenced by trkf::CCTrackMaker::FillEndMatch(), trkf::CCTrackMaker::PlnMatch(), geo::WireReadoutGeom::ThirdPlaneSlope(), and geo::WireReadoutGeom::WireEndPoints().

   {
     CheckIndependentPlanesOnSameTPC(pid1, pid2, "ThirdPlaneSlope()");
 
     // We need the "wire coordinate direction" for each plane.  This is perpendicular to
     // the wire orientation.  PlaneGeo::PhiZ() defines the right orientation too.
     return ComputeThirdPlaneSlope(
       Plane(pid1).PhiZ(), slope1, Plane(pid2).PhiZ(), slope2, Plane(output_plane).PhiZ());
   }

double geo::WireReadoutGeom::ThirdPlaneSlope	(	PlaneID const &	pid1,
		double	slope1,
		PlaneID const &	pid2,
		double	slope2
	)		const

inherited

Returns the slope on the third plane, given it in the other two.

Parameters

pid1	ID of the plane of the first slope
slope1	slope as seen on the first plane
pid2	ID of the plane of the second slope
slope2	slope as seen on the second plane

Returns: the slope on the third plane, or -999. if slope would be infinity

Exceptions

cet::exception	(category: "GeometryCore") if different TPC
cet::exception	(category: "GeometryCore") if same plane
cet::exception	(category: "GeometryCore") if other than 3 planes

Given a slope as projected in two planes, returns the slope as projected in the third plane. This function is a shortcut assuming exactly three wire planes in the TPC, in which case the output plane is chosen as the one that is neither of the input planes.

Definition at line 654 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ThirdPlane(), and geo::WireReadoutGeom::ThirdPlaneSlope().

   {
     return ThirdPlaneSlope(pid1, slope1, pid2, slope2, ThirdPlane(pid1, pid2));
   }

double geo::WireReadoutGeom::ThirdPlaneSlope	(	PlaneID::PlaneID_t	plane1,
		double	slope1,
		PlaneID::PlaneID_t	plane2,
		double	slope2,
		TPCID const &	tpcid
	)		const

inlineinherited

Returns the slope on the third plane, given it in the other two.

Parameters

plane1	index of the plane of the first slope
slope1	slope as seen on the first plane
plane2	index of the plane of the second slope
slope2	slope as seen on the second plane
tpcid	TPC where the two planes belong

Returns: the slope on the third plane, or -999. if slope would be infinity

Exceptions

cet::exception	(category: "GeometryCore") if different TPC
cet::exception	(category: "GeometryCore") if same plane
cet::exception	(category: "GeometryCore") if other than 3 planes

Given a slope as projected in two planes, returns the slope as projected in the third plane.

Definition at line 451 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::ChannelsInTPCs(), geo::WireReadoutGeom::ComputeThirdPlane_dTdW(), geo::WireReadoutGeom::ComputeThirdPlaneSlope(), geo::WireReadoutGeom::Nchannels(), geo::WireReadoutGeom::Plane0Pitch(), geo::WireReadoutGeom::PlanePitch(), geo::WireReadoutGeom::ThirdPlane_dTdW(), and geo::WireReadoutGeom::ThirdPlaneSlope().

     {
       return ThirdPlaneSlope(PlaneID(tpcid, plane1), slope1, PlaneID(tpcid, plane2), slope2);
     }

std::vector< TPCID > geo::WireReadoutStandardGeom::TPCsetToTPCs ( readout::TPCsetID const & tpcsetid ) const

overridevirtual

Returns a list of ID of TPCs belonging to the specified TPC set.

Parameters

tpcsetid ID of the TPC set to convert into TPC IDs

Returns: the list of TPCs, empty if TPC set is invalid

In this mapping, TPC sets and TPCs are mapped one-to-one. The returned list contains always one entry, unless the specified TPC set ID is invalid, in which case the list is empty. Note that the check is performed on the validity of the TPC set ID, that does not necessarily imply that the TPC set specified by the ID actually exists. Check the existence of the TPC set first (HasTPCset()). Behaviour on valid, non-existent TPC set IDs is undefined.

Implements geo::WireReadoutGeom.

Definition at line 354 of file WireReadoutStandardGeom.cxx.

References geo::CryostatID::isValid.

   {
     std::vector<TPCID> IDs;
     if (tpcsetid.isValid) IDs.emplace_back(ConvertTPCsetToTPC(tpcsetid));
     return IDs;
   }

readout::TPCsetID geo::WireReadoutStandardGeom::TPCtoTPCset ( TPCID const & tpcid ) const

overridevirtual

Returns the ID of the TPC set the specified TPC belongs to.

Parameters

tpcid ID of the TPC

Returns: the ID of the corresponding TPC set, or invalid ID when tpcid is

In this mapping, TPC sets and TPCs are mapped one-to-one. The returned value mirrors the TPC ID in the readout space. If the TPC ID is not valid, an invalid TPC set ID is returned. Note that this check is performed on the validity of the TPC ID, that does not necessarily imply that the TPC specified by the ID actually exists.

Implements geo::WireReadoutGeom.

Definition at line 348 of file WireReadoutStandardGeom.cxx.

   {
     return ConvertTPCtoTPCset(tpcid);
   }

View_t geo::WireReadoutGeom::View ( raw::ChannelID_t const channel ) const

inherited

Returns the view (wire orientation) on the specified TPC channel.

Parameters

channel TPC channel ID

Returns: the type of signal on the specified channel, or geo::kUnknown

The view of the readout plane channel belongs to is returned, as in View(readout::ROPID const&) const.

Definition at line 518 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::ChannelToROP(), raw::InvalidChannelID, and geo::kUnknown.

Referenced by nnet::WaveformDenoiseTest::analyze(), nnet::NoiseWaveformDump::analyze(), nnet::RawWaveformDump::analyze(), nnet::RawWaveformClnSigDump::analyze(), apa::APAGeometryAlg::APAChannelsIntersect(), apa::APAGeometryAlg::APAGeometryAlg(), apa::APAGeometryAlg::APAView(), tca::FillWireHitRange(), apa::APAGeometryAlg::FirstChannelInView(), quad::GetPts2D(), geo::WireReadoutGeom::HasChannel(), apa::DisambigAlg::MakeCloseHits(), apa::APAGeometryAlg::NearestWireIDOnChan(), and cheat::BackTracker::TrackIdToSimIDEs_Ps().

   {
     return (channel == raw::InvalidChannelID) ? kUnknown : View(ChannelToROP(channel));
   }

View_t geo::WireReadoutGeom::View ( readout::ROPID const & ropid ) const

inherited

Returns the view of the channels in the specified readout plane.

Parameters

ropid readout plane ID

Returns: the type of signal on the specified ROP

Returns the view (wire orientation) on the channels of specified readout plane. If ropid is an invalid ID, geo::kUnknown is returned. If ropid is a valid ID (i.e. an ID whose isValid flag is set) that points to a non-existent readout plane, the result is undefined. Use HasROP() to check if the readout plane actually exists.

Definition at line 512 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::FirstWirePlaneInROP(), geo::kUnknown, geo::WireReadoutGeom::Plane(), and geo::PlaneGeo::View().

   {
     auto const pid = FirstWirePlaneInROP(ropid);
     return pid ? Plane(pid).View() : kUnknown;
   }

std::set<View_t> const& geo::WireReadoutGeom::Views ( ) const

inlineinherited

Returns a list of possible views in the detector.

Returns: the set of views

Definition at line 123 of file WireReadoutGeom.h.

References geo::WireReadoutGeom::allViews, geo::WireReadoutGeom::FirstPlane(), and geo::WireReadoutGeom::Plane().

Referenced by pma::PMAlgTracker::PMAlgTracker(), and geo::WireReadoutGeom::WireReadoutGeom().

123 { return allViews; }

geo::WireReadoutGeom::allViews

std::set< View_t > allViews

All views in the detector.

Definition: WireReadoutGeom.h:1090

std::set< View_t > geo::WireReadoutGeom::Views ( TPCID const & tpcid ) const

inherited

Definition at line 241 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::fPlanes.

   {
     auto it = fPlanes.find(tpcid);
     if (it == fPlanes.cend()) { return {}; }
     return ::Views(it->second);
   }

WireGeo const & geo::WireReadoutGeom::Wire ( WireID const & wireid ) const

inherited

Returns the specified wire.

Parameters

wireid ID of the wire

Returns: a constant reference to the specified wire

Exceptions

cet::exception if not found

The GetElement() method is overloaded and its return depends on the type of ID.

Definition at line 270 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::Plane(), and geo::PlaneGeo::Wire().

Referenced by lar_cluster3d::SnippetHit3DBuilder::DistanceFromPointToHitWire(), geo::WireReadoutGeom::GetElement(), geo::WireReadoutGeom::GetElementPtr(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire_(), lar_cluster3d::PrincipalComponentsAlg::PCAAnalysis_2D(), lar_cluster3d::PrincipalComponentsAlg::PCAAnalysis_calc2DDocas(), trkf::TrackKalmanFitter::setupInputStates(), apa::APAGeometryAlg::ThreeChanPos(), geo::WireReadoutGeom::WireEndPoints(), lar_cluster3d::SnippetHit3DBuilder::WireIDsIntersect(), and geo::WireReadoutGeom::WireIDsIntersect().

   {
     return Plane(wireid).Wire(wireid);
   }

double geo::WireReadoutGeom::WireAngleToVertical	(	View_t	view,
		TPCID const &	tpcid
	)		const

inherited

Returns the angle of the wires in the specified view from vertical.

Parameters

view	the view
tpcid	ID of the TPC

Returns: the angle [radians]

Exceptions

cet::exception ("GeometryCore" category) if no such view

The angle is defined as in WireGeo::ThetaZ().

This method assumes all wires in the view have the same angle (it queries for the first).

Deprecated:: This does not feel APA-ready

Definition at line 376 of file WireReadoutGeom.cxx.

References geo::PlaneGeo::ViewName().

Referenced by ShowerRecoTools::ShowerUnidirectiondEdx::CalculateElement(), util::GeometryUtilities::CalculatePitch(), util::GeometryUtilities::CalculatePitchPolar(), util::GeometryUtilities::Get2DangleFrom3D(), shower::EMShowerAlg::MakeShower(), util::GeometryUtilities::PitchInView(), lar_pandora::detector_functions::WireAngle(), and geo::WireReadoutGeom::WireEndPoints().

   {
     for (PlaneGeo const& plane : Iterate<PlaneGeo>(tpcid)) {
       if (plane.View() == view) return plane.ThetaZ();
     } // for
     throw cet::exception("WireReadoutGeom")
       << "WireAngleToVertical(): no view \"" << PlaneGeo::ViewName(view) << "\" (#" << ((int)view)
       << ") in " << std::string(tpcid);
   }

double geo::WireReadoutStandardGeom::WireCoordinate	(	double	YPos,
		double	ZPos,
		PlaneID const &	planeID
	)		const

overridevirtual

Returns the index of the wire nearest to the specified position.

Parameters

YPos	y coordinate on the wire plane
ZPos	z coordinate on the wire plane
planeID	ID of the plane

Returns: an index interpolation between the two nearest wires

See also: NearestWireID()

Respect to NearestWireID(), this method returns a real number, representing a continuous coordinate in the wire axis, with the round values corresponding to the actual wires.

The plane is required to be valid and exist in the detector. Otherwise, the behaviour is undefined.

Implements geo::WireReadoutGeom.

Definition at line 220 of file WireReadoutStandardGeom.cxx.

References geo::WireReadoutGeom::AccessElement(), fFirstWireProj, fOrthVectorsY, and fOrthVectorsZ.

Referenced by NearestWireID().

   {
     return YPos * AccessElement(fOrthVectorsY, planeID) +
            ZPos * AccessElement(fOrthVectorsZ, planeID) - AccessElement(fFirstWireProj, planeID);
   }

unsigned int geo::WireReadoutStandardGeom::WireCount ( PlaneID const & id ) const

inlineprivate

Retrieved the wire cound for the specified plane ID.

Definition at line 213 of file WireReadoutStandardGeom.h.

References geo::WireReadoutGeom::AccessElement(), and MaxTPCs().

Referenced by Nchannels(), and NearestWireID().

213 { return AccessElement(fWireCounts, id); }

geo::WireReadoutStandardGeom::fWireCounts

PlaneInfoMap_t< float > fWireCounts

Definition: WireReadoutStandardGeom.h:200

geo::WireReadoutGeom::AccessElement

T const & AccessElement(TPCInfoMap_t< T > const &map, TPCID const &id) const

Returns the specified element of the TPC map.

Definition: WireReadoutGeom.h:994

void geo::WireReadoutGeom::WireEndPoints	(	WireID const &	wireid,
		double *	xyzStart,
		double *	xyzEnd
	)		const

inherited

Fills two arrays with the coordinates of the wire end points.

Parameters

wireid	ID of the wire
xyzStart	(output) an array with the start coordinate
xyzEnd	(output) an array with the end coordinate

Exceptions

cet::exception wire not present

The starting point is the wire end with lower z coordinate.

Deprecated:: use the wire ID interface instead (but note that it does not sort the ends)

Definition at line 276 of file WireReadoutGeom.cxx.

References util::abs(), geo::Iterable< details::ReadoutIterationPolicy, details::ToReadoutGeometryElement >::end(), and geo::WireReadoutGeom::start().

Referenced by lar_cluster3d::StandardHit3DBuilder::DistanceFromPointToHitWire(), geo::WireReadoutGeom::GetElement(), quad::GetPts2D(), apa::DisambigAlg::TrivialDisambig(), and geo::WireReadoutGeom::WireIDsIntersect().

   {
     auto const [start, end] = WireEndPoints(wireid);
 
     xyzStart[0] = start.X();
     xyzStart[1] = start.Y();
     xyzStart[2] = start.Z();
     xyzEnd[0] = end.X();
     xyzEnd[1] = end.Y();
     xyzEnd[2] = end.Z();
 
     if (xyzEnd[2] < xyzStart[2]) {
       // Ensure that "End" has higher z-value than "Start"
       std::swap(xyzStart[0], xyzEnd[0]);
       std::swap(xyzStart[1], xyzEnd[1]);
       std::swap(xyzStart[2], xyzEnd[2]);
     }
     if (xyzEnd[1] < xyzStart[1] && std::abs(xyzEnd[2] - xyzStart[2]) < 0.01) {
       // If wire is vertical ensure that "End" has higher y-value than "Start"
       std::swap(xyzStart[0], xyzEnd[0]);
       std::swap(xyzStart[1], xyzEnd[1]);
       std::swap(xyzStart[2], xyzEnd[2]);
     }
   }

Segment geo::WireReadoutGeom::WireEndPoints ( WireID const & wireid ) const

inlineinherited

Returns a segment whose ends are the wire end points.

Parameters

wireid ID of the wire

Returns: a segment whose ends are the wire end points

Exceptions

cet::exception wire not present

The start and end are assigned as returned from the geo::WireGeo object. The rules for this assignment are documented in that class.

Deprecated:: use the wire ID interface instead (but note that it does not sort the ends)

Definition at line 260 of file WireReadoutGeom.h.

References geo::WireGeo::GetEnd(), geo::WireGeo::GetStart(), geo::WireReadoutGeom::ThirdPlane(), geo::WireReadoutGeom::ThirdPlaneSlope(), geo::WireReadoutGeom::Wire(), geo::WireReadoutGeom::WireAngleToVertical(), and geo::WireReadoutGeom::WireIDsIntersect().

     {
       WireGeo const& wire = Wire(wireid);
       return {wire.GetStart(), wire.GetEnd()};
     }

bool geo::WireReadoutGeom::WireIDsIntersect	(	WireID const &	wid1,
		WireID const &	wid2,
		Point_t &	intersection
	)		const

inherited

Computes the intersection between two wires.

Parameters

	wid1	ID of the first wire
	wid2	ID of the other wire
[out]	intersection	the intersection point (global coordinates)

Returns: whether an intersection was found inside the TPC the wires belong

See also: geo::WiresIntersection(), geo::LineClosestPoint()

The wires identified by wid1 and wid2 are intersected, and the 3D intersection point is written into the intersection parameter. The "intersection" point is actually the point belonging to the first wire (wid2) which is the closest (in Euclidean 3D metric) to the second wire.

The intersection is computed only if the wires belong to different planes of the same TPC. If that is not the case (i.e. they belong to different TPC or cryostat, or if they belong to the same plane), false is returned and intersection is set with all components to infinity (std::numeric_limits<>::infinity()).

When the intersection is computed, it is always stored in the intersection output parameter. Return value is true if this intersection lies within the physical boundaries first wire, while it is instead false if it lies on the extrapolation of the wire direction, but not within the wire physical extension.

To test that the result is not infinity (nor NaN), use geo::vect::isfinite(intersection) etc.

Note: If geo::WireGeo objects are already available, using instead the free function geo::WiresIntersection() or the method geo::WireGeo::IntersectionWith() is faster (and recommended). For purely geometric intersection, geo::LineClosestPoint() is also available.

Definition at line 344 of file WireReadoutGeom.cxx.

References util::abs(), geo::WireGeo::HalfL(), geo::IntersectionPointAndOffsets< Point >::offset1, geo::IntersectionPointAndOffsets< Point >::offset2, geo::IntersectionPointAndOffsets< Point >::point, geo::WireReadoutGeom::Wire(), geo::WireReadoutGeom::WireIDIntersectionCheck(), and geo::WiresIntersectionAndOffsets().

Referenced by pma::ProjectionMatchingAlg::addEndpointRef_(), apa::APAGeometryAlg::APAChannelsIntersect(), geo::WireReadoutGeom::ChannelsIntersect(), apa::DisambigAlg::CompareViews(), shower::EMShowerAlg::Construct3DPoint_(), trkf::CCTrackMaker::FillEndMatch(), tca::FillWireIntersections(), tca::Find3DVertices(), reco3d::IntersectionCache::ISect(), lar_cluster3d::StandardHit3DBuilder::makeDeadChannelPair(), lar_cluster3d::SnippetHit3DBuilder::makeDeadChannelPair(), lar_cluster3d::StandardHit3DBuilder::makeHitPair(), tca::Match2Planes(), trkf::CCTrackMaker::PlnMatch(), trkf::TrackTrajectoryAlg::ShortTrackTrajectory(), apa::APAGeometryAlg::ThreeChanPos(), and geo::WireReadoutGeom::WireEndPoints().

   {
     // This is not a real 3D intersection: the wires do not cross, since they are required
     // to belong to two different planes.
     //
     // We take the point on the first wire which is closest to the other one.
     static_assert(std::numeric_limits<decltype(intersection.X())>::has_infinity,
                   "the vector coordinate type can't represent infinity!");
     constexpr auto infinity = std::numeric_limits<decltype(intersection.X())>::infinity();
 
     if (!WireIDIntersectionCheck(wid1, wid2)) {
       intersection = {infinity, infinity, infinity};
       return false;
     }
 
     WireGeo const& wire1 = Wire(wid1);
     WireGeo const& wire2 = Wire(wid2);
 
     // Distance of the intersection point from the center of the two wires:
     IntersectionPointAndOffsets<Point_t> intersectionAndOffset =
       WiresIntersectionAndOffsets(wire1, wire2);
     intersection = intersectionAndOffset.point;
 
     return std::abs(intersectionAndOffset.offset1) <= wire1.HalfL() &&
            std::abs(intersectionAndOffset.offset2) <= wire2.HalfL();
   }

std::optional< WireIDIntersection > geo::WireReadoutGeom::WireIDsIntersect	(	WireID const &	wid1,
		WireID const &	wid2
	)		const

inherited

Computes the intersection between two wires.

Parameters

wid1	ID of the first wire
wid2	ID of the other wire

Returns: whether an intersection was found within the TPC

The "intersection" refers to the projection of the wires into the same $ x = 0 $ plane. Wires are assumed to have at most one intersection. If wires are parallel, widIntersect will have the two components set to infinity (std::numeric_limits<>::infinity()) and the TPC number set to invalid (geo::TPCID::InvalidID). Also, false is returned. If the intersection is outside the TPC, false is also returned, but the widIntersect will contain the coordinates of that intersection. The TPC number is still set to invalid, although the intersection might belong to a valid TPC somewhere else.

Deprecated:: This method uses arbitrary assumptions and should not be used. Use the interface returning a full vector instead.

Definition at line 302 of file WireReadoutGeom.cxx.

References geo::vect::cross(), geo::WireReadoutGeom::finish(), geo::IntersectLines(), lar::util::PointWithinSegments(), geo::WireReadoutGeom::start(), geo::TPCID::TPC, geo::WireIDIntersection::TPC, geo::WireReadoutGeom::WireEndPoints(), geo::WireReadoutGeom::WireIDIntersectionCheck(), geo::Y, geo::WireIDIntersection::y, geo::Z, and geo::WireIDIntersection::z.

   {
     if (!WireIDIntersectionCheck(wid1, wid2)) { return std::nullopt; }
 
     // get the endpoints to see if wires intersect
     Segment const w1 = WireEndPoints(wid1);
     Segment const w2 = WireEndPoints(wid2);
 
     // TODO extract the coordinates in the right way;
     // is it any worth, since then the result is in (y, z), whatever it means?
     WireIDIntersection result;
     bool const cross = IntersectLines(start(w1).Y(),
                                       start(w1).Z(),
                                       finish(w1).Y(),
                                       finish(w1).Z(),
                                       start(w2).Y(),
                                       start(w2).Z(),
                                       finish(w2).Y(),
                                       finish(w2).Z(),
                                       result.y,
                                       result.z);
     if (!cross) { return std::nullopt; }
 
     bool const within = lar::util::PointWithinSegments(start(w1).Y(),
                                                        start(w1).Z(),
                                                        finish(w1).Y(),
                                                        finish(w1).Z(),
                                                        start(w2).Y(),
                                                        start(w2).Z(),
                                                        finish(w2).Y(),
                                                        finish(w2).Z(),
                                                        result.y,
                                                        result.z);
 
     if (!within) { return std::nullopt; }
 
     result.TPC = wid1.TPC;
     return result;
   }

readout::ROPID geo::WireReadoutStandardGeom::WirePlaneToROP ( PlaneID const & planeid ) const

overridevirtual

Returns the ID of the ROP planeid belongs to, or invalid if none.

Parameters

planeid ID of the plane

Returns: the ID of the corresponding ROP, or invalid ID when planeid is

In this mapping, readout planes and wire planes are mapped one-to-one. The returned value mirrors the plane ID in the readout space. If the plane ID is not valid, an invalid readout plane ID is returned. Note that this check is performed on the validity of the plane ID, that does not necessarily imply that the plane specified by the ID actually exists.

Implements geo::WireReadoutGeom.

Definition at line 400 of file WireReadoutStandardGeom.cxx.

Referenced by ChannelToROP().

   {
     return ConvertWirePlaneToROP(planeid);
   }

WireGeo const * geo::WireReadoutGeom::WirePtr ( WireID const & wireid ) const

inherited

Returns the specified wire.

Parameters

wireid wire ID

Returns: a constant pointer to the specified wire, or nullptr if none

The GetElementPtr() method is overloaded and its return depends on the type of ID.

Definition at line 263 of file WireReadoutGeom.cxx.

References geo::WireReadoutGeom::PlanePtr(), and geo::PlaneGeo::WirePtr().

Referenced by geo::details::getElementPtr(), geo::WireReadoutGeom::GetElementPtr(), geo::WireReadoutGeom::HasElement(), and lar_cluster3d::MinSpanTreeAlg::MinSpanTreeAlg().

   {
     PlaneGeo const* pPlane = PlanePtr(wireid);
     return pPlane ? pPlane->WirePtr(wireid) : nullptr;
   }

Member Data Documentation

PlaneInfoMap_t<raw::ChannelID_t> geo::WireReadoutGeom::fFirstChannelInNextPlane

protectedinherited

Definition at line 1071 of file WireReadoutGeom.h.

Referenced by ChannelToWire(), geo::WireReadoutGeom::FirstChannelInThisPlane(), SignalTypeForChannelImpl(), and WireReadoutStandardGeom().

PlaneInfoMap_t<raw::ChannelID_t> geo::WireReadoutGeom::fFirstChannelInThisPlane

protectedinherited

Definition at line 1070 of file WireReadoutGeom.h.

Referenced by ChannelToWire(), geo::WireReadoutGeom::FirstChannelInNextPlane(), SignalTypeForChannelImpl(), and WireReadoutStandardGeom().

PlaneInfoMap_t<float> geo::WireReadoutStandardGeom::fFirstWireProj

private

Distance (0,0,0) to first wire along orth vector per plane per TPC

Definition at line 194 of file WireReadoutStandardGeom.h.

Referenced by WireCoordinate(), and WireReadoutStandardGeom().

unsigned int geo::WireReadoutStandardGeom::fNchannels

private

number of channels in the detector

Definition at line 189 of file WireReadoutStandardGeom.h.

Referenced by Nchannels(), SignalTypeForChannelImpl(), and WireReadoutStandardGeom().

unsigned int geo::WireReadoutStandardGeom::fNcryostat

private

number of cryostats in the detector

Definition at line 188 of file WireReadoutStandardGeom.h.

Referenced by ChannelToWire(), and WireReadoutStandardGeom().

TPCInfoMap_t<unsigned int> geo::WireReadoutStandardGeom::fNPlanes

private

Number of planes in each TPC - for range checking after calculation

Definition at line 202 of file WireReadoutStandardGeom.h.

Referenced by MaxROPs(), NROPs(), SignalTypeForChannelImpl(), and WireReadoutStandardGeom().

std::vector<unsigned int> geo::WireReadoutStandardGeom::fNTPC

private

number of TPCs in each cryostat

Definition at line 191 of file WireReadoutStandardGeom.h.

Referenced by ChannelToWire(), MaxTPCs(), NTPCsets(), SignalTypeForChannelImpl(), and WireReadoutStandardGeom().

PlaneInfoMap_t<float> geo::WireReadoutStandardGeom::fOrthVectorsY

private

Unit vectors orthogonal to wires in.

Definition at line 196 of file WireReadoutStandardGeom.h.

Referenced by WireCoordinate(), and WireReadoutStandardGeom().

PlaneInfoMap_t<float> geo::WireReadoutStandardGeom::fOrthVectorsZ

private

each plane - stored as 2 components to avoid having to invoke any bulky TObjects / CLHEP vectors etc

Definition at line 197 of file WireReadoutStandardGeom.h.

Referenced by WireCoordinate(), and WireReadoutStandardGeom().

PlaneInfoMap_t<unsigned int> geo::WireReadoutStandardGeom::fPlaneBaselines

private

The number of wires in all the tpcs and planes up to this one in the heirachy

Definition at line 204 of file WireReadoutStandardGeom.h.

Referenced by FirstChannelInROP(), PlaneWireToChannel(), and WireReadoutStandardGeom().

std::set<PlaneID> geo::WireReadoutStandardGeom::fPlaneIDs

private

vector of the PlaneIDs present in the detector

Definition at line 193 of file WireReadoutStandardGeom.h.

Referenced by PlaneIDs(), and WireReadoutStandardGeom().

raw::ChannelID_t geo::WireReadoutStandardGeom::fTopChannel

private

book keeping highest channel #

Definition at line 190 of file WireReadoutStandardGeom.h.

Referenced by ChannelToWire(), and WireReadoutStandardGeom().

std::set<View_t> geo::WireReadoutStandardGeom::fViews

private

vector of the views present in the detector

Definition at line 192 of file WireReadoutStandardGeom.h.

PlaneInfoMap_t<float> geo::WireReadoutStandardGeom::fWireCounts

private

Number of wires in each plane - for range checking after calculation

Definition at line 200 of file WireReadoutStandardGeom.h.

Referenced by WireReadoutStandardGeom().

PlaneInfoMap_t<unsigned int> geo::WireReadoutStandardGeom::fWiresPerPlane

private

The number of wires in this plane in the heirachy

Definition at line 207 of file WireReadoutStandardGeom.h.

Referenced by WireReadoutStandardGeom().

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

larcorealg/v10_00_02/source/larcorealg/Geometry/WireReadoutStandardGeom.h
larcorealg/v10_00_02/source/larcorealg/Geometry/WireReadoutStandardGeom.cxx

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation