sbnd::TPCGeoUtil Namespace Reference

Functions
int	DetectedInTPC (std::vector< art::Ptr< recob::Hit >> hits)
std::pair< double, double >	XLimitsFromHits (const geo::GeometryCore *GeometryService, std::vector< art::Ptr< recob::Hit >> hits)
int	DriftDirectionFromHits (const geo::GeometryCore *GeometryService, std::vector< art::Ptr< recob::Hit >> hits)
bool	InsideTPC (geo::Point_t point, const geo::TPCGeo &tpc, double buffer)

Function Documentation

int sbnd::TPCGeoUtil::DetectedInTPC

(

std::vector< art::Ptr< recob::Hit >>

hits

)

Definition at line 5 of file sbndcode/sbndcode/CRT/CRTUtils/TPCGeoUtil.cc.

                                                   {
  // Return tpc of hit collection or -1 if in multiple
  if(hits.size() == 0) return -1;
  int tpc = hits[0]->WireID().TPC;
  for(size_t i = 0; i < hits.size(); i++){
    if((int)hits[i]->WireID().TPC != tpc) return -1;
  }
  return tpc;
}

int sbnd::TPCGeoUtil::DriftDirectionFromHits	(	const geo::GeometryCore *	GeometryService,
		std::vector< art::Ptr< recob::Hit >>	hits
	)

Definition at line 28 of file sbndcode/sbndcode/CRT/CRTUtils/TPCGeoUtil.cc.

                                                                                                    {
  // If there are no hits then return 0
  if(hits.size() == 0) return 0;
  
  // If the track is stitched (in multiple TPCs) return 0
  if(DetectedInTPC(hits) == -1) return 0;

  // Work out the drift direction
  geo::TPCID tpcID = hits[0]->WireID().asTPCID();
  const geo::TPCGeo& tpcGeo = GeometryService->GetElement(tpcID);
  double driftDirection = tpcGeo.DetectDriftDirection();
  if(std::abs(driftDirection) != 1) driftDirection = 0;
  return driftDirection;
}

bool sbnd::TPCGeoUtil::InsideTPC	(	geo::Point_t	point,
		const geo::TPCGeo &	tpc,
		double	buffer
	)

Definition at line 44 of file sbndcode/sbndcode/CRT/CRTUtils/TPCGeoUtil.cc.

                                                                     {
  if(point.X() < (tpc.MinX()-buffer) || point.X() > (tpc.MaxX()+buffer)
      || point.Y() < (tpc.MinY()-buffer) || point.Y() > (tpc.MaxY()+buffer)
      || point.Z() < (tpc.MinZ()-buffer) || point.Z() > (tpc.MaxZ()+buffer)) return false;
  return true;
}

std::pair< double, double > sbnd::TPCGeoUtil::XLimitsFromHits	(	const geo::GeometryCore *	GeometryService,
		std::vector< art::Ptr< recob::Hit >>	hits
	)

Definition at line 15 of file sbndcode/sbndcode/CRT/CRTUtils/TPCGeoUtil.cc.

                                                                                                                 {
  // If there are no hits then return 0
  if(hits.size() == 0) return std::make_pair(0, 0);
  
  // If the track is stitched (in multiple TPCs) return 0
  if(DetectedInTPC(hits) == -1) return std::make_pair(0, 0);

  // Work out the drift direction
  geo::TPCID tpcID = hits[0]->WireID().asTPCID();
  const geo::TPCGeo& tpcGeo = GeometryService->GetElement(tpcID);
  return std::make_pair(tpcGeo.MinX(), tpcGeo.MaxX());
}