#include <AuxDetChannelMapAlg.h>

Public Member Functions
virtual	~AuxDetChannelMapAlg ()=default

virtual void	Initialize (AuxDetGeometryData_t &geodata)=0

virtual void	Uninitialize ()=0

virtual size_t	NearestAuxDet (const double *point, std::vector< geo::AuxDetGeo > const &auxDets, double tolerance=0) const

virtual size_t	NearestSensitiveAuxDet (const double *point, std::vector< geo::AuxDetGeo > const &auxDets, size_t &ad, double tolerance=0) const

virtual size_t	ChannelToAuxDet (std::vector< geo::AuxDetGeo > const &auxDets, std::string const &detName, uint32_t const &channel) const

virtual std::pair< size_t, size_t >	ChannelToSensitiveAuxDet (std::vector< geo::AuxDetGeo > const &auxDets, std::string const &detName, uint32_t const &channel) const

virtual uint32_t	PositionToAuxDetChannel (double const worldLoc[3], std::vector< geo::AuxDetGeo > const &auxDets, size_t &ad, size_t &sv) const =0

virtual const TVector3	AuxDetChannelToPosition (uint32_t const &channel, std::string const &auxDetName, std::vector< geo::AuxDetGeo > const &auxDets) const =0

Protected Attributes
std::map< size_t, std::string >	fADGeoToName
	map the AuxDetGeo index to the name More...

std::map< std::string, size_t >	fNameToADGeo
	map the names to the AuxDetGeo index More...

std::map< size_t, std::vector< chanAndSV > >	fADGeoToChannelAndSV

Private Types
typedef std::pair< uint32_t, size_t >	chanAndSV

Detailed Description

Definition at line 28 of file AuxDetChannelMapAlg.h.

Member Typedef Documentation

typedef std::pair<uint32_t, size_t> geo::AuxDetChannelMapAlg::chanAndSV

private

Definition at line 30 of file AuxDetChannelMapAlg.h.

Constructor & Destructor Documentation

virtual geo::AuxDetChannelMapAlg::~AuxDetChannelMapAlg ( )

virtualdefault

Member Function Documentation

virtual const TVector3 geo::AuxDetChannelMapAlg::AuxDetChannelToPosition	(	uint32_t const &	channel,
		std::string const &	auxDetName,
		std::vector< geo::AuxDetGeo > const &	auxDets
	)		const

pure virtual

size_t geo::AuxDetChannelMapAlg::ChannelToAuxDet	(	std::vector< geo::AuxDetGeo > const &	auxDets,
		std::string const &	detName,
		uint32_t const &	channel
	)		const

virtual

Definition at line 94 of file AuxDetChannelMapAlg.cxx.

   {
     // loop over the map of AuxDet names to Geo object numbers to determine which auxdet
     // we have.  If no name in the map matches the provided string, throw an exception;
     // the list of AuxDetGeo passed as argument is ignored!
     // Note that fADGeoToName must have been updated by a derived class.
     for(auto itr : fADGeoToName)
       if( itr.second.compare(detName) == 0 ) return itr.first;
 
 
     throw cet::exception("Geometry") << "No AuxDetGeo matching name: " << detName;
 
     return UINT_MAX;
   }

std::pair< size_t, size_t > geo::AuxDetChannelMapAlg::ChannelToSensitiveAuxDet	(	std::vector< geo::AuxDetGeo > const &	auxDets,
		std::string const &	detName,
		uint32_t const &	channel
	)		const

virtual

Definition at line 114 of file AuxDetChannelMapAlg.cxx.

   {
     size_t adGeoIdx = this->ChannelToAuxDet(auxDets, detName, channel);
 
     // look for the index of the sensitive volume for the given channel
     if( fADGeoToChannelAndSV.count(adGeoIdx) > 0 ){
 
       auto itr = fADGeoToChannelAndSV.find(adGeoIdx);
 
       // get the vector of channels to AuxDetSensitiveGeo index
       if( channel < itr->second.size() )
         return std::make_pair(adGeoIdx, itr->second[channel].second);
 
       throw cet::exception("Geometry") << "Given AuxDetSensitive channel, " << channel
                                        << ", cannot be found in vector associated to AuxDetGeo index: "
                                        << adGeoIdx << ". Vector has size " << itr->second.size();
     }
 
     throw cet::exception("Geometry") << "Given AuxDetGeo with index " << adGeoIdx
                                      << " does not correspond to any vector of sensitive volumes";
 
     return std::make_pair(adGeoIdx, UINT_MAX);
   }

virtual void geo::AuxDetChannelMapAlg::Initialize ( AuxDetGeometryData_t & geodata )

pure virtual

size_t geo::AuxDetChannelMapAlg::NearestAuxDet	(	const double *	point,
		std::vector< geo::AuxDetGeo > const &	auxDets,
		double	tolerance = `0`
	)		const

virtual

Definition at line 20 of file AuxDetChannelMapAlg.cxx.

   {
     double HalfCenterWidth = 0.;
     double localPoint[3] = {0.};
 
     for(size_t a = 0; a < auxDets.size(); ++a) {
 
       auxDets[a].WorldToLocal(point, localPoint);
 
       HalfCenterWidth = 0.5 * (auxDets[a].HalfWidth1() + auxDets[a].HalfWidth2());
 
       if(localPoint[2] >= (- auxDets[a].Length()/2 - tolerance)      &&
          localPoint[2] <= (  auxDets[a].Length()/2 + tolerance)      &&
          localPoint[1] >= (- auxDets[a].HalfHeight() - tolerance)    &&
          localPoint[1] <=   (auxDets[a].HalfHeight() + tolerance)    &&
          // if AuxDet a is a box, then HalfSmallWidth = HalfWidth
          localPoint[0] >= (- HalfCenterWidth + localPoint[2]*(HalfCenterWidth - auxDets[a].HalfWidth2())/(0.5 * auxDets[a].Length()) - tolerance) &&
          localPoint[0] <= (  HalfCenterWidth - localPoint[2]*(HalfCenterWidth - auxDets[a].HalfWidth2())/(0.5 * auxDets[a].Length()) + tolerance)
          ) return a;
 
     }// for loop over AudDet a
 
     // throw an exception because we couldn't find the sensitive volume
     throw cet::exception("AuxDetChannelMapAlg") << "Can't find AuxDet for position ("
                                              << point[0] << ","
                                              << point[1] << ","
                                              << point[2] << ")\n";
 
     return UINT_MAX;
 
   }

size_t geo::AuxDetChannelMapAlg::NearestSensitiveAuxDet	(	const double *	point,
		std::vector< geo::AuxDetGeo > const &	auxDets,
		size_t &	ad,
		double	tolerance = `0`
	)		const

virtual

Definition at line 54 of file AuxDetChannelMapAlg.cxx.

   {
     double HalfCenterWidth = 0.;
     double localPoint[3] = {0.};
 
     ad = this->NearestAuxDet(point, auxDets, tolerance);
 
     geo::AuxDetGeo const& adg = auxDets[ad];
 
     for(size_t a = 0; a < adg.NSensitiveVolume(); ++a) {
 
       geo::AuxDetSensitiveGeo const& adsg = adg.SensitiveVolume(a);
       adsg.WorldToLocal(point, localPoint);
 
       HalfCenterWidth = 0.5 * (adsg.HalfWidth1() + adsg.HalfWidth2());
 
 
       if(localPoint[2] >= - (adsg.Length()/2 + tolerance)      &&
          localPoint[2] <=   adsg.Length()/2 + tolerance      &&
          localPoint[1] >= - (adsg.HalfHeight() + tolerance)    &&
          localPoint[1] <=   adsg.HalfHeight() + tolerance    &&
          // if AuxDet a is a box, then HalfSmallWidth = HalfWidth
          localPoint[0] >= - HalfCenterWidth + localPoint[2]*(HalfCenterWidth - adsg.HalfWidth2())/(0.5 * adsg.Length()) - tolerance &&
          localPoint[0] <=   HalfCenterWidth - localPoint[2]*(HalfCenterWidth - adsg.HalfWidth2())/(0.5 * adsg.Length()) + tolerance
          ) return a;
     }// for loop over AuxDetSensitive a
 
     // throw an exception because we couldn't find the sensitive volume
     throw cet::exception("Geometry") << "Can't find AuxDetSensitive for position ("
                                      << point[0] << ","
                                      << point[1] << ","
                                      << point[2] << ")\n";
 
     return UINT_MAX;
   }

virtual uint32_t geo::AuxDetChannelMapAlg::PositionToAuxDetChannel	(	double const	worldLoc[3],
		std::vector< geo::AuxDetGeo > const &	auxDets,
		size_t &	ad,
		size_t &	sv
	)		const

pure virtual

virtual void geo::AuxDetChannelMapAlg::Uninitialize ( )

pure virtual

Member Data Documentation

std::map<size_t, std::vector<chanAndSV> > geo::AuxDetChannelMapAlg::fADGeoToChannelAndSV

protected

map the AuxDetGeo index to a vector of pairs corresponding to the channel and AuxDetSensitiveGeo index

Definition at line 71 of file AuxDetChannelMapAlg.h.

std::map<size_t, std::string> geo::AuxDetChannelMapAlg::fADGeoToName

protected

map the AuxDetGeo index to the name

Definition at line 69 of file AuxDetChannelMapAlg.h.

std::map<std::string, size_t> geo::AuxDetChannelMapAlg::fNameToADGeo

protected

map the names to the AuxDetGeo index

Definition at line 70 of file AuxDetChannelMapAlg.h.

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

larcorealg/larcorealg/Geometry/AuxDetChannelMapAlg.h
larcorealg/larcorealg/Geometry/AuxDetChannelMapAlg.cxx

Public Member Functions

Protected Attributes

Private Types

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation