LArSoft: geo::GeometryCore Class Reference

geo::GeometryData_t::auxDets

GeometryData_t fGeoData

The detector description data.

Definition: GeometryCore.h:5391

AuxDetList_t auxDets

The auxiliary detectors.

Definition: GeometryCore.h:1392

AuxDetList_t const& geo::GeometryCore::AuxDets ( ) const

inlineprotected

Return the interfal auxiliary detectors list.

Definition at line 5347 of file GeometryCore.h.

5347 { return fGeoData.auxDets; }

geo::GeometryData_t::auxDets

GeometryData_t fGeoData

The detector description data.

Definition: GeometryCore.h:5391

AuxDetList_t auxDets

The auxiliary detectors.

Definition: GeometryCore.h:1392

cryostat_iterator geo::GeometryCore::begin_cryostat ( ) const

inline

Returns an iterator pointing to the first cryostat.

Definition at line 2124 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), and IterateCryostatIDs().

2125 { return cryostat_iterator(this, cryostat_iterator::begin_pos); }

geo::GeometryCore::cryostat_iterator

geo::cryostat_iterator cryostat_iterator

Forward-iterator browsing all cryostats in the detector.

Definition: GeometryCore.h:1531

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

cryostat_id_iterator geo::GeometryCore::begin_cryostat_id ( ) const

inline

Returns an iterator pointing to the first cryostat ID.

Definition at line 2116 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by end_cryostat().

2117 { return cryostat_id_iterator(this, cryostat_id_iterator::begin_pos); }

geo::GeometryCore::cryostat_id_iterator

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::cryostat_id_iterator cryostat_id_iterator

Forward-iterator browsing all cryostat IDs in the detector.

Definition: GeometryCore.h:1510

plane_iterator geo::GeometryCore::begin_plane ( ) const

inline

Returns an iterator pointing to the first plane in the detector.

Definition at line 3035 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by IteratePlaneIDs(), and IteratePlanes().

3036 { return plane_iterator(this, plane_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

plane_iterator geo::GeometryCore::begin_plane ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing to the first plane in the specified cryostat.

Definition at line 3044 of file GeometryCore.h.

3045 { return plane_iterator(this, GetBeginPlaneID(ID)); }

geo::PlaneID GetBeginPlaneID(geo::CryostatID const &id) const

Returns the ID of the first plane of the specified cryostat.

Definition: GeometryCore.h:2989

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

plane_iterator geo::GeometryCore::begin_plane ( geo::TPCID const & ID ) const

inline

Returns an iterator pointing to the first plane in the specified TPC.

Definition at line 3053 of file GeometryCore.h.

3054 { return plane_iterator(this, GetBeginPlaneID(ID)); }

geo::PlaneID GetBeginPlaneID(geo::CryostatID const &id) const

Returns the ID of the first plane of the specified cryostat.

Definition: GeometryCore.h:2989

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

plane_id_iterator geo::GeometryCore::begin_plane_id ( ) const

inline

Returns an iterator pointing to the first plane ID in the detector.

Definition at line 3007 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by end_plane(), and IteratePlaneIDs().

3008 { return plane_id_iterator(this, plane_id_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

plane_id_iterator geo::GeometryCore::begin_plane_id ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing to the first plane ID in the specified cryostat.

Definition at line 3016 of file GeometryCore.h.

3017 { return plane_id_iterator(this, GetBeginPlaneID(ID)); }

geo::PlaneID GetBeginPlaneID(geo::CryostatID const &id) const

Returns the ID of the first plane of the specified cryostat.

Definition: GeometryCore.h:2989

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

plane_id_iterator geo::GeometryCore::begin_plane_id ( geo::TPCID const & ID ) const

inline

Returns an iterator pointing to the first plane ID in the specified TPC.

Definition at line 3026 of file GeometryCore.h.

3027 { return plane_id_iterator(this, GetBeginPlaneID(ID)); }

geo::PlaneID GetBeginPlaneID(geo::CryostatID const &id) const

Returns the ID of the first plane of the specified cryostat.

Definition: GeometryCore.h:2989

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

ROP_id_iterator geo::GeometryCore::begin_ROP_id ( ) const

inline

Returns an iterator pointing to the first ROP ID in the detector.

Definition at line 5086 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by end_ROP_id(), and IterateROPIDs().

5087 { return ROP_id_iterator(this, ROP_id_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

ROP_id_iterator geo::GeometryCore::begin_ROP_id ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing to the first readout plane ID in the specified cryostat.

Definition at line 5095 of file GeometryCore.h.

5096 { return ROP_id_iterator(this, GetBeginROPID(ID)); }

geo::GeometryCore::GetBeginROPID

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

readout::ROPID GetBeginROPID(geo::CryostatID const &id) const

Returns the ID of the first readout plane of the specified cryostat.

Definition: GeometryCore.h:5068

ROP_id_iterator geo::GeometryCore::begin_ROP_id ( readout::TPCsetID const & ID ) const

inline

Returns an iterator pointing to the first readout plane ID in the specified TPC set.

Definition at line 5105 of file GeometryCore.h.

5106 { return ROP_id_iterator(this, GetBeginROPID(ID)); }

geo::GeometryCore::GetBeginROPID

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

readout::ROPID GetBeginROPID(geo::CryostatID const &id) const

Returns the ID of the first readout plane of the specified cryostat.

Definition: GeometryCore.h:5068

TPC_iterator geo::GeometryCore::begin_TPC ( ) const

inline

Returns an iterator pointing to the first TPC in the detector.

Definition at line 2534 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), IterateTPCIDs(), IterateTPCs(), and sim::MCRecoPart::MCRecoPart().

2535 { return TPC_iterator(this, TPC_iterator::begin_pos); }

geo::TPC_iterator TPC_iterator

Forward-iterator browsing all TPCs in the detector.

Definition: GeometryCore.h:1577

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

TPC_iterator geo::GeometryCore::begin_TPC ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing to the first TPC in the detector.

Definition at line 2542 of file GeometryCore.h.

2543 { return TPC_iterator(this, GetBeginTPCID(cid)); }

geo::GeometryCore::GetBeginTPCID

geo::TPC_iterator TPC_iterator

Forward-iterator browsing all TPCs in the detector.

Definition: GeometryCore.h:1577

geo::TPCID GetBeginTPCID(geo::CryostatID const &id) const

Returns the ID of the first TPC in the specified cryostat.

Definition: GeometryCore.h:2506

TPC_id_iterator geo::GeometryCore::begin_TPC_id ( ) const

inline

Returns an iterator pointing to the first TPC ID in the detector.

Definition at line 2516 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by pma::PMAlgTracker::build(), end_TPC(), IterateTPCIDs(), and tss::TrackShowerHits::produce().

2517 { return TPC_id_iterator(this, TPC_id_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::TPC_id_iterator TPC_id_iterator

Forward-iterator browsing all TPC IDs in the detector.

Definition: GeometryCore.h:1555

TPC_id_iterator geo::GeometryCore::begin_TPC_id ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing to the first TPC ID in the specified cryostat.

Definition at line 2525 of file GeometryCore.h.

2526 { return TPC_id_iterator(this, GetBeginTPCID(cid)); }

geo::GeometryCore::GetBeginTPCID

geo::TPC_id_iterator TPC_id_iterator

Forward-iterator browsing all TPC IDs in the detector.

Definition: GeometryCore.h:1555

geo::TPCID GetBeginTPCID(geo::CryostatID const &id) const

Returns the ID of the first TPC in the specified cryostat.

Definition: GeometryCore.h:2506

TPCset_id_iterator geo::GeometryCore::begin_TPCset_id ( ) const

inline

Returns an iterator pointing to the first TPC set ID in the detector.

Definition at line 4849 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by end_TPCset_id(), and IterateTPCsetIDs().

4850 { return TPCset_id_iterator(this, TPCset_id_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::GeometryCore::GetBeginTPCsetID

details::TPCset_id_iterator_base< readout::TPCsetID > TPCset_id_iterator

Forward iterator browsing all TPC sets in the detector.

Definition: GeometryCore.h:1288

TPCset_id_iterator geo::GeometryCore::begin_TPCset_id ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing to the first TPC set ID in the specified cryostat.

Definition at line 4858 of file GeometryCore.h.

4859 { return TPCset_id_iterator(this, GetBeginTPCsetID(cid)); }

readout::TPCsetID GetBeginTPCsetID(geo::CryostatID const &id) const

Returns the ID of the first TPC set in the specified cryostat.

Definition: GeometryCore.h:4839

details::TPCset_id_iterator_base< readout::TPCsetID > TPCset_id_iterator

Forward iterator browsing all TPC sets in the detector.

Definition: GeometryCore.h:1288

wire_iterator geo::GeometryCore::begin_wire ( ) const

inline

Returns an iterator pointing to the first wire in the detector.

Definition at line 3455 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by IterateWireIDs(), and IterateWires().

3456 { return wire_iterator(this, wire_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::GeometryCore::begin_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::begin_wire ( geo::CryostatID const & id ) const

inline

Returns an iterator pointing to the first wire in specified cryostat.

Definition at line 3463 of file GeometryCore.h.

3464 { return wire_iterator(begin_wire_id(id)); }

wire_id_iterator begin_wire_id() const

Returns an iterator pointing to the first wire ID in the detector.

Definition: GeometryCore.h:3422

geo::GeometryCore::begin_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::begin_wire ( geo::TPCID const & id ) const

inline

Returns an iterator pointing to the first wire in specified TPC.

Definition at line 3471 of file GeometryCore.h.

3472 { return wire_iterator(begin_wire_id(id)); }

wire_id_iterator begin_wire_id() const

Returns an iterator pointing to the first wire ID in the detector.

Definition: GeometryCore.h:3422

geo::GeometryCore::begin_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::begin_wire ( geo::PlaneID const & id ) const

inline

Returns an iterator pointing to the first wire in specified plane.

Definition at line 3479 of file GeometryCore.h.

3480 { return wire_iterator(begin_wire_id(id)); }

wire_id_iterator begin_wire_id() const

Returns an iterator pointing to the first wire ID in the detector.

Definition: GeometryCore.h:3422

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_id_iterator geo::GeometryCore::begin_wire_id ( ) const

inline

Returns an iterator pointing to the first wire ID in the detector.

Definition at line 3422 of file GeometryCore.h.

References geo::details::geometry_iterator_types::begin_pos.

Referenced by end_wire(), IterateWireIDs(), lariov::SIOVElectronicsCalibProvider::Reconfigure(), and lariov::DetPedestalRetrievalAlg::Reconfigure().

3423 { return wire_id_iterator(this, wire_id_iterator::begin_pos); }

static constexpr BeginPos_t begin_pos

Definition: GeometryCore.h:119

geo::GeometryCore::GetBeginWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::begin_wire_id ( geo::CryostatID const & id ) const

inline

Returns an iterator pointing to the first wire ID in specified cryostat.

Definition at line 3430 of file GeometryCore.h.

3431 { return wire_id_iterator(this, GetBeginWireID(id)); }

geo::WireID GetBeginWireID(geo::CryostatID const &id) const

Returns the ID of the first wire in the specified cryostat.

Definition: GeometryCore.h:3395

geo::GeometryCore::GetBeginWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::begin_wire_id ( geo::TPCID const & id ) const

inline

Returns an iterator pointing to the first wire ID in specified TPC.

Definition at line 3439 of file GeometryCore.h.

3440 { return wire_id_iterator(this, GetBeginWireID(id)); }

geo::WireID GetBeginWireID(geo::CryostatID const &id) const

Returns the ID of the first wire in the specified cryostat.

Definition: GeometryCore.h:3395

geo::GeometryCore::GetBeginWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::begin_wire_id ( geo::PlaneID const & id ) const

inline

Returns an iterator pointing to the first wire ID in specified plane.

Definition at line 3447 of file GeometryCore.h.

3448 { return wire_id_iterator(this, GetBeginWireID(id)); }

geo::WireID GetBeginWireID(geo::CryostatID const &id) const

Returns the ID of the first wire in the specified cryostat.

Definition: GeometryCore.h:3395

geo::GeometryCore::fChannelMapAlg

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

geo::ChannelMapAlg const* geo::GeometryCore::ChannelMap ( ) const

inlineprotected

Returns the object handling the channel map.

Definition at line 5335 of file GeometryCore.h.

Referenced by geo::Geometry::InitializeChannelMap().

5336 { return fChannelMapAlg.get(); }

std::shared_ptr< const geo::ChannelMapAlg > fChannelMapAlg

Object containing the channel to wire mapping.

Definition: GeometryCore.h:5400

bool geo::GeometryCore::ChannelsIntersect	(	raw::ChannelID_t	c1,
		raw::ChannelID_t	c2,
		double &	y,
		double &	z
	)		const

Returns an intersection point of two channels.

Parameters

c1	one channel ID
c2	the other channel ID
y	(output) y coordinate of the intersection
z	(output) z coordinate of the intersection

Returns: whether a intersection point was found

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 1335 of file GeometryCore.cxx.

References ChannelToWire(), WireIDsIntersect(), geo::WireIDIntersection::y, and geo::WireIDIntersection::z.

Referenced by vertex::FeatureVertexFinder::Find3dVtxFrom2dClusterVtxCand(), evd::TWQMultiTPCProjectionView::FindEndPoint(), evd::TWQMultiTPCProjectionView::FindLineLength(), cmtool::CFAlgoStartTimeCompat::Float(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), util::GeometryUtilities::GetProjectedPoint(), util::GeometryUtilities::GetYZ(), vertex::VertexFinder2D::produce(), and shwf::ShowerReco::produce().

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

const AuxDetGeo & geo::GeometryCore::ChannelToAuxDet	(	std::string const &	auxDetName,
		uint32_t const &	channel
	)		const

Returns the number of auxiliary detectors.

Todo:: use a AutDetID_t instead of unsigned int?

This method returns the total number of scintillator paddles (Auxiliary Detectors aka AuxDet) outside of the cryostat

Todo:: Change return type to size_t

Definition at line 514 of file GeometryCore.cxx.

References AuxDet(), AuxDets(), and fChannelMapAlg.

   {
     size_t adIdx = fChannelMapAlg->ChannelToAuxDet(AuxDets(), auxDetName, channel);
     return this->AuxDet(adIdx);
   }

const AuxDetSensitiveGeo & geo::GeometryCore::ChannelToAuxDetSensitive	(	std::string const &	auxDetName,
		uint32_t const &	channel
	)		const

Returns the number of auxiliary detectors.

Todo:: use a AutDetID_t instead of unsigned int?

This method returns the total number of scintillator paddles (Auxiliary Detectors aka AuxDet) outside of the cryostat

Todo:: Change return type to size_t

Definition at line 522 of file GeometryCore.cxx.

References AuxDet(), AuxDets(), fChannelMapAlg, and geo::AuxDetGeo::SensitiveVolume().

   {
     auto idx = fChannelMapAlg->ChannelToSensitiveAuxDet(AuxDets(), auxDetName, channel);
     return this->AuxDet(idx.first).SensitiveVolume(idx.second);
   }

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

Returns the ID of the ROP the channel belongs to

Exceptions

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

Definition at line 1174 of file GeometryCore.cxx.

References fChannelMapAlg, and WireCoordinate().

Referenced by reco3d::TripletFinder::FillBadMap(), reco3d::TripletFinder::FillHitMap(), and View().

   {
     return fChannelMapAlg->ChannelToROP(channel);
   } // GeometryCore::ChannelToROP()

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

Returns a list of wires connected to the specified TPC channel.

Parameters

channel TPC channel ID

Returns: vector containing the ID of all the connected wires

Exceptions

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

Definition at line 1168 of file GeometryCore.cxx.

References fChannelMapAlg.

Referenced by hit::MCHitAnaExample::analyze(), MCBTDemo::analyze(), detsim::WienerFilterAna::analyze(), apa::APAGeometryAlg::APAChannelsIntersect(), lar_cluster3d::StandardHit3DBuilder::BuildChannelStatusVec(), ChannelsIntersect(), apa::APAGeometryAlg::ChanSegsPerSide(), apa::DisambigAlg::CompareViews(), lar_cluster3d::MinSpanTreeAlg::configure(), reco3d::TripletFinder::FillHitMap(), evd::RawDataDrawer::FillQHisto(), evd::RecoBaseDrawer::FillQHisto(), lar_pandora::PFParticleHitDumper::FillRecoWires(), evd::RecoBaseDrawer::FillTQHisto(), evd::RecoBaseDrawer::FillTQHistoDP(), reco_tool::CandHitStandard::findHitCandidates(), reco_tool::CandHitDerivative::findHitCandidates(), reco_tool::CandHitMorphological::findHitCandidates(), hit::HitCheater::FindHitsOnChannel(), cluster::SmallClusterFinderAlg::GetPlaneAndTPC(), corner::CornerFinderAlg::GrabWires(), apa::APAGeometryAlg::Init(), hit::DisambigCheater::InitHitToWids(), reco3d::IntersectionCache::ISect(), apa::APAGeometryAlg::LineSegChanIntersect(), apa::DisambigAlg::MakeCloseHits(), sim::MCRecoEdep::MakeMCEdep(), apa::APAGeometryAlg::NearestWireIDOnChan(), btutil::MCBTAlg::ProcessSimChannel(), hit::TTHitFinder::produce(), hit::FFTHitFinder::produce(), hit::RawHitFinder::produce(), hit::GausHitFinder::produce(), hit::DPRawHitFinder::produce(), hit::RFFHitFinderAlg::Run(), hit::CCHitFinderAlg::RunCCHitFinder(), evd::RawDataDrawer::RunOperation(), img::DataProviderAlg::setWireDriftData(), apa::APAGeometryAlg::ThreeChanPos(), apa::DisambigAlg::TrivialDisambig(), and evd::RecoBaseDrawer::Wire2D().

   {
     return fChannelMapAlg->ChannelToWire(channel);
   }

void geo::GeometryCore::CheckIndependentPlanesOnSameTPC	(	geo::PlaneID const &	pid1,
		geo::PlaneID const &	pid2,
		const char *	caller
	)

staticprivate

Throws an exception ("GeometryCore" category) unless pid1 and pid2 are on different planes of the same TPC (ID validity is not checked)

Definition at line 1606 of file GeometryCore.cxx.

References geo::TPCID::asTPCID().

Referenced by ThirdPlane(), ThirdPlane_dTdW(), and ThirdPlaneSlope().

   {
     if(pid1.asTPCID() != pid2.asTPCID()) {
       throw cet::exception("GeometryCore")
         << caller << " needs two planes on the same TPC (got "
         << std::string(pid1) << " and " << std::string(pid2) << ")\n";
     }
     if(pid1 == pid2) { // was: return 999;
       throw cet::exception("GeometryCore")
         << caller << " needs two different planes, got "
         << std::string(pid1) << " twice\n";
     }
   } // GeometryCore::CheckIndependentPlanesOnSameTPC()

void geo::GeometryCore::ClearGeometry ( )

private

Deletes the detector geometry structures.

Definition at line 133 of file GeometryCore.cxx.

References AuxDets(), evd::details::begin(), Cryostats(), and evd::details::end().

Referenced by LoadGeometryFile(), and ~GeometryCore().

                                    {
     
     Cryostats().clear();
     
     // auxiliary detectors
     std::for_each(AuxDets().begin(), AuxDets().end(),
       std::default_delete<AuxDetGeo>());
     AuxDets().clear();
     
   } // GeometryCore::ClearGeometry()

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

static

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 1724 of file GeometryCore.cxx.

References ComputeThirdPlaneSlope().

Referenced by ThirdPlane_dTdW().

   {
     // 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);
   } // GeometryCore::ComputeThirdPlane_dTdW()

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

static

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 1699 of file GeometryCore.cxx.

Referenced by ComputeThirdPlane_dTdW(), ThirdPlane_dTdW(), and 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;
   } // GeometryCore::ComputeThirdPlaneSlope()

CryostatGeo const & geo::GeometryCore::Cryostat ( geo::CryostatID const & cryoid ) const

Returns the specified cryostat.

Parameters

cstat	number of cryostat
cryoid	cryostat ID

Returns: a constant reference to the specified cryostat

Exceptions

cet::exception (GeometryCore category) if cryostat not present

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

Todo:: Make the cryostat number mandatory (as CryostatID)

Definition at line 280 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, and CryostatPtr().

Referenced by trkf::SpacePointAna::bookHistograms(), pma::ProjectionMatchingAlg::buildShowerSeg(), lar_pandora::LArPandoraOutput::CalculateT0(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadoutGeometry::Construct(), ems::MultiEMShowers::convCluster(), lar_pandora::LArPandoraInput::CreatePandoraHits2D(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), CryostatBoundaries(), CryostatHalfHeight(), CryostatHalfWidth(), CryostatLength(), hit::DisambigCheater::DisambigCheater(), evd::RecoBaseDrawer::DrawTrack2D(), cluster::EndPointAlg::EndPoint(), cluster::HoughBaseAlg::FastTransform(), trkf::SpacePointAlg::fillComplexSpacePoint(), cosmic::BeamFlashTrackMatchTaggerAlg::FillFlashProperties(), opdet::FlashHypothesisAnaAlg::FillOpDetPositions(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), vertex::FeatureVertexFinder::Find3dVtxFrom2dClusterVtxCand(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), corner::CornerFinderAlg::get_feature_points_fast(), GetCryostatVolumeName(), trkf::BezierTrack::GetProjectedPointUVWT(), trkf::BezierTrack::GetProjectedPointUVWX(), trkf::BezierTrack::GetTrackPitch(), lar_pandora::LArPandoraInput::GetTrueX0(), lar_pandora::PFParticleHitDumper::GetUVW(), cluster::fuzzyClusterAlg::InitFuzzy(), ems::EMShower3D::LinkCandidates(), trkf::SpacePointAlg::makeSpacePoints(), NOpDets(), OpDetFromCryo(), OpDetGeoFromOpDet(), OpDetGeoName(), PositionToCryostat(), cluster::ClusterCheater::produce(), vertex::HarrisVertexFinder::produce(), vertex::FeatureVertexFinder::produce(), larg4::LArG4::produce(), evd::SimulationDrawer::SimulationDrawer(), apa::APAGeometryAlg::ThreeChanPos(), cluster::HoughBaseAlg::Transform(), apa::DisambigAlg::TrivialDisambig(), DUNE::NeutrinoTrackingEff::truthLength(), trkf::SpacePointAlg::update(), apa::DisambigAlg::UseEndPts(), and trkf::VertexFitAlg::VertexFit().

                                                                           {
     CryostatGeo const* pCryo = CryostatPtr(cryoid);
     if(!pCryo) {
       throw cet::exception("GeometryCore") << "Cryostat #"
                                            << cryoid.Cryostat
                                            << " does not exist\n";
     }
     return *pCryo;
   } // GeometryCore::Cryostat(CryostatID)

CryostatGeo const& geo::GeometryCore::Cryostat ( unsigned int const cstat = 0 ) const

inline

Returns the specified cryostat.

Parameters

cstat	number of cryostat
cryoid	cryostat ID

Returns: a constant reference to the specified cryostat

Exceptions

cet::exception (GeometryCore category) if cryostat not present

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

Todo:: Make the cryostat number mandatory (as CryostatID)

Definition at line 2001 of file GeometryCore.h.

2002 { return Cryostat(geo::CryostatID(cstat)); }

geo::GeometryCore::Cryostat

CryostatGeo const & Cryostat(geo::CryostatID const &cryoid) const

Returns the specified cryostat.

Definition: GeometryCore.cxx:280

geo::GeometryCore::CryostatBoundaries

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

void geo::GeometryCore::CryostatBoundaries	(	double *	boundaries,
		geo::CryostatID const &	cid
	)		const

Returns the boundaries of the specified cryostat.

Parameters

boundaries	(output) pointer to an area of 6 doubles for boundaries
cid	cryostat ID

Exceptions

cet::exception ("GeometryCore" category) if cryostat not present

See also: CryostatGeo::Boundaries()

The boundaries array is filled with: [0] lower x coordinate [1] upper x coordinate [2] lower y coordinate [3] upper y coordinate [4] lower z coordinate [5] upper z coordinate

Deprecated:: Use CryostatGeo::Boundaries() (from Cryostat(cid)).

Todo:: What happen on invalid cryostat?

Definition at line 750 of file GeometryCore.cxx.

References geo::CryostatGeo::Boundaries(), and Cryostat().

Referenced by evgen::GaisserParam::beginJob(), ems::MergeEMShower3D::collectshowers(), CryostatLength(), evgen::CORSIKAGen::populateNShowers(), evgen::CosmicsGen::produce(), evgen::CORSIKAGen::produce(), and phot::PhotonVisibilityService::reconfigure().

   {
     geo::CryostatGeo const& cryo = Cryostat(cid);
     cryo.Boundaries(boundaries);
   } // GeometryCore::CryostatBoundaries()

void geo::GeometryCore::CryostatBoundaries	(	double *	boundaries,
		unsigned int	cstat = `0`
	)		const

inline

Returns the boundaries of the specified cryostat.

Parameters

boundaries	(output) pointer to an area of 6 doubles for boundaries
cstat	number of cryostat

Exceptions

cet::exception ("GeometryCore" category) if cryostat not present

See also: CryostatGeo::Boundaries()

The boundaries array is filled with: [0] lower x coordinate [1] upper x coordinate [2] lower y coordinate [3] upper y coordinate [4] lower z coordinate [5] upper z coordinate

Deprecated:: Use CryostatBoundaries(double*, geo::CryostatID const&) or (recommended) CryostatGeo::Boundaries() from Cryostat(cid) instead

Definition at line 2230 of file GeometryCore.h.

2231 { CryostatBoundaries(boundaries, geo::CryostatID(cstat)); }

void CryostatBoundaries(double *boundaries, geo::CryostatID const &cid) const

Returns the boundaries of the specified cryostat.

Definition: GeometryCore.cxx:750

geo::GeometryCore::CryostatHalfHeight

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

geo::Length_t geo::GeometryCore::CryostatHalfHeight ( geo::CryostatID const & cid ) const

Returns the height of the cryostat (y direction)

Definition at line 737 of file GeometryCore.cxx.

References Cryostat(), and geo::CryostatGeo::HalfHeight().

Referenced by CryostatHalfWidth(), simfilter::FilterNoDirtNeutrinos::filter(), and evgen::FileMuons::ReadEvents().

   {
     return Cryostat(cid).HalfHeight();
   }

geo::Length_t geo::GeometryCore::CryostatHalfHeight ( unsigned int cstat = 0 ) const

inline

Returns the height of the cryostat (y direction)

Definition at line 2184 of file GeometryCore.h.

2185 { return CryostatHalfHeight(geo::CryostatID(cstat)); }

geo::Length_t CryostatHalfHeight(geo::CryostatID const &cid) const

Returns the height of the cryostat (y direction)

Definition: GeometryCore.cxx:737

geo::GeometryCore::CryostatHalfWidth

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

geo::Length_t geo::GeometryCore::CryostatHalfWidth ( geo::CryostatID const & cid ) const

Returns the half width of the cryostat (x direction)

Definition at line 730 of file GeometryCore.cxx.

References Cryostat(), CryostatHalfHeight(), and geo::CryostatGeo::HalfWidth().

Referenced by DetLength(), simfilter::FilterNoDirtNeutrinos::filter(), and evgen::FileMuons::ReadEvents().

   {
     return Cryostat(cid).HalfWidth();
   }

geo::Length_t geo::GeometryCore::CryostatHalfWidth ( unsigned int cstat = 0 ) const

inline

Returns the half width of the cryostat (x direction)

Definition at line 2177 of file GeometryCore.h.

2178 { return CryostatHalfWidth(geo::CryostatID(cstat)); }

geo::Length_t CryostatHalfWidth(geo::CryostatID const &cid) const

Returns the half width of the cryostat (x direction)

Definition: GeometryCore.cxx:730

geo::GeometryCore::CryostatLength

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

geo::Length_t geo::GeometryCore::CryostatLength ( geo::CryostatID const & cid ) const

Returns the length of the cryostat (z direction)

Definition at line 743 of file GeometryCore.cxx.

References Cryostat(), CryostatBoundaries(), and geo::CryostatGeo::Length().

Referenced by evgen::FileMuons::ReadEvents().

   {
     return Cryostat(cid).Length();
   }

geo::Length_t geo::GeometryCore::CryostatLength ( unsigned int cstat = 0 ) const

inline

Returns the length of the cryostat (z direction)

Definition at line 2191 of file GeometryCore.h.

2192 { return CryostatLength(geo::CryostatID(cstat)); }

geo::Length_t CryostatLength(geo::CryostatID const &cid) const

Returns the length of the cryostat (z direction)

Definition: GeometryCore.cxx:743

geo::GeometryCore::Cryostats

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

CryostatGeo const* geo::GeometryCore::CryostatPtr ( geo::CryostatID const & cryoid ) const

inline

Returns the specified cryostat.

Parameters

cryoid cryostat ID

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

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

Definition at line 2016 of file GeometryCore.h.

References geo::CryostatID::Cryostat.

Referenced by Cryostat().

2017 { return HasCryostat(cryoid)? &(Cryostats()[cryoid.Cryostat]): nullptr; }

CryostatList_t & Cryostats()

Return the internal cryostat list.

Definition: GeometryCore.h:5340

geo::GeometryCore::HasCryostat

bool HasCryostat(geo::CryostatID const &cryoid) const

Returns whether we have the specified cryostat.

Definition: GeometryCore.h:1981

CryostatList_t& geo::GeometryCore::Cryostats ( )

inlineprotected

Return the internal cryostat list.

Definition at line 5340 of file GeometryCore.h.

Referenced by ClearGeometry(), MakeCryostat(), MaxPlanes(), MaxTPCs(), MaxWires(), SortGeometry(), TotalNTPC(), and UpdateAfterSorting().

5340 { return fGeoData.cryostats; }

geo::GeometryData_t::cryostats

GeometryData_t fGeoData

The detector description data.

Definition: GeometryCore.h:5391

CryostatList_t cryostats

The detector cryostats.

Definition: GeometryCore.h:1391

CryostatList_t const& geo::GeometryCore::Cryostats ( ) const

inlineprotected

Return the internal cryostat list.

Definition at line 5341 of file GeometryCore.h.

5341 { return fGeoData.cryostats; }

geo::GeometryData_t::cryostats

GeometryData_t fGeoData

The detector description data.

Definition: GeometryCore.h:5391

CryostatList_t cryostats

The detector cryostats.

Definition: GeometryCore.h:1391

double geo::GeometryCore::DefaultWiggle ( ) const

inline

Returns the tolerance used in looking for positions.

Returns: the tolerance value

This parameter is used as tolerance ("wiggle") for methods that require it (e.g. geo::CryostatGeo::FindTPCAtPosition()). Typically, it's a additional fraction of tolerance: 0 means no tolerance, 0.1 means 10% tolerance.

Todo:: Confirm the definition of wiggle: this one is taken from other doc

Definition at line 1705 of file GeometryCore.h.

1705 { return fPositionWiggle; }

geo::GeometryCore::fPositionWiggle

double fPositionWiggle

accounting for rounding errors when testing positions

Definition: GeometryCore.h:5399

geo::BoxBoundedGeo geo::GeometryCore::DetectorEnclosureBox ( std::string const & name = "volDetEnclosure" ) const

Returns the absolute coordinates of the detector enclosure volume [cm].

Parameters

name	name of the volume to be sought (default: `volDetEnclosure`)

Exceptions

cet::exception if the specified volume is not found

Definition at line 588 of file GeometryCore.cxx.

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

Referenced by GetWorldVolumeName().

   {
     auto const& path = FindDetectorEnclosure(name);
     if (path.empty()) {
       throw cet::exception("GeometryCore")
         << "DetectorEnclosureBox(): can't find enclosure volume '" << name << "'\n";
     }
     
     TGeoVolume const* pEncl = path.back()->GetVolume();
     auto const* pBox = dynamic_cast<TGeoBBox const*>(pEncl->GetShape());
     
     // check that this is indeed a box
     if (!pBox) {
       // at initialisation time we don't know yet our real ID
       throw cet::exception("GeometryCore") << "Detector enclosure '"
         << name << "' is not a box! (it is a " << pEncl->GetShape()->IsA()->GetName()
         << ")\n";
     }
 
     geo::LocalTransformation<TGeoHMatrix> trans(path, path.size() - 1);
     // get the half width, height, etc of the cryostat
     const double halfwidth  = pBox->GetDX();
     const double halfheight = pBox->GetDY();
     const double halflength = pBox->GetDZ();
     
     return {
       trans.LocalToWorld(geo::Point_t{ -halfwidth, -halfheight, -halflength }),
       trans.LocalToWorld(geo::Point_t{ +halfwidth, +halfheight, +halflength })
       };
   } // geo::GeometryCore::DetectorEnclosureBox()

std::string geo::GeometryCore::DetectorName ( ) const

inline

Returns a string with the name of the detector, as configured.

Definition at line 1735 of file GeometryCore.h.

1735 { return fDetectorName; }

geo::GeometryCore::fDetectorName

std::string fDetectorName

Name of the detector.

Definition: GeometryCore.h:5394

geo::Length_t geo::GeometryCore::DetHalfHeight ( geo::TPCID const & tpcid ) const

Returns the half height of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the half height of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveHalfHeight()

See geo::TPCGeo::ActiveHalfHeight() for more details.

Todo:

deprecate this function

rename the function

Definition at line 717 of file GeometryCore.cxx.

References geo::TPCGeo::ActiveHalfHeight(), and TPC().

Referenced by lar::example::AnalysisExample::analyze(), evgen::NUANCEGen::beginJob(), evgen::NDKGen::beginJob(), evgen::GENIEGen::beginJob(), trkf::SpacePointAna::bookHistograms(), LArStackingAction::ClassifyNewTrack(), evd_tool::MicroBooNEDrawer::DetOutline3D(), evd::Display3DPad::Draw(), lar_pandora::PFParticleCosmicAna::FillRecoTree(), tca::FillWireHitRange(), simfilter::FilterStoppingMuon::filter(), simfilter::FilterNoDirtNeutrinos::filter(), cmtool::CFAlgoVolumeOverlap::Float(), cmtool::CFAlgoWireOverlap::Float(), util::GeometryUtilities::Get2DangleFrom3D(), trk::TrackContainmentAlg::GetCosmicTagID(), ToyOneShowerGen::GetXYZPosition(), cosmic::BeamFlashTrackMatchTaggerAlg::InDetector(), tca::InsideTPC(), trk::TrackContainmentAlg::IsContained(), trkf::SeedAna::MCHists::MCHists(), trkf::TrackAna::MCHists::MCHists(), evd::Ortho3DPad::Ortho3DPad(), trk::TrackContainmentAlg::ProcessTracks(), trkf::SpacePts::produce(), trkf::Track3Dreco::produce(), evgen::NDKGen::produce(), vertex::VertexFinder2D::produce(), trkf::Track3DKalmanSPS::produce(), trkf::SeedAna::RecoHists::RecoHists(), trkf::TrackAna::RecoHists::RecoHists(), and evd::TWQProjectionView::SetSeeds().

   {
     return TPC(tpcid).ActiveHalfHeight();
   }

geo::Length_t geo::GeometryCore::DetHalfHeight	(	unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the half height of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the half height of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveHalfHeight()

See geo::TPCGeo::ActiveHalfHeight() for more details.

Todo:

deprecate this function

rename the function

Definition at line 2708 of file GeometryCore.h.

2709 { return DetHalfHeight(geo::TPCID(cstat, tpc)); }

geo::GeometryCore::DetHalfHeight

geo::Length_t DetHalfHeight(geo::TPCID const &tpcid) const

Returns the half height of the active volume of the specified TPC.

Definition: GeometryCore.cxx:717

Referenced by lar::example::AnalysisExample::analyze(), evgen::NUANCEGen::beginJob(), evgen::NDKGen::beginJob(), evgen::GENIEGen::beginJob(), trkf::SpacePointAna::bookHistograms(), LArStackingAction::ClassifyNewTrack(), evd_tool::MicroBooNEDrawer::DetOutline3D(), evd::Display3DPad::Draw(), lar_pandora::PFParticleCosmicAna::FillRecoTree(), tca::FillWireHitRange(), simfilter::FilterStoppingMuon::filter(), simfilter::FilterNoDirtNeutrinos::filter(), util::GeometryUtilities::Get2DangleFrom3D(), trk::TrackContainmentAlg::GetCosmicTagID(), ToyOneShowerGen::GetXYZPosition(), cosmic::BeamFlashTrackMatchTaggerAlg::InDetector(), cosmic::BeamFlashTrackMatchTaggerAlg::InDriftWindow(), tca::InsideTPC(), trk::TrackContainmentAlg::IsContained(), trkf::SeedAna::MCHists::MCHists(), trkf::TrackAna::MCHists::MCHists(), evd::RecoBaseDrawer::OpFlashOrtho(), evd::Ortho3DPad::Ortho3DPad(), trk::TrackContainmentAlg::ProcessTracks(), evgen::NUANCEGen::produce(), evgen::NDKGen::produce(), trkf::Track3DKalmanSPS::produce(), trkf::SeedAna::RecoHists::RecoHists(), trkf::TrackAna::RecoHists::RecoHists(), evd::TWQProjectionView::SetSeeds(), and evd::RecoBaseDrawer::Vertex2D().

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

geo::Length_t geo::GeometryCore::DetHalfWidth ( geo::TPCID const & tpcid ) const

Returns the half width of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the half width of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveHalfWidth()

Todo:

deprecate this function

rename the function

Definition at line 711 of file GeometryCore.cxx.

References geo::TPCGeo::ActiveHalfWidth(), and TPC().

   {
     return TPC(tpcid).ActiveHalfWidth();
   }

geo::Length_t geo::GeometryCore::DetHalfWidth	(	unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the half width of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the half width of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveHalfWidth()

Todo:

deprecate this function

rename the function

Definition at line 2686 of file GeometryCore.h.

2687 { return DetHalfWidth(geo::TPCID(cstat, tpc)); }

geo::GeometryCore::DetHalfWidth

geo::Length_t DetHalfWidth(geo::TPCID const &tpcid) const

Returns the half width of the active volume of the specified TPC.

Definition: GeometryCore.cxx:711

Referenced by lar::example::AnalysisExample::analyze(), evgen::NUANCEGen::beginJob(), evgen::NDKGen::beginJob(), evgen::GENIEGen::beginJob(), trkf::SpacePointAna::bookHistograms(), LArStackingAction::ClassifyNewTrack(), evd_tool::MicroBooNEDrawer::DetOutline3D(), evd::Display3DPad::Draw(), lar_pandora::PFParticleCosmicAna::FillRecoTree(), tca::FillWireHitRange(), simfilter::FilterNoDirtNeutrinos::filter(), simfilter::FilterStoppingMuon::filter(), cmtool::CFAlgoVolumeOverlap::Float(), cmtool::CFAlgoWireOverlap::Float(), util::GeometryUtilities::Get2DangleFrom3D(), trk::TrackContainmentAlg::GetCosmicTagID(), ToyOneShowerGen::GetXYZPosition(), cosmic::BeamFlashTrackMatchTaggerAlg::InDetector(), tca::InsideTPC(), trk::TrackContainmentAlg::IsContained(), trkf::SeedAna::MCHists::MCHists(), trkf::TrackAna::MCHists::MCHists(), lar_cluster3d::StandardHit3DBuilder::NearestWireID(), evd::Ortho3DPad::Ortho3DPad(), trk::TrackContainmentAlg::ProcessTracks(), evgen::CosmicsGen::produce(), evgen::NUANCEGen::produce(), evgen::NDKGen::produce(), trkf::Track3DKalmanSPS::produce(), trkf::SeedAna::RecoHists::RecoHists(), trkf::TrackAna::RecoHists::RecoHists(), evd::TWQProjectionView::SetSeeds(), and evd::RecoBaseDrawer::Vertex2D().

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

geo::Length_t geo::GeometryCore::DetLength ( geo::TPCID const & tpcid ) const

Returns the length of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the length of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveLength()

See geo::TPCGeo::ActiveLength() for more details.

Todo:

deprecate this function

rename the function

Definition at line 723 of file GeometryCore.cxx.

References geo::TPCGeo::ActiveLength(), CryostatHalfWidth(), and TPC().

   { 
     return TPC(tpcid).ActiveLength();
   }

geo::Length_t geo::GeometryCore::DetLength	(	unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the length of the active volume of the specified TPC.

Parameters

tpcid	ID of the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

Returns: the value of the length of the specified TPC

Exceptions

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

See also: geo::TPCGeo::ActiveLength()

See geo::TPCGeo::ActiveLength() for more details.

Todo:

deprecate this function

rename the function

Definition at line 2729 of file GeometryCore.h.

2730 { return DetLength(geo::TPCID(cstat, tpc)); }

geo::GeometryCore::DetLength

geo::Length_t DetLength(geo::TPCID const &tpcid) const

Returns the length of the active volume of the specified TPC.

Definition: GeometryCore.cxx:723

geo::GeometryCore::cryostat_iterator

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

cryostat_iterator geo::GeometryCore::end_cryostat ( ) const

inline

Returns an iterator pointing after the last cryostat.

Definition at line 2128 of file GeometryCore.h.

References begin_cryostat_id(), end_cryostat_id(), and geo::details::geometry_iterator_types::end_pos.

Referenced by evd_tool::ICARUSDrawer::DetOutline3D(), and IterateCryostatIDs().

2129 { return cryostat_iterator(this, cryostat_iterator::end_pos); }

geo::cryostat_iterator cryostat_iterator

Forward-iterator browsing all cryostats in the detector.

Definition: GeometryCore.h:1531

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

cryostat_id_iterator geo::GeometryCore::end_cryostat_id ( ) const

inline

Returns an iterator pointing after the last cryostat ID.

Definition at line 2120 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by end_cryostat().

2121 { return cryostat_id_iterator(this, cryostat_id_iterator::end_pos); }

geo::GeometryCore::cryostat_id_iterator

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::cryostat_id_iterator cryostat_id_iterator

Forward-iterator browsing all cryostat IDs in the detector.

Definition: GeometryCore.h:1510

plane_iterator geo::GeometryCore::end_plane ( ) const

inline

Returns an iterator pointing after the last plane in the detector.

Definition at line 3039 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by IteratePlaneIDs(), and IteratePlanes().

3040 { return plane_iterator(this, plane_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

plane_iterator geo::GeometryCore::end_plane ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing after the last plane in the specified cryostat.

Definition at line 3049 of file GeometryCore.h.

3050 { return plane_iterator(this, GetEndPlaneID(ID)); }

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

geo::PlaneID GetEndPlaneID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2994

plane_iterator geo::GeometryCore::end_plane ( geo::TPCID const & ID ) const

inline

Returns an iterator pointing after the last plane in the specified TPC.

Definition at line 3057 of file GeometryCore.h.

References begin_plane_id(), and end_plane_id().

3058 { return plane_iterator(this, GetEndPlaneID(ID)); }

geo::plane_iterator plane_iterator

Forward-iterator browsing all planes in the detector.

Definition: GeometryCore.h:1623

geo::PlaneID GetEndPlaneID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2994

plane_id_iterator geo::GeometryCore::end_plane_id ( ) const

inline

Returns an iterator pointing after the last plane ID in the detector.

Definition at line 3011 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by end_plane(), and IteratePlaneIDs().

3012 { return plane_id_iterator(this, plane_id_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

plane_id_iterator geo::GeometryCore::end_plane_id ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing after the last plane ID in the specified cryostat.

Definition at line 3021 of file GeometryCore.h.

3022 { return plane_id_iterator(this, GetEndPlaneID(ID)); }

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

geo::PlaneID GetEndPlaneID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2994

plane_id_iterator geo::GeometryCore::end_plane_id ( geo::TPCID const & ID ) const

inline

Returns an iterator pointing after the last plane ID in the specified TPC.

Definition at line 3031 of file GeometryCore.h.

3032 { return plane_id_iterator(this, GetEndPlaneID(ID)); }

geo::plane_id_iterator plane_id_iterator

Forward-iterator browsing all plane IDs in the detector.

Definition: GeometryCore.h:1601

geo::PlaneID GetEndPlaneID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2994

ROP_id_iterator geo::GeometryCore::end_ROP_id ( ) const

inline

Returns an iterator pointing after the last ROP ID in the detector.

Definition at line 5090 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by end_ROP_id(), and IterateROPIDs().

5091 { return ROP_id_iterator(this, ROP_id_iterator::end_pos); }

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

geo::GeometryCore::GetEndROPID

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

ROP_id_iterator geo::GeometryCore::end_ROP_id ( geo::CryostatID const & ID ) const

inline

Returns an iterator pointing after the last readout plane ID in the specified cryostat.

Definition at line 5100 of file GeometryCore.h.

5101 { return ROP_id_iterator(this, GetEndROPID(ID)); }

readout::ROPID GetEndROPID(geo::CryostatID const &id) const

Definition: GeometryCore.h:5073

geo::GeometryCore::GetEndROPID

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

ROP_id_iterator geo::GeometryCore::end_ROP_id ( readout::TPCsetID const & ID ) const

inline

Returns an iterator pointing after the last readout plane ID in the specified TPC set.

Definition at line 5110 of file GeometryCore.h.

References begin_ROP_id(), and end_ROP_id().

5111 { return ROP_id_iterator(this, GetEndROPID(ID)); }

readout::ROPID GetEndROPID(geo::CryostatID const &id) const

Definition: GeometryCore.h:5073

details::ROP_id_iterator_base< readout::ROPID > ROP_id_iterator

Forward iterator browsing all readout planes in the detector.

Definition: GeometryCore.h:1309

TPC_iterator geo::GeometryCore::end_TPC ( ) const

inline

Returns an iterator pointing after the last TPC in the detector.

Definition at line 2538 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by lar::example::SpacePointIsolationAlg::fillAlgConfigFromGeometry(), IterateTPCIDs(), IterateTPCs(), and sim::MCRecoPart::MCRecoPart().

2539 { return TPC_iterator(this, TPC_iterator::end_pos); }

geo::TPC_iterator TPC_iterator

Forward-iterator browsing all TPCs in the detector.

Definition: GeometryCore.h:1577

geo::GeometryCore::GetEndTPCID

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

TPC_iterator geo::GeometryCore::end_TPC ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing after the last TPC in the detector.

Definition at line 2546 of file GeometryCore.h.

References begin_TPC_id(), and end_TPC_id().

2547 { return TPC_iterator(this, GetEndTPCID(cid)); }

geo::TPCID GetEndTPCID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2511

geo::TPC_iterator TPC_iterator

Forward-iterator browsing all TPCs in the detector.

Definition: GeometryCore.h:1577

TPC_id_iterator geo::GeometryCore::end_TPC_id ( ) const

inline

Returns an iterator pointing after the last TPC ID in the detector.

Definition at line 2520 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by pma::PMAlgTracker::build(), end_TPC(), IterateTPCIDs(), and tss::TrackShowerHits::produce().

2521 { return TPC_id_iterator(this, TPC_id_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::GeometryCore::GetEndTPCID

geo::TPC_id_iterator TPC_id_iterator

Forward-iterator browsing all TPC IDs in the detector.

Definition: GeometryCore.h:1555

TPC_id_iterator geo::GeometryCore::end_TPC_id ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing after the last TPC ID in the specified cryostat.

Definition at line 2530 of file GeometryCore.h.

2531 { return TPC_id_iterator(this, GetEndTPCID(cid)); }

geo::TPCID GetEndTPCID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2511

geo::TPC_id_iterator TPC_id_iterator

Forward-iterator browsing all TPC IDs in the detector.

Definition: GeometryCore.h:1555

TPCset_id_iterator geo::GeometryCore::end_TPCset_id ( ) const

inline

Returns an iterator pointing after the last TPC set ID in the detector.

Definition at line 4853 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by end_TPCset_id(), and IterateTPCsetIDs().

4854 { return TPCset_id_iterator(this, TPCset_id_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::GeometryCore::GetEndTPCsetID

details::TPCset_id_iterator_base< readout::TPCsetID > TPCset_id_iterator

Forward iterator browsing all TPC sets in the detector.

Definition: GeometryCore.h:1288

TPCset_id_iterator geo::GeometryCore::end_TPCset_id ( geo::CryostatID const & cid ) const

inline

Returns an iterator pointing after the last TPC set ID in the specified cryostat.

Definition at line 4863 of file GeometryCore.h.

References begin_TPCset_id(), and end_TPCset_id().

4864 { return TPCset_id_iterator(this, GetEndTPCsetID(cid)); }

readout::TPCsetID GetEndTPCsetID(geo::CryostatID const &id) const

Definition: GeometryCore.h:4844

details::TPCset_id_iterator_base< readout::TPCsetID > TPCset_id_iterator

Forward iterator browsing all TPC sets in the detector.

Definition: GeometryCore.h:1288

wire_iterator geo::GeometryCore::end_wire ( ) const

inline

Returns an iterator pointing after the last wire in the detector.

Definition at line 3459 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by IterateWireIDs(), and IterateWires().

3460 { return wire_iterator(this, wire_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::GeometryCore::end_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::end_wire ( geo::CryostatID const & id ) const

inline

Returns an iterator pointing after the last wire in specified cryostat.

Definition at line 3467 of file GeometryCore.h.

3468 { return wire_iterator(end_wire_id(id)); }

wire_id_iterator end_wire_id() const

Returns an iterator pointing after the last wire ID in the detector.

Definition: GeometryCore.h:3426

geo::GeometryCore::end_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::end_wire ( geo::TPCID const & id ) const

inline

Returns an iterator pointing after the last wire in specified TPC.

Definition at line 3475 of file GeometryCore.h.

3476 { return wire_iterator(end_wire_id(id)); }

wire_id_iterator end_wire_id() const

Returns an iterator pointing after the last wire ID in the detector.

Definition: GeometryCore.h:3426

geo::GeometryCore::end_wire_id

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_iterator geo::GeometryCore::end_wire ( geo::PlaneID const & id ) const

inline

Returns an iterator pointing after the last wire in specified plane.

Definition at line 3483 of file GeometryCore.h.

References begin_wire_id(), and end_wire_id().

3484 { return wire_iterator(end_wire_id(id)); }

wire_id_iterator end_wire_id() const

Returns an iterator pointing after the last wire ID in the detector.

Definition: GeometryCore.h:3426

geo::wire_iterator wire_iterator

Forward-iterator browsing all wires in the detector.

Definition: GeometryCore.h:1669

wire_id_iterator geo::GeometryCore::end_wire_id ( ) const

inline

Returns an iterator pointing after the last wire ID in the detector.

Definition at line 3426 of file GeometryCore.h.

References geo::details::geometry_iterator_types::end_pos.

Referenced by end_wire(), IterateWireIDs(), lariov::SIOVElectronicsCalibProvider::Reconfigure(), and lariov::DetPedestalRetrievalAlg::Reconfigure().

3427 { return wire_id_iterator(this, wire_id_iterator::end_pos); }

static constexpr EndPos_t end_pos

Definition: GeometryCore.h:120

geo::GeometryCore::GetEndWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::end_wire_id ( geo::CryostatID const & id ) const

inline

Returns an iterator pointing after the last wire ID in specified cryostat.

Definition at line 3435 of file GeometryCore.h.

3436 { return wire_id_iterator(this, GetEndWireID(id)); }

geo::WireID GetEndWireID(geo::CryostatID const &id) const

Definition: GeometryCore.h:3400

geo::GeometryCore::GetEndWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::end_wire_id ( geo::TPCID const & id ) const

inline

Returns an iterator pointing after the last wire ID in specified TPC.

Definition at line 3443 of file GeometryCore.h.

3444 { return wire_id_iterator(this, GetEndWireID(id)); }

geo::WireID GetEndWireID(geo::CryostatID const &id) const

Definition: GeometryCore.h:3400

geo::GeometryCore::GetEndWireID

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

wire_id_iterator geo::GeometryCore::end_wire_id ( geo::PlaneID const & id ) const

inline

Returns an iterator pointing after the last wire ID in specified plane.

Definition at line 3451 of file GeometryCore.h.

3452 { return wire_id_iterator(this, GetEndWireID(id)); }

geo::WireID GetEndWireID(geo::CryostatID const &id) const

Definition: GeometryCore.h:3400

geo::GeometryCore::FindAuxDetAtPosition

geo::wire_id_iterator wire_id_iterator

Forward-iterator browsing all wire IDs in the detector.

Definition: GeometryCore.h:1647

std::vector< std::vector< TGeoNode const * > > geo::GeometryCore::FindAllVolumePaths ( std::set< std::string > const & vol_names ) const

Returns paths of all nodes with volumes with the specified names.

Parameters

vol_names list of names of volumes

Returns: list paths of the found nodes

All the nodes in the geometry are checked, and the path of all the ones that contain a volume with a name among the ones specified in vol_names is saved in the collection and returned. A node path is a ordered list of all nodes leading to the final one, starting from thetop level (root) down. The node at the back() of the path is the one with name in vol_names. No empty paths are returned.

Definition at line 682 of file GeometryCore.cxx.

References geo::CollectPathsByName::paths, and ROOTGeoManager().

Referenced by FindAllVolumes().

   {
     CollectPathsByName path_collector(vol_names);
     
     ROOTGeoNodeForwardIterator iNode(ROOTGeoManager()->GetTopNode());
     
     while (*iNode) {
       path_collector(iNode);
       ++iNode;
     } // while
     
     return path_collector.paths;
   } // GeometryCore::FindAllVolumePaths()

std::vector< TGeoNode const * > geo::GeometryCore::FindAllVolumes ( std::set< std::string > const & vol_names ) const

Returns all the nodes with volumes with any of the specified names.

Parameters

vol_names list of names of volumes

Returns: list of nodes found

All the nodes in the geometry are checked, and all the ones that contain a volume with a name among the ones specified in vol_names are saved in the collection and returned.

Definition at line 665 of file GeometryCore.cxx.

References FindAllVolumePaths(), geo::CollectNodesByName::nodes, and ROOTGeoManager().

   {
     CollectNodesByName node_collector(vol_names);
     
     ROOTGeoNodeForwardIterator iNode(ROOTGeoManager()->GetTopNode());
     TGeoNode const* pCurrentNode;
     
     while ((pCurrentNode = *iNode)) {
       node_collector(*pCurrentNode);
       ++iNode;
     } // while
     
     return node_collector.nodes;
   } // GeometryCore::FindAllVolumes()

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

private

Definition at line 1077 of file GeometryCore.cxx.

References MakeAuxDet().

Referenced by LoadGeometryFile().

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

unsigned int geo::GeometryCore::FindAuxDetAtPosition ( double const worldLoc[3] ) const

Returns the index of the auxiliary detector at specified location.

Parameters

worldLoc 3D coordinates of the point (world reference frame)

Returns: the index of the detector, or std::numeric_limits<unsigned int>max() if no detector is there

Bug:: Actually, an exception is thrown.

Deprecated:: Use the version with geo::Point_t.

Definition at line 464 of file GeometryCore.cxx.

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

Referenced by FindAuxDetSensitiveAtPosition(), and PositionToAuxDet().

465 { return FindAuxDetAtPosition(geo::vect::makePointFromCoords(worldPos)); }

unsigned int FindAuxDetAtPosition(double const worldLoc[3]) const

Returns the index of the auxiliary detector at specified location.

Definition: GeometryCore.cxx:464

geo::GeometryCore::FindAuxDetSensitiveAtPosition

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

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

Definition: geo_vectors_utils.h:1307

unsigned int geo::GeometryCore::FindAuxDetAtPosition ( geo::Point_t const & point ) const

Returns the index of the auxiliary detector at specified location.

Parameters

point location to be tested

Returns: the index of the detector, or std::numeric_limits<unsigned int>max() if no detector is there

Bug:: Actually, an exception is thrown.

Definition at line 455 of file GeometryCore.cxx.

References AuxDets(), and fChannelMapAlg.

   {
     // BUG the double brace syntax is required to work around clang bug 21629
     // (https://bugs.llvm.org/show_bug.cgi?id=21629)
     std::array<double, 3U> worldPos = {{ point.X(), point.Y(), point.Z() }};
     return fChannelMapAlg->NearestAuxDet(worldPos.data(), AuxDets());
   } // GeometryCore::FindAuxDetAtPosition()

void geo::GeometryCore::FindAuxDetSensitiveAtPosition	(	geo::Point_t const &	point,
		std::size_t &	adg,
		std::size_t &	sv
	)		const

Fills the indices of the sensitive auxiliary detector at location.

Parameters

point	location to be tested
adg	_(output)_ auxiliary detector index
sv	_(output)_ sensitive volume index

Definition at line 484 of file GeometryCore.cxx.

References AuxDets(), fChannelMapAlg, and FindAuxDetAtPosition().

Referenced by larg4::AuxDetReadoutGeometry::FindAndMakeAuxDetSensitive(), FindAuxDetSensitiveAtPosition(), PositionToAuxDet(), and PositionToAuxDetSensitive().

   {
     adg = FindAuxDetAtPosition(point);
     // BUG the double brace syntax is required to work around clang bug 21629
     // (https://bugs.llvm.org/show_bug.cgi?id=21629)
     std::array<double, 3U> const worldPos = {{ point.X(), point.Y(), point.Z() }};
     sv  = fChannelMapAlg->NearestSensitiveAuxDet(worldPos.data(), AuxDets());
   } // GeometryCore::FindAuxDetAtPosition()

void geo::GeometryCore::FindAuxDetSensitiveAtPosition	(	double const	worldLoc[3],
		size_t &	adg,
		size_t &	sv
	)		const

Fills the indices of the sensitive auxiliary detector at location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
adg	(output) auxiliary detector index
sv	(output) sensitive volume index

Deprecated:: Use the version with geo::Point_t.

Definition at line 495 of file GeometryCore.cxx.

References FindAuxDetSensitiveAtPosition(), geo::vect::makePointFromCoords(), and PositionToAuxDetSensitive().

496 { return FindAuxDetSensitiveAtPosition(geo::vect::makePointFromCoords(worldPos), adg, sv); }

void FindAuxDetSensitiveAtPosition(geo::Point_t const &point, std::size_t &adg, std::size_t &sv) const

Fills the indices of the sensitive auxiliary detector at location.

Definition: GeometryCore.cxx:484

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), ems::MultiEMShowers::convCluster(), hit::HitCheater::FindHitsOnChannel(), shower::EMShowerAlg::FindInitialTrackHits(), FindTPCsetAtPosition(), opdet::GetHitGeometryInfo(), nnet::TrainingDataAlg::getProjection(), lar_pandora::LArPandoraInput::GetTrueStartAndEndPoints(), hit::DisambigCheater::InitHitToWids(), evd::SimulationDrawer::MCTruthOrtho(), shower::EMShowerAlg::OrderShowerHits(), calo::Calorimetry::produce(), shower::EMShowerAlg::Project3DPointOntoPlane(), filt::LArG4ParticleFilter::StartInTPCCheck(), filt::LArG4ParticleFilter::StopInTPCCheck(), filt::LArG4ParticleFilter::TPCTrajLengthCheck(), DUNE::NeutrinoTrackingEff::truthLength(), and DUNE::MuonTrackingEff::truthLength().

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

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

Definition: geo_vectors_utils.h:1307

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

private

Definition at line 1046 of file GeometryCore.cxx.

References MakeCryostat().

Referenced by LoadGeometryFile().

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

geo::CryostatID::CryostatID_t geo::GeometryCore::FindCryostatAtPosition ( geo::Point_t const & worldLoc ) const

Returns the index of the cryostat at specified location.

Parameters

worldLoc 3D coordinates of the point (world reference frame)

Returns: the index of the cryostat, or UINT_MAX if no cryostat is there

Deprecated:: Use PositionToCryostatID() instead

Definition at line 352 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, geo::CryostatGeo::ID(), geo::CryostatID::InvalidID, and PositionToCryostatPtr().

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), ems::MultiEMShowers::convCluster(), FindCryostatAtPosition(), evd::SimulationDrawer::MCTruthOrtho(), PositionToCryostat(), and PositionToCryostatID().

   {
     geo::CryostatGeo const* cryo = PositionToCryostatPtr(worldLoc);
     return cryo? cryo->ID().Cryostat: geo::CryostatID::InvalidID;
   } // GeometryCore::FindCryostatAtPosition(Point)

geo::CryostatID::CryostatID_t geo::GeometryCore::FindCryostatAtPosition ( double const worldLoc[3] ) const

Returns the index of the cryostat at specified location.

Parameters

worldLoc 3D coordinates of the point (world reference frame)

Returns: the index of the cryostat, or UINT_MAX if no cryostat is there

Deprecated:: Use PositionToCryostatID() instead

Definition at line 361 of file GeometryCore.cxx.

References FindCryostatAtPosition(), geo::vect::makePointFromCoords(), and PositionToTPCptr().

   {
     return FindCryostatAtPosition(geo::vect::makePointFromCoords(worldLoc));
   } // GeometryCore::FindCryostatAtPosition(double[])

std::vector< TGeoNode const * > geo::GeometryCore::FindDetectorEnclosure ( std::string const & name = "volDetEnclosure" ) const

private

Definition at line 957 of file GeometryCore.cxx.

References FindFirstVolume(), and ROOTGeoManager().

Referenced by DetectorEnclosureBox(), and MaterialName().

   {
      std::vector<TGeoNode const*> path { ROOTGeoManager()->GetTopNode() };
      if (!FindFirstVolume(name, path)) path.clear();
      return path;
   } // FindDetectorEnclosure()

bool geo::GeometryCore::FindFirstVolume	(	std::string const &	name,
		std::vector< const TGeoNode * > &	path
	)		const

private

Definition at line 966 of file GeometryCore.cxx.

Referenced by FindDetectorEnclosure().

   {
     assert(!path.empty());
 
     auto const* pCurrent = path.back();
 
     // first check the current layer
     if (strncmp(name.c_str(), pCurrent->GetName(), name.length()) == 0)
       return true;
     
     //explore the next layer down
     auto const* pCurrentVolume = pCurrent->GetVolume();
     unsigned int nd = pCurrentVolume->GetNdaughters();
     for(unsigned int i = 0; i < nd; ++i) {
       path.push_back(pCurrentVolume->GetNode(i));
       if (FindFirstVolume(name, path)) return true;
       path.pop_back();
     } // for
     return false;
   } // FindFirstVolume()

geo::TPCID geo::GeometryCore::FindTPCAtPosition ( double const worldLoc[3] ) const

inline

Returns the ID of the TPC at specified location.

Parameters

worldLoc 3D coordinates of the point (world reference frame) [cm]

Returns: the TPC ID, or an invalid one if no TPC is there

Definition at line 2424 of file GeometryCore.h.

References geo::vect::makePointFromCoords().

2425 { return FindTPCAtPosition(geo::vect::makePointFromCoords(worldLoc)); }

geo::GeometryCore::FindTPCAtPosition

geo::TPCID FindTPCAtPosition(double const worldLoc[3]) const

Returns the ID of the TPC at specified location.

Definition: GeometryCore.h:2424

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

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

Definition: geo_vectors_utils.h:1307

geo::TPCID geo::GeometryCore::FindTPCAtPosition ( geo::Point_t const & point ) const

Returns the ID of the TPC at specified location.

Parameters

worldLoc 3D point (world reference frame, centimeters)

Returns: the TPC ID, or an invalid one if no TPC is there

Definition at line 311 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, fPositionWiggle, geo::CryostatGeo::ID(), geo::CryostatID::markInvalid(), PositionToCryostatPtr(), and geo::CryostatGeo::PositionToTPCID().

                                                                         {
     
     // first find the cryostat
     geo::CryostatGeo const* cryo = PositionToCryostatPtr(point);
     if (!cryo) return {};
     
     // then ask it about the TPC
     geo::TPCID tpcid = cryo->PositionToTPCID(point, 1. + fPositionWiggle);
     if (tpcid) return tpcid;
     
     // return an invalid TPC ID with cryostat information set:
     tpcid.Cryostat = cryo->ID().Cryostat;
     tpcid.markInvalid();
     return tpcid;
     
   } // GeometryCore::FindTPCAtPosition()

geo::TPCID geo::GeometryCore::FindTPCAtPosition ( TVector3 const & point ) const

inline

Returns the ID of the TPC at specified location.

Parameters

worldLoc 3D point (world reference frame, centimeters)

Returns: the TPC ID, or an invalid one if no TPC is there

Definition at line 2434 of file GeometryCore.h.

References geo::vect::toPoint().

2435 { return FindTPCAtPosition(geo::vect::toPoint(point)); }

geo::GeometryCore::FindTPCAtPosition

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

geo::TPCID FindTPCAtPosition(double const worldLoc[3]) const

Returns the ID of the TPC at specified location.

Definition: GeometryCore.h:2424

readout::TPCsetID geo::GeometryCore::FindTPCsetAtPosition ( double const worldLoc[3] ) const

Returns the ID of the TPC set at specified location.

Parameters

worldLoc 3D coordinates of the point (world reference frame)

Returns: the TPC set ID, or an invalid one if no TPC set is there

Definition at line 1783 of file GeometryCore.cxx.

References FindTPCAtPosition(), and TPCtoTPCset().

Referenced by HasTPCset().

   {
     return TPCtoTPCset(FindTPCAtPosition(worldLoc));
   } // GeometryCore::FindTPCsetAtPosition()

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

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.

Definition at line 1850 of file GeometryCore.cxx.

References fChannelMapAlg.

Referenced by ROPtoTPCs().

   {
     return fChannelMapAlg->FirstChannelInROP(ropid);
   } // GeometryCore::FirstChannelInROP()

std::string geo::GeometryCore::GDMLFile ( ) const

inline

Returns the full directory path to the GDML file source.

Returns: the full directory path to the GDML file source

This is the full path of the source of the detector geometry handed to the detector simulation (GEANT).

Definition at line 1723 of file GeometryCore.h.

Referenced by larg4::LArG4::beginJob(), and evd::Landed::Landed().

1723 { return fGDMLfile; }

geo::GeometryCore::fGDMLfile

std::string fGDMLfile

path to geometry file used for Geant4 simulation

Definition: GeometryCore.h:5395

template<typename GeoID >

GeoID geo::GeometryCore::GetBeginID ( ) const

inline

Returns the ID of the first element of the detector.

Template Parameters

GeoID type of the ID to be returned

Returns: ID of the first subelement in the detector

Definition at line 1903 of file GeometryCore.h.

References GetBeginID().

Referenced by GetBeginID().

1903 { GeoID id; GetBeginID(id); return id; }

geo::GeometryCore::HasCryostat

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

template<typename GeoID , typename ContextID >

GeoID geo::GeometryCore::GetBeginID ( ContextID const & id ) const

Returns the ID of the first subelement of the specified element.

Template Parameters

GeoID	type of the ID to be returned
ContextID	type of the ID of the containing element

Parameters

id	ID of the containing element

Returns: ID of the first subelement in the specified geometry element

void geo::GeometryCore::GetBeginID ( geo::CryostatID & id ) const

inline

Initializes the specified ID with the ID of the first cryostat.

Definition at line 2104 of file GeometryCore.h.

2105 { id = geo::CryostatID(0, HasCryostat(geo::CryostatID(0))); }

readout::CryostatID

geo::CryostatID CryostatID

Definition: readout_types.h:30

bool HasCryostat(geo::CryostatID const &cryoid) const

Returns whether we have the specified cryostat.

Definition: GeometryCore.h:1981

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

void geo::GeometryCore::GetBeginID ( geo::TPCID & id ) const

inline

Initializes the specified ID with the ID of the first TPC.

iterators

Definition at line 2494 of file GeometryCore.h.

2495 { GetBeginID(id.asCryostatID()); id.TPC = 0; }

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

void geo::GeometryCore::GetBeginID ( geo::PlaneID & id ) const

inline

Initializes the specified ID with the ID of the first plane.

Definition at line 2977 of file GeometryCore.h.

2978 { GetBeginID(id.asTPCID()); id.Plane = 0; }

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

void geo::GeometryCore::GetBeginID ( geo::WireID & id ) const

inline

Initializes the specified ID with the ID of the first wire.

Definition at line 3383 of file GeometryCore.h.

3384 { GetBeginID(id.asPlaneID()); id.Wire = 0; }

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

void geo::GeometryCore::GetBeginID ( readout::TPCsetID & id ) const

inline

Initializes the specified ID with the ID of the first TPC set.

iterators

Definition at line 4827 of file GeometryCore.h.

4828 { GetBeginID(id.asCryostatID()); id.TPCset = 0; }

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

void geo::GeometryCore::GetBeginID ( readout::ROPID & id ) const

inline

Initializes the specified ID with the ID of the first readout plane.

iterators

Definition at line 5056 of file GeometryCore.h.

5057 { GetBeginID(id.asTPCsetID()); id.ROP = 0; }

geo::GeometryCore::GetBeginTPCID

GeoID GetBeginID() const

Returns the ID of the first element of the detector.

Definition: GeometryCore.h:1903

geo::PlaneID geo::GeometryCore::GetBeginPlaneID ( geo::CryostatID const & id ) const

inline

Returns the ID of the first plane of the specified cryostat.

Definition at line 2989 of file GeometryCore.h.

2990 { return { GetBeginTPCID(id), 0 }; }

geo::TPCID GetBeginTPCID(geo::CryostatID const &id) const

Returns the ID of the first TPC in the specified cryostat.

Definition: GeometryCore.h:2506

geo::PlaneID geo::GeometryCore::GetBeginPlaneID ( geo::TPCID const & id ) const

inline

Returns the ID of the first plane of the specified TPC.

Definition at line 2998 of file GeometryCore.h.

2999 { return { id, 0 }; }

readout::ROPID geo::GeometryCore::GetBeginROPID ( geo::CryostatID const & id ) const

inline

Returns the ID of the first readout plane of the specified cryostat.

Definition at line 5068 of file GeometryCore.h.

5069 { return { GetBeginTPCsetID(id), 0 }; }

geo::GeometryCore::GetBeginTPCsetID

readout::TPCsetID GetBeginTPCsetID(geo::CryostatID const &id) const

Returns the ID of the first TPC set in the specified cryostat.

Definition: GeometryCore.h:4839

readout::ROPID geo::GeometryCore::GetBeginROPID ( readout::TPCsetID const & id ) const

inline

Returns the ID of the first readout plane of the specified TPC set.

Definition at line 5077 of file GeometryCore.h.

5078 { return { id, 0 }; }

geo::TPCID geo::GeometryCore::GetBeginTPCID ( geo::CryostatID const & id ) const

inline

Returns the ID of the first TPC in the specified cryostat.

Definition at line 2506 of file GeometryCore.h.

2507 { return { id, 0 }; }

readout::TPCsetID geo::GeometryCore::GetBeginTPCsetID ( geo::CryostatID const & id ) const

inline

Returns the ID of the first TPC set in the specified cryostat.

Definition at line 4839 of file GeometryCore.h.

4840 { return { id, 0 }; }

geo::WireID geo::GeometryCore::GetBeginWireID ( geo::CryostatID const & id ) const

inline

Returns the ID of the first wire in the specified cryostat.

Definition at line 3395 of file GeometryCore.h.

3396 { return { GetBeginPlaneID(id), 0 }; }

geo::PlaneID GetBeginPlaneID(geo::CryostatID const &id) const

Returns the ID of the first plane of the specified cryostat.

Definition: GeometryCore.h:2989

geo::WireID geo::GeometryCore::GetBeginWireID ( geo::TPCID const & id ) const

inline

Returns the ID of the first wire of the specified TPC.

Definition at line 3404 of file GeometryCore.h.

3405 { return { geo::PlaneID(id, 0), 0 }; }

geo::GeometryCore::GetClosestOpDet

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID geo::GeometryCore::GetBeginWireID ( geo::PlaneID const & id ) const

inline

Returns the ID of the first wire of the specified wire plane.

Definition at line 3413 of file GeometryCore.h.

3414 { return { id, 0 }; }

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

Find the nearest OpChannel to some point.

Parameters

xyz	point to be queried, in world coordinates

Returns: the nearest OpChannel to the point, or std::numeric_limits<unsigned int>max() if invalid point

Deprecated:: This method does not tell in which cryostat the detector is; use geo::CryostatGeo::GetClosestOpDet() instead (find the cryostat with PositionToCryostatPtr()).

Definition at line 1959 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, geo::CryostatGeo::GetClosestOpDet(), geo::CryostatGeo::ID(), max, OpDetFromCryo(), and PositionToCryostatPtr().

Referenced by GetClosestOpDet().

   {
     geo::CryostatGeo const* cryo = PositionToCryostatPtr(point);
     if (!cryo) return std::numeric_limits<unsigned int>::max();
     int o = cryo->GetClosestOpDet(point);
     return OpDetFromCryo(o, cryo->ID().Cryostat);
   }

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

Find the nearest OpChannel to some point.

Parameters

xyz	point to be queried, in world coordinates

Returns: the nearest OpChannel to the point, or std::numeric_limits<unsigned int>max() if invalid point

Deprecated:: This method does not tell in which cryostat the detector is; use geo::CryostatGeo::GetClosestOpDet() instead (find the cryostat with PositionToCryostatPtr()).

Definition at line 1970 of file GeometryCore.cxx.

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

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

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

Find the nearest OpChannel to some point.

Definition: GeometryCore.cxx:1959

geo::GeometryCore::GetCryostatVolumeName

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

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

Definition: geo_vectors_utils.h:1307

std::string geo::GeometryCore::GetCryostatVolumeName ( geo::CryostatID const & cid ) const

Return the name of LAr TPC volume.

Parameters

cstat index of the cryostat

Returns: the name of the specified TPC

This information is used in the event display.

Todo:

Use a cryostat ID instead

What if it does not exist?

Definition at line 705 of file GeometryCore.cxx.

References Cryostat().

   {
     return std::string(Cryostat(cid).Volume()->GetName());
   }

std::string geo::GeometryCore::GetCryostatVolumeName ( unsigned int const cstat = 0 ) const

inline

Return the name of LAr TPC volume.

Parameters

cstat index of the cryostat

Returns: the name of the specified TPC

This information is used in the event display.

Todo:

Use a cryostat ID instead

What if it does not exist?

Definition at line 2250 of file GeometryCore.h.

2251 { return GetCryostatVolumeName(geo::CryostatID(cstat)); }

std::string GetCryostatVolumeName(geo::CryostatID const &cid) const

Return the name of LAr TPC volume.

Definition: GeometryCore.cxx:705

geo::GeometryCore::Cryostat

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

CryostatGeo const& geo::GeometryCore::GetElement ( geo::CryostatID const & cryoid ) const

inline

Returns the specified cryostat.

Parameters

cstat	number of cryostat
cryoid	cryostat ID

Returns: a constant reference to the specified cryostat

Exceptions

cet::exception (GeometryCore category) if cryostat not present

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

Todo:: Make the cryostat number mandatory (as CryostatID)

Definition at line 2003 of file GeometryCore.h.

Referenced by evd::SimulationDrawer::MCTruthOrtho().

2004 { return Cryostat(cryoid); }

CryostatGeo const & Cryostat(geo::CryostatID const &cryoid) const

Returns the specified cryostat.

Definition: GeometryCore.cxx:280

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

inline

Returns the specified TPC.

Parameters

tpcid	ID of the tpc
tpc	tpc number within the cryostat
cstat	number of cryostat

Returns: a constant reference to the specified TPC

Exceptions

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

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

Todo:: remove the version with integers

Definition at line 2397 of file GeometryCore.h.

2398 { return TPC(tpcid); }

geo::GeometryCore::CryostatPtr

TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified TPC.

Definition: GeometryCore.h:2393

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

inline

Returns the specified wire.

Parameters

planeid	ID of the plane
p	plane number within the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

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.

Todo:: remove the version with integers

Definition at line 2950 of file GeometryCore.h.

2951 { return Plane(planeid); }

geo::GeometryCore::Plane

PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified wire.

Definition: GeometryCore.h:2945

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

inline

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 3374 of file GeometryCore.h.

3375 { return Wire(wireid); }

geo::GeometryCore::Wire

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

Returns the specified wire.

Definition: GeometryCore.h:3370

CryostatGeo const* geo::GeometryCore::GetElementPtr ( geo::CryostatID const & cryoid ) const

inline

Returns the specified cryostat.

Parameters

cryoid cryostat ID

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

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

Definition at line 2018 of file GeometryCore.h.

2019 { return CryostatPtr(cryoid); }

CryostatGeo const * CryostatPtr(geo::CryostatID const &cryoid) const

Returns the specified cryostat.

Definition: GeometryCore.h:2016

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

inline

Returns the specified TPC.

Parameters

tpcid TPC ID

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

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

Definition at line 2415 of file GeometryCore.h.

2416 { return TPCPtr(tpcid); }

geo::GeometryCore::TPCPtr

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

Returns the specified TPC.

Definition: GeometryCore.h:2410

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

inline

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 2968 of file GeometryCore.h.

2969 { return PlanePtr(planeid); }

geo::GeometryCore::PlanePtr

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

Returns the specified plane.

Definition: GeometryCore.h:2963

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

inline

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 3356 of file GeometryCore.h.

3357 { return WirePtr(wireid); }

geo::GeometryCore::WirePtr

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

Returns the specified wire.

Definition: GeometryCore.h:3351

template<typename GeoID >

GeoID geo::GeometryCore::GetEndID ( ) const

inline

Returns the (possibly invalid) ID after the last subelement of the detector.

Template Parameters

GeoID type of the ID to be returned

Returns: ID after the last subelement in the specified geometry element

Definition at line 1922 of file GeometryCore.h.

References GetEndID().

Referenced by GetEndID().

1922 { GeoID id; GetEndID(id); return id; }

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

template<typename GeoID , typename ContextID >

GeoID geo::GeometryCore::GetEndID ( ContextID const & id ) const

Returns the (possibly invalid) ID after the last subelement of the specified element.

Template Parameters

GeoID	type of the ID to be returned
ContextID	type of the ID of the containing element

Parameters

id	ID of the containing element

Returns: ID (possibly invalid) after the last subelement in the specified geometry element

void geo::GeometryCore::GetEndID ( geo::CryostatID & id ) const

inline

Initializes the specified ID with the invalid ID after the last cryostat.

Definition at line 2108 of file GeometryCore.h.

2109 { id = geo::CryostatID(Ncryostats(), false); }

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

readout::CryostatID

geo::CryostatID CryostatID

Definition: readout_types.h:30

void geo::GeometryCore::GetEndID ( geo::TPCID & id ) const

inline

Initializes the specified ID with the invalid ID after the last TPC.

Definition at line 2498 of file GeometryCore.h.

2499 { GetEndID(id.asCryostatID()); id.TPC = 0; }

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

void geo::GeometryCore::GetEndID ( geo::PlaneID & id ) const

inline

Initializes the specified ID with the invalid ID after the last plane.

Definition at line 2981 of file GeometryCore.h.

2982 { GetEndID(id.asTPCID()); id.Plane = 0; }

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

void geo::GeometryCore::GetEndID ( geo::WireID & id ) const

inline

Initializes the specified ID with the invalid ID after the last wire.

Definition at line 3387 of file GeometryCore.h.

3388 { GetEndID(id.asPlaneID()); id.Wire = 0; }

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

void geo::GeometryCore::GetEndID ( readout::TPCsetID & id ) const

inline

Initializes the specified ID with the invalid ID after the last TPC set.

Definition at line 4831 of file GeometryCore.h.

4832 { GetEndID(id.asCryostatID()); id.TPCset = 0; }

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

void geo::GeometryCore::GetEndID ( readout::ROPID & id ) const

inline

Initializes the specified ID with the invalid ID after the last ROP.

Definition at line 5060 of file GeometryCore.h.

5061 { GetEndID(id.asTPCsetID()); id.ROP = 0; }

geo::GeometryCore::GetEndTPCID

GeoID GetEndID() const

Returns the (possibly invalid) ID after the last subelement of the detector.

Definition: GeometryCore.h:1922

geo::PlaneID geo::GeometryCore::GetEndPlaneID ( geo::CryostatID const & id ) const

inline

Returns the (possibly invalid) ID after the last plane of the specified cryostat.

Definition at line 2994 of file GeometryCore.h.

2995 { return { GetEndTPCID(id), 0 }; }

geo::TPCID GetEndTPCID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2511

geo::PlaneID geo::GeometryCore::GetEndPlaneID ( geo::TPCID const & id ) const

inline

Returns the (possibly invalid) ID after the last plane of the specified TPC.

Definition at line 3003 of file GeometryCore.h.

3004 { return { GetNextID(id), 0 }; }

geo::GeometryCore::GetEndTPCsetID

GeoID GetNextID(GeoID const &id) const

Returns the ID next to the specified one.

Definition: GeometryCore.h:1912

readout::ROPID geo::GeometryCore::GetEndROPID ( geo::CryostatID const & id ) const

inline

Returns the (possibly invalid) ID after the last readout plane of the specified cryostat.

Definition at line 5073 of file GeometryCore.h.

5074 { return { GetEndTPCsetID(id), 0 }; }

readout::TPCsetID GetEndTPCsetID(geo::CryostatID const &id) const

Definition: GeometryCore.h:4844

readout::ROPID geo::GeometryCore::GetEndROPID ( readout::TPCsetID const & id ) const

inline

Returns the (possibly invalid) ID after the last readout plane of the specified TPC set.

Definition at line 5082 of file GeometryCore.h.

5083 { return { GetNextID(id), 0 }; }

GeoID GetNextID(GeoID const &id) const

Returns the ID next to the specified one.

Definition: GeometryCore.h:1912

geo::TPCID geo::GeometryCore::GetEndTPCID ( geo::CryostatID const & id ) const

inline

Returns the (possibly invalid) ID after the last TPC of the specified cryostat.

Definition at line 2511 of file GeometryCore.h.

References geo::CryostatID::Cryostat.

2512 { return { id.Cryostat + 1, 0 }; }

readout::TPCsetID geo::GeometryCore::GetEndTPCsetID ( geo::CryostatID const & id ) const

inline

Returns the (possibly invalid) ID after the last TPC set of the specified cryostat.

Definition at line 4844 of file GeometryCore.h.

References geo::CryostatID::Cryostat.

4845 { return { id.Cryostat + 1, 0 }; }

geo::WireID geo::GeometryCore::GetEndWireID ( geo::CryostatID const & id ) const

inline

Returns the (possibly invalid) ID after the last wire in the specified cryostat.

Definition at line 3400 of file GeometryCore.h.

3401 { return { GetEndPlaneID(id), 0 }; }

geo::PlaneID GetEndPlaneID(geo::CryostatID const &id) const

Definition: GeometryCore.h:2994

geo::WireID geo::GeometryCore::GetEndWireID ( geo::TPCID const & id ) const

inline

Returns the (possibly invalid) ID after the last wire of the specified TPC.

Definition at line 3409 of file GeometryCore.h.

3410 { return { geo::PlaneID(GetNextID(id), 0), 0 }; }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

GeoID GetNextID(GeoID const &id) const

Returns the ID next to the specified one.

Definition: GeometryCore.h:1912

geo::WireID geo::GeometryCore::GetEndWireID ( geo::PlaneID const & id ) const

inline

Returns the (possibly invalid) ID after the last wire of the specified wire plane.

Definition at line 3418 of file GeometryCore.h.

3419 { return { GetNextID(id), 0 }; }

GeoID GetNextID(GeoID const &id) const

Returns the ID next to the specified one.

Definition: GeometryCore.h:1912

std::string geo::GeometryCore::GetLArTPCVolumeName ( geo::TPCID const & tpcid ) const

Return the name of specified LAr TPC volume.

Parameters

tpcid	ID of the TPC
tpc	index of TPC in the cryostat
cstat	index of the cryostat

Returns: the name of the specified TPC

This information is used by Geant4 simulation

Todo:

Use a TPCID instead

What if it does not exist?

Definition at line 699 of file GeometryCore.cxx.

References TPC().

Referenced by LArStackingAction::ClassifyNewTrack(), trkf::SpacePts::produce(), trkf::Track3Dreco::produce(), vertex::VertexFinder2D::produce(), trkf::Track3DKalman::produce(), and trkf::Track3DKalmanSPS::produce().

   {
     return std::string(TPC(tpcid).ActiveVolume()->GetName());
   }

std::string geo::GeometryCore::GetLArTPCVolumeName	(	unsigned int const	tpc = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Return the name of specified LAr TPC volume.

Parameters

tpcid	ID of the TPC
tpc	index of TPC in the cryostat
cstat	index of the cryostat

Returns: the name of the specified TPC

This information is used by Geant4 simulation

Todo:

Use a TPCID instead

What if it does not exist?

Definition at line 2775 of file GeometryCore.h.

2776 { return GetLArTPCVolumeName(geo::TPCID(cstat, tpc)); }

geo::GeometryCore::GetLArTPCVolumeName

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

std::string GetLArTPCVolumeName(geo::TPCID const &tpcid) const

Return the name of specified LAr TPC volume.

Definition: GeometryCore.cxx:699

template<typename GeoID >

GeoID geo::GeometryCore::GetNextID ( GeoID const & id ) const

inline

Returns the ID next to the specified one.

Template Parameters

GeoID type of the ID to be returned

Parameters

id	the element ID to be incremented

Returns: ID of the next subelement after id

Definition at line 1912 of file GeometryCore.h.

1913 { auto nextID(id); IncrementID(nextID); return nextID; }

geo::GeometryCore::IncrementID

bool IncrementID(geo::CryostatID &id) const

Definition: GeometryCore.h:5475

template<typename Point = DefaultPoint_t>

Point geo::GeometryCore::GetTPCFrontFaceCenter ( geo::TPCID const & tpcid ) const

inline

Returns the center of side of the detector facing the beam.

Parameters

tpcid	ID of the TPC
tpc	tpc number within the cryostat
cstat	number of cryostat

Returns: the position of center of TPC face toward the beam

Effectively, this is the center of the side of TPC active volume which faces the negative z direction, the first that a beam following

Definition at line 2747 of file GeometryCore.h.

2748 { return TPC(tpcid).GetFrontFaceCenter<Point>(); }

dcel2d::Point

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

Definitions used by the VoronoiDiagram algorithm.

Definition: DCEL.h:34

geo::TPCGeo::GetFrontFaceCenter

Point GetFrontFaceCenter() const

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

Definition: TPCGeo.h:258

TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified TPC.

Definition: GeometryCore.h:2393

template<typename Point = DefaultPoint_t>

Point geo::GeometryCore::GetTPCFrontFaceCenter	(	unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the center of side of the detector facing the beam.

Parameters

tpcid	ID of the TPC
tpc	tpc number within the cryostat
cstat	number of cryostat

Returns: the position of center of TPC face toward the beam

Effectively, this is the center of the side of TPC active volume which faces the negative z direction, the first that a beam following

Definition at line 2751 of file GeometryCore.h.

2752 { return GetTPCFrontFaceCenter<Point>(geo::TPCID(cstat, tpc)); }

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

const std::string geo::GeometryCore::GetWorldVolumeName ( ) const

Return the name of the world volume (needed by Geant4 simulation)

Definition at line 578 of file GeometryCore.cxx.

References DetectorEnclosureBox().

Referenced by evgen::CosmicsGen::reconfigure(), WorldBox(), and WorldVolume().

   {
     // For now, and possibly forever, this is a constant (given the
     // definition of "nodeNames" above).
     return std::string("volWorld");
   }

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

Convert unique channel to hardware channel.

Definition at line 242 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->HardwareChannelFromOpChannel(opChannel);
   }

bool geo::GeometryCore::HasChannel ( raw::ChannelID_t channel ) const

Returns whether the specified channel exists and is valid.

Parameters

channel the ID of the channel

Returns: whether the specified channel exists

A channel is defined as existing and valid if its ID is not invalid and if the channel is physical.

Definition at line 566 of file GeometryCore.cxx.

References fChannelMapAlg.

                                                             {
     return fChannelMapAlg->HasChannel(channel);
   } // GeometryCore::HasChannel()

bool geo::GeometryCore::HasCryostat ( geo::CryostatID const & cryoid ) const

inline

Returns whether we have the specified cryostat.

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

Definition at line 1981 of file GeometryCore.h.

References geo::CryostatID::Cryostat.

1982 { return cryoid.Cryostat < Ncryostats(); }

geo::GeometryCore::HasCryostat

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

bool geo::GeometryCore::HasElement ( geo::CryostatID const & cryoid ) const

inline

Returns whether we have the specified cryostat.

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

Definition at line 1983 of file GeometryCore.h.

1984 { return HasCryostat(cryoid); }

bool HasCryostat(geo::CryostatID const &cryoid) const

Returns whether we have the specified cryostat.

Definition: GeometryCore.h:1981

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

inline

Returns whether we have the specified TPC.

Definition at line 2374 of file GeometryCore.h.

2374 { return HasTPC(tpcid); }

geo::GeometryCore::HasTPC

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

Returns whether we have the specified TPC.

Definition: GeometryCore.h:2367

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

inline

Returns whether we have the specified plane.

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

Definition at line 2923 of file GeometryCore.h.

2924 { return HasPlane(planeid); }

geo::GeometryCore::HasPlane

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

Returns whether we have the specified plane.

Definition: GeometryCore.h:2918

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

inline

Returns whether we have the specified wire.

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

Definition at line 3339 of file GeometryCore.h.

3339 { return HasWire(wireid); }

geo::GeometryCore::HasWire

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

Returns whether we have the specified wire.

Definition: GeometryCore.h:3334

bool geo::GeometryCore::HasElement ( readout::TPCsetID const & tpcsetid ) const

inline

Returns whether we have the specified TPC set.

Definition at line 4791 of file GeometryCore.h.

4792 { return HasTPCset(tpcsetid); }

geo::GeometryCore::HasTPCset

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

Definition: GeometryCore.cxx:1776

bool geo::GeometryCore::HasElement ( readout::ROPID const & ropid ) const

inline

Returns whether we have the specified readout plane

Returns: whether the readout plane is valid and exists

Definition at line 4995 of file GeometryCore.h.

4995 { return HasROP(ropid); }

geo::GeometryCore::HasROP

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

Definition: GeometryCore.cxx:1820

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

inline

Returns whether we have the specified plane.

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

Definition at line 2918 of file GeometryCore.h.

References geo::TPCGeo::HasPlane().

       {
         geo::TPCGeo const* pTPC = TPCPtr(planeid);
         return pTPC? pTPC->HasPlane(planeid): false;
       }

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

Returns whether we have the specified readout plane

Returns: whether the readout plane is valid and exists

Definition at line 1820 of file GeometryCore.cxx.

References fChannelMapAlg.

                                                            {
     return fChannelMapAlg->HasROP(ropid);
   } // GeometryCore::HasROP()

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

inline

Returns whether we have the specified TPC.

Definition at line 2367 of file GeometryCore.h.

References geo::CryostatGeo::HasTPC().

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), ems::MultiEMShowers::convCluster(), and evd::TWQProjectionView::SelectTPC().

       {
         CryostatGeo const* pCryo = CryostatPtr(tpcid);
         return pCryo? pCryo->HasTPC(tpcid): false;
       }

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

Returns whether we have the specified TPC set

Returns: whether the TPC set is valid and exists

Definition at line 1776 of file GeometryCore.cxx.

References fChannelMapAlg, and FindTPCsetAtPosition().

                                                                     {
     return fChannelMapAlg->HasTPCset(tpcsetid);
   } // GeometryCore::HasTPCset()

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

inline

Returns whether we have the specified wire.

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

Definition at line 3334 of file GeometryCore.h.

References geo::PlaneGeo::HasWire().

Referenced by tca::FillmAllTraj(), tca::Find3DVertices(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), pma::ProjectionMatchingAlg::validate_on_adc_test(), and WireIDIntersectionCheck().

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

bool geo::GeometryCore::IncrementID ( geo::CryostatID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5475 of file GeometryCore.h.

                                                                 {
   ++id.Cryostat;
   if (id) id.isValid = HasCryostat(id); // if invalid already, it stays so
   return bool(id);
 } // geo::GeometryCore::IncrementID(geo::CryostatID)

bool geo::GeometryCore::IncrementID ( geo::TPCID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5481 of file GeometryCore.h.

                                                            {
   unsigned int const nTPCsInCryo = NTPC(id);
   if (++id.TPC < nTPCsInCryo) return bool(id); // if was invalid, it stays so
   // no more TPCs in this cryostat
   id.TPC = 0;
   return IncrementID(id.asCryostatID()); // also sets validity
 } // geo::GeometryCore::IncrementID(geo::TPCID)

bool geo::GeometryCore::IncrementID ( geo::PlaneID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5489 of file GeometryCore.h.

                                                              {
   // this implementation is non-optimal, in that the cryostat lookup is
   // performed both here and, potentially, in IncrementID(TPCID)
   unsigned int const nPlanesInTPC = Nplanes(id);
   if (++id.Plane < nPlanesInTPC) return bool(id); // if was invalid, stays so
   // no more planes in this TPCs
   id.Plane = 0;
   return IncrementID(id.asTPCID()); // also sets validity
 } // geo::GeometryCore::IncrementID(geo::PlaneID)

bool geo::GeometryCore::IncrementID ( geo::WireID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5499 of file GeometryCore.h.

                                                             {
   // this implementation is non-optimal, in that the TPC lookup is
   // performed both here and, potentially, in IncrementID(PlaneID)
   unsigned int const nWiresInPlane = Nwires(id);
   if (++id.Wire < nWiresInPlane) return bool(id); // if was invalid, stays so
   // no more wires in this plane
   id.Wire = 0;
   return IncrementID(id.asPlaneID()); // also sets validity
 } // geo::GeometryCore::IncrementID(geo::WireID)

bool geo::GeometryCore::IncrementID ( readout::TPCsetID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5509 of file GeometryCore.h.

                                                                   {
   unsigned int const nTPCsetsInCryo = NTPCsets(id);
   if (++id.TPCset < nTPCsetsInCryo)
     return bool(id); // if was invalid, it stays so
   // no more TPC sets in this cryostat
   id.TPCset = 0;
   return IncrementID(id.asCryostatID()); // also sets validity
 } // geo::GeometryCore::IncrementID(readout::TPCsetID)

bool geo::GeometryCore::IncrementID ( readout::ROPID & id ) const

inline

Sets the ID to the ID after the specified one.

Returns: whether the ID is actually valid (validity flag is also set)

Definition at line 5518 of file GeometryCore.h.

References WireEndPoints().

                                                                {
   // this implementation is non-optimal, in that the cryostat lookup is
   // performed both here and, potentially, in IncrementID(TPCsetID)
   unsigned int const nROPinTPC = NROPs(id);
   if (++id.ROP < nROPinTPC) return bool(id); // if was invalid, stays so
   // no more readout planes in this TPC set
   id.ROP = 0;
   return IncrementID(id.asTPCsetID()); // also sets validity
 } // geo::GeometryCore::IncrementID(readout::ROPID)

bool geo::GeometryCore::IntersectionPoint	(	geo::WireID const &	wid1,
		geo::WireID const &	wid2,
		double &	y,
		double &	z
	)		const

Returns the intersection point of two wires.

Parameters

wid1	ID of the first wire
wid2	ID of the other wire
y	(output) y coordinate of the intersection point
z	(output) z coordinate of the intersection point

Returns: whether an intersection was found within the TPC

See also: WireIDsIntersect()

The behaviour of this method reflects the one of WireIDsIntersect(), which supersedes this one.

To test if the result is infinity, use e.g. std::isfinite(y).

Deprecated:: This method uses arbitrary assumptions and should not be used. Use WireIDsIntersect() returning a vector, instead.

Definition at line 1749 of file GeometryCore.cxx.

References WireIDsIntersect(), geo::WireIDIntersection::y, and geo::WireIDIntersection::z.

Referenced by pma::ProjectionMatchingAlg::addEndpointRef(), sppt::SpacePointAlg_TimeSort::fillCoordinatesArrays(), filter::MuonFilter::filter(), tca::Find3DVertices(), cmtool::CFAlgoStartPointCompat::Float(), cmtool::CFAlgoStartPointMatch::Float(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), trkf::BezierTrackerAlgorithm::MakeTracks(), cluster::ClusterMatchAlg::Match_RoughZ(), and trkf::TrackTrajectoryAlg::ShortTrackTrajectory().

   {
     geo::WireIDIntersection widIntersect;
     bool const found = WireIDsIntersect(wid1, wid2, widIntersect);
     y = widIntersect.y;
     z = widIntersect.z;
     return found;
   }

bool geo::GeometryCore::IntersectionPoint	(	unsigned int	wire1,
		unsigned int	wire2,
		unsigned int	plane1,
		unsigned int	plane2,
		unsigned int	cstat,
		unsigned int	tpc,
		double &	y,
		double &	z
	)		const

inline

Returns the intersection point of two wires.

Parameters

wire1	wire index of the first wire
wire2	wire index of the other wire
plane1	plane index of the first wire
plane2	plane index of the other wire
cstat	cryostat number
tpc	tpc number within the cryostat where the planes belong
y	(output) y coordinate of the intersection point
z	(output) z coordinate of the intersection point

Returns: whether an intersection was found

No check is performed, not any information provided, about the validity of the result.

Deprecated:: This method uses arbitrary assumptions and should not be used. Use WireIDsIntersect() returning a vector, instead.

Definition at line 4188 of file GeometryCore.h.

       {
         return IntersectionPoint(
           geo::WireID(cstat, tpc, plane1, wire1),
           geo::WireID(cstat, tpc, plane2, wire2),
           y, z
           );
       }

bool geo::GeometryCore::IntersectLines	(	double	A_start_x,
		double	A_start_y,
		double	A_end_x,
		double	A_end_y,
		double	B_start_x,
		double	B_start_y,
		double	B_end_x,
		double	B_end_y,
		double &	x,
		double &	y
	)		const

Computes the intersection between two lines on a plane.

Parameters

A_start_x	x coordinate of one point of the first segment
A_start_y	y coordinate of one point of the first segment
A_end_x	x coordinate of another point of the first segment
A_end_y	y coordinate of another point of the first segment
B_start_x	x coordinate of one point of the second segment
B_start_y	y coordinate of one point of the second segment
B_end_x	x coordinate of another point of the second segment
B_end_y	y coordinate of another point of the second segment
x	_(output)_ variable to store the x coordinate of intersection
y	_(output)_ variable to store the y coordinate of intersection

Returns: whether intersection exists

The order of the ends is not relevant. The return value is false if the two segments are parallel. In that case, x and y variables are not changed. Otherwise, they hold the intersection coordinate, even if the intersection point is beyond one or both the segments.

Definition at line 1384 of file GeometryCore.cxx.

References B, coordIs, and lar::util::RealComparisons< RealType >::zero().

Referenced by IntersectSegments().

           {
     
     // Equation from http://en.wikipedia.org/wiki/Line%E2%80%93line_intersection
     // T.Yang Nov, 2014
     // Notation: x => coordinate orthogonal to the drift direction and to the beam direction
     //           y => direction orthogonal to the previous and to beam direction
     
     double const denom = (A_start_x - A_end_x)*(B_start_y - B_end_y)
       - (A_start_y - A_end_y)*(B_start_x - B_end_x);
     
     if (coordIs.zero(denom)) return false;
     
     double const A = (A_start_x * A_end_y - A_start_y * A_end_x) / denom;
     double const B = (B_start_x * B_end_y - B_start_y * B_end_x) / denom;
     
     x = (B_start_x - B_end_x) * A - (A_start_x - A_end_x) * B;
     y = (B_start_y - B_end_y) * A - (A_start_y - A_end_y) * B;
     
     return true;
     
   } // GeometryCore::IntersectLines()

bool geo::GeometryCore::IntersectSegments	(	double	A_start_x,
		double	A_start_y,
		double	A_end_x,
		double	A_end_y,
		double	B_start_x,
		double	B_start_y,
		double	B_end_x,
		double	B_end_y,
		double &	x,
		double &	y
	)		const

Computes the intersection between two segments on a plane.

Parameters

A_start_x	x coordinate of the start of the first segment
A_start_y	y coordinate of the start of the first segment
A_end_x	x coordinate of the end of the first segment
A_end_y	y coordinate of the end of the first segment
B_start_x	x coordinate of the start of the second segment
B_start_y	y coordinate of the start of the second segment
B_end_x	x coordinate of the end of the second segment
B_end_y	y coordinate of the end of the second segment
x	_(output)_ variable to store the x coordinate of intersection
y	_(output)_ variable to store the y coordinate of intersection

Returns: whether intersection exists and is on both segments

The order of the ends is not relevant. The return value is false if the two segments are parallel, or if their intersection point is not on both the segments. If the segments are parallel, x and y variables are not changed. Otherwise, they hold the intersection coordinate, even if the intersection point is beyond one or both the segments.

Definition at line 1411 of file GeometryCore.cxx.

References c1, c2, geo::vect::cross(), geo::WireGeo::Direction(), geo::vect::dot(), geo::GeometryCore::Segment< Point >::end(), geo::WireGeo::GetCenter(), geo::WireGeo::HalfL(), IntersectLines(), geo::TPCID::InvalidID, PointWithinSegments(), geo::sqr(), geo::GeometryCore::Segment< Point >::start(), ThirdPlane(), geo::vect::toTVector3(), geo::TPCID::TPC, geo::WireIDIntersection::TPC, Wire(), WireEndPoints(), WireIDIntersectionCheck(), WireIDsIntersect(), geo::WireIDIntersection::y, and geo::WireIDIntersection::z.

           {
     
     bool bCross = IntersectLines(
       A_start_x, A_start_y, A_end_x, A_end_y,
       B_start_x, B_start_y, B_end_x, B_end_y,
       x, y
       );
     
     if (bCross) {
       mf::LogWarning("IntersectSegments") << "The segments are parallel!";
       return false;
     }
     
     return PointWithinSegments(
       A_start_x, A_start_y, A_end_x, A_end_y,
       B_start_x, B_start_y, B_end_x, B_end_y,
       x, y
       );
     
   } // GeometryCore::IntersectSegments()

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

Is this a valid OpChannel number?

Definition at line 249 of file GeometryCore.cxx.

References fChannelMapAlg, and NOpDets().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison(), lariov::SIOVPmtGainProvider::Reconfigure(), and opdet::RunHitFinder().

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

IteratorBox< cryostat_id_iterator, &GeometryCore::begin_cryostat_id, &GeometryCore::end_cryostat_id > geo::GeometryCore::IterateCryostatIDs ( ) const

inline

Enables ranged-for loops on all cryostat IDs of the detector.

Returns: an object suitable for ranged-for loops on all cryostat IDs

Example of usage:

for (geo::CryostatID const& cID: geom->IterateCryostatIDs()) {
  geo::CryostatGeo const& Cryo = geom->Cryostat(cID);

  // useful code here

} // for all cryostats

Definition at line 2149 of file GeometryCore.h.

References begin_cryostat(), and end_cryostat().

2149 { return { this }; }

IteratorBox< cryostat_iterator, &GeometryCore::begin_cryostat, &GeometryCore::end_cryostat > geo::GeometryCore::IterateCryostats ( ) const

inline

Enables ranged-for loops on all cryostats of the detector.

Returns: an object suitable for ranged-for loops on all cryostats

Example of usage:

for (geo::CryostatGeo const& Cryo: geom->IterateCryostats()) {

  // useful code here

} // for all cryostats

Definition at line 2168 of file GeometryCore.h.

Referenced by PositionToCryostatPtr().

2168 { return { this }; }

IteratorBox< plane_id_iterator, &GeometryCore::begin_plane_id, &GeometryCore::end_plane_id > geo::GeometryCore::IteratePlaneIDs ( ) const

inline

Enables ranged-for loops on all plane IDs of the detector.

Returns: an object suitable for ranged-for loops on all plane IDs

Example of usage:

for (geo::PlaneID const& pID: geom->IteratePlaneIDs()) {
  geo::PlaneGeo const& Plane = geom->Plane(pID);
  // useful code here
} // for all plane IDs

Definition at line 3078 of file GeometryCore.h.

References begin_plane_id(), and end_plane_id().

Referenced by vertex::FeatureVertexFinderAna::analyze(), corner::CornerFinderAlg::create_smaller_histos(), sim::details::createPlaneIndexMap(), corner::CornerFinderAlg::get_feature_points(), corner::CornerFinderAlg::get_feature_points_LineIntegralScore(), cluster::DBScan3DAlg::init(), vertex::VertexMatch::produce(), and vertex::HarrisVertexFinder::produce().

3078 { return { this }; }

LocalIteratorBox< plane_id_iterator, geo::CryostatID, &GeometryCore::begin_plane_id, &GeometryCore::end_plane_id > geo::GeometryCore::IteratePlaneIDs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all plane IDs of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the plane IDs of

Returns: an object suitable for ranged-for loops on plane IDs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::PlaneID const& pID: geom->IteratePlaneIDs(cid)) {
  geo::PlaneGeo const& plane = geom->Plane(pID);
  // useful code here
} // for all planes in cryostat #1

Definition at line 3103 of file GeometryCore.h.

References begin_plane_id(), and end_plane_id().

3103 { return { this, cid }; }

LocalIteratorBox< plane_id_iterator, geo::TPCID, &GeometryCore::begin_plane_id, &GeometryCore::end_plane_id > geo::GeometryCore::IteratePlaneIDs ( geo::TPCID const & tid ) const

inline

Enables ranged-for loops on all plane IDs of the specified TPC.

Parameters

tid	the ID of the TPC to loop the plane IDs of

Returns: an object suitable for ranged-for loops on plane IDs

If the TPC ID is invalid, the effect is undefined.

Example of usage:

geo::TPCID tid{ 0, 1 }; // C:0 T:1 (hope it exists!)
for (geo::PlaneID const& pID: geom->IteratePlaneIDs(tid)) {
  geo::PlaneGeo const& plane = geom->Plane(pID);
  // useful code here
} // for all planes in C:0 T:1

Definition at line 3127 of file GeometryCore.h.

References begin_plane(), and end_plane().

3127 { return { this, tid }; }

void geo::GeometryCore::IteratePlaneIDs ( geo::PlaneID const & pid ) const

delete

IteratePlaneIDs() is not supported on plane IDs.

void geo::GeometryCore::IteratePlaneIDs ( geo::WireID const & pid ) const

delete

IteratePlaneIDs() is not supported on wire IDs.

void geo::GeometryCore::IteratePlaneIDs ( readout::TPCsetID const & ) const

delete

IteratePlaneIDs() is not supported on readout IDs.

void geo::GeometryCore::IteratePlaneIDs ( readout::ROPID const & ) const

delete

IteratePlaneIDs() is not supported on readout IDs.

IteratorBox< plane_iterator, &GeometryCore::begin_plane, &GeometryCore::end_plane > geo::GeometryCore::IteratePlanes ( ) const

inline

Enables ranged-for loops on all planes of the detector.

Returns: an object suitable for ranged-for loops on all planes

Example of usage:

for (geo::PlaneGeo const& Plane: geom->IteratePlanes()) {
  // useful code here
} // for all planes

Definition at line 3160 of file GeometryCore.h.

References begin_plane(), and end_plane().

Referenced by mvapid::MVAAlg::GetWireNormals().

3160 { return { this }; }

LocalIteratorBox< plane_iterator, geo::CryostatID, &GeometryCore::begin_plane, &GeometryCore::end_plane > geo::GeometryCore::IteratePlanes ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all planes of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the planes of

Returns: an object suitable for ranged-for loops on planes

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::PlaneGeo const& plane: geom->IteratePlanes(cid)) {
  // useful code here
} // for planes in cryostat 1

Definition at line 3183 of file GeometryCore.h.

References begin_plane(), and end_plane().

3183 { return { this, cid }; }

LocalIteratorBox< plane_iterator, geo::TPCID, &GeometryCore::begin_plane, &GeometryCore::end_plane > geo::GeometryCore::IteratePlanes ( geo::TPCID const & tid ) const

inline

Enables ranged-for loops on all planes of the specified TPC.

Parameters

tid	the ID of the TPC to loop the planes of

Returns: an object suitable for ranged-for loops on planes

If the TPC ID is invalid, the effect is undefined.

Example of usage:

geo::TPCID tid{ 0, 1 }; // C:0 T:1 (hope it exists!)
for (geo::PlaneGeo const& plane: geom->IteratePlanes(tid)) {
  // useful code here
} // for planes in C:0 T:1

Definition at line 3206 of file GeometryCore.h.

3206 { return { this, tid }; }

void geo::GeometryCore::IteratePlanes ( geo::PlaneID const & pid ) const

delete

IteratePlanes() is not supported on plane IDs.

void geo::GeometryCore::IteratePlanes ( geo::WireID const & pid ) const

delete

IteratePlanes() is not supported on wire IDs.

void geo::GeometryCore::IteratePlanes ( readout::TPCsetID const & ) const

delete

IteratePlanes() is not supported on readout IDs.

void geo::GeometryCore::IteratePlanes ( readout::ROPID const & ) const

delete

IteratePlanes() is not supported on readout IDs.

IteratorBox< ROP_id_iterator, &GeometryCore::begin_ROP_id, &GeometryCore::end_ROP_id > geo::GeometryCore::IterateROPIDs ( ) const

inline

Enables ranged-for loops on all readout plane IDs of the detector.

Returns: an object suitable for ranged-for loops on all ROP IDs

Example of usage:

for (readout::ROPID const& rID: geom->IterateROPIDs()) {
  // useful code here
} // for all ROPs

Definition at line 5131 of file GeometryCore.h.

References begin_ROP_id(), and end_ROP_id().

5131 { return { this }; }

LocalIteratorBox< ROP_id_iterator, geo::CryostatID, &GeometryCore::begin_ROP_id, &GeometryCore::end_ROP_id > geo::GeometryCore::IterateROPIDs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all readout plane IDs of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the readout plane IDs of

Returns: an object suitable for ranged-for loops on readout plane IDs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (readout::ROPID const& rID: geom->IterateROPIDs(cid)) {
  // useful code here
} // for all readout planes in cryostat #1

Definition at line 5155 of file GeometryCore.h.

References begin_ROP_id(), and end_ROP_id().

5155 { return { this, cid }; }

LocalIteratorBox< ROP_id_iterator, readout::TPCsetID, &GeometryCore::begin_ROP_id, &GeometryCore::end_ROP_id > geo::GeometryCore::IterateROPIDs ( readout::TPCsetID const & sid ) const

inline

Enables ranged-for loops on all readout plane IDs of the specified TPC set.

Parameters

sid	the ID of the TPC set to loop the readout plane IDs of

Returns: an object suitable for ranged-for loops on readout plane IDs

If the TPC set ID is invalid, the effect is undefined.

Example of usage:

readout::TPCsetID sid{ 0, 1 }; // C:0 S:1 (hope it exists!)
for (readout::ROPID const& rID: geom->IterateROPIDs(sid)) {
  // useful code here
} // for all readout planes in C:0 S:1

Definition at line 5179 of file GeometryCore.h.

References max, min, and fhicl::detail::atom::value().

5179 { return { this, sid }; }

IteratorBox< TPC_id_iterator, &GeometryCore::begin_TPC_id, &GeometryCore::end_TPC_id > geo::GeometryCore::IterateTPCIDs ( ) const

inline

Enables ranged-for loops on all TPC IDs of the detector.

Returns: an object suitable for ranged-for loops on all TPC IDs

Example of usage:

for (geo::TPCID const& tID: geom->IterateTPCIDs()) {
  geo::TPCGeo const& TPC = geom->TPC(tID);
  // useful code here
} // for all TPC

Definition at line 2567 of file GeometryCore.h.

References begin_TPC_id(), and end_TPC_id().

Referenced by tca::InsideTPC(), and tca::TruthMatcher::MatchAndSum().

2567 { return { this }; }

LocalIteratorBox< TPC_id_iterator, geo::CryostatID, &GeometryCore::begin_TPC_id, &GeometryCore::end_TPC_id > geo::GeometryCore::IterateTPCIDs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all TPC IDs of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the TPC IDs of

Returns: an object suitable for ranged-for loops on TPC IDs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::TPCID const& tID: geom->IterateTPCIDs(cid)) {
  geo::TPCGeo const& TPC = geom->TPC(tID);
  // useful code here
} // for all TPC in cryostat #1

Definition at line 2591 of file GeometryCore.h.

References begin_TPC(), and end_TPC().

2591 { return { this, cid }; }

void geo::GeometryCore::IterateTPCIDs ( geo::TPCID const & pid ) const

delete

IterateTPCIDs() is not supported on TPC IDs.

void geo::GeometryCore::IterateTPCIDs ( geo::PlaneID const & pid ) const

delete

IterateTPCIDs() is not supported on plane IDs.

void geo::GeometryCore::IterateTPCIDs ( geo::WireID const & pid ) const

delete

IterateTPCIDs() is not supported on wire IDs.

void geo::GeometryCore::IterateTPCIDs ( readout::TPCsetID const & ) const

delete

IterateTPCIDs() is not supported on readout IDs.

void geo::GeometryCore::IterateTPCIDs ( readout::ROPID const & ) const

delete

IterateTPCIDs() is not supported on readout IDs.

IteratorBox<TPC_iterator, &GeometryCore::begin_TPC, &GeometryCore::end_TPC> geo::GeometryCore::IterateTPCs ( ) const

inline

Enables ranged-for loops on all TPCs of the detector.

Returns: an object suitable for ranged-for loops on all TPCs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

for (geo::TPCGeo const& TPC: geom->IterateTPCs()) {
  // useful code here
} // for TPCs

Definition at line 2625 of file GeometryCore.h.

References begin_TPC(), and end_TPC().

Referenced by ems::MCinfo::Findtpcborders(), and UpdateAfterSorting().

2625 { return { this }; }

LocalIteratorBox< TPC_iterator, geo::CryostatID, &GeometryCore::begin_TPC, &GeometryCore::end_TPC > geo::GeometryCore::IterateTPCs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all TPCs of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the TPCs of

Returns: an object suitable for ranged-for loops on TPCs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::TPCGeo const& TPC: geom->IterateTPCs(cid)) {
  // useful code here
} // for TPCs in cryostat 1

Definition at line 2648 of file GeometryCore.h.

2648 { return { this, cid }; }

void geo::GeometryCore::IterateTPCs ( geo::TPCID const & pid ) const

delete

IterateTPCs() is not supported on TPC IDs.

void geo::GeometryCore::IterateTPCs ( geo::PlaneID const & pid ) const

delete

IterateTPCs() is not supported on plane IDs.

void geo::GeometryCore::IterateTPCs ( geo::WireID const & pid ) const

delete

IterateTPCs() is not supported on wire IDs.

void geo::GeometryCore::IterateTPCs ( readout::TPCsetID const & ) const

delete

IterateTPCs() is not supported on readout IDs.

void geo::GeometryCore::IterateTPCs ( readout::ROPID const & ) const

delete

IterateTPCs() is not supported on readout IDs.

IteratorBox< TPCset_id_iterator, &GeometryCore::begin_TPCset_id, &GeometryCore::end_TPCset_id > geo::GeometryCore::IterateTPCsetIDs ( ) const

inline

Enables ranged-for loops on all TPC set IDs of the detector.

Returns: an object suitable for ranged-for loops on all TPC set IDs

Example of usage:

for (readout::TPCsetID const& sID: geom->IterateTPCsetIDs()) {
  // useful code here
} // for all TPC sets

Definition at line 4884 of file GeometryCore.h.

References begin_TPCset_id(), and end_TPCset_id().

4884 { return { this }; }

LocalIteratorBox< TPCset_id_iterator, geo::CryostatID, &GeometryCore::begin_TPCset_id, &GeometryCore::end_TPCset_id > geo::GeometryCore::IterateTPCsetIDs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all TPC set IDs of the specified cryostat.

Parameters

cid	the ID of the cryostat to loop the TPC set IDs of

Returns: an object suitable for ranged-for loops on TPC set IDs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (readout::TPCsetID const& tID: geom->IterateTPCsetIDs(cid)) {
  // useful code here
} // for all TPC sets in cryostat #1

Definition at line 4908 of file GeometryCore.h.

4908 { return { this, cid }; }

IteratorBox< wire_id_iterator, &GeometryCore::begin_wire_id, &GeometryCore::end_wire_id > geo::GeometryCore::IterateWireIDs ( ) const

inline

Enables ranged-for loops on all wire IDs of the detector.

Returns: an object suitable for ranged-for loops on all wire IDs

Example of usage:

for (geo::WireID const& wID: geom->IterateWireIDs()) {
  geo::WireGeo const& Wire = geom->Wire(wID);
  // useful code here
} // for all wires

Definition at line 3505 of file GeometryCore.h.

References begin_wire_id(), and end_wire_id().

Referenced by cluster::DBScan3DAlg::init().

3505 { return { this }; }

LocalIteratorBox< wire_id_iterator, geo::CryostatID, &GeometryCore::begin_wire_id, &GeometryCore::end_wire_id > geo::GeometryCore::IterateWireIDs ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all wire IDs of specified cryostat.

Parameters

cid	the ID of the cryostat to loop the wires of

Returns: an object suitable for ranged-for loops on cryostat wire IDs

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::WireID const& wID: geom->IterateWireIDs(cid)) {
  geo::WireGeo const& Wire = geom->Wire(wID);
  // useful code here
} // for all wires

Definition at line 3530 of file GeometryCore.h.

References begin_wire_id(), and end_wire_id().

3530 { return { this, cid }; }

LocalIteratorBox< wire_id_iterator, geo::TPCID, &GeometryCore::begin_wire_id, &GeometryCore::end_wire_id > geo::GeometryCore::IterateWireIDs ( geo::TPCID const & tid ) const

inline

Enables ranged-for loops on all wire IDs of specified TPC.

Parameters

tid	the ID of the TPC to loop the wires of

Returns: an object suitable for ranged-for loops on TPC wire IDs

If the TPC ID is invalid, the effect is undefined.

Example of usage:

geo::TPCID tid{0, 1}; // C:0 T:1 (hope it exists!)
for (geo::WireID const& wID: geom->IterateWireIDs(tid)) {
  geo::WireGeo const& Wire = geom->Wire(wID);
  // useful code here
} // for all wires in C:0 T:1

Definition at line 3555 of file GeometryCore.h.

References begin_wire_id(), and end_wire_id().

3555 { return { this, tid }; }

LocalIteratorBox< wire_id_iterator, geo::PlaneID, &GeometryCore::begin_wire_id, &GeometryCore::end_wire_id > geo::GeometryCore::IterateWireIDs ( geo::PlaneID const & pid ) const

inline

Enables ranged-for loops on all wire IDs of specified wire plane.

Parameters

pid	the ID of the wire plane to loop the wires of

Returns: an object suitable for ranged-for loops on plane wire IDs

If the wire plane ID is invalid, the effect is undefined.

Example of usage:

geo::PlaneID pid{0, 0, 1}; // C:0 T:0 P:1
for (geo::WireID const& wID: geom->IterateWireIDs(pid)) {
  geo::WireGeo const& Wire = geom->Wire(wID);
  // useful code here
} // for all wires in C:0 T:0 P:1

Definition at line 3580 of file GeometryCore.h.

References begin_wire(), and end_wire().

3580 { return { this, pid }; }

pid

Definition: Chi2ParticleID_module.cc:41

void geo::GeometryCore::IterateWireIDs ( geo::WireID const & pid ) const

delete

IterateWireIDs() is not supported on wire IDs.

void geo::GeometryCore::IterateWireIDs ( readout::TPCsetID const & ) const

delete

IterateWireIDs() is not supported on readout IDs.

void geo::GeometryCore::IterateWireIDs ( readout::ROPID const & ) const

delete

IterateWireIDs() is not supported on readout IDs.

IteratorBox< wire_iterator, &GeometryCore::begin_wire, &GeometryCore::end_wire > geo::GeometryCore::IterateWires ( ) const

inline

Enables ranged-for loops on all wires of the detector.

Returns: an object suitable for ranged-for loops on all wires

Example of usage:

for (geo::WireGeo const& Wire: geom->IterateWires()) {
  // useful code here
} // for all wires

Definition at line 3609 of file GeometryCore.h.

References begin_wire(), and end_wire().

3609 { return { this }; }

LocalIteratorBox< wire_iterator, geo::CryostatID, &GeometryCore::begin_wire, &GeometryCore::end_wire > geo::GeometryCore::IterateWires ( geo::CryostatID const & cid ) const

inline

Enables ranged-for loops on all wires of specified cryostat.

Parameters

cid	the ID of the cryostat to loop the wires of

Returns: an object suitable for ranged-for loops on cryostat wires

If the cryostat ID is invalid, the effect is undefined.

Example of usage:

geo::CryostatID cid{1}; // cryostat #1 (hope it exists!)
for (geo::WireID const& Wire: geom->IterateWires(cid)) {
  // useful code here
} // for all wires

Definition at line 3633 of file GeometryCore.h.

References begin_wire(), and end_wire().

3633 { return { this, cid }; }

LocalIteratorBox< wire_iterator, geo::TPCID, &GeometryCore::begin_wire, &GeometryCore::end_wire > geo::GeometryCore::IterateWires ( geo::TPCID const & tid ) const

inline

Enables ranged-for loops on all wires of specified TPC.

Parameters

tid	the ID of the TPC to loop the wires of

Returns: an object suitable for ranged-for loops on TPC wires

If the TPC ID is invalid, the effect is undefined.

Example of usage:

geo::TPCID tid{0, 1}; // C:0 T:1 (hope it exists!)
for (geo::WireID const& Wire: geom->IterateWires(tid)) {
  // useful code here
} // for all wires in C:0 T:1

Definition at line 3657 of file GeometryCore.h.

References begin_wire(), and end_wire().

3657 { return { this, tid }; }

LocalIteratorBox< wire_iterator, geo::PlaneID, &GeometryCore::begin_wire, &GeometryCore::end_wire > geo::GeometryCore::IterateWires ( geo::PlaneID const & tid ) const

inline

Enables ranged-for loops on all wires of specified wire plane.

Parameters

pid	the ID of the wire plane to loop the wires of

Returns: an object suitable for ranged-for loops on plane wires

If the wire plane ID is invalid, the effect is undefined.

Example of usage:

geo::PlaneID pid{0, 1}; // C:0 T:0 P:1
for (geo::WireID const& Wire: geom->IterateWires(pid)) {
  // useful code here
} // for all wires in C:0 T:0 T:1

Definition at line 3681 of file GeometryCore.h.

3681 { return { this, tid }; }

void geo::GeometryCore::IterateWires ( geo::WireID const & pid ) const

delete

IterateWires() is not supported on wire IDs.

void geo::GeometryCore::IterateWires ( readout::TPCsetID const & ) const

delete

IterateWires() is not supported on readout IDs.

void geo::GeometryCore::IterateWires ( readout::ROPID const & ) const

delete

IterateWires() is not supported on readout IDs.

void geo::GeometryCore::LoadGeometryFile	(	std::string	gdmlfile,
		std::string	rootfile,
		bool	bForceReload = `false`
	)

Loads the geometry information from the specified files.

Parameters

gdmlfile	path to file to be used for Geant4 simulation
rootfile	path to file for internal geometry representation
bForceReload	reload even if there is already a valid geometry

See also: ApplyChannelMap()

Both paths must directly resolve to an available file, as no search is performed for them.

The gdmlfile parameter does not have to necessarily be in GDML format, as long as it's something supported by Geant4. This file is not used by the geometry, but its path is provided on request by the simulation modules (see LArSoft LArG4 module). The rootfile also does not need to be a ROOT file, but just anything that TGeoManager::Import() supports. This file is parsed immediately and the internal geometry representation is built out of it.

Note: After calling this method, the detector geometry information can be considered complete, but the geometry service provider is not fully initialized yet, since it's still necessary to provide or update the channel mapping.

Definition at line 92 of file GeometryCore.cxx.

References ClearGeometry(), fGDMLfile, FindAuxDet(), FindCryostat(), fROOTfile, and MaxWireDepthInGDML.

Referenced by geo::Geometry::LoadNewGeometry().

     {
     
     if (gdmlfile.empty()) {
       throw cet::exception("GeometryCore")
         << "No GDML Geometry file specified!\n";
     }
     
     if (rootfile.empty()) {
       throw cet::exception("GeometryCore")
         << "No ROOT Geometry file specified!\n";
     }
     
     ClearGeometry();
     
     // Open the GDML file, and convert it into ROOT TGeoManager format.
     // Then lock the gGeoManager to prevent future imports, for example
     // in AuxDetGeometry
     if( !gGeoManager || bForceReload ){
       if (gGeoManager) TGeoManager::UnlockGeometry();
       TGeoManager::Import(rootfile.c_str());
       gGeoManager->LockGeometry();
     }
 
     std::vector<const TGeoNode*> path(MaxWireDepthInGDML);
     path[0] = gGeoManager->GetTopNode();
     FindCryostat(path, 0);
     FindAuxDet(path, 0);
     
     fGDMLfile = gdmlfile;
     fROOTfile = rootfile;
     
     mf::LogInfo("GeometryCore") << "New detector geometry loaded from "
                                 << "\n\t" << fROOTfile 
                                 << "\n\t" << fGDMLfile << "\n";
     
   } // GeometryCore::LoadGeometryFile()

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

private

Definition at line 1102 of file GeometryCore.cxx.

References AuxDets().

Referenced by FindAuxDet().

   {
     AuxDets().push_back(new AuxDetGeo(path, depth));
   }

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

private

Definition at line 1071 of file GeometryCore.cxx.

References Cryostats().

Referenced by FindCryostat().

   {
     Cryostats().emplace_back(path, depth);
   }

template<typename T >

geo::PlaneDataContainer<T> geo::GeometryCore::makePlaneData ( ) const

inline

Returns a container with one entry per wire plane.

Template Parameters

T	type of data in the container

Returns: a container with one default-constructed T per plane

See also: geo::PlaneDataContainer

The working assumption is that all cryostats have the same number of TPCs, and all TPCs have the same number of planes. It is always guaranteed that all existing planes have an entry in the container, although if the previous working assumption is not satisfied there will be entries in the containers which are not associated to a valid plane.

The interface of the container is detailed in the documentation of the container itself, geo::PlaneDataContainer. Example of usage:

auto const* geom = lar::providerFrom<geo::GeometryCore>();
auto hitsPerPlane
  = geom->makePlaneData<std::vector<recob::Hit const*>>();
for (recob::Hit const& hit: hits) {
  if (hit.WireID()) hitsPerPlane[hit.WireID()].push_back(&hit);
} // for

where the container will be filled with pointers to all hits on the given wire plane (wire IDs are implicitly converted into plane IDs in the index operator[] call).

Definition at line 2833 of file GeometryCore.h.

2834 { return { Ncryostats(), MaxTPCs(), MaxPlanes() }; }

geo::GeometryCore::MaxPlanes

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int MaxPlanes() const

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

Definition: GeometryCore.cxx:837

unsigned int MaxTPCs() const

Returns the largest number of TPCs a cryostat in the detector has.

Definition: GeometryCore.cxx:817

template<typename T >

geo::PlaneDataContainer<T> geo::GeometryCore::makePlaneData ( T const & defValue ) const

inline

Returns a container with one entry per wire plane.

Template Parameters

T	type of data in the container

Parameters

defValue the initial value of all elements in the container

Returns: a container with one default-constructed T per plane

See also: geo::PlaneDataContainer

This function operates as makePlaneData() const, except that copies the specified value into all the entries of the container. Example:

auto const* geom = lar::providerFrom<geo::GeometryCore>();
auto nHitsPerPlane = geom->makePlaneData(0U);
for (recob::Hit const& hit: hits) {
  if (hit.WireID()) ++(hitsPerPlane[hit.WireID()]);
} // for

Definition at line 2855 of file GeometryCore.h.

2856 { return { Ncryostats(), MaxTPCs(), MaxPlanes(), defValue }; }

geo::GeometryCore::MaxPlanes

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int MaxPlanes() const

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

Definition: GeometryCore.cxx:837

unsigned int MaxTPCs() const

Returns the largest number of TPCs a cryostat in the detector has.

Definition: GeometryCore.cxx:817

template<typename T >

geo::TPCDataContainer<T> geo::GeometryCore::makeTPCData ( ) const

inline

Returns a container with one entry per TPC.

Template Parameters

T	type of data in the container

Returns: a container with one default-constructed T per TPC

See also: geo::TPCDataContainer

The working assumption is that all cryostats have the same number of TPCs. It is always guaranteed that all existing TPCs have an entry in the container, although if the previous working assumption is not satisfied there will be entries in the containers which are not associated to a valid TPC.

The interface of the container is detailed in the documentation of the container itself, geo::TPCDataContainer. Example of usage:

auto const* geom = lar::providerFrom<geo::GeometryCore>();
auto tracksPerTPC
  = geom->makeTPCData<std::vector<recob::Track const*>>();
for (recob::Track const& track: tracks) {
  geo::TPCGeo const* tpc = geom->PositionToTPCptr(track.Start());
  if (tpc) tracksPerTPC[tpc->ID()].push_back(&track);
} // for

where the container will be filled with pointers to all tracks starting from a given TPC (tracks reconstructed as starting outside the TPCs will be not saved in the container).

Definition at line 2312 of file GeometryCore.h.

2313 { return { Ncryostats(), MaxTPCs() }; }

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int MaxTPCs() const

Returns the largest number of TPCs a cryostat in the detector has.

Definition: GeometryCore.cxx:817

template<typename T >

geo::TPCDataContainer<T> geo::GeometryCore::makeTPCData ( T const & defValue ) const

inline

Returns a container with one entry per TPC.

Template Parameters

T	type of data in the container

Parameters

defValue the initial value of all elements in the container

Returns: a container with a value defValue per each TPC

See also: geo::TPCDataContainer

This function operates as makeTPCData() const, except that copies the specified value into all the entries of the container. Example:

auto const* geom = lar::providerFrom<geo::GeometryCore>();
auto nTracksPerTPC = geom->makeTPCData(0U);
for (recob::Track const& track: tracks) {
  geo::TPCGeo const* tpc = geom->PositionToTPCptr(track.Start());
  if (tpc) ++(tracksPerTPC[tpc->ID()]);
} // for

Definition at line 2335 of file GeometryCore.h.

2336 { return { Ncryostats(), MaxTPCs(), defValue }; }

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int MaxTPCs() const

Returns the largest number of TPCs a cryostat in the detector has.

Definition: GeometryCore.cxx:817

double geo::GeometryCore::MassBetweenPoints	(	geo::Point_t const &	p1,
		geo::Point_t const &	p2
	)		const

Returns the column density between two points.

Parameters

p1	the first point
p2	the second point

Returns: the column density [kg / cm²]

The column density is defined as $\int_{\vec{p}_{1}}^{\vec{p}_{2}} \rho(\vec{p}) d\vec{p}$ where $\rho(\vec{p})$ is the density at point $\vec{p}$ , which the integral leads from p1 to p2 in a straight line.

Both points are specified in world coordinates.

Definition at line 1126 of file GeometryCore.cxx.

References dir, geo::vect::makePointFromCoords(), and R.

Referenced by TotalMass().

   {
 
     //The purpose of this method is to determine the column density
     //between the two points given.  Do that by starting at p1 and 
     //stepping until you get to the node of p2.  calculate the distance
     //between the point just inside that node and p2 to get the last
     //bit of column density
     double columnD = 0.;
 
     //first initialize a track - get the direction cosines
     geo::Vector_t const dir = (p2 - p1).Unit();
     
     double const dxyz[3] = { dir.X(), dir.Y(), dir.Z() };
     double const cp1[3] = { p1.X(), p1.Y(), p1.Z() };
     gGeoManager->InitTrack(cp1, dxyz);
 
     //might be helpful to have a point to a TGeoNode
     TGeoNode *node = gGeoManager->GetCurrentNode();
 
     //check that the points are not in the same volume already.  
     //if they are in different volumes, keep stepping until you 
     //are in the same volume as the second point
     while(!gGeoManager->IsSameLocation(p2.X(), p2.Y(), p2.Z())){
       gGeoManager->FindNextBoundary();
       columnD += gGeoManager->GetStep()*node->GetMedium()->GetMaterial()->GetDensity();
     
       //the act of stepping puts you in the next node and returns that node
       node = gGeoManager->Step();
     }//end loop to get to volume of second point
 
     //now you are in the same volume as the last point, but not at that point.
     //get the distance between the current point and the last one
     geo::Point_t const last
       = geo::vect::makePointFromCoords(gGeoManager->GetCurrentPoint());
     double const lastStep = (p2 - last).R();
     columnD += lastStep*node->GetMedium()->GetMaterial()->GetDensity();
 
     return columnD;
   }

double geo::GeometryCore::MassBetweenPoints	(	double *	p1,
		double *	p2
	)		const

Returns the column density between two points.

Parameters

p1	the first point
p2	the second point

Returns: the column density [kg / cm²]

The column density is defined as $\int_{\vec{p}_{1}}^{\vec{p}_{2}} \rho(\vec{p}) d\vec{p}$ where $\rho(\vec{p})$ is the density at point $\vec{p}$ , which the integral leads from p1 to p2 in a straight line.

Both points are specified in world coordinates.

TGeoMaterial const * geo::GeometryCore::Material ( geo::Point_t const & point ) const

Returns the material at the specified position.

Definition at line 925 of file GeometryCore.cxx.

Referenced by MaterialName().

                                                                           {
     auto const pNode = gGeoManager->FindNode(point.X(), point.Y(), point.Z());
     if (!pNode) return nullptr;
     auto const pMedium = pNode->GetMedium();
     return pMedium? pMedium->GetMaterial(): nullptr;
   }

std::string geo::GeometryCore::MaterialName ( TVector3 const & point ) const

inline

Name of the deepest material containing the point xyz.

Returns: material of the origin by default

Definition at line 1857 of file GeometryCore.h.

References art::detail::indent(), pyG4Element::Print(), and geo::vect::toPoint().

Referenced by larg4::LArG4Ana::analyze().

1858 { return MaterialName(geo::vect::toPoint(point)); }

geo::GeometryCore::MaterialName

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

std::string MaterialName(TVector3 const &point) const

Name of the deepest material containing the point xyz.

Definition: GeometryCore.h:1857

std::string geo::GeometryCore::MaterialName ( geo::Point_t const & point ) const

Name of the deepest material containing the point xyz.

Returns: material of the origin by default

Definition at line 933 of file GeometryCore.cxx.

References geo::BoxBoundedGeo::ContainsPosition(), FindDetectorEnclosure(), Material(), geo::BoxBoundedGeo::Max(), geo::BoxBoundedGeo::Min(), and WorldBox().

   {
     // check that the given point is in the World volume at least
     geo::BoxBoundedGeo worldBox = WorldBox();
     if (!worldBox.ContainsPosition(point)) {
       mf::LogWarning("GeometryCoreBadInputPoint")
         << "point " << point << " is not inside the world volume "
         << worldBox.Min() << " -- " << worldBox.Max()
         << "; returning unknown material name";
       return { "unknownMaterial" };
     }
     auto const pMaterial = Material(point);
     if (!pMaterial) {
       mf::LogWarning("GeometryCoreBadInputPoint")
         << "material for point " << point
         << " not found! returning unknown material name";
       return { "unknownMaterial" };
     }
     return pMaterial->GetName();
   }

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

Largest optical channel number.

Definition at line 218 of file GeometryCore.cxx.

References fChannelMapAlg, and NOpDets().

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

   {
     return fChannelMapAlg->MaxOpChannel(this->NOpDets());
   }

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

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

Definition at line 837 of file GeometryCore.cxx.

References Cryostats().

Referenced by pma::PMAlgTracker::build(), shower::EMShowerAlg::CheckIsolatedHits(), img::DataProviderAlg::DataProviderAlg(), trkf::TrackKalmanFitter::doFitWork(), lar_pandora::LArPandoraGeometry::LoadGeometry(), shower::EMShowerAlg::MakeShower(), shower::EMShowerAlg::MakeSpacePoints(), Nviews(), shower::EMShowerAlg::OrderShowerHits(), pma::PMAlgTracker::PMAlgTracker(), trkf::PMAlgTrackMaker::PMAlgTrackMaker(), lar_pandora::LArPandoraTrackCreation::produce(), shower::EMShowerAlg::RelativeWireWidth(), trkf::TrackKalmanFitter::sortOutput(), and shower::EMShowerAlg::WorstPlane().

                                              {
     unsigned int maxPlanes = 0;
     for (geo::CryostatGeo const& cryo: Cryostats()) {
       unsigned int maxPlanesInCryo = cryo.MaxPlanes();
       if (maxPlanesInCryo > maxPlanes) maxPlanes = maxPlanesInCryo;
     } // for
     return maxPlanes;
   } // GeometryCore::MaxPlanes()

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

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

Definition at line 1814 of file GeometryCore.cxx.

References fChannelMapAlg.

                                            {
     return fChannelMapAlg->MaxROPs();
   } // GeometryCore::MaxROPs()

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

Returns the largest number of TPCs a cryostat in the detector has.

Definition at line 817 of file GeometryCore.cxx.

References Cryostats().

Referenced by evd::TWQProjectionView::SetUpTPCselection().

                                            {
     unsigned int maxTPCs = 0;
     for (geo::CryostatGeo const& cryo: Cryostats()) {
       unsigned int maxTPCsInCryo = cryo.NTPC();
       if (maxTPCsInCryo > maxTPCs) maxTPCs = maxTPCsInCryo;
     } // for
     return maxTPCs;
   } // GeometryCore::MaxTPCs()

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

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

Definition at line 1770 of file GeometryCore.cxx.

References fChannelMapAlg.

                                               {
     return fChannelMapAlg->MaxTPCsets();
   } // GeometryCore::MaxTPCsets()

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

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

Definition at line 847 of file GeometryCore.cxx.

References Cryostats().

Referenced by cluster::BlurredClusteringAlg::ConvertRecobHitsToVector().

                                             {
     unsigned int maxWires = 0;
     for (geo::CryostatGeo const& cryo: Cryostats()) {
       unsigned int maxWiresInCryo = cryo.MaxWires();
       if (maxWiresInCryo > maxWires) maxWires = maxWiresInCryo;
     } // for
     return maxWires;
   } // GeometryCore::MaxWires()

unsigned int geo::GeometryCore::NAuxDets ( ) const

inline

Returns the number of auxiliary detectors.

Todo:: use a AutDetID_t instead of unsigned int?

This method returns the total number of scintillator paddles (Auxiliary Detectors aka AuxDet) outside of the cryostat

Todo:: Change return type to size_t

Definition at line 4458 of file GeometryCore.h.

Referenced by AuxDet(), NAuxDetSensitive(), and larg4::LArG4::produce().

4458 { return AuxDets().size(); }

geo::GeometryCore::AuxDets

AuxDetList_t & AuxDets()

Return the interfal auxiliary detectors list.

Definition: GeometryCore.h:5346

unsigned int geo::GeometryCore::NAuxDetSensitive ( size_t const & aid ) const

Returns the number of sensitive components of auxiliary detector.

Parameters

aid	ID of the auxiliary detector

Returns: number of sensitive components in the auxiliary detector aid cet::exception (category "Geometry") if aid does not exist

Definition at line 255 of file GeometryCore.cxx.

References AuxDets(), and NAuxDets().

   {
     if( aid > NAuxDets() - 1)
       throw cet::exception("Geometry") << "Requested AuxDet index " << aid 
                                        << " is out of range: " << NAuxDets();
 
     return AuxDets()[aid]->NSensitiveVolume();
   }

unsigned int geo::GeometryCore::Nchannels ( ) const

Returns the number of TPC readout channels in the detector.

Definition at line 191 of file GeometryCore.cxx.

References fChannelMapAlg.

Referenced by larg4::LArG4Ana::beginJob(), lar_cluster3d::StandardHit3DBuilder::BuildChannelStatusVec(), trkf::SeedFinderAlgorithm::CalculateGeometricalElements(), trkf::BezierTrackerAlgorithm::MakeTracks(), hit::MCHitFinder::produce(), and detsim::SimWire::produce().

   {
     return fChannelMapAlg->Nchannels();
   }

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

Returns the number of channels in the specified ROP.

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

Definition at line 197 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->Nchannels(ropid);
   } // GeometryCore::Nchannels(ROPID)

unsigned int geo::GeometryCore::Ncryostats ( ) const

inline

Returns the number of cryostats in the detector.

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 1964 of file GeometryCore.h.

Referenced by detsim::WienerFilterAna::beginJob(), evgen::GaisserParam::beginJob(), detsim::SimDriftElectrons::beginJob(), trkf::SpacePointAna::bookHistograms(), lar_cluster3d::StandardHit3DBuilder::BuildHitPairMap(), util::DetectorPropertiesServiceArgoNeuT::CalculateXTicksParams(), detinfo::DetectorPropertiesStandard::CalculateXTicksParams(), larg4::LArVoxelReadout::ClearSimChannels(), lar_cluster3d::StandardHit3DBuilder::CollectArtHits(), larg4::LArVoxelReadoutGeometry::Construct(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), hit::DisambigCheater::DisambigCheater(), detsim::WienerFilterAna::endJob(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), vertex::FeatureVertexFinder::Find3dVtxFrom2dClusterVtxCand(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), corner::CornerFinderAlg::get_feature_points_fast(), lar_pandora::LArPandoraInput::GetTrueStartAndEndPoints(), lar_cluster3d::SpacePointHit3DBuilder::Hit3DBuilder(), lar_pandora::LArPandoraGeometry::LoadGeometry(), trkf::SpacePointAlg::makeSpacePoints(), NOpDets(), OpDetFromCryo(), OpDetGeoFromOpDet(), evgen::CORSIKAGen::populateNShowers(), lar_pandora::LArPandoraTrackCreation::produce(), evgen::CosmicsGen::produce(), evgen::CORSIKAGen::produce(), vertex::FeatureVertexFinder::produce(), larg4::LArG4::produce(), evd::TWQProjectionView::SelectTPC(), evd::TWQProjectionView::SetUpTPCselection(), evd::SimulationDrawer::SimulationDrawer(), cheat::BackTracker::SpacePointHitsToWeightedXYZ(), trkf::SpacePointAlg::update(), and UpdateAfterSorting().

1964 { return Cryostats().size(); }

geo::GeometryCore::Cryostats

CryostatList_t & Cryostats()

Return the internal cryostat list.

Definition: GeometryCore.h:5340

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

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 1271 of file GeometryCore.cxx.

References raw::InvalidChannelID, NearestWireID(), and PlaneWireToChannel().

Referenced by larg4::LArVoxelReadout::DriftIonizationElectrons(), NearestChannel(), NearestWireID(), detsim::SimDriftElectrons::produce(), vertex::FeatureVertexFinder::produce(), 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.
     geo::WireID const wireID = NearestWireID(worldPos, planeid);
     return wireID? PlaneWireToChannel(wireID): raw::InvalidChannelID;
   } // GeometryCore::NearestChannel()

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	const double	worldLoc[3],
		geo::PlaneID const &	planeid
	)		const

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 1253 of file GeometryCore.cxx.

References NearestChannel().

   {
     return NearestChannel(TVector3(worldPos), planeid);
   }

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	std::vector< double > const &	worldLoc,
		geo::PlaneID const &	planeid
	)		const

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 1260 of file GeometryCore.cxx.

References NearestChannel().

   {
     if(worldPos.size() > 3) throw cet::exception("GeometryCore") << "bad size vector for "
                                                              << "worldPos: " 
                                                              << worldPos.size() << "\n";
     TVector3 wp(&(worldPos[0]));
     return NearestChannel(wp, planeid);
   }

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	const TVector3 &	worldLoc,
		geo::PlaneID const &	planeid
	)		const

inline

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 4714 of file GeometryCore.h.

References geo::vect::toPoint().

4715 { return NearestChannel(geo::vect::toPoint(worldLoc), planeid); }

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

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

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

Definition: GeometryCore.cxx:1271

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	const double	worldLoc[3],
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 4716 of file GeometryCore.h.

4720 { return NearestChannel(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

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

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

Definition: GeometryCore.cxx:1271

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	std::vector< double > const &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 4721 of file GeometryCore.h.

4725 { return NearestChannel(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

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

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

Definition: GeometryCore.cxx:1271

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	const TVector3 &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 4726 of file GeometryCore.h.

4730 { return NearestChannel(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

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

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

Definition: GeometryCore.cxx:1271

raw::ChannelID_t geo::GeometryCore::NearestChannel	(	geo::Point_t const &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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

Parameters

worldLoc	3D coordinates of the point (world reference frame)
PlaneNo	the number of plane
TPCNo	the number of TPC
cstat	the number of cryostat

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

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

The different versions allow different way to provide the position.

Todo:: remove the integers version

Definition at line 4731 of file GeometryCore.h.

References c1, and c2.

4735 { return NearestChannel(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

Referenced by vertex::FeatureVertexFinderAna::analyze(), trkf::TrackAna::analyze(), evd::RecoBaseDrawer::DrawTrack2D(), evd::TWQMultiTPCProjectionView::FindEndPoint(), evd::TWQProjectionView::FindEndPoint(), hit::HitCheater::FindHitsOnChannel(), cmtool::CFAlgoStartPointCompat::Float(), cmtool::CFAlgoStartTimeCompat::Float(), util::GeometryUtilities::Get2DPointProjection(), util::GeometryUtilities::Get2DPointProjectionCM(), trkf::FeatureTracker::GetConnectionMap(), opdet::GetHitGeometryInfo(), trkf::BezierTrack::GetProjectedPointUVWT(), trkf::BezierTrack::GetProjectedPointUVWX(), nnet::TrainingDataAlg::getProjection(), apa::APAGeometryAlg::LineSegChanIntersect(), NearestWire(), cluster::ClusterCheater::produce(), shwf::ShowerReco::produce(), vertex::FeatureVertexFinder::produce(), evd::RecoBaseDrawer::Seed2D(), and WireCoordinate().

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

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

Definition: GeometryCore.cxx:1271

unsigned int geo::GeometryCore::NearestWire	(	geo::Point_t const &	point,
		geo::PlaneID const &	planeid
	)		const

Returns the index of wire closest to position in the specified TPC.

Parameters

point	the point to be tested [cm]
planeid	ID of the plane

Returns: the index of the wire, or geo::WireID::InvalidID on failure

Bug:: Actually, on failure an exception geo::InvalidWireError is thrown

Deprecated:: Use NearestWireID() instead.

Definition at line 1202 of file GeometryCore.cxx.

References NearestWireID(), and geo::WireID::Wire.

   {
     return NearestWireID(point, planeid).Wire;
     // return fChannelMapAlg->NearestWire(worldPos, planeid);
   }

unsigned int geo::GeometryCore::NearestWire	(	const double	worldLoc[3],
		geo::PlaneID const &	planeid
	)		const

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 1210 of file GeometryCore.cxx.

References NearestWire().

   {
     return NearestWire(TVector3(worldPos), planeid);
   }

unsigned int geo::GeometryCore::NearestWire	(	std::vector< double > const &	worldLoc,
		geo::PlaneID const &	planeid
	)		const

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 1217 of file GeometryCore.cxx.

References NearestWire(), and NearestWireID().

   {
     if(worldPos.size() > 3) throw cet::exception("GeometryCore") << "bad size vector for "
                                                              << "worldPos: " 
                                                              << worldPos.size() << "\n";
     TVector3 wp(&(worldPos[0]));
     return NearestWire(wp, planeid);
   }

unsigned int geo::GeometryCore::NearestWire	(	const TVector3 &	worldLoc,
		geo::PlaneID const &	planeid
	)		const

inline

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 3951 of file GeometryCore.h.

References geo::vect::toPoint().

3952 { return NearestWire(geo::vect::toPoint(worldLoc), planeid); }

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

geo::WireID::WireID_t NearestWire(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the index of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1202

unsigned int geo::GeometryCore::NearestWire	(	const double	worldLoc[3],
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 3953 of file GeometryCore.h.

3957 { return NearestWire(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID::WireID_t NearestWire(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the index of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1202

unsigned int geo::GeometryCore::NearestWire	(	std::vector< double > const &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 3958 of file GeometryCore.h.

3962 { return NearestWire(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID::WireID_t NearestWire(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the index of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1202

unsigned int geo::GeometryCore::NearestWire	(	const TVector3 &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 3963 of file GeometryCore.h.

3967 { return NearestWire(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID::WireID_t NearestWire(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the index of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1202

unsigned int geo::GeometryCore::NearestWire	(	geo::Point_t const &	worldLoc,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the index of wire closest to position in the specified TPC.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the index of the wire

The different versions allow different way to provide the position.

Deprecated:: Use NearestWireID() instead.

Todo:

remove the integers version

what happens when no wire is found?

Definition at line 3968 of file GeometryCore.h.

3972 { return NearestWire(worldLoc, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID::WireID_t NearestWire(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the index of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1202

geo::WireID geo::GeometryCore::NearestWireID	(	geo::Point_t const &	point,
		geo::PlaneID const &	planeid
	)		const

Returns the ID of wire closest to position in the specified TPC.

Parameters

point	the point to be tested [cm]
planeid	ID of the plane

Returns: the ID of the wire, or an invalid wire ID

See also: geo::PlaneGeo::ClosestWireID()

Bug:: Instead of returning an invalid wire ID, an exception is thrown!

If the nearest wire is not closer than half a wire pitch, the result is marked invalid. The returned (invalid) ID will contain the non-existing wire that would be the nearest, if it existed.

If the wire ID is invalid and the existing closest wire is desired, a possible solution is (when the BUG will be solved):

geo::WireID wireID = geom->NearestWireID(point, planeID);

if (!wireID) wireID = geom->Plane(planeID).ClosestWireID(wireID);

Note however that this will execute plane lookup twice, and a more efficient approach would be to ask the plane everything directly:

geo::PlaneGeo const& plane = geom->Plane(planeID);
geo::WireID wireID = plane.NearestWireID(point);
if (!wireID) wireID = plane.ClosestWireID(wireID);

Until the BUG is fixed, the actual working code is:

geo::PlaneGeo const& plane = geom->Plane(planeID);
geo::WireID wireID;
try {
  wireID = plane.NearestWireID(point);
}
catch (geo::InvalidWireError const& e) {
  if (!e.hasSuggestedWire()) throw;
  wireID = plane.ClosestWireID(e.suggestedWireID());
}

Definition at line 1228 of file GeometryCore.cxx.

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

Referenced by hit::DisambigCheater::InitHitToWids(), NearestChannel(), NearestWire(), NearestWireID(), apa::APAGeometryAlg::NearestWireIDOnChan(), evd::RecoBaseDrawer::OpFlash2D(), trk::TrackContainmentAlg::ProcessTracks(), and shower::EMShowerAlg::Project3DPointOntoPlane().

   {
     return Plane(planeid).NearestWireID(worldPos);
   }

geo::WireID geo::GeometryCore::NearestWireID	(	const double	point[3],
		geo::PlaneID const &	planeid
	)		const

Returns the ID of wire closest to position in the specified TPC.

Parameters

point	the point to be tested [cm]
planeid	ID of the plane
PlaneNo	plane number within the TPC
TPCNo	tpc number within the cryostat
cstat	cryostat number

Returns: the ID of the wire, or an invalid wire ID

Bug:: Instead of returning an invalid wire ID, an exception is thrown!

The different versions allow different way to provide the position.

Deprecated:: Use the version with a geo::Point_t and PlaneID arguments

Todo:: remove the integers version

Definition at line 1246 of file GeometryCore.cxx.

References NearestChannel(), and NearestWireID().

   {
     return NearestWireID(TVector3(worldPos), planeid);
   }

geo::WireID geo::GeometryCore::NearestWireID	(	std::vector< double > const &	point,
		geo::PlaneID const &	planeid
	)		const

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 1235 of file GeometryCore.cxx.

References NearestWireID().

   {
     if(worldPos.size() > 3) throw cet::exception("GeometryCore") << "bad size vector for "
                                                              << "worldPos: " 
                                                              << worldPos.size() << "\n";
     TVector3 wp(&(worldPos[0]));
     return NearestWireID(wp, planeid);
   }

geo::WireID geo::GeometryCore::NearestWireID	(	const TVector3 &	point,
		geo::PlaneID const &	planeid
	)		const

inline

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 3894 of file GeometryCore.h.

References geo::vect::toPoint().

3895 { return NearestWireID(geo::vect::toPoint(point), planeid); }

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

geo::WireID NearestWireID(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the ID of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1228

geo::WireID geo::GeometryCore::NearestWireID	(	const double	point[3],
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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 3896 of file GeometryCore.h.

3900 { return NearestWireID(point, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID NearestWireID(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the ID of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1228

geo::WireID geo::GeometryCore::NearestWireID	(	std::vector< double > const &	point,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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 3901 of file GeometryCore.h.

3905 { return NearestWireID(point, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID NearestWireID(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the ID of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1228

geo::WireID geo::GeometryCore::NearestWireID	(	const TVector3 &	point,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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 3906 of file GeometryCore.h.

3910 { return NearestWireID(point, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID NearestWireID(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the ID of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1228

geo::WireID geo::GeometryCore::NearestWireID	(	geo::Point_t const &	point,
		unsigned int const	PlaneNo,
		unsigned int const	TPCNo = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

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 3911 of file GeometryCore.h.

3915 { return NearestWireID(point, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::WireID NearestWireID(geo::Point_t const &point, geo::PlaneID const &planeid) const

Returns the ID of wire closest to position in the specified TPC.

Definition: GeometryCore.cxx:1228

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

inline

Returns the number of cryostats in the detector.

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 1965 of file GeometryCore.h.

1965 { return Ncryostats(); }

geo::GeometryCore::Nplanes

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int geo::GeometryCore::NElements ( geo::CryostatID const & cryoid ) const

inline

Returns the total number of TPCs in the specified cryostat.

Parameters

cstat cryostat number

Todo:

Make the cryostat number mandatory (as CryostatID)

Change return type to size_t

what happens if it does not exist?

Definition at line 2356 of file GeometryCore.h.

2357 { return NTPC(cryoid); }

geo::GeometryCore::NTPC

unsigned int NTPC(unsigned int cstat=0) const

Returns the total number of TPCs in the specified cryostat.

Definition: GeometryCore.h:2273

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

inline

Returns the total number of wire planes in the specified TPC.

Parameters

tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as TPCID)

Change return type to size_t

what happens if TPC does not exist?

Definition at line 2875 of file GeometryCore.h.

2876 { return Nplanes(tpcid); }

unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wire planes in the specified TPC.

Definition: GeometryCore.h:2799

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

inline

Returns the total number of wires in the specified plane.

Parameters

p	plane number within the TPC
tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as PlaneID)

Change return type to size_t

what happens if it does not exist?

Definition at line 3312 of file GeometryCore.h.

3313 { return Nwires(planeid); }

geo::GeometryCore::Nwires

unsigned int Nwires(unsigned int p, unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wires in the specified plane.

Definition: GeometryCore.h:3293

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

Number of electronics channels for all the optical detectors.

Todo:: add explanation of the different IDs

Definition at line 212 of file GeometryCore.cxx.

References fChannelMapAlg, and NOpDets().

Referenced by opdet::OpHitAna::analyze(), opdet::OpFlashAna::analyze(), opdet::SimPhotonCounter::analyze(), opdet::SimPhotonCounter::beginJob(), opdet::OpDetResponseInterface::doNOpChannels(), opdet::OpDigiProperties::FillGainArray(), opdet::OpDigiProperties::FillPedMeanArray(), opdet::FlashClusterMatch::produce(), opdet::TrackTimeAssoc::produce(), opdet::BeamFlashCompatabilityCheck::produce(), opdet::OptDetDigitizer::produce(), and opdet::OpDigiProperties::reconfigure().

   {
     return fChannelMapAlg->NOpChannels(this->NOpDets());
   }

unsigned int geo::GeometryCore::NOpDets ( ) const

Number of OpDets in the whole detector.

Definition at line 203 of file GeometryCore.cxx.

References Cryostat(), Ncryostats(), and geo::CryostatGeo::NOpDet().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), phot::CreateHybridLibrary::CreateHybridLibrary(), opdet::FlashHypothesisAnaAlg::FillOpDetPositions(), larg4::OpDetLookup::FindClosestOpDet(), opdet::FlashClusterMatch::GetLightHypothesis(), opdet::TrackTimeAssoc::GetMIPHypotheses(), opdet::BeamFlashCompatabilityCheck::GetMIPHypotheses(), IsValidOpChannel(), phot::PhotonVisibilityService::LoadLibrary(), MaxOpChannel(), NOpChannels(), cosmic::BeamFlashTrackMatchTaggerAlg::NormalizeLightHypothesis(), phot::PhotonLibraryHybrid::PhotonLibraryHybrid(), cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison(), wc::CellTree::processOpFlash(), opdet::FlashClusterMatch::produce(), opdet::TrackTimeAssoc::produce(), opdet::BeamFlashCompatabilityCheck::produce(), larg4::LArG4::produce(), lariov::SIOVPmtGainProvider::Reconfigure(), opdet::TrackTimeAssoc::ScanMIPHypotheses(), and opdet::FlashHypothesisAnaAlg::SetOutputObjects().

   {
     int N=0;
     for(size_t cstat=0; cstat!=Ncryostats(); cstat++)
       N += this->Cryostat(cstat).NOpDet();
     return N;
   }

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

Number of electronics channels for all the optical detectors.

Todo:: add explanation of the different IDs

Definition at line 224 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->NOpHardwareChannels(opDet);
   }

unsigned int geo::GeometryCore::Nplanes	(	unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the total number of wire planes in the specified TPC.

Parameters

tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as TPCID)

Change return type to size_t

what happens if TPC does not exist?

Definition at line 2799 of file GeometryCore.h.

Referenced by hit::MCHitAnaExample::analyze(), cluster::ClusterAna::analyze(), pfpf::PFPAna::analyze(), hit::MagDriftAna::analyze(), cluster::DBclusterAna::analyze(), trkf::TrackAna::analyze(), tca::AnalyzeHits(), cluster::SmallClusterFilter::beginJob(), cluster::SmallClusterFinder::beginJob(), detsim::WienerFilterAna::beginJob(), shwf::ShowerReco::beginJob(), lar_cluster3d::StandardHit3DBuilder::BuildChannelStatusVec(), btutil::MCMatchAlg::BuildMap(), mvapid::MVAAlg::CalcSegmentdEdxDist(), util::GeometryUtilities::CalculatePitch(), util::GeometryUtilities::CalculatePitchPolar(), tca::TruthMatcher::CanReconstruct(), cluster::ClusterMatchTQ::ClusterMatch(), opdet::ConstructFlash(), evd::AnalysisBaseDrawer::DrawDeDx(), evd::AnalysisBaseDrawer::DrawKineticEnergy(), detsim::WienerFilterAna::endJob(), cluster::HoughBaseAlg::FastTransform(), hit::GausHitFinder::FillOutHitParameterVector(), hit::DPRawHitFinder::FillOutHitParameterVector(), trkf::BezierTrack::FillTrajectoryVectors(), tca::FillWireHitRange(), filter::MuonFilter::filter(), calo::GeneralCalorimetry::GeneralCalorimetry(), opdet::GetHitGeometryInfo(), mvapid::MVAAlg::GetWireNormals(), corner::CornerFinderAlg::GrabWires(), evd::GraphCluster::GraphCluster(), evd::GraphClusterAlg::GraphClusterAlg(), evd::InfoTransfer::InfoTransfer(), corner::CornerFinderAlg::InitializeGeometry(), tca::TruthMatcher::MatchAndSum(), hit::MCHitAnaExample::MCHitAnaExample(), cluster::ClusterMatchAlg::PrepareDetParams(), btutil::MCBTAlg::ProcessSimChannel(), trk::TrackContainmentAlg::ProcessTracks(), calo::Calorimetry::produce(), shwf::ShowerReco::produce(), evd::RawDataDrawer::RawDataDrawer(), evd::InfoTransfer::Rebuild(), evd::RecoBaseDrawer::RecoBaseDrawer(), evd::InfoTransfer::reconfigure(), util::GeometryUtilities::Reconfigure(), hit::RFFHitFinderAlg::SetFitterParamsVectors(), showerreco::ShowerRecoManager::ShowerRecoManager(), cheat::BackTracker::SpacePointHitsToWeightedXYZ(), ThirdPlane(), evd::TQPad::TQPad(), cluster::HoughBaseAlg::Transform(), evd::TWireProjPad::TWireProjPad(), evd::TWQMultiTPCProjectionView::TWQMultiTPCProjectionView(), and evd::TWQProjectionView::TWQProjectionView().

2800 { return Nplanes(geo::TPCID(cstat, tpc)); }

geo::GeometryCore::Nplanes

unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wire planes in the specified TPC.

Definition: GeometryCore.h:2799

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

unsigned int geo::GeometryCore::Nplanes ( geo::TPCID const & tpcid ) const

inline

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 2870 of file GeometryCore.h.

References geo::TPCGeo::NElements().

       {
         TPCGeo const* pTPC = GetElementPtr(tpcid);
         return pTPC? pTPC->NElements(): 0;
       }

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

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

Parameters

tpcsetid TPC set ID

Returns: number of readout planes in the TPC set, or 0 if no TPC set found

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

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

Definition at line 1808 of file GeometryCore.cxx.

References fChannelMapAlg.

                                                                         {
     return fChannelMapAlg->NROPs(tpcsetid);
   } // GeometryCore::NROPs()

unsigned int geo::GeometryCore::NSiblingElements ( geo::CryostatID const & ) const

inline

Returns the number of cryostats in the detector.

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 1966 of file GeometryCore.h.

1967 { return Ncryostats(); }

geo::GeometryCore::Nplanes

unsigned int Ncryostats() const

Returns the number of cryostats in the detector.

Definition: GeometryCore.h:1964

unsigned int geo::GeometryCore::NSiblingElements ( geo::TPCID const & tpcid ) const

inline

Returns the total number of TPCs in the specified cryostat.

Parameters

cstat cryostat number

Todo:

Make the cryostat number mandatory (as CryostatID)

Change return type to size_t

what happens if it does not exist?

Definition at line 2358 of file GeometryCore.h.

2359 { return NTPC(tpcid); }

geo::GeometryCore::NTPC

unsigned int NTPC(unsigned int cstat=0) const

Returns the total number of TPCs in the specified cryostat.

Definition: GeometryCore.h:2273

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

inline

Returns the total number of wire planes in the specified TPC.

Parameters

tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as TPCID)

Change return type to size_t

what happens if TPC does not exist?

Definition at line 2877 of file GeometryCore.h.

2878 { return Nplanes(planeid); }

unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wire planes in the specified TPC.

Definition: GeometryCore.h:2799

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

inline

Returns the total number of wires in the specified plane.

Parameters

p	plane number within the TPC
tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as PlaneID)

Change return type to size_t

what happens if it does not exist?

Definition at line 3314 of file GeometryCore.h.

3315 { return Nwires(wireid); }

geo::GeometryCore::Nwires

unsigned int Nwires(unsigned int p, unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wires in the specified plane.

Definition: GeometryCore.h:3293

unsigned int geo::GeometryCore::NSiblingElements ( readout::TPCsetID const & tpcsetid ) const

inline

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

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

Definition at line 4775 of file GeometryCore.h.

4776 { return NTPCsets(tpcsetid); }

geo::GeometryCore::NTPCsets

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

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

Definition: GeometryCore.cxx:1764

unsigned int geo::GeometryCore::NSiblingElements ( readout::ROPID const & ropid ) const

inline

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

Parameters

tpcsetid TPC set ID

Returns: number of readout planes in the TPC set, or 0 if no TPC set found

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

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

Definition at line 4978 of file GeometryCore.h.

4979 { return NROPs(ropid); }

geo::GeometryCore::NROPs

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

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

Definition: GeometryCore.cxx:1808

unsigned int geo::GeometryCore::NTPC ( unsigned int cstat = 0 ) const

inline

Returns the total number of TPCs in the specified cryostat.

Parameters

cstat cryostat number

Todo:

Make the cryostat number mandatory (as CryostatID)

Change return type to size_t

what happens if it does not exist?

Definition at line 2273 of file GeometryCore.h.

Referenced by detsim::WienerFilterAna::beginJob(), detsim::SimDriftElectrons::beginJob(), lar_cluster3d::StandardHit3DBuilder::BuildHitPairMap(), mvapid::MVAAlg::CalcSegmentdEdxDist(), larg4::LArVoxelReadout::ClearSimChannels(), lar_cluster3d::StandardHit3DBuilder::CollectArtHits(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), hit::DisambigCheater::DisambigCheater(), detsim::WienerFilterAna::endJob(), lar_pandora::LArPandoraInput::GetTrueStartAndEndPoints(), mvapid::MVAAlg::GetWireNormals(), lar_cluster3d::SpacePointHit3DBuilder::Hit3DBuilder(), apa::APAGeometryAlg::Init(), lar_pandora::LArPandoraGeometry::LoadGeometry(), evd::RecoBaseDrawer::OpFlashOrtho(), evd::Ortho3DPad::Ortho3DPad(), lar_pandora::LArPandoraTrackCreation::produce(), cluster::TrajCluster::produce(), evgen::NucleonDecay::produce(), evgen::NDKGen::produce(), evgen::NeutronOsc::produce(), evd::RecoBaseDrawer::RecoBaseDrawer(), evgen::ActiveVolumeVertexSampler::reconfigure(), evd::TWQProjectionView::SelectTPC(), evd::TWQProjectionView::SetUpTPCselection(), cheat::BackTracker::SpacePointHitsToWeightedXYZ(), tca::StitchPFPs(), apa::DisambigAlg::TrivialDisambig(), and evd::TWQMultiTPCProjectionView::TWQMultiTPCProjectionView().

2274 { return NTPC(geo::CryostatID(cstat)); }

geo::GeometryCore::NTPC

unsigned int NTPC(unsigned int cstat=0) const

Returns the total number of TPCs in the specified cryostat.

Definition: GeometryCore.h:2273

Referenced by lar_cluster3d::StandardHit3DBuilder::BuildChannelStatusVec(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), noisefilteralg::WireCellNoiseFilter::DoNoiseFilter(), evd_tool::ICARUSDrawer::DrawBadChannels(), evd_tool::MicroBooNEDrawer::DrawBadChannels(), sppt::SpacePointAlg_TimeSort::fillCoordinatesArrays(), tca::FillWireHitRange(), vertex::FeatureVertexFinder::Find2dClusterVertexCandidates(), evd::RawDataDrawer::RoIextractorClass::Finish(), cmtool::CFAlgoVolumeOverlap::Float(), cmtool::CFAlgoWireOverlap::Float(), evd::RawDataDrawer::GetRegionOfInterest(), evd::RecoBaseDrawer::GetRegionOfInterest(), util::GeometryUtilities::GetYZ(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire(), corner::CornerFinderAlg::GrabWires(), corner::CornerFinderAlg::InitializeGeometry(), lar_cluster3d::StandardHit3DBuilder::NearestWireID(), cluster::fuzzyClusterAlg::run_fuzzy_cluster(), util::GeometryUtilities::SelectPolygonHitList(), evd::TWQProjectionView::SetWire(), img::DataProviderAlg::setWireDriftData(), cluster::HoughBaseAlg::Transform(), evd::TWireProjPad::TWireProjPad(), evd::TWQMultiTPCProjectionView::TWQMultiTPCProjectionView(), evd::TWQProjectionView::TWQProjectionView(), evd::RecoBaseDrawer::Wire2D(), WireIDIntersectionCheck(), evd::TWQMultiTPCProjectionView::ZoomInterest(), and evd::TWQProjectionView::ZoomInterest().

The data type to uniquely identify a cryostat.

Definition: geo_types.h:120

unsigned int geo::GeometryCore::NTPC ( geo::CryostatID const & cryoid ) const

inline

Returns the total number of TPCs in the specified cryostat.

Parameters

cryoid cryostat number

Returns: number of TPCs in specified cryostat, or 0 if no cryostat 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 2351 of file GeometryCore.h.

References geo::CryostatGeo::NElements().

       {
         CryostatGeo const* pCryo = GetElementPtr(cryoid);
         return pCryo? pCryo->NElements(): 0;
       }

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

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

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

Definition at line 1764 of file GeometryCore.cxx.

References fChannelMapAlg.

                                                                            {
     return fChannelMapAlg->NTPCsets(cryoid);
   } // GeometryCore::NTPCsets()

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

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 266 of file GeometryCore.cxx.

References MaxPlanes().

Referenced by evd_tool::ICARUSDrawer::DrawBadChannels(), and evd_tool::MicroBooNEDrawer::DrawBadChannels().

   {
     return MaxPlanes();
   }

unsigned int geo::GeometryCore::Nwires	(	unsigned int	p,
		unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the total number of wires in the specified plane.

Parameters

p	plane number within the TPC
tpc	tpc number within the cryostat
cstat	cryostat number

Todo:

Make all the arguments mandatory (as PlaneID)

Change return type to size_t

what happens if it does not exist?

Definition at line 3293 of file GeometryCore.h.

3294 { return Nwires(geo::PlaneID(cstat, tpc, p)); }

geo::GeometryCore::Nwires

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

unsigned int Nwires(unsigned int p, unsigned int tpc=0, unsigned int cstat=0) const

Returns the total number of wires in the specified plane.

Definition: GeometryCore.h:3293

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

inline

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 3307 of file GeometryCore.h.

References geo::PlaneGeo::NElements().

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

unsigned int geo::GeometryCore::OpChannel	(	int	detNum,
		int	hardwareChannel
	)		const

Convert detector number and hardware channel to unique channel.

Definition at line 230 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->OpChannel(detNum, hardwareChannel);
   }

unsigned int geo::GeometryCore::OpDetFromCryo	(	unsigned int	o,
		unsigned int	c
	)		const

Get unique opdet number from cryo and internal count.

Definition at line 1880 of file GeometryCore.cxx.

References Cryostat(), Ncryostats(), and geo::CryostatGeo::NOpDet().

Referenced by GetClosestOpDet().

   {
     static bool Loaded=false;
     static std::vector<unsigned int> LowestID;
     static unsigned int NCryo;
     // If not yet loaded static parameters, do it
     if(Loaded == false){
         
       Loaded = true;
       
       // Store the lowest ID for each cryostat
       NCryo=Ncryostats();
       LowestID.resize(NCryo + 1);
       LowestID.at(0)=0;        
       for(size_t cryo=0; cryo!=NCryo; ++cryo){
         LowestID.at(cryo+1)=LowestID.at(cryo)+Cryostat(c).NOpDet();
       }
       
     }
 
     if( (c < NCryo) && (o < Cryostat(c).NOpDet())){
       return LowestID.at(c)+o;
     }
     else{
       throw cet::exception("OpDetCryoToOpID Error") << "Coordinates c=" << c 
                                                     << ", o=" << o 
                                                     << " out of range. Abort\n";
     }
     
     // if all is well, we never get to this point in the method
     // but still a good idea to be sure to always return something.
 
     return INT_MAX;  
   }

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

Convert unique channel to detector number.

Definition at line 236 of file GeometryCore.cxx.

References fChannelMapAlg.

Referenced by cheat::PhotonBackTracker::ChannelToTrackSDPs(), cosmic::BeamFlashTrackMatchTaggerAlg::CheckCompatibility(), opdet::OpDetResponseInterface::doReadoutToGeoChannel(), OpDetGeoFromOpChannel(), cheat::PhotonBackTracker::OpHitToSDPs(), cheat::PhotonBackTracker::OpHitToSimSDPs_Ps(), cheat::PhotonBackTracker::OpHitToTrackSDPs(), cosmic::BeamFlashTrackMatchTaggerAlg::PrintHypothesisFlashComparison(), opdet::FlashClusterMatch::produce(), opdet::TrackTimeAssoc::produce(), opdet::BeamFlashCompatabilityCheck::produce(), and cheat::PhotonBackTracker::TrackIdsToOpHits_Ps().

   {
     return fChannelMapAlg->OpDetFromOpChannel(opChannel);
   }

const OpDetGeo & geo::GeometryCore::OpDetGeoFromOpChannel ( unsigned int OpChannel ) const

Access the OpDetGeo object by OpDet or Channel Number.

Definition at line 1916 of file GeometryCore.cxx.

References OpDetFromOpChannel(), and OpDetGeoFromOpDet().

Referenced by opdet::GetHitGeometryInfo().

   {
     return this->OpDetGeoFromOpDet(this->OpDetFromOpChannel(OpChannel));
   }

const OpDetGeo & geo::GeometryCore::OpDetGeoFromOpDet ( unsigned int OpDet ) const

Number of OpDets in the whole detector.

Definition at line 1922 of file GeometryCore.cxx.

References Cryostat(), Ncryostats(), geo::CryostatGeo::NOpDet(), and geo::CryostatGeo::OpDet().

Referenced by phot::CreateHybridLibrary::CreateHybridLibrary(), phot::PhotonVisibilityService::DistanceToOpDet(), larg4::OpDetLookup::FindClosestOpDet(), phot::PhotonLibraryHybrid::GetCount(), OpDetGeoFromOpChannel(), larg4::OpFastScintillation::propagation_time(), and phot::PhotonVisibilityService::SolidAngleFactor().

   {
     static bool Loaded=false;
     static std::vector<unsigned int> LowestID;
     static size_t NCryo;
     // If not yet loaded static parameters, do it
     if(Loaded == false){
       
       Loaded = true;
 
       // Store the lowest ID for each cryostat
       NCryo=Ncryostats();
       LowestID.resize(NCryo + 1);
       LowestID[0] = 0;        
       for(size_t cryo = 0; cryo != NCryo; ++cryo){
         LowestID[cryo+1] = LowestID[cryo] + Cryostat(cryo).NOpDet();
       }
         
     }
 
     for(size_t i=0; i!=NCryo; ++i){
       if( (OpDet >= LowestID[i]) && (OpDet < LowestID[i+1]) ){
         int c = i;
         int o = OpDet-LowestID[i]; 
         return this->Cryostat(c).OpDet(o);
       }
     }
     // If we made it here, we didn't find the right combination. abort
     throw cet::exception("OpID To OpDetCryo error")<<"OpID out of range, "<< OpDet << "\n";
 
     // Will not reach due to exception
     return this->Cryostat(0).OpDet(0);
   }

std::string geo::GeometryCore::OpDetGeoName ( unsigned int c = 0 ) const

Returns gdml string which gives sensitive opdet name.

Parameters

c	ID of the cryostat the detector is in

This name is defined in the geometry (GDML) description.

Todo:: Change to use CryostatID

Definition at line 1873 of file GeometryCore.cxx.

References Cryostat(), and geo::CryostatGeo::OpDetGeoName().

Referenced by larg4::LArG4::beginJob().

   {
     return Cryostat(c).OpDetGeoName();
   }

GeometryCore& geo::GeometryCore::operator= ( GeometryCore const & )

delete

GeometryCore& geo::GeometryCore::operator= ( GeometryCore && )

delete

PlaneGeo const& geo::GeometryCore::Plane	(	unsigned int const	p,
		unsigned int const	tpc = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the specified wire.

Parameters

planeid	ID of the plane
p	plane number within the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

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.

Todo:: remove the version with integers

Definition at line 2945 of file GeometryCore.h.

2947 { return Plane(geo::PlaneID(cstat, tpc, p)); }

geo::GeometryCore::Plane

PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified wire.

Definition: GeometryCore.h:2945

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

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

inline

Returns the specified wire.

Parameters

planeid	ID of the plane
p	plane number within the TPC
tpc	TPC number within the cryostat
cstat	number of cryostat

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.

Todo:: remove the version with integers

Definition at line 2948 of file GeometryCore.h.

2949 { return TPC(planeid).Plane(planeid); }

Referenced by hit::MCHitAnaExample::analyze(), cluster::BlurredClusteringAlg::ConvertRecobHitsToVector(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), evd_tool::MicroBooNEDrawer::DrawBadChannels(), evd_tool::ICARUSDrawer::DrawBadChannels(), evd::RawDataDrawer::FillTQHisto(), tca::FillWireHitRange(), vertex::FeatureVertexFinder::Find3dVtxFrom2dClusterVtxCand(), evd::TWQMultiTPCProjectionView::FindEndPoint(), evd::TWQMultiTPCProjectionView::FindLineLength(), cmtool::CFAlgoStartTimeCompat::Float(), vertex::FeatureVertexFinder::Get3dVertexCandidates(), tca::GetOrigin(), util::GeometryUtilities::GetProjectedPoint(), util::GeometryUtilities::GetYZ(), cluster::DBScan3DAlg::init(), apa::DisambigAlg::MakeCloseHits(), NearestChannel(), wc::MergeWireCell::produce(), calo::Calorimetry::produce(), vertex::VertexFinder2D::produce(), shwf::ShowerReco::produce(), lariov::SIOVElectronicsCalibProvider::Reconfigure(), lariov::DetPedestalRetrievalAlg::Reconfigure(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), pma::ProjectionMatchingAlg::validate_on_adc_test(), and evd::RecoBaseDrawer::Wire2D().

TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified TPC.

Definition: GeometryCore.h:2393

geo::TPCGeo::Plane

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

Return the plane in the tpc with View_t view.

Definition: TPCGeo.cxx:299

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

Returns a list of possible PlaneIDs in the detector.

Returns: a constant reference to the set of plane IDs

Deprecated:: use IteratePlaneIDs() instead plane IDs of DUNE FD? probably better to use iterators instead

Definition at line 572 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->PlaneIDs();
   }

double geo::GeometryCore::PlanePitch	(	geo::TPCID const &	tpcid,
		geo::PlaneID::PlaneID_t	p1 = `0`,
		geo::PlaneID::PlaneID_t	p2 = `1`
	)		const

Returns the distance between two planes.

Parameters

p1	index of the first plane
p2	index of the second plane
tpc	tpc number within the cryostat
cstat	cryostat number

Returns: distance between the planes

Todo:

add a version with plane IDs

deprecate this function

add a default version for a given TPCID

add a version with two plane indices for a given TPCID

return the absolute value of the distance (makes the order unimportant)

document what will happen (in the future methods) with planes on different TPCs

Definition at line 759 of file GeometryCore.cxx.

References geo::TPCGeo::PlanePitch(), and TPC().

Referenced by PlanePitch(), trkf::SpacePts::produce(), and trkf::Track3Dreco::produce().

   {
     return TPC(tpcid).PlanePitch(p1, p2);
   }

double geo::GeometryCore::PlanePitch	(	geo::PlaneID const &	pid1,
		geo::PlaneID const &	pid2
	)		const

Returns the distance between two planes.

Parameters

p1	index of the first plane
p2	index of the second plane
tpc	tpc number within the cryostat
cstat	cryostat number

Returns: distance between the planes

Todo:

add a version with plane IDs

deprecate this function

add a default version for a given TPCID

add a version with two plane indices for a given TPCID

return the absolute value of the distance (makes the order unimportant)

document what will happen (in the future methods) with planes on different TPCs

Definition at line 768 of file GeometryCore.cxx.

References geo::TPCID::asTPCID(), geo::PlaneID::Plane, and PlanePitch().

   {
     return PlanePitch(pid1.asTPCID(), pid1.Plane, pid2.Plane);
   }

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

Returns the distance between two planes.

Parameters

p1	index of the first plane
p2	index of the second plane
tpc	tpc number within the cryostat
cstat	cryostat number

Returns: distance between the planes

Todo:

add a version with plane IDs

deprecate this function

add a default version for a given TPCID

add a version with two plane indices for a given TPCID

return the absolute value of the distance (makes the order unimportant)

document what will happen (in the future methods) with planes on different TPCs

Definition at line 773 of file GeometryCore.cxx.

References PlanePitch().

   { 
     return PlanePitch(geo::TPCID(cstat, tpc), p1, p2);
   }

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

inline

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 2963 of file GeometryCore.h.

References geo::TPCGeo::PlanePtr().

       {
         geo::TPCGeo const* pTPC = TPCPtr(planeid);
         return pTPC? pTPC->PlanePtr(planeid): nullptr;
       } // PlanePtr()

raw::ChannelID_t geo::GeometryCore::PlaneWireToChannel ( WireID const & wireid ) const

Returns the ID of the TPC channel connected to the specified wire.

Parameters

plane	the number of plane
wire	the number of wire
tpc	the number of TPC
cryostat	the number of cryostat
wireid	the ID of the wire

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

Todo:

Verify the raw::InvalidChannelID part

remove the integers version

Definition at line 1290 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->PlaneWireToChannel(wireid);
   }

raw::ChannelID_t geo::GeometryCore::PlaneWireToChannel	(	unsigned int const	plane,
		unsigned int const	wire,
		unsigned int const	tpc = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the number of TPC readout channels in the detector.

Definition at line 4633 of file GeometryCore.h.

4637 { return PlaneWireToChannel(geo::WireID(cstat, tpc, plane, wire)); }

geo::WireID

Definition: geo_types.h:305

geo::GeometryCore::PlaneWireToChannel

raw::ChannelID_t PlaneWireToChannel(WireID const &wireid) const

Returns the ID of the TPC channel connected to the specified wire.

Definition: GeometryCore.cxx:1290

bool geo::GeometryCore::PointWithinSegments	(	double	A_start_x,
		double	A_start_y,
		double	A_end_x,
		double	A_end_y,
		double	B_start_x,
		double	B_start_y,
		double	B_end_x,
		double	B_end_y,
		double	x,
		double	y
	)

staticprivate

Returns whether x and y are within both specified ranges (A and B).

Definition at line 2005 of file GeometryCore.cxx.

References geo::details::geometry_iterator_types::begin_pos, coordIs, geo::details::geometry_iterator_types::end_pos, geo::details::geometry_iterator_types::undefined_pos, and lar::util::RealComparisons< RealType >::withinSorted().

Referenced by IntersectSegments().

     {
     return coordIs.withinSorted(x, A_start_x, A_end_x)
       && coordIs.withinSorted(y, A_start_y, A_end_y)
       && coordIs.withinSorted(x, B_start_x, B_end_x)
       && coordIs.withinSorted(y, B_start_y, B_end_y)
       ;
     
   } // GeometryCore::PointWithinSegments()

const AuxDetGeo & geo::GeometryCore::PositionToAuxDet	(	geo::Point_t const &	point,
		unsigned int &	ad
	)		const

Returns the auxiliary detector at specified location.

Parameters

point	location to be tested
ad	_(output)_ the auxiliary detector index

Returns: constant reference to AuxDetGeo object of the auxiliary detector

Todo:: what happens if it does not exist?

Definition at line 470 of file GeometryCore.cxx.

References AuxDet(), and FindAuxDetAtPosition().

Referenced by larg4::AuxDetReadoutGeometry::FindAndMakeAuxDet(), FindAuxDetAtPosition(), and PositionToAuxDet().

   {
     // locate the desired Auxiliary Detector
     ad = FindAuxDetAtPosition(point);
     return AuxDet(ad);
   }

const AuxDetGeo & geo::GeometryCore::PositionToAuxDet	(	double const	worldLoc[3],
		unsigned int &	ad
	)		const

Returns the auxiliary detector at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
ad	(output) the auxiliary detector index

Returns: constant reference to AuxDetGeo object of the auxiliary detector

Deprecated:: Use the version with geo::Point_t.

Todo:: what happens if it does not exist?

Definition at line 479 of file GeometryCore.cxx.

References FindAuxDetSensitiveAtPosition(), geo::vect::makePointFromCoords(), and PositionToAuxDet().

480 { return PositionToAuxDet(geo::vect::makePointFromCoords(worldLoc), ad); }

geo::GeometryCore::PositionToAuxDet

AuxDetGeo const & PositionToAuxDet(geo::Point_t const &point, unsigned int &ad) const

Returns the auxiliary detector at specified location.

Definition: GeometryCore.cxx:470

geo::GeometryCore::PositionToAuxDetSensitive

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

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

Definition: geo_vectors_utils.h:1307

const AuxDetSensitiveGeo & geo::GeometryCore::PositionToAuxDetSensitive	(	geo::Point_t const &	point,
		size_t &	ad,
		size_t &	sv
	)		const

Returns the auxiliary detector at specified location.

Parameters

point	location to be tested
ad	_(output)_ the auxiliary detector index
sv	_(output)_ the auxiliary detector sensitive volume index

Returns: reference to AuxDetSensitiveGeo object of the auxiliary detector

Todo:: what happens if it does not exist?

Definition at line 501 of file GeometryCore.cxx.

References AuxDet(), FindAuxDetSensitiveAtPosition(), and geo::AuxDetGeo::SensitiveVolume().

Referenced by FindAuxDetSensitiveAtPosition(), and PositionToAuxDetSensitive().

   {
     // locate the desired Auxiliary Detector
     FindAuxDetSensitiveAtPosition(point, ad, sv);
     return AuxDet(ad).SensitiveVolume(sv);
   }

const AuxDetSensitiveGeo & geo::GeometryCore::PositionToAuxDetSensitive	(	double const	worldLoc[3],
		size_t &	ad,
		size_t &	sv
	)		const

Returns the auxiliary detector at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
ad	(output) the auxiliary detector index
sv	(output) the auxiliary detector sensitive volume index

Returns: reference to AuxDetSensitiveGeo object of the auxiliary detector

Todo:: what happens if it does not exist?

Deprecated:: Use the version with geo::Point_t.

Definition at line 510 of file GeometryCore.cxx.

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

511 { return PositionToAuxDetSensitive(geo::vect::makePointFromCoords(worldLoc), ad, sv); }

const AuxDetSensitiveGeo & PositionToAuxDetSensitive(geo::Point_t const &point, size_t &ad, size_t &sv) const

Returns the auxiliary detector at specified location.

Definition: GeometryCore.cxx:501

geo::GeometryCore::PositionToCryostat

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

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

Definition: geo_vectors_utils.h:1307

geo::CryostatGeo const & geo::GeometryCore::PositionToCryostat ( geo::Point_t const & point ) const

Returns the cryostat at specified location.

Parameters

point the location [cm]

Returns: a constant reference to the geo::CryostatGeo containing point

Exceptions

cet::exception ("Geometry" category) if no cryostat matches

The tolerance used here is the one returned by DefaultWiggle().

Definition at line 419 of file GeometryCore.cxx.

References PositionToCryostatPtr().

Referenced by larg4::LArVoxelReadoutGeometry::FindNestedVolume(), PositionToCryostat(), PositionToTPCID(), and detsim::SimDriftElectrons::produce().

   {
     geo::CryostatGeo const* cstat = PositionToCryostatPtr(point);
     if (!cstat) {
       throw cet::exception("GeometryCore")
         << "Can't find any cryostat at position " << point << "\n";
     }
     return *cstat;
   } // GeometryCore::PositionToCryostat()

CryostatGeo const& geo::GeometryCore::PositionToCryostat ( double const point[3] ) const

inline

Returns the cryostat at specified location.

Parameters

point the location [cm]

Returns: a constant reference to the geo::CryostatGeo containing point

Exceptions

cet::exception ("Geometry" category) if no cryostat matches

The tolerance used here is the one returned by DefaultWiggle().

Definition at line 2067 of file GeometryCore.h.

References geo::vect::makePointFromCoords().

2068 { return PositionToCryostat(geo::vect::makePointFromCoords(point)); }

CryostatGeo const & PositionToCryostat(geo::Point_t const &point) const

Returns the cryostat at specified location.

Definition: GeometryCore.cxx:419

Referenced by larg4::LArVoxelReadoutGeometry::FindNestedVolume(), lar_pandora::PFParticleMonitoring::GetStartAndEndPoints(), lar_pandora::LArPandoraInput::GetTrueX0(), apa::APAGeometryAlg::LineSegChanIntersect(), PositionToTPC(), PositionToTPCptr(), detsim::SimDriftElectrons::produce(), cheat::BackTracker::SpacePointHitsToWeightedXYZ(), and opreco::OpticalRecoAna::update_MC_particle_time().

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

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

Definition: geo_vectors_utils.h:1307

const CryostatGeo & geo::GeometryCore::PositionToCryostat	(	double const	worldLoc[3],
		geo::CryostatID &	cid
	)		const

Returns the cryostat at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
cid	(output) cryostat ID

Returns: a constant reference to the CryostatGeo object of the cryostat

Exceptions

cet::exception ("Geometry" category) if no cryostat matches

The tolerance used here is the one returned by DefaultWiggle().

Deprecated:: Use PositionToCryostat(geo::Point_t const&) instead.

Definition at line 432 of file GeometryCore.cxx.

References Cryostat(), FindCryostatAtPosition(), geo::CryostatID::InvalidID, and PositionToCryostat().

   {
     geo::CryostatID::CryostatID_t cstat = FindCryostatAtPosition(worldLoc);
     
     if(cstat == geo::CryostatID::InvalidID)
       throw cet::exception("GeometryCore") << "Can't find Cryostat for position (" 
                                        << worldLoc[0] << ","
                                        << worldLoc[1] << "," 
                                        << worldLoc[2] << ")\n";
     cid = geo::CryostatID(cstat);
     return Cryostat(cid);
   } // GeometryCore::PositionToCryostat(double[3], CryostatID)

const CryostatGeo & geo::GeometryCore::PositionToCryostat	(	double const	worldLoc[3],
		unsigned int &	cstat
	)		const

Returns the cryostat at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
cstat	(output) number of cryostat

Returns: a constant reference to the CryostatGeo object of the cryostat

Exceptions

cet::exception ("Geometry" category) if no cryostat matches

The tolerance used here is the one returned by DefaultWiggle().

Deprecated:: Use PositionToCryostat(geo::Point_t const&) instead.

Definition at line 446 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, and PositionToCryostat().

   {
     geo::CryostatID cid;
     geo::CryostatGeo const& cryo = PositionToCryostat(worldLoc, cid);
     cstat = cid.Cryostat;
     return cryo;
   } // GeometryCore::PositionToCryostat(double[3], unsigned int)

geo::CryostatID geo::GeometryCore::PositionToCryostatID ( geo::Point_t const & point ) const

Returns the ID of the cryostat at specified location.

Parameters

point the location [cm]

Returns: ID of the cryostat including point (invalid if none)

The tolerance used here is the one returned by DefaultWiggle().

Definition at line 343 of file GeometryCore.cxx.

References FindCryostatAtPosition(), geo::CryostatGeo::ID(), and PositionToCryostatPtr().

Referenced by PositionToCryostatPtr().

   {
     geo::CryostatGeo const* cryo = PositionToCryostatPtr(point);
     return cryo? cryo->ID(): geo::CryostatID{};
   } // GeometryCore::PositionToCryostatID()

geo::CryostatGeo const * geo::GeometryCore::PositionToCryostatPtr ( geo::Point_t const & point ) const

Returns the cryostat at specified location.

Parameters

point the location [cm]

Returns: pointer to the geo::CryostatGeo including point, or nullptr

The tolerance used here is the one returned by DefaultWiggle().

Definition at line 331 of file GeometryCore.cxx.

References fPositionWiggle, IterateCryostats(), and PositionToCryostatID().

Referenced by FindCryostatAtPosition(), FindTPCAtPosition(), GetClosestOpDet(), PositionToCryostat(), PositionToCryostatID(), and PositionToTPCptr().

   {
     for (geo::CryostatGeo const& cryostat: IterateCryostats()) {
       if (cryostat.ContainsPosition(point, 1.0 + fPositionWiggle))
         return &cryostat;
     }
     return nullptr;
   } // GeometryCore::PositionToCryostatPtr()

geo::TPCGeo const & geo::GeometryCore::PositionToTPC ( geo::Point_t const & point ) const

Returns the TPC at specified location.

Parameters

point the location [cm]

Returns: a constant reference to the geo::TPCGeo including point

Exceptions

cet::exception ("Geometry" category) if no TPC matches

Definition at line 378 of file GeometryCore.cxx.

References PositionToTPCptr().

   {
     geo::TPCGeo const* tpc = PositionToTPCptr(point);
     if (!tpc) {
       throw cet::exception("GeometryCore")
         << "Can't find any TPC at position " << point << "\n";
     }
     return *tpc;
   } // GeometryCore::PositionToTPC()

TPCGeo const& geo::GeometryCore::PositionToTPC ( double const point[3] ) const

inline

Returns the TPC at specified location.

Parameters

point the location [cm]

Returns: a constant reference to the geo::TPCGeo including point

Exceptions

cet::exception ("Geometry" category) if no TPC matches

Definition at line 2454 of file GeometryCore.h.

References geo::vect::makePointFromCoords().

2455 { return PositionToTPC(geo::vect::makePointFromCoords(point)); }

geo::GeometryCore::PositionToTPC

geo::TPCGeo const & PositionToTPC(geo::Point_t const &point) const

Returns the TPC at specified location.

Definition: GeometryCore.cxx:378

References geo::PlaneGeo::ID(), art::detail::indent(), geo::AuxDetSensitiveGeo::MaxVerbosity, geo::OpDetGeo::MaxVerbosity, geo::AuxDetGeo::MaxVerbosity, geo::WireGeo::MaxVerbosity, geo::TPCGeo::MaxVerbosity, geo::PlaneGeo::MaxVerbosity, geo::TPCGeo::Nplanes(), geo::AuxDetGeo::NSensitiveVolume(), geo::PlaneGeo::Nwires(), geo::TPCGeo::Plane(), geo::AuxDetSensitiveGeo::PrintAuxDetInfo(), geo::AuxDetGeo::PrintAuxDetInfo(), geo::OpDetGeo::PrintOpDetInfo(), geo::PlaneGeo::PrintPlaneInfo(), geo::TPCGeo::PrintTPCInfo(), geo::WireGeo::PrintWireInfo(), geo::AuxDetGeo::SensitiveVolume(), w, and geo::PlaneGeo::Wire().

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

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

Definition: geo_vectors_utils.h:1307

TPCGeo const & geo::GeometryCore::PositionToTPC	(	double const	worldLoc[3],
		unsigned int &	tpc,
		unsigned int &	cstat
	)		const

Returns the TPC at specified location.

Parameters

point	the location [cm]
tpc	_(output)_ where to store the number of TPC
cstat	_(output)_ where to store the number of cryostat

Returns: a constant reference to the geo::TPCGeo including point

Exceptions

cet::exception ("Geometry" category) if no TPC matches

Deprecated:: Use PositionToTPCID() or PositionToTPC().ID()

Definition at line 401 of file GeometryCore.cxx.

References geo::CryostatID::Cryostat, geo::TPCGeo::ID(), PositionToTPC(), geo::TPCID::TPC, and TPC().

   {
     geo::TPCGeo const& TPC = PositionToTPC(worldLoc);
     cstat = TPC.ID().Cryostat;
     tpc = TPC.ID().TPC;
     return TPC;
   } // GeometryCore::PositionToTPC(double*, TPCID&)

TPCGeo const & geo::GeometryCore::PositionToTPC	(	double const	worldLoc[3],
		TPCID &	tpcid
	)		const

Returns the TPC at specified location.

Parameters

point	the location [cm]
tpcid	_(output)_ where to store the TPC ID

Returns: a constant reference to the geo::TPCGeo including point

Exceptions

cet::exception ("Geometry" category) if no TPC matches

Deprecated:: Use PositionToTPCID() or PositionToTPC().ID()

Definition at line 391 of file GeometryCore.cxx.

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

   {
     geo::TPCGeo const& TPC = PositionToTPC(worldLoc);
     tpcid = TPC.ID();
     return TPC;
   } // GeometryCore::PositionToTPC(double*, TPCID&)

geo::TPCID geo::GeometryCore::PositionToTPCID ( geo::Point_t const & point ) const

Returns the ID of the TPC at specified location.

Parameters

point the location [cm]

Returns: ID of the TPC at specified location, invalid if none

See also: PositionToTPC()

Definition at line 411 of file GeometryCore.cxx.

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

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

geo::TPCGeo const * geo::GeometryCore::PositionToTPCptr ( geo::Point_t const & point ) const

Returns the TPC at specified location.

Parameters

point the location [cm]

Returns: the geo::TPCGeo including point, or nullptr if none

Definition at line 369 of file GeometryCore.cxx.

References fPositionWiggle, PositionToCryostatPtr(), PositionToTPC(), and geo::CryostatGeo::PositionToTPCptr().

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

   {
     geo::CryostatGeo const* cryo = PositionToCryostatPtr(point);
     return cryo? cryo->PositionToTPCptr(point, 1. + fPositionWiggle): nullptr;
   } // GeometryCore::PositionToTPCptr()

template<typename Stream >

void geo::GeometryCore::Print	(	Stream &&	out,
		std::string	indent = `" "`
	)		const

Prints geometry information with maximum verbosity.

Definition at line 5546 of file GeometryCore.h.

Referenced by geo::DumpGeometry::dumpGeometryCore(), and WireEndPoints().

 {
   
   out << "Detector " << DetectorName() << " has "
     << Ncryostats() << " cryostats and "
     << NAuxDets() << " auxiliary detectors:";
  
   auto const& detEnclosureBox = DetectorEnclosureBox();
   out << "\n" << indent << "Detector enclosure: "
     << detEnclosureBox.Min() << " -- " << detEnclosureBox.Max()
     << " cm => ( " << detEnclosureBox.SizeX() << " x "
     << detEnclosureBox.SizeY() << " x "
     << detEnclosureBox.SizeZ() << " ) cm^3"
     ;
   
   for (geo::CryostatGeo const& cryostat: IterateCryostats()) {
     out << "\n" << indent;
     cryostat.PrintCryostatInfo
       (std::forward<Stream>(out), indent + "  ", cryostat.MaxVerbosity);
     
     const unsigned int nTPCs = cryostat.NTPC();
     for(unsigned int t = 0;  t < nTPCs; ++t) {
       const geo::TPCGeo& tpc = cryostat.TPC(t);
       
       out << "\n" << indent << "  ";
       tpc.PrintTPCInfo
         (std::forward<Stream>(out), indent + "    ", tpc.MaxVerbosity);
       
       const unsigned int nPlanes = tpc.Nplanes();
       for(unsigned int p = 0; p < nPlanes; ++p) {
         const geo::PlaneGeo& plane = tpc.Plane(p);
         const unsigned int nWires = plane.Nwires();
         
         out << "\n" << indent << "    ";
         plane.PrintPlaneInfo
           (std::forward<Stream>(out), indent + "      ", plane.MaxVerbosity);
         
         for(unsigned int w = 0;  w < nWires; ++w) {
           const geo::WireGeo& wire = plane.Wire(w);
           geo::WireID wireID(plane.ID(), w);
           
           // the wire should be aligned on z axis, half on each side of 0,
           // in its local frame
           out << "\n" << indent << "      " << wireID << " ";
           wire.PrintWireInfo
             (std::forward<Stream>(out), indent + "      ", wire.MaxVerbosity);
         } // for wire
       } // for plane
     } // for TPC
     
     unsigned int nOpDets = cryostat.NOpDet();
     for (unsigned int iOpDet = 0; iOpDet < nOpDets; ++iOpDet) {
       geo::OpDetGeo const& opDet = cryostat.OpDet(iOpDet);
       out << "\n" << indent << "  [OpDet #" << iOpDet << "] ";
       opDet.PrintOpDetInfo
         (std::forward<Stream>(out), indent + "  ", opDet.MaxVerbosity);
     } // for
   } // for cryostat
   
   unsigned int const nAuxDets = NAuxDets();
   for (unsigned int iDet = 0; iDet < nAuxDets; ++iDet) {
     geo::AuxDetGeo const& auxDet = AuxDet(iDet);
     
     out << "\n" << indent << "[#" << iDet << "] ";
     auxDet.PrintAuxDetInfo
       (std::forward<Stream>(out), indent + "  ", auxDet.MaxVerbosity);
     
     unsigned int const nSensitive = auxDet.NSensitiveVolume();
     switch (nSensitive) {
       case 0: break;
       case 1: {
         geo::AuxDetSensitiveGeo const& auxDetS = auxDet.SensitiveVolume(0U);
         out << "\n" << indent << "  ";
         auxDetS.PrintAuxDetInfo
           (std::forward<Stream>(out), indent + "    ", auxDetS.MaxVerbosity);
         break;
       }
       default:
         for (unsigned int iSens = 0; iSens < nSensitive; ++iSens) {
           out << "\n" << indent << "[#" << iSens << "] ";
           geo::AuxDetSensitiveGeo const& auxDetS
             = auxDet.SensitiveVolume(iSens);
           auxDetS.PrintAuxDetInfo
             (std::forward<Stream>(out), indent + "    ", auxDetS.MaxVerbosity);
         } // for
         break;
     } // if sensitive detectors
   
   } // for auxiliary detector
   
   out << '\n';
   
 } // geo::GeometryCore::Print()

std::string geo::GeometryCore::ROOTFile ( ) const

inline

Returns the full directory path to the geometry file source.

Returns: the full directory path to the geometry file source

This is the full path of the source of the detector geometry GeometryCore relies on.

Definition at line 1714 of file GeometryCore.h.

Referenced by geo::DumpGeometry::dumpGeometryCore(), and evgen::GENIEGen::GENIEGen().

1714 { return fROOTfile; }

geo::GeometryCore::fROOTfile

std::string fROOTfile

path to geometry file for geometry in GeometryCore

Definition: GeometryCore.h:5396

TGeoManager * geo::GeometryCore::ROOTGeoManager ( ) const

Access to the ROOT geometry description manager.

Definition at line 185 of file GeometryCore.cxx.

Referenced by FindAllVolumePaths(), FindAllVolumes(), FindDetectorEnclosure(), evgen::GENIEGen::GENIEGen(), and evgen::RadioGen::SampleOne().

   {
     return gGeoManager;
   }

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

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

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.

Definition at line 1842 of file GeometryCore.cxx.

References fChannelMapAlg, and FirstChannelInROP().

Referenced by reco3d::TripletFinder::FillBadMap(), reco3d::TripletFinder::FillHitMap(), and ROPtoWirePlanes().

   {
     return fChannelMapAlg->ROPtoTPCs(ropid);
   } // GeometryCore::ROPtoTPCs()

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

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

Parameters

ropid ID of the readout plane

Returns: list of ID of wire planes belonging to the specified ROP

If ropid is an invalid ID, an empty list 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 1834 of file GeometryCore.cxx.

References fChannelMapAlg, and ROPtoTPCs().

Referenced by WirePlaneToROP().

   {
     return fChannelMapAlg->ROPtoWirePlanes(ropid);
   } // GeometryCore::ROPtoWirePlanes()

void geo::GeometryCore::SetDetectorName ( std::string new_name )

inlineprotected

Sets the detector name.

Definition at line 5332 of file GeometryCore.h.

Referenced by geo::Geometry::preBeginRun().

5332 { fDetectorName = new_name; }

geo::GeometryCore::fDetectorName

std::string fDetectorName

Name of the detector.

Definition: GeometryCore.h:5394

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

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 537 of file GeometryCore.cxx.

References WirePlaneToROP().

Referenced by caldata::CalWireAna::analyze(), evd::Landed::analyze(), lar::example::AnalysisExample::analyze(), cluster::HoughBaseAlg::FastTransform(), reco3d::TripletFinder::FillHitMap(), hit::HitCheater::FindHitsOnChannel(), calo::GeneralCalorimetry::GeneralCalorimetry(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire(), shwf::ShowerReco::LongTransEnergy(), trkf::SpacePointAlg::makeSpacePoints(), evd::TWQMultiTPCProjectionView::PrintCharge(), evd::TWQProjectionView::PrintCharge(), hit::TTHitFinder::produce(), hit::FFTHitFinder::produce(), wc::MergeWireCell::produce(), caldata::CalWire::produce(), cluster::DBcluster::produce(), cluster::fuzzyCluster::produce(), trkf::SpacePts::produce(), hit::RawHitFinder::produce(), reco3d::SpacePointSolver::produce(), detsim::SimWire::produce(), evd::RawDataDrawer::QueueDrawingBoxes(), lariov::DetPedestalRetrievalAlg::Reconfigure(), hit::RFFHitFinderAlg::Run(), cluster::fuzzyClusterAlg::run_fuzzy_cluster(), cluster::HoughBaseAlg::Transform(), evd::TWireProjPad::TWireProjPad(), trkf::SpacePointAlg::update(), and evd::RecoBaseDrawer::Wire2D().

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

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

Returns the type of signal on the specified TPC channel.

Parameters

channel TPC channel ID

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

Todo:: verify that kMysteryType is returned on invalid channel

Definition at line 531 of file GeometryCore.cxx.

References fChannelMapAlg.

   {
     return fChannelMapAlg->SignalTypeForChannel(channel);
   }

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

Returns the type of signal of channels in specified readout plane.

Parameters

ropid readout plane ID

Returns: the type of signal on the specified ROP

Assumes that all the channels on the readout plane have the same signal type. If ropid is an invalid ID, geo::kMysteryType 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 1863 of file GeometryCore.cxx.

References fChannelMapAlg.

                                                                        {
     return fChannelMapAlg->SignalTypeForROPID(ropid);
   } // GeometryCore::SignalType(ROPID)

void geo::GeometryCore::SortGeometry ( geo::GeoObjectSorter const & sorter )

private

Runs the sorting of geometry with the sorter provided by channel mapping.

Definition at line 146 of file GeometryCore.cxx.

References AuxDets(), Cryostats(), geo::GeoObjectSorter::SortAuxDets(), util::SortByPointers(), and geo::GeoObjectSorter::SortCryostats().

Referenced by ApplyChannelMap().

                                                                   {
     
     mf::LogInfo("GeometryCore") << "Sorting volumes...";
     
     sorter.SortAuxDets(AuxDets());
     
     //
     // cryostats
     //
     util::SortByPointers(Cryostats(),
       [&sorter](auto& coll){ sorter.SortCryostats(coll); });
     
     geo::CryostatID::CryostatID_t c = 0;
     for (geo::CryostatGeo& cryo: Cryostats())
     {
       cryo.SortSubVolumes(sorter);
       cryo.UpdateAfterSorting(geo::CryostatID(c));
       ++c;
     } // for
     
   } // GeometryCore::SortGeometry()

geo::Length_t geo::GeometryCore::SurfaceY ( ) const

inline

The position of the detector respect to earth surface.

Returns: typical y position at surface in units of cm

This is the depth (y) of the surface (where earth meets air) for this detector site. The number is expressed in world coordinates and in centimetres, and it represents the y coordinate of earth surface. A negative value means that the origin of coordinates, typically matching the detector centre, is above surface.

Todo:: check that this is actually how it is used

Definition at line 1786 of file GeometryCore.h.

Referenced by evgen::CosmicsGen::produce().

1786 { return fSurfaceY; }

geo::GeometryCore::fSurfaceY

double fSurfaceY

The point where air meets earth for this detector.

Definition: GeometryCore.h:5393

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

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 1574 of file GeometryCore.cxx.

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

Referenced by IntersectSegments(), ThirdPlane_dTdW(), and ThirdPlaneSlope().

   {
     // how many planes in the TPC pid1 belongs to:
     const unsigned int nPlanes = Nplanes(pid1);
     if(nPlanes != 3) {
       throw cet::exception("GeometryCore")
         << "ThirdPlane() supports only TPCs with 3 planes, and I see "
         << nPlanes << " instead\n";
     }
     
     geo::PlaneID::PlaneID_t target_plane = nPlanes;
     for (geo::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 geo::PlaneID(pid1, target_plane);
   } // GeometryCore::ThirdPlane()

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

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 1651 of file GeometryCore.cxx.

References CheckIndependentPlanesOnSameTPC(), ComputeThirdPlane_dTdW(), geo::PlaneGeo::PhiZ(), geo::TPCGeo::Plane(), TPC(), and geo::PlaneGeo::WirePitch().

Referenced by ThirdPlane_dTdW().

           {
     
     CheckIndependentPlanesOnSameTPC(pid1, pid2, "ThirdPlane_dTdW()");
     
     geo::TPCGeo const& TPC = this->TPC(pid1);
     
     double angle[3], pitch[3];
     geo::PlaneGeo const* const planes[3]
       = { &TPC.Plane(pid1), &TPC.Plane(pid2), &TPC.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();
     } // for
 
     return ComputeThirdPlane_dTdW(
       angle[0], pitch[0], slope1,
       angle[1], pitch[1], slope2,
       angle[2], pitch[2]
       );
     
   } // GeometryCore::ThirdPlane_dTdW()

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

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 1682 of file GeometryCore.cxx.

References ComputeThirdPlaneSlope(), ThirdPlane(), and ThirdPlane_dTdW().

           {
     geo::PlaneID target_plane = ThirdPlane(pid1, pid2);
     return ThirdPlane_dTdW(pid1, slope1, pid2, slope2, target_plane);
   } // ThirdPlane_dTdW()

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

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 1621 of file GeometryCore.cxx.

References CheckIndependentPlanesOnSameTPC(), ComputeThirdPlaneSlope(), geo::PlaneGeo::PhiZ(), geo::TPCGeo::Plane(), and TPC().

Referenced by ThirdPlaneSlope().

           {
     
     CheckIndependentPlanesOnSameTPC(pid1, pid2, "ThirdPlaneSlope()");
     
     geo::TPCGeo const& TPC = this->TPC(pid1);
 
     // 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(
       TPC.Plane(pid1).PhiZ(), slope1,
       TPC.Plane(pid2).PhiZ(), slope2,
       TPC.Plane(output_plane).PhiZ()
       );
   } // ThirdPlaneSlope()

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

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 1642 of file GeometryCore.cxx.

References ThirdPlane(), and ThirdPlaneSlope().

           {
     geo::PlaneID target_plane = ThirdPlane(pid1, pid2);
     return ThirdPlaneSlope(pid1, slope1, pid2, slope2, target_plane);
   } // ThirdPlaneSlope()

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

inline

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 4278 of file GeometryCore.h.

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

double geo::GeometryCore::ThirdPlaneSlope	(	unsigned int	plane1,
		double	slope1,
		unsigned int	plane2,
		double	slope2,
		unsigned int	tpc,
		unsigned int	cstat
	)		const

inline

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 4287 of file GeometryCore.h.

       {
         return ThirdPlaneSlope
           (plane1, slope1, plane2, slope2, geo::TPCID(cstat, tpc));
       }

double geo::GeometryCore::TotalMass ( ) const

inline

Returns the total mass [kg] of the specified volume (default: world).

Definition at line 1865 of file GeometryCore.h.

Referenced by evgen::GENIEGen::GENIEGen().

1865 { return TotalMass(GetWorldVolumeName()); }

geo::GeometryCore::GetWorldVolumeName

const std::string GetWorldVolumeName() const

Return the name of the world volume (needed by Geant4 simulation)

Definition: GeometryCore.cxx:578

geo::GeometryCore::TotalMass

double TotalMass() const

Returns the total mass [kg] of the specified volume (default: world).

Definition: GeometryCore.h:1865

double geo::GeometryCore::TotalMass ( std::string vol ) const

Returns the total mass [kg] of the specified volume (default: world).

Definition at line 1112 of file GeometryCore.cxx.

References MassBetweenPoints().

   {
     //the TGeoNode::GetVolume() returns the TGeoVolume of the detector outline
     //and ROOT calculates the mass in kg for you
     TGeoVolume *gvol = gGeoManager->FindVolumeFast(vol.c_str());
     if(gvol) return gvol->Weight();
 
     throw cet::exception("GeometryCore") << "could not find specified volume '" 
                                          << vol
                                          << " 'to determine total mass\n"; 
   }

unsigned int geo::GeometryCore::TotalNTPC ( ) const

Returns the total number of TPCs in the detector.

Definition at line 827 of file GeometryCore.cxx.

References evd::details::begin(), Cryostats(), evd::details::end(), and geo::CryostatGeo::NTPC().

Referenced by mvapid::MVAAlg::GetDetectorEdges().

                                              {
     // it looks like C++11 lambdas have made STL algorithms easier to use,
     // but only so much:
     return std::accumulate(Cryostats().begin(), Cryostats().end(), 0U,
       [](unsigned int sum, geo::CryostatGeo const& cryo)
         { return sum + cryo.NTPC(); }
       );
   } // GeometryCore::TotalNTPC()

TPCGeo const& geo::GeometryCore::TPC	(	unsigned int const	tpc = `0`,
		unsigned int const	cstat = `0`
	)		const

inline

Returns the specified TPC.

Parameters

tpcid	ID of the tpc
tpc	tpc number within the cryostat
cstat	number of cryostat

Returns: a constant reference to the specified TPC

Exceptions

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

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

Todo:: remove the version with integers

Definition at line 2393 of file GeometryCore.h.

2394 { return TPC(geo::TPCID(cstat, tpc)); }

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

geo::GeometryCore::Cryostat

TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified TPC.

Definition: GeometryCore.h:2393

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

inline

Returns the specified TPC.

Parameters

tpcid	ID of the tpc
tpc	tpc number within the cryostat
cstat	number of cryostat

Returns: a constant reference to the specified TPC

Exceptions

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

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

Todo:: remove the version with integers

Definition at line 2395 of file GeometryCore.h.

References geo::CryostatGeo::TPC().

2396 { return Cryostat(tpcid).TPC(tpcid); }

CryostatGeo const & Cryostat(geo::CryostatID const &cryoid) const

Returns the specified cryostat.

Definition: GeometryCore.cxx:280

geo::CryostatGeo::TPC

const TPCGeo & TPC(unsigned int itpc) const

Return the itpc'th TPC in the cryostat.

Definition: CryostatGeo.cxx:174

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

inline

Returns the specified TPC.

Parameters

tpcid TPC ID

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

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

Definition at line 2410 of file GeometryCore.h.

References geo::CryostatGeo::TPCPtr().

       {
         CryostatGeo const* pCryo = CryostatPtr(tpcid);
         return pCryo? pCryo->TPCPtr(tpcid): nullptr;
       } // TPCPtr()

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

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

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.

Definition at line 1798 of file GeometryCore.cxx.

References fChannelMapAlg.

Referenced by TPCtoTPCset().

   {
     return fChannelMapAlg->TPCsetToTPCs(tpcsetid);
   } // GeometryCore::TPCsetToTPCs()

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

Returns the ID of the TPC set tpcid belongs to.

Definition at line 1790 of file GeometryCore.cxx.

References fChannelMapAlg, and TPCsetToTPCs().

Referenced by FindTPCsetAtPosition().

   {
     return fChannelMapAlg->TPCtoTPCset(tpcid);
   } // GeometryCore::TPCtoTPCset()

void geo::GeometryCore::UpdateAfterSorting ( )

private

Performs all the updates needed after sorting.

Definition at line 170 of file GeometryCore.cxx.

References allViews, Cryostats(), IterateTPCs(), and Ncryostats().

Referenced by ApplyChannelMap().

                                         {
     
     for (size_t c = 0; c < Ncryostats(); ++c)
       Cryostats()[c].UpdateAfterSorting(geo::CryostatID(c));
     
     allViews.clear();
     for (geo::TPCGeo const& tpc: IterateTPCs()) {
       auto const& TPCviews = tpc.Views();
       allViews.insert(TPCviews.cbegin(), TPCviews.cend());
     }
     
   } // GeometryCore::UpdateAfterSorting()

bool geo::GeometryCore::ValueInRange	(	double	value,
		double	min,
		double	max
	)		const

Returns whether a value is within the specified range.

Parameters

value	the value to be tested
min	the lower boundary
max	the upper boundary

Returns: whether the value is within range

If min is larger than max, they are swapped. A tolerance of 10^-6 (absolute) is used.

Todo:

Use wiggle instead of 10^-6

resort source code for a bit of speed up

Definition at line 1298 of file GeometryCore.cxx.

References e, max, and WireEndPoints().

Referenced by apa::APAGeometryAlg::LineSegChanIntersect().

   {
     if(min>max) std::swap(min,max);//protect against funny business due to wire angles
     if (std::abs(value-min)<1e-6||std::abs(value-max)<1e-6) return true;
     return (value>=min) && (value<=max);
   }

View_t geo::GeometryCore::View ( geo::PlaneID const & pid ) const

Returns the view (wire orientation) on the channels of specified TPC plane.

Parameters

plane TPC plane ID

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

Definition at line 560 of file GeometryCore.cxx.

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

Referenced by calo::TrackCalorimetryAlg::AnalyzeHit(), apa::APAGeometryAlg::APAChannelsIntersect(), apa::APAGeometryAlg::APAView(), tca::FillWireHitRange(), hit::HitCheater::FindHitsOnChannel(), apa::APAGeometryAlg::FirstChannelInView(), calo::GeneralCalorimetry::GeneralCalorimetry(), corner::CornerFinderAlg::get_feature_points(), corner::CornerFinderAlg::get_feature_points_LineIntegralScore(), apa::APAGeometryAlg::Init(), apa::DisambigAlg::MakeCloseHits(), apa::APAGeometryAlg::NearestWireIDOnChan(), wc::MergeWireCell::produce(), cluster::TrajCluster::produce(), vertex::VertexMatch::produce(), reco3d::SpacePointSolver::produce(), vertex::HarrisVertexFinder::produce(), vertex::FeatureVertexFinder::produce(), shower::TCShowerTemplateMaker::showerProfileTrue(), cheat::BackTracker::TrackIdToSimIDEs_Ps(), and View().

   {
     return pid? Plane(pid).View(): geo::kUnknown;
   }

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

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 554 of file GeometryCore.cxx.

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

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

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

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 1857 of file GeometryCore.cxx.

References fChannelMapAlg, and View().

                                                               {
     return View(fChannelMapAlg->FirstWirePlaneInROP(ropid));
   } // GeometryCore::View()

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

inline

Returns a list of possible views in the detector.

Returns: the set of views

Definition at line 4602 of file GeometryCore.h.

Referenced by cluster::HoughLineFinderAna::analyze(), vertex::FeatureVertexFinderAna::analyze(), detinfo::DetectorPropertiesStandard::CheckTimeOffsetConfigurationAfterSetup(), cluster::EndPointAlg::EndPoint(), cluster::HoughBaseAlg::FastTransform(), pma::PMAlgTracker::PMAlgTracker(), and vertex::VertexFinder2D::produce().

4602 { return allViews; }

geo::GeometryCore::allViews

std::set< geo::View_t > allViews

All views in the detector.

Definition: GeometryCore.h:5404

std::string geo::GeometryCore::VolumeName ( geo::Point_t const & point ) const

Returns the name of the deepest volume containing specified point.

Parameters

point the location to query, in world coordinates

Returns: name of the volume containing the point

Todo:

what happens if none?

Unify the coordinates type

Definition at line 899 of file GeometryCore.cxx.

References art::productstatus::unknown(), and WorldVolume().

Referenced by larg4::LArG4Ana::analyze(), and LArStackingAction::InsideTPC().

   {
     // check that the given point is in the World volume at least
     TGeoVolume const*volWorld = WorldVolume();
     double halflength = ((TGeoBBox*)volWorld->GetShape())->GetDZ();
     double halfheight = ((TGeoBBox*)volWorld->GetShape())->GetDY();
     double halfwidth  = ((TGeoBBox*)volWorld->GetShape())->GetDX();
     if(std::abs(point.x()) > halfwidth  ||
        std::abs(point.y()) > halfheight ||
        std::abs(point.z()) > halflength
        ){
       mf::LogWarning("GeometryCoreBadInputPoint") << "point (" << point.x() << ","
                                               << point.y() << "," << point.z() << ") "
                                               << "is not inside the world volume "
                                               << " half width = " << halfwidth
                                               << " half height = " << halfheight
                                               << " half length = " << halflength
                                               << " returning unknown volume name";
       const std::string unknown("unknownVolume");
       return unknown;
     }
     
     return gGeoManager->FindNode(point.X(), point.Y(), point.Z())->GetName();
   }

std::string geo::GeometryCore::VolumeName ( TVector3 const & point ) const

inline

Returns a string with the name of the detector, as configured.

Definition at line 1816 of file GeometryCore.h.

References geo::vect::toPoint().

1817 { return VolumeName(geo::vect::toPoint(point)); }

geo::GeometryCore::VolumeName

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

std::string VolumeName(geo::Point_t const &point) const

Returns the name of the deepest volume containing specified point.

Definition: GeometryCore.cxx:899

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

inline

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 3370 of file GeometryCore.h.

Referenced by IntersectSegments(), and trkf::Track3DKalmanHitAlg::makeSeed().

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

geo::GeometryCore::Plane

PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified wire.

Definition: GeometryCore.h:2945

geo::PlaneGeo::Wire

WireGeo const & Wire(unsigned int iwire) const

Definition: PlaneGeo.cxx:506

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

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

Parameters

view	the view
TPC	the index of the TPC in the specified cryostat
Cryo	the cryostat
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 802 of file GeometryCore.cxx.

References geo::TPCGeo::Nplanes(), geo::TPCGeo::Plane(), geo::PlaneGeo::ThetaZ(), TPC(), geo::PlaneGeo::View(), and geo::PlaneGeo::ViewName().

Referenced by util::GeometryUtilities::CalculatePitch(), util::GeometryUtilities::CalculatePitchPolar(), util::GeometryUtilities::Get2DangleFrom3D(), lar_pandora::LArPandoraGeometry::LoadGeometry(), shower::EMShowerAlg::MakeShower(), util::GeometryUtilities::PitchInView(), calo::Calorimetry::produce(), WirePitch(), lar_pandora::PFParticleHitDumper::YZtoU(), and lar_pandora::PFParticleHitDumper::YZtoV().

   {
     // loop over the planes in cryostat 0, tpc 0 to find the plane with the 
     // specified view
     geo::TPCGeo const& TPC = this->TPC(tpcid);
     for (unsigned int p = 0; p < TPC.Nplanes(); ++p) {
       geo::PlaneGeo const& plane = TPC.Plane(p);
       if (plane.View() == view) return plane.ThetaZ();
     } // for
     throw cet::exception("GeometryCore") << "WireAngleToVertical(): no view \""
       << geo::PlaneGeo::ViewName(view) << "\" (#" << ((int) view)
       << ") in " << std::string(tpcid);
   } // GeometryCore::WireAngleToVertical()

double geo::GeometryCore::WireAngleToVertical	(	geo::View_t	view,
		int	TPC = `0`,
		int	Cryo = `0`
	)		const

inline

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

Parameters

view	the view
TPC	the index of the TPC in the specified cryostat
Cryo	the cryostat
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 3753 of file GeometryCore.h.

3754 { return WireAngleToVertical(view, geo::TPCID(Cryo, TPC)); }

The data type to uniquely identify a TPC.

Definition: geo_types.h:195

geo::GeometryCore::WireAngleToVertical

TPCGeo const & TPC(unsigned int const tpc=0, unsigned int const cstat=0) const

Returns the specified TPC.

Definition: GeometryCore.h:2393

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

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

Definition: GeometryCore.cxx:802

geo::Length_t geo::GeometryCore::WireCoordinate	(	double	YPos,
		double	ZPos,
		geo::PlaneID const &	planeid
	)		const

Returns the index of the nearest wire 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

Deprecated:: Use WireCoordinate(geo::Point_t const&, geo::PlaneID const&) instead

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.

Todo:: Unify (y, z) coordinate

Definition at line 1189 of file GeometryCore.cxx.

References fChannelMapAlg, and NearestWire().

Referenced by ChannelToROP(), tca::ChgFracBetween(), tca::ChgFracNearEnd(), evd::RecoBaseDrawer::DrawProng2D(), evd::RecoBaseDrawer::DrawTrackVertexAssns2D(), tca::Find3DVertices(), tca::FindCompleteness(), tca::FindXMatches(), tca::FitTp3(), tca::FitTp3s(), tca::FollowTp3s(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire(), ems::Hit2D::Hit2D(), lar_cluster3d::PCASeedFinderAlg::LineFit2DHits(), lar_cluster3d::HoughSeedFinderAlg::LineFit2DHits(), tca::MakeBareTP(), tca::MakeTp3(), shower::EMShowerAlg::OrderShowerHits(), tca::PosInPlane(), evd::RecoBaseDrawer::Prong2D(), evd::RecoBaseDrawer::Slice2D(), trkf::TrackLineFitAlg::TrkLineFit(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), pma::ProjectionMatchingAlg::validate_on_adc_test(), evd::RecoBaseDrawer::Vertex2D(), trkf::VertexFitAlg::VertexFit(), and WireCoordinate().

   {
     return fChannelMapAlg->WireCoordinate(YPos, ZPos, planeid);
   }

geo::Length_t geo::GeometryCore::WireCoordinate	(	double	YPos,
		double	ZPos,
		unsigned int	PlaneNo,
		unsigned int	TPCNo,
		unsigned int	cstat
	)		const

inline

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

Parameters

YPos	y coordinate on the wire plane
ZPos	z coordinate on the wire plane
PlaneNo	number of plane
TPCNo	number of TPC
cstat	number of cryostat

Returns: an index interpolation between the two nearest wires

See also: ChannelMapAlg::WireCoordinate()

Deprecated:: Use the version with plane ID instead

Definition at line 4007 of file GeometryCore.h.

4011 { return WireCoordinate(YPos, ZPos, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

geo::GeometryCore::WireCoordinate

geo::Length_t WireCoordinate(double YPos, double ZPos, geo::PlaneID const &planeid) const

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

Definition: GeometryCore.cxx:1189

geo::GeometryCore::WireCoordinate

The data type to uniquely identify a Plane.

Definition: geo_types.h:250

geo::Length_t geo::GeometryCore::WireCoordinate	(	geo::Point_t const &	pos,
		geo::PlaneID const &	planeid
	)		const

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

Parameters

pos	world coordinates of the position (it will be projected)
planeid	ID of the plane

Returns: an index interpolation between the two nearest wires

See also: ChannelMapAlg::WireCoordinate()

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.

Definition at line 1182 of file GeometryCore.cxx.

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

   {
     return Plane(planeid).WireCoordinate(pos);
   }

geo::Length_t geo::GeometryCore::WireCoordinate	(	TVector3 const &	pos,
		geo::PlaneID const &	planeid
	)		const

inline

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

Parameters

pos	world coordinates of the position (it will be projected)
planeid	ID of the plane

Returns: an index interpolation between the two nearest wires

See also: ChannelMapAlg::WireCoordinate()

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.

Definition at line 4028 of file GeometryCore.h.

References geo::vect::toPoint(), x, y, and z.

4029 { return WireCoordinate(geo::vect::toPoint(pos), planeid); }

geo::Length_t WireCoordinate(double YPos, double ZPos, geo::PlaneID const &planeid) const

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

Definition: GeometryCore.cxx:1189

Referenced by trkf::BezierTrackerAlgorithm::CalculateGeometricalElements(), trkf::SeedFinderAlgorithm::CalculateGeometricalElements(), cmtool::CFAlgoWireOverlap::Float(), cmtool::CFAlgoVolumeOverlap::Float(), trkf::BezierTrack::GetClosestApproach(), trkf::BezierTrack::GetClosestApproaches(), trkf::SeedFinderAlgorithm::GetHitDistAndProj(), trkf::BezierTrack::GetTrackPitch(), IncrementID(), IntersectSegments(), apa::DisambigAlg::TrivialDisambig(), and ValueInRange().

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

Convert the specified point into a geo::Point_t.

Definition: geo_vectors_utils.h:1281

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

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 1307 of file GeometryCore.cxx.

References geo::GeometryCore::Segment< Point >::end(), and geo::GeometryCore::Segment< Point >::start().

   {
     Segment_t result = WireEndPoints(wireid);
     
     xyzStart[0] = result.start().X();
     xyzStart[1] = result.start().Y();
     xyzStart[2] = result.start().Z();
     xyzEnd[0]   = result.end().X();
     xyzEnd[1]   = result.end().Y();
     xyzEnd[2]   = result.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]);
     }
     
   } // GeometryCore::WireEndPoints(WireID, 2x double*)

void geo::GeometryCore::WireEndPoints	(	unsigned int	cstat,
		unsigned int	tpc,
		unsigned int	plane,
		unsigned int	wire,
		double *	xyzStart,
		double *	xyzEnd
	)		const

inline

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

Parameters

cstat	cryostat number
tpc	tpc number within the cryostat
plane	plane number within the TPC
wire	wire number within the plane
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 3807 of file GeometryCore.h.

3811 { WireEndPoints(geo::WireID(cstat, tpc, plane, wire), xyzStart, xyzEnd); }

geo::WireID

Definition: geo_types.h:305

geo::GeometryCore::WireEndPoints

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

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

Definition: GeometryCore.cxx:1307

template<typename Point = DefaultPoint_t>

Segment<Point> geo::GeometryCore::WireEndPoints ( geo::WireID const & wireID ) const

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)

template<typename Point >

geo::GeometryCore::Segment<Point> geo::GeometryCore::WireEndPoints ( geo::WireID const & wireid ) const

Definition at line 5536 of file GeometryCore.h.

References geo::WireGeo::GetEnd(), geo::WireGeo::GetStart(), and Print().

 {
   geo::WireGeo const& wire = Wire(wireid);
   return { wire.GetStart<Point>(), wire.GetEnd<Point>() };
 } // geo::GeometryCore::WireEndPoints(WireID)

bool geo::GeometryCore::WireIDIntersectionCheck	(	const geo::WireID &	wid1,
		const geo::WireID &	wid2
	)		const

private

Wire ID check for WireIDsIntersect methods.

Definition at line 1976 of file GeometryCore.cxx.

References geo::TPCID::asTPCID(), HasWire(), Nwires(), geo::PlaneID::Plane, and geo::WireID::planeID().

Referenced by GetClosestOpDet(), and IntersectSegments().

   {
     if (wid1.asTPCID() != wid2) {
       mf::LogError("WireIDIntersectionCheck")
         << "Comparing two wires on different TPCs: return failure.";
       return false;
     }
     if (wid1.Plane == wid2.Plane) {
       mf::LogError("WireIDIntersectionCheck")
         << "Comparing two wires in the same plane: return failure";
       return false;
     }
     if (!HasWire(wid1)) {
       mf::LogError("WireIDIntersectionCheck")
         << "1st wire " << wid1 << " does not exist (max wire number: "
         << Nwires(wid1.planeID()) << ")";
       return false;
     }
     if (!HasWire(wid2)) {
       mf::LogError("WireIDIntersectionCheck")
         << "2nd wire " << wid2 << " does not exist (max wire number: "
         << Nwires(wid2.planeID()) << ")";
       return false;
     }
     return true;
   } // GeometryCore::WireIDIntersectionCheck()

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

Computes the intersection between two wires.

Parameters

wid1	ID of the first wire
wid2	ID of the other wire
intersection	(output) the intersection point (global coordinates)

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

The "intersection" refers to the projection of the wires into the same wire plane. The coordinate along the drift direction is arbitrarily set to the one of the first wire. Wires are assumed to have at most one intersection. If wires are parallel, intersection will have all components set to infinity (std::numeric_limits<>::infinity()) and false is returned. If the intersection is outside the TPC, false is also returned, but the intersection point will contain that intersection.

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

Referenced by apa::APAGeometryAlg::APAChannelsIntersect(), ChannelsIntersect(), apa::DisambigAlg::CompareViews(), shower::EMShowerAlg::Construct3DPoint(), evd::TWQProjectionView::FindEndPoint(), evd::TWQProjectionView::FindLineLength(), IntersectionPoint(), IntersectSegments(), reco3d::IntersectionCache::ISect(), lar_cluster3d::StandardHit3DBuilder::makeDeadChannelPair(), lar_cluster3d::StandardHit3DBuilder::makeHitPair(), and apa::APAGeometryAlg::ThreeChanPos().

bool geo::GeometryCore::WireIDsIntersect	(	WireID const &	wid1,
		WireID const &	wid2,
		TVector3 &	intersection
	)		const

Computes the intersection between two wires.

Parameters

wid1	ID of the first wire
wid2	ID of the other wire
intersection	(output) the intersection point (global coordinates)

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

The "intersection" refers to the projection of the wires into the same wire plane. The coordinate along the drift direction is arbitrarily set to the one of the first wire. Wires are assumed to have at most one intersection. If wires are parallel, intersection will have all components set to infinity (std::numeric_limits<>::infinity()) and false is returned. If the intersection is outside the TPC, false is also returned, but the intersection point will contain that intersection.

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

bool geo::GeometryCore::WireIDsIntersect	(	WireID const &	wid1,
		WireID const &	wid2,
		WireIDIntersection &	widIntersect
	)		const

Computes the intersection between two wires.

Parameters

wid1	ID of the first wire
wid2	ID of the other wire
widIntersect	(output) the coordinate of the intersection point

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.

WireGeo const& geo::GeometryCore::WireIDToWireGeo ( geo::WireID const & wireid ) const

inline

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 3372 of file GeometryCore.h.

Referenced by trkf::SpacePointAlg::compatible(), trkf::SpacePointAlg::fillComplexSpacePoint(), trkf::SpacePointAlg::fillSpacePoint(), lar_cluster3d::PrincipalComponentsAlg::getHit2DPocaToAxis(), cluster::MergeClusterAlg::GlobalWire(), cluster::BlurredClusteringAlg::GlobalWire(), shower::EMShowerAlg::GlobalWire(), ems::Hit2D::Hit2D(), trkf::SpacePointAlg::makeSpacePoints(), lar_cluster3d::PrincipalComponentsAlg::PCAAnalysis_2D(), lar_cluster3d::PrincipalComponentsAlg::PCAAnalysis_calc2DDocas(), trkf::SpacePointAlg::separation(), trkf::TrackKalmanFitter::setupInputStates(), trkf::SurfWireLine::SurfWireLine(), trkf::SurfWireX::SurfWireX(), and apa::APAGeometryAlg::ThreeChanPos().

3373 { return Wire(wireid); }

geo::GeometryCore::Wire

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

Returns the specified wire.

Definition: GeometryCore.h:3370

geo::Length_t geo::GeometryCore::WirePitch ( geo::PlaneID const & planeid ) const

Returns the distance between two consecutive wires.

Parameters

p	plane number within the TPC
tpc	tpc number within the cryostat
cstat	cryostat number

Returns: the distance between the two wires

Note: The current geometry assumptions imply that wire pitch is constant between all wires on the same wire plane. This is an assumption non-trivial to remove.

Todo:

add a version with wire IDs

deprecate this function

document what will happen (in the future methods) with wires on different planes

Definition at line 783 of file GeometryCore.cxx.

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

Referenced by vertex::FeatureVertexFinderAna::analyze(), shwf::ShowerReco::beginJob(), shwf::ShowerReco::beginRun(), mvapid::MVAAlg::CalcSegmentdEdxDist(), trkf::BezierTrackerAlgorithm::CalculateGeometricalElements(), trkf::SeedFinderAlgorithm::CalculateGeometricalElements(), util::GeometryUtilities::CalculatePitch(), util::GeometryUtilities::CalculatePitchPolar(), evd::RecoBaseDrawer::Cluster2D(), cluster::ClusterMatchTQ::ClusterMatch(), lar_cluster3d::StandardHit3DBuilder::configure(), shower::HitPosition::ConvertWireTickToCm(), lar_pandora::LArPandoraInput::CreatePandoraHits2D(), lar_pandora::LArPandoraInput::CreatePandoraReadoutGaps(), trkf::Track3DKalmanSPS::dQdxCalc(), evd::RecoBaseDrawer::DrawTrack2D(), cluster::EndPointAlg::EndPoint(), cluster::HoughBaseAlg::FastTransform(), trkf::SpacePointAlg::fillComplexSpacePoint(), tca::FilldEdx(), trkf::SpacePointAlg::fillSpacePoint(), tca::FillWireHitRange(), apa::DisambigAlg::FindChanTimeEndPts(), lar_cluster3d::HoughSeedFinderAlg::findHoughClusters(), trkf::BezierTrack::GetCalorimetryObject(), lar_pandora::LArPandoraInput::GetMips(), calo::Calorimetry::GetPitch(), trkf::BezierTrack::GetTrackPitch(), shower::TCShowerAlg::goodHit(), evd::RecoBaseDrawer::Hit2D(), shower::EMShowerAlg::HitPosition(), apa::APAGeometryAlg::Init(), lar_cluster3d::PCASeedFinderAlg::LineFit2DHits(), lar_cluster3d::HoughSeedFinderAlg::LineFit2DHits(), lar_pandora::LArPandoraGeometry::LoadGeometry(), shower::EMShowerAlg::MakeShower(), trkf::SpacePointAlg::makeSpacePoints(), util::GeometryUtilities::PitchInView(), lar_pandora::LArPandoraTrackCreation::produce(), cluster::DBCluster3D::produce(), trkf::SpacePts::produce(), trkf::Track3Dreco::produce(), calo::Calorimetry::produce(), vertex::VertexFinder2D::produce(), util::GeometryUtilities::Reconfigure(), trkf::KHitWireLine::subpredict(), trkf::KHitWireX::subpredict(), cluster::HoughBaseAlg::Transform(), trkf::TrackLineFitAlg::TrkLineFit(), evd::details::ADCCorrectorClass::update(), and trkf::VertexFitAlg::VertexFit().

   {
     return Plane(planeid).WirePitch();
   }

geo::Length_t geo::GeometryCore::WirePitch	(	unsigned int	plane = `0`,
		unsigned int	tpc = `0`,
		unsigned int	cstat = `0`
	)		const

inline

Returns the distance between two consecutive wires.

Parameters

p	plane number within the TPC
tpc	tpc number within the cryostat
cstat	cryostat number

Returns: the distance between the two wires

Note: The current geometry assumptions imply that wire pitch is constant between all wires on the same wire plane. This is an assumption non-trivial to remove.

Todo:

add a version with wire IDs

deprecate this function

document what will happen (in the future methods) with wires on different planes

Definition at line 3714 of file GeometryCore.h.

3717 { return WirePitch(geo::PlaneID(cstat, tpc, plane)); }