geo::ChannelMapCRMAlg Class Reference

#include <ChannelMapCRMAlg.h>

Inheritance diagram for geo::ChannelMapCRMAlg:

Public Member Functions
	ChannelMapCRMAlg (fhicl::ParameterSet const &p)
void	Initialize (GeometryData_t const &geodata) override
	Geometry initialisation. More...
void	Uninitialize () override
	Deconfiguration: prepare for a following call of Initialize() More...
std::vector< WireID >	ChannelToWire (raw::ChannelID_t channel) const override
unsigned int	Nchannels () const override
	Returns the total number of channels present (not necessarily contiguous) More...
virtual unsigned int	Nchannels (readout::ROPID const &ropid) const override
	Returns the number of channels in the specified ROP. More...
virtual SigType_t	SignalTypeForChannelImpl (raw::ChannelID_t const channel) const override
	Return the signal type of the specified channel. More...
virtual std::set< PlaneID > const &	PlaneIDs () const override
	Returns a list of the plane IDs in the whole detector. More...
virtual geo::GeoObjectSorter const &	Sorter () const override
	Returns the object to sort geometry with. More...
virtual double	WireCoordinate (double YPos, double ZPos, geo::PlaneID const &planeID) const override
	Returns the index of the wire nearest to the specified position. More...
virtual double	WireCoordinate (double YPos, double ZPos, unsigned int PlaneNo, unsigned int TPCNo, unsigned int cstat) const override
	Returns the index of the wire nearest to the specified position. More...
virtual WireID	NearestWireID (const TVector3 &worldPos, geo::PlaneID const &planeID) const override
	Returns the ID of the wire nearest to the specified position. More...
virtual WireID	NearestWireID (const TVector3 &worldPos, unsigned int PlaneNo, unsigned int TPCNo, unsigned int cstat) const override
	Returns the ID of the wire nearest to the specified position. More...
virtual raw::ChannelID_t	PlaneWireToChannel (geo::WireID const &wireID) const override
	Returns the channel ID a wire is connected to. More...
virtual raw::ChannelID_t	PlaneWireToChannel (unsigned int plane, unsigned int wire, unsigned int tpc, unsigned int cstat) const override
	Returns the channel ID a wire is connected to. More...
TPC set mapping
virtual unsigned int	NTPCsets (readout::CryostatID const &cryoid) const override
	Returns the total number of TPC sets in the specified cryostat. More...
virtual unsigned int	MaxTPCsets () const override
	Returns the largest number of TPC sets any cryostat in the detector has. More...
virtual bool	HasTPCset (readout::TPCsetID const &tpcsetid) const override
virtual readout::TPCsetID	TPCtoTPCset (geo::TPCID const &tpcid) const override
	Returns the ID of the TPC set tpcid belongs to. More...
virtual std::vector< geo::TPCID >	TPCsetToTPCs (readout::TPCsetID const &tpcsetid) const override
	Returns a list of ID of TPCs belonging to the specified TPC set. More...
virtual geo::TPCID	FirstTPCinTPCset (readout::TPCsetID const &tpcsetid) const override
	Returns the ID of the first TPC belonging to the specified TPC set. More...
Readout plane mapping
virtual unsigned int	NROPs (readout::TPCsetID const &tpcsetid) const override
	Returns the total number of ROP in the specified TPC set. More...
virtual unsigned int	MaxROPs () const override
	Returns the largest number of ROPs a TPC set in the detector has. More...
virtual bool	HasROP (readout::ROPID const &ropid) const override
virtual readout::ROPID	WirePlaneToROP (geo::PlaneID const &planeid) const override
	Returns the ID of the ROP planeid belongs to. More...
virtual std::vector< geo::PlaneID >	ROPtoWirePlanes (readout::ROPID const &ropid) const override
	Returns a list of ID of planes belonging to the specified ROP. More...
virtual geo::PlaneID	FirstWirePlaneInROP (readout::ROPID const &ropid) const override
	Returns the ID of the first plane belonging to the specified ROP. More...
virtual std::vector< geo::TPCID >	ROPtoTPCs (readout::ROPID const &ropid) const override
	Returns a list of ID of TPCs the specified ROP spans. More...
virtual readout::ROPID	ChannelToROP (raw::ChannelID_t channel) const override
virtual raw::ChannelID_t	FirstChannelInROP (readout::ROPID const &ropid) const override
	Returns the ID of the first channel in the specified readout plane. More...
Public Member Functions inherited from geo::ChannelMapAlg
virtual	~ChannelMapAlg ()=default
	Virtual destructor. More...
virtual unsigned int	NOpChannels (unsigned int NOpDets) const
	Returns the number of optical channels contained in some detectors. More...
virtual unsigned int	MaxOpChannel (unsigned int NOpDets) const
	Returns the number of optical channels contained in some detectors. More...
virtual unsigned int	NOpHardwareChannels (unsigned int opDet) const
	Returns the number of channels in the specified optical detectors. More...
virtual bool	IsValidOpChannel (unsigned int opChannel, unsigned int NOpDets) const
	Returns whether the ID identifies a valid optical detector channel. More...
virtual unsigned int	OpChannel (unsigned int detNum, unsigned int hwchannel=0) const
	Returns the channel ID of the specified hardware channel. More...
virtual unsigned int	OpDetFromOpChannel (unsigned int opChannel) const
	Returns the optical detector the specified optical channel belongs. More...
virtual unsigned int	HardwareChannelFromOpChannel (unsigned int opChannel) const
	Returns the hardware channel number of specified optical channel. More...
unsigned int	NearestWire (const TVector3 &worldPos, geo::PlaneID const &planeID) const
	Returns the index of the wire nearest to the specified position. More...
unsigned int	NearestWire (const TVector3 &worldPos, unsigned int PlaneNo, unsigned int TPCNo, unsigned int cstat) const
	Returns the index of the wire nearest to the specified position. More...
virtual size_t	NearestAuxDet (const double *point, std::vector< geo::AuxDetGeo > const &auxDets, double tolerance=0) const
	Returns the auxiliary detector closest to the specified point. More...
virtual size_t	NearestSensitiveAuxDet (const double *point, std::vector< geo::AuxDetGeo > const &auxDets, double tolerance=0) const
	Returns sensitive auxiliary detector closest to specified point. More...
virtual size_t	ChannelToAuxDet (std::vector< geo::AuxDetGeo > const &auxDets, std::string const &detName, uint32_t const &channel) const
	Returns the index of the detector containing the specified channel. More...
virtual std::pair< size_t, size_t >	ChannelToSensitiveAuxDet (std::vector< geo::AuxDetGeo > const &auxDets, std::string const &detName, uint32_t const &channel) const
	Returns the index of the sensitive detector containing the channel. More...
const std::vector< std::vector< std::vector< raw::ChannelID_t > > >	FirstChannelInNextPlane () const
	Retrieve the private fFirstChannelInNextPlane vector for testing. More...
const std::vector< std::vector< std::vector< raw::ChannelID_t > > >	FirstChannelInThisPlane () const
	Retrieve the private fFirstChannelInThisPlane vector for testing. More...
virtual bool	HasChannel (raw::ChannelID_t channel) const
	Returns whether the specified channel is valid This default implementation assumes all channels up to Nchannels() valid. More...
geo::SigType_t	SignalTypeForChannel (raw::ChannelID_t const channel) const
	Return the signal type of the specified channel. More...
geo::SigType_t	SignalTypeForROPID (readout::ROPID const &ropid) const
	Return the signal type on the specified readout plane. More...

Private Member Functions
unsigned int	WireCount (geo::PlaneID const &id) const
	Retrieved the wire cound for the specified plane ID. More...
unsigned int	MaxTPCs () const
	Returns the largest number of TPCs in a single cryostat. More...

Static Private Member Functions
static readout::TPCsetID	ConvertTPCtoTPCset (geo::TPCID const &tpcid)
	Converts a TPC ID into a TPC set ID using the same numerical indices. More...
static geo::TPCID	ConvertTPCsetToTPC (readout::TPCsetID const &tpcsetid)
	Converts a TPC set ID into a TPC ID using the same numerical indices. More...
static readout::ROPID	ConvertWirePlaneToROP (geo::PlaneID const &planeid)
	Converts a ROP ID into a wire plane ID using the same numerical indices. More...
static geo::PlaneID	ConvertROPtoWirePlane (readout::ROPID const &ropid)
	Converts a wire plane ID into a ROP ID using the same numerical indices. More...

Private Attributes
unsigned int	fNcryostat
	number of cryostats in the detector More...
unsigned int	fNchannels
	number of channels in the detector More...
raw::ChannelID_t	fTopChannel
	book keeping highest channel # More...
std::vector< unsigned int >	fNTPC
	number of TPCs in each cryostat More...
std::set< View_t >	fViews
	vector of the views present in the detector More...
std::set< PlaneID >	fPlaneIDs
	vector of the PlaneIDs present in the detector More...
PlaneInfoMap_t< float >	fFirstWireProj
PlaneInfoMap_t< float >	fOrthVectorsY
	Unit vectors orthogonal to wires in. More...
PlaneInfoMap_t< float >	fOrthVectorsZ
PlaneInfoMap_t< float >	fWireCounts
TPCInfoMap_t< unsigned int >	fNPlanes
PlaneInfoMap_t< unsigned int >	fPlaneBaselines
PlaneInfoMap_t< unsigned int >	fWiresPerPlane
geo::GeoObjectSorterCRM	fSorter
	class to sort geo objects More...

Additional Inherited Members
Protected Types inherited from geo::ChannelMapAlg
template<typename T >
using	TPCInfoMap_t = std::vector< std::vector< T >>
	Data type for per-TPC information. More...
template<typename T >
using	PlaneInfoMap_t = TPCInfoMap_t< std::vector< T >>
	Data type for per-plane information. More...
Protected Member Functions inherited from geo::ChannelMapAlg
template<typename T >
T const &	AccessElement (PlaneInfoMap_t< T > const &map, geo::PlaneID const &id) const
	Returns the specified element of the plane map. More...
template<typename T >
size_t	AccessElementSize (PlaneInfoMap_t< T > const &map, geo::TPCID const &id) const
	Returns the number of elements in the specified TPC of the plane map. More...
template<typename T >
T const *	GetElementPtr (PlaneInfoMap_t< T > const &map, geo::PlaneID const &id) const
	Returns a pointer to the specified element, or nullptr if invalid. More...
template<typename T >
T const &	AccessElement (TPCInfoMap_t< T > const &map, geo::TPCID const &id) const
	Returns the specified element of the TPC map. More...
template<typename T >
size_t	AccessElementSize (TPCInfoMap_t< T > const &map, geo::CryostatID const &id) const
	Returns the number of elements in the specified cryostat of the TPC map. More...
template<typename T >
bool	isValidElement (TPCInfoMap_t< T > const &map, geo::CryostatID const &id) const
	Returns whether the ID specifies a valid entry. More...
template<typename T >
bool	isValidElement (TPCInfoMap_t< T > const &map, geo::TPCID const &id) const
template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, geo::CryostatID const &id) const
	Returns whether the ID specifies a valid entry. More...
template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, geo::TPCID const &id) const
template<typename T >
bool	isValidElement (PlaneInfoMap_t< T > const &map, geo::PlaneID const &id) const
virtual geo::SigType_t	SignalTypeForROPIDImpl (readout::ROPID const &ropid) const
	Return the signal type on the specified readout plane. More...
Protected Attributes inherited from geo::ChannelMapAlg
PlaneInfoMap_t< raw::ChannelID_t >	fFirstChannelInThisPlane
PlaneInfoMap_t< raw::ChannelID_t >	fFirstChannelInNextPlane
std::map< std::string, size_t >	fADNameToGeo
	map the names of the dets to the AuxDetGeo objects More...
std::map< size_t, std::vector< size_t > >	fADChannelToSensitiveGeo

Detailed Description

Definition at line 23 of file ChannelMapCRMAlg.h.

Constructor & Destructor Documentation

geo::ChannelMapCRMAlg::ChannelMapCRMAlg

(

fhicl::ParameterSet const &

p

)

Definition at line 23 of file ChannelMapCRMAlg.cxx.

    : fSorter(geo::GeoObjectSorterCRM(p))
  {
  }

Member Function Documentation

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

overridevirtual

Returns the ID of the ROP the channel belongs to

Exceptions

cet::exception

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

Implements geo::ChannelMapAlg.

Definition at line 525 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

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

Exceptions

cet::exception

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

Implements geo::ChannelMapAlg.

Definition at line 162 of file ChannelMapCRMAlg.cxx.

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

    geo::WireID CodeWire(cstat, tpc, plane, wire);
    
    AllSegments.push_back(CodeWire);

    return AllSegments;
  }

geo::PlaneID geo::ChannelMapCRMAlg::ConvertROPtoWirePlane ( readout::ROPID const &  ropid )

staticprivate

Converts a wire plane ID into a ROP ID using the same numerical indices.

Definition at line 600 of file ChannelMapCRMAlg.cxx.

  {
    if (!ropid.isValid) return {};
    return {
      (geo::CryostatID::CryostatID_t) ropid.Cryostat,
      (geo::TPCID::TPCID_t) ropid.TPCset,
      (geo::PlaneID::PlaneID_t) ropid.ROP
      };
  } // ChannelMapCRMAlg::ConvertROPtoWirePlane()

geo::TPCID geo::ChannelMapCRMAlg::ConvertTPCsetToTPC ( readout::TPCsetID const &  tpcsetid )

staticprivate

Converts a TPC set ID into a TPC ID using the same numerical indices.

Definition at line 574 of file ChannelMapCRMAlg.cxx.

  {
    if (!tpcsetid.isValid) return {};
    return {
      (geo::CryostatID::CryostatID_t) tpcsetid.Cryostat,
      (geo::TPCID::TPCID_t) tpcsetid.TPCset
      };
  } // ChannelMapCRMAlg::ConvertTPCsetToTPC()

readout::TPCsetID geo::ChannelMapCRMAlg::ConvertTPCtoTPCset ( geo::TPCID const &  tpcid )

staticprivate

Converts a TPC ID into a TPC set ID using the same numerical indices.

Definition at line 562 of file ChannelMapCRMAlg.cxx.

  {
    if (!tpcid.isValid) return {}; // invalid ID, default-constructed
    return {
      (readout::CryostatID::CryostatID_t) tpcid.Cryostat,
      (readout::TPCsetID::TPCsetID_t) tpcid.TPC
      };
  } // ChannelMapCRMAlg::ConvertTPCtoTPCset()

readout::ROPID geo::ChannelMapCRMAlg::ConvertWirePlaneToROP ( geo::PlaneID const &  planeid )

staticprivate

Converts a ROP ID into a wire plane ID using the same numerical indices.

Definition at line 586 of file ChannelMapCRMAlg.cxx.

  {
    if (!planeid.isValid) return {}; // invalid ID, default-constructed
    return {
      (readout::CryostatID::CryostatID_t) planeid.Cryostat,
      (readout::TPCsetID::TPCsetID_t) planeid.TPC,
      (readout::ROPID::ROPID_t) planeid.Plane
      };
    
  } // ChannelMapCRMAlg::ConvertWirePlaneToROP()

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

overridevirtual

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

Parameters

ropid

ID of the readout plane

Returns

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

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

Implements geo::ChannelMapAlg.

Definition at line 552 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 423 of file ChannelMapCRMAlg.cxx.

  {
    return ConvertTPCsetToTPC(tpcsetid);
  } // ChannelMapCRMAlg::FirstTPCinTPCset()

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 495 of file ChannelMapCRMAlg.cxx.

  {
    return ConvertROPtoWirePlane(ropid);
  } // ChannelMapCRMAlg::FirstWirePlaneInROP()

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

overridevirtual

Returns whether we have the specified ROP

Returns

whether the readout plane is valid and exists

Implements geo::ChannelMapAlg.

Definition at line 467 of file ChannelMapCRMAlg.cxx.

                                                               {
    return ropid.ROP < NROPs(ropid);
  } // ChannelMapCRMAlg::HasROP()

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

overridevirtual

Returns whether we have the specified TPC set

Returns

whether the TPC set is valid and exists

Implements geo::ChannelMapAlg.

Definition at line 386 of file ChannelMapCRMAlg.cxx.

                                                                        {
    return tpcsetid.TPCset < NTPCsets(tpcsetid);
  } // ChannelMapCRMAlg::HasTPCset()

void geo::ChannelMapCRMAlg::Initialize ( GeometryData_t const &  geodata )

overridevirtual

Geometry initialisation.

Implements geo::ChannelMapAlg.

Definition at line 29 of file ChannelMapCRMAlg.cxx.

  {
    // start over:
    Uninitialize();
    
    std::vector<geo::CryostatGeo> const& cgeo = geodata.cryostats;
    
    fNcryostat = cgeo.size();
    
    mf::LogInfo("ChannelMapCRMAlg") << "Initializing CRM ChannelMap...";

    fNTPC.resize(fNcryostat);
    fWireCounts.resize(fNcryostat);
    fNPlanes.resize(fNcryostat);
    fFirstWireProj.resize(fNcryostat);
    fOrthVectorsY.resize(fNcryostat);
    fOrthVectorsZ.resize(fNcryostat);
    fPlaneBaselines.resize(fNcryostat);
    fWiresPerPlane.resize(fNcryostat);
    fFirstChannelInNextPlane.resize(fNcryostat);
    fFirstChannelInThisPlane.resize(fNcryostat);
    fPlaneIDs.clear();
    fTopChannel = 0;

    int RunningTotal = 0;

    for(unsigned int cs = 0; cs != fNcryostat; ++cs){
      
      fNTPC[cs] = cgeo[cs].NTPC();
      
      // Size up all the vectors 
      fWireCounts[cs]             .resize(fNTPC[cs]);
      fFirstWireProj[cs]          .resize(fNTPC[cs]);
      fOrthVectorsY[cs]           .resize(fNTPC[cs]);
      fOrthVectorsZ[cs]           .resize(fNTPC[cs]);
      fPlaneBaselines[cs]         .resize(fNTPC[cs]);
      fWiresPerPlane[cs]          .resize(fNTPC[cs]);
      fNPlanes[cs]                .resize(fNTPC[cs]);
      fFirstChannelInThisPlane[cs].resize(fNTPC[cs]);
      fFirstChannelInNextPlane[cs].resize(fNTPC[cs]);

      for(unsigned int TPCCount = 0; TPCCount != fNTPC[cs]; ++TPCCount)
        {
          unsigned int PlanesThisTPC = cgeo[cs].TPC(TPCCount).Nplanes();
          fWireCounts[cs][TPCCount]   .resize(PlanesThisTPC);
          fFirstWireProj[cs][TPCCount].resize(PlanesThisTPC);
          fOrthVectorsY[cs][TPCCount] .resize(PlanesThisTPC);
          fOrthVectorsZ[cs][TPCCount] .resize(PlanesThisTPC);
          fNPlanes[cs][TPCCount]=PlanesThisTPC;
          
          for(unsigned int PlaneCount = 0; PlaneCount != PlanesThisTPC; ++PlaneCount){
            
            if(PlaneCount >= 2) //should only find two views for dual-phase
              throw cet::exception("Geometry") <<"CANNOT HAVE more than two readout planes for dual-phase" 
                                               << "\n";
            /*
            if(cgeo[cs].TPC(TPCCount).Plane(PlaneCount).View() == geo::kU)
              std::cout<<" View is kU "<<std::endl;
            else if(cgeo[cs].TPC(TPCCount).Plane(PlaneCount).View() == geo::kV)
              std::cout<<" View is kV "<<std::endl;
            else if(cgeo[cs].TPC(TPCCount).Plane(PlaneCount).View() == geo::kZ)
              std::cout<<" View is kZ "<<std::endl;

            if(cgeo[cs].TPC(TPCCount).Plane(PlaneCount).SignalType() == geo::kInduction)
              std::cout<<" View is kInduction "<<std::endl;
            if(cgeo[cs].TPC(TPCCount).Plane(PlaneCount).SignalType() == geo::kCollection)
              std::cout<<" View is kCollection "<<std::endl;
            */

            fPlaneIDs.emplace(PlaneID(cs, TPCCount, PlaneCount));
            double ThisWirePitch = cgeo[cs].TPC(TPCCount).WirePitch(PlaneCount);
            fWireCounts[cs][TPCCount][PlaneCount] = cgeo[cs].TPC(TPCCount).Plane(PlaneCount).Nwires();
            
            double  WireCentre1[3] = {0.,0.,0.};
            double  WireCentre2[3] = {0.,0.,0.};
          
            const geo::WireGeo& firstWire = cgeo[cs].TPC(TPCCount).Plane(PlaneCount).Wire(0);
            const double sth = firstWire.SinThetaZ(), cth = firstWire.CosThetaZ();

            firstWire.GetCenter(WireCentre1,0);
            cgeo[cs].TPC(TPCCount).Plane(PlaneCount).Wire(1).GetCenter(WireCentre2,0);
            
            // figure out if we need to flip the orthogonal vector 
            // (should point from wire n -> n+1)
            double OrthY = cth, OrthZ = -sth;
            if(((WireCentre2[1] - WireCentre1[1])*OrthY 
                + (WireCentre2[2] - WireCentre1[2])*OrthZ) < 0){
              OrthZ *= -1;
              OrthY *= -1;
            }
          
          // Overall we are trying to build an expression that looks like
          //   int NearestWireNumber = round((worldPos.OrthVector - FirstWire.OrthVector)/WirePitch);     
          // That runs as fast as humanly possible.
          // We predivide everything by the wire pitch so we don't do this in the loop.
          //
          // Putting this together into the useful constants we will use later per plane and tpc:
          fOrthVectorsY[cs][TPCCount][PlaneCount] = OrthY / ThisWirePitch;
          fOrthVectorsZ[cs][TPCCount][PlaneCount] = OrthZ / ThisWirePitch;
          
          fFirstWireProj[cs][TPCCount][PlaneCount]  = WireCentre1[1]*OrthY + WireCentre1[2]*OrthZ;
          fFirstWireProj[cs][TPCCount][PlaneCount] /= ThisWirePitch;
          
          // now to count up wires in each plane and get first channel in each plane
          int WiresThisPlane = cgeo[cs].TPC(TPCCount).Plane(PlaneCount).Nwires();        
          fWiresPerPlane[cs] .at(TPCCount).push_back(WiresThisPlane);
          fPlaneBaselines[cs].at(TPCCount).push_back(RunningTotal);
          
          RunningTotal += WiresThisPlane;

          fFirstChannelInThisPlane[cs].at(TPCCount).push_back(fTopChannel);
          fTopChannel += WiresThisPlane;
          fFirstChannelInNextPlane[cs].at(TPCCount).push_back(fTopChannel);

        }// end loop over planes
      }// end loop over TPCs
    }// end loop over cryostats

    // calculate the total number of channels in the detector
    fNchannels = fTopChannel;

    MF_LOG_DEBUG("ChannelMapCRM") << "# of channels is " << fNchannels;


    return;
  }

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 455 of file ChannelMapCRMAlg.cxx.

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

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

private

Returns the largest number of TPCs in a single cryostat.

Definition at line 430 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 378 of file ChannelMapCRMAlg.cxx.

                                                  {
    return MaxTPCs();
  } // ChannelMapCRMAlg::MaxTPCsets()

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 205 of file ChannelMapCRMAlg.cxx.

  {
    return fNchannels;
  }

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

overridevirtual

Returns the number of channels in the specified ROP.

Returns

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

Todo:

Needs to be implemented

Implements geo::ChannelMapAlg.

Definition at line 211 of file ChannelMapCRMAlg.cxx.

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

WireID geo::ChannelMapCRMAlg::NearestWireID ( const TVector3 &  worldPos, geo::PlaneID const &  planeID  ) const

overridevirtual

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

Parameters

worldPos	position to be tested
planeID	plane containing the wire

Returns

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

Exceptions

InvalidWireIDError

the ID found is not present in the detector

See also

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

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

Reimplemented from geo::ChannelMapAlg.

Definition at line 233 of file ChannelMapCRMAlg.cxx.

  {

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

    return geo::WireID(planeID, (geo::WireID::WireID_t) NearestWireNumber);
  }

virtual WireID geo::ChannelMapCRMAlg::NearestWireID ( const TVector3 &  worldPos, unsigned int  PlaneNo, unsigned int  TPCNo, unsigned int  cstat  ) const

inlineoverridevirtual

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

Parameters

worldPos	position to be tested
PlaneNo	number of plane containing the wire
TPCNo	number of TPC containing the wire
cstat	number of cryostat containing the wire

Returns

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

See also

NearestWireID(const TVector3&, geo::PlaneID const&)

Deprecated:

Use the version with geo::PlaneID instead

Implements geo::ChannelMapAlg.

Definition at line 54 of file ChannelMapCRMAlg.h.

      { return NearestWireID(worldPos, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

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

overridevirtual

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.

Implements geo::ChannelMapAlg.

Definition at line 446 of file ChannelMapCRMAlg.cxx.

  {
    if (!HasTPCset(tpcsetid)) return 0;
    return AccessElement(fNPlanes, FirstTPCinTPCset(tpcsetid));
  } // ChannelMapCRMAlg::NROPs()

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

overridevirtual

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

Parameters

cryoid

cryostat ID

Returns

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

Implements geo::ChannelMapAlg.

Definition at line 368 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 355 of file ChannelMapCRMAlg.cxx.

  {
    return fPlaneIDs;
  }

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

overridevirtual

Returns the channel ID a wire is connected to.

Parameters

wireID

ID of the wire

Returns

the ID of the channel

See also

PlaneWireToChannel(geo::WireID const&)

Behaviour on an invalid or not present wires is undefined.

Deprecated:

Use the version with geo::WireID

Reimplemented from geo::ChannelMapAlg.

Definition at line 278 of file ChannelMapCRMAlg.cxx.

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

  }

virtual raw::ChannelID_t geo::ChannelMapCRMAlg::PlaneWireToChannel ( unsigned int  plane, unsigned int  wire, unsigned int  tpc, unsigned int  cstat  ) const

inlineoverridevirtual

Returns the channel ID a wire is connected to.

Parameters

plane	number of plane
wire	number of wire
tpc	number of TPC
cstat	number of cryostat

Returns

the ID of the channel

See also

PlaneWireToChannel(geo::WireID const&)

Deprecated:

Use the version with geo::WireID

Implements geo::ChannelMapAlg.

Definition at line 64 of file ChannelMapCRMAlg.h.

      { return PlaneWireToChannel(geo::WireID(cstat, tpc, plane, wire)); }

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

overridevirtual

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

Parameters

ropid

ID of the readout plane

Returns

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

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

Implements geo::ChannelMapAlg.

Definition at line 513 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

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

Implements geo::ChannelMapAlg.

Definition at line 484 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

Return the signal type of the specified channel.

Parameters

channel

ID of the channel

Returns

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

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

Implements geo::ChannelMapAlg.

Definition at line 302 of file ChannelMapCRMAlg.cxx.

  {

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

    return sigt;
  }

virtual geo::GeoObjectSorter const& geo::ChannelMapCRMAlg::Sorter ( ) const

inlineoverridevirtual

Returns the object to sort geometry with.

Implements geo::ChannelMapAlg.

Definition at line 184 of file ChannelMapCRMAlg.h.

      { return fSorter; }

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

overridevirtual

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

Parameters

tpcsetid

ID of the TPC set to convert into TPC IDs

Returns

the list of TPCs, empty if TPC set is invalid

Note that the check is performed on the validity of the TPC set ID, that does not necessarily imply that the TPC set specified by the ID actually exists. Check the existence of the TPC set first (HasTPCset()). Behaviour on valid, non-existent TPC set IDs is undefined.

Implements geo::ChannelMapAlg.

Definition at line 412 of file ChannelMapCRMAlg.cxx.

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

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

overridevirtual

Returns the ID of the TPC set tpcid belongs to.

Implements geo::ChannelMapAlg.

Definition at line 394 of file ChannelMapCRMAlg.cxx.

  {
    return ConvertTPCtoTPCset(tpcid);
  } // ChannelMapCRMAlg::TPCtoTPCset()

void geo::ChannelMapCRMAlg::Uninitialize ( )

overridevirtual

Deconfiguration: prepare for a following call of Initialize()

Implements geo::ChannelMapAlg.

Definition at line 157 of file ChannelMapCRMAlg.cxx.

  {
  }

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

overridevirtual

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

Parameters

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

Returns

an index interpolation between the two nearest wires

See also

NearestWireID()

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

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

Reimplemented from geo::ChannelMapAlg.

Definition at line 220 of file ChannelMapCRMAlg.cxx.

  {
    // Returns the wire number corresponding to a (Y,Z) position in PlaneNo 
    // with float precision.
    // B. Baller August 2014
    return YPos*AccessElement(fOrthVectorsY, planeID) 
         + ZPos*AccessElement(fOrthVectorsZ, planeID)
         - AccessElement(fFirstWireProj, planeID);
  }

virtual double geo::ChannelMapCRMAlg::WireCoordinate ( double  YPos, double  ZPos, unsigned int  PlaneNo, unsigned int  TPCNo, unsigned int  cstat  ) const

inlineoverridevirtual

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

Parameters

YPos	y coordinate on the wire plane
ZPos	z coordinate on the wire plane
PlaneNo	number of plane
TPCNo	number of TPC
cstat	number of cryostat

Returns

an index interpolation between the two nearest wires

See also

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

Deprecated:

Use the version with geo::PlaneID instead

Implements geo::ChannelMapAlg.

Definition at line 44 of file ChannelMapCRMAlg.h.

      { return WireCoordinate(YPos, ZPos, geo::PlaneID(cstat, TPCNo, PlaneNo)); }

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

inlineprivate

Retrieved the wire cound for the specified plane ID.

Definition at line 214 of file ChannelMapCRMAlg.h.

    { return AccessElement(fWireCounts, id); }

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

overridevirtual

Returns the ID of the ROP planeid belongs to.

Implements geo::ChannelMapAlg.

Definition at line 475 of file ChannelMapCRMAlg.cxx.

  {
    return ConvertWirePlaneToROP(planeid);
  } // ChannelMapCRMAlg::WirePlaneToROP()

Member Data Documentation

PlaneInfoMap_t<float> geo::ChannelMapCRMAlg::fFirstWireProj

private

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

Definition at line 195 of file ChannelMapCRMAlg.h.

unsigned int geo::ChannelMapCRMAlg::fNchannels

private

number of channels in the detector

Definition at line 190 of file ChannelMapCRMAlg.h.

unsigned int geo::ChannelMapCRMAlg::fNcryostat

private

number of cryostats in the detector

Definition at line 189 of file ChannelMapCRMAlg.h.

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

private

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

Definition at line 203 of file ChannelMapCRMAlg.h.

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

private

number of TPCs in each cryostat

Definition at line 192 of file ChannelMapCRMAlg.h.

PlaneInfoMap_t<float> geo::ChannelMapCRMAlg::fOrthVectorsY

private

Unit vectors orthogonal to wires in.

Definition at line 197 of file ChannelMapCRMAlg.h.

PlaneInfoMap_t<float> geo::ChannelMapCRMAlg::fOrthVectorsZ

private

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

Definition at line 198 of file ChannelMapCRMAlg.h.

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

private

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

Definition at line 205 of file ChannelMapCRMAlg.h.

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

private

vector of the PlaneIDs present in the detector

Definition at line 194 of file ChannelMapCRMAlg.h.

geo::GeoObjectSorterCRM geo::ChannelMapCRMAlg::fSorter

private

class to sort geo objects

Definition at line 210 of file ChannelMapCRMAlg.h.

raw::ChannelID_t geo::ChannelMapCRMAlg::fTopChannel

private

book keeping highest channel #

Definition at line 191 of file ChannelMapCRMAlg.h.

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

private

vector of the views present in the detector

Definition at line 193 of file ChannelMapCRMAlg.h.

PlaneInfoMap_t<float> geo::ChannelMapCRMAlg::fWireCounts

private

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

Definition at line 201 of file ChannelMapCRMAlg.h.

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

private

The number of wires in this plane in the heirachy

Definition at line 208 of file ChannelMapCRMAlg.h.

---

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

• dunecore/dunecore/Geometry/ChannelMapCRMAlg.h
• dunecore/dunecore/Geometry/ChannelMapCRMAlg.cxx