WireSelector Class Reference

#include <WireSelector.h>

Classes
struct	WireInfo
class	WireSummary

Public Types
using	Index = unsigned int
using	IndexVector = std::vector< Index >
using	GeometryCore = geo::GeometryCore
using	View = geo::View_t
using	PlaneID = geo::PlaneID
using	PlaneIDVector = std::vector< PlaneID >
using	WireInfoVector = std::vector< WireInfo >
using	WireInfoMap = std::multimap< Index, const WireInfo * >

Public Member Functions
	WireSelector (Index icry=0)
const GeometryCore *	geometry () const
const IndexVector &	cryostats () const
View	view () const
double	wireAngle () const
double	wireAngleTolerance () const
double	driftMin () const
double	driftMax () const
const PlaneIDVector &	planeIDs () const
bool	haveData () const
const WireInfoVector &	data () const
const WireInfoMap &	dataMap () const
const WireSummary &	wireSummary () const
const GeometryCore *	geometry ()
void	selectCryostats (const IndexVector &a_icrys)
void	selectView (View view)
void	selectWireAngle (double wireAngle, double tol=0.001)
void	selectDrift (double dmin, double dmax=1.e20)
void	selectTpcSet (Index itps)
void	selectTpcSets (const IndexVector &itpss)
void	selectPlanes ()
const WireInfoVector &	fillData ()
const WireInfoMap &	fillDataMap ()
const WireSummary &	fillWireSummary ()
void	clearData ()

Private Attributes
const GeometryCore *	m_pgeo =nullptr
IndexVector	m_icrys
View	m_view =geo::kUnknown
double	m_wireAngle =999.
double	m_wireAngleTolerance =999.
double	m_driftMin = 0.0
double	m_driftMax = 1.e20
PlaneIDVector	m_pids
IndexVector	m_tpcSets
bool	m_haveData =false
WireInfoVector	m_data
WireInfoMap	m_datamap
WireSummary	m_wireSummary

Detailed Description

Definition at line 29 of file WireSelector.h.

Member Typedef Documentation

using WireSelector::GeometryCore = geo::GeometryCore

Definition at line 35 of file WireSelector.h.

using WireSelector::Index = unsigned int

Definition at line 33 of file WireSelector.h.

using WireSelector::IndexVector = std::vector<Index>

Definition at line 34 of file WireSelector.h.

using WireSelector::PlaneID = geo::PlaneID

Definition at line 37 of file WireSelector.h.

using WireSelector::PlaneIDVector = std::vector<PlaneID>

Definition at line 38 of file WireSelector.h.

using WireSelector::View = geo::View_t

Definition at line 36 of file WireSelector.h.

using WireSelector::WireInfoMap = std::multimap<Index, const WireInfo*>

Definition at line 98 of file WireSelector.h.

using WireSelector::WireInfoVector = std::vector<WireInfo>

Definition at line 97 of file WireSelector.h.

Constructor & Destructor Documentation

WireSelector::WireSelector ( Index  icry = 0 )

explicit

Definition at line 15 of file WireSelector.cxx.

                                     {
  selectCryostats({icry});
}

Member Function Documentation

void WireSelector::clearData

(

)

Definition at line 193 of file WireSelector.cxx.

                             {
  m_data.clear();
  m_datamap.clear();
  m_wireSummary.clear();
  m_haveData = false;
}

const IndexVector& WireSelector::cryostats ( ) const

inline

Definition at line 109 of file WireSelector.h.

{ return m_icrys; }

const WireInfoVector& WireSelector::data ( ) const

inline

Definition at line 127 of file WireSelector.h.

{ return m_data; }

const WireInfoMap& WireSelector::dataMap ( ) const

inline

Definition at line 128 of file WireSelector.h.

{ return m_datamap; }

double WireSelector::driftMax ( ) const

inline

Definition at line 120 of file WireSelector.h.

{ return m_driftMax; }

double WireSelector::driftMin ( ) const

inline

Definition at line 119 of file WireSelector.h.

{ return m_driftMin; }

const WireSelector::WireInfoVector & WireSelector::fillData

(

)

Definition at line 117 of file WireSelector.cxx.

                                                         {
  const string myname = "WireSelector::fillData: ";
  const double piOver2 = 0.5*acos(-1.0);
  if ( haveData() ) return m_data;
  for ( geo::PlaneID pid : planeIDs() ) {
    const geo::PlaneGeo& gpla = geometry()->Plane(pid);
    const geo::TPCGeo& gtpc = geometry()->TPC(pid);
    const geo::PlaneGeo& gplaLast = gtpc.LastPlane();  // Plane furthest from TPC center.
    double xLastPlane = gplaLast.MiddleWire().GetCenter<TVector3>().x();
    double xThisPlane = gpla.MiddleWire().GetCenter<TVector3>().x();
    double driftOffset = fabs(xThisPlane - xLastPlane);
    for ( Index iwir=0; iwir<gpla.Nwires(); ++iwir ) {
      geo::WireID wid(pid, iwir);
      const geo::WireGeo* pgwir = gpla.WirePtr(wid);
      double wireAngle = gpla.ThetaZ() - piOver2;
      TVector3 xyzWire = pgwir->GetCenter<TVector3>();
      xyzWire.RotateX(wireAngle);
      double driftSign = gpla.GetNormalDirection().x();
      if ( fabs(fabs(driftSign) - 1.0) > 0.001 ) {
        cout << myname << "ERROR: Plane normal is not along x." << endl;
        continue;
      }
      double driftDist = gtpc.DriftDistance() - driftOffset;
      m_data.emplace_back(xyzWire.x(), xyzWire.y(), xyzWire.z(),
                          driftSign*driftDist,
                          pgwir->Length(),
                          gpla.WirePitch(),
                          geometry()->PlaneWireToChannel(wid));
    }
  }
  m_haveData = true;
  return m_data;
}

const WireSelector::WireInfoMap & WireSelector::fillDataMap

(

)

Definition at line 153 of file WireSelector.cxx.

                                                         {
  if ( haveData() && m_datamap.size() == m_data.size() ) return m_datamap;
  for ( const WireInfo& dat : fillData() ) {
    m_datamap.emplace(dat.channel, &dat);
  }
  return m_datamap;
}

const WireSelector::WireSummary & WireSelector::fillWireSummary

(

)

Definition at line 163 of file WireSelector.cxx.

                                                             {
  if ( haveData() && m_wireSummary.size() == m_data.size() ) return m_wireSummary;
  WireSummary& ws = m_wireSummary;
  Index nwir = data().size();
  ws.xmin = 1000.0;
  ws.ymin = 1000.0;
  ws.zmin = 1000.0;
  ws.xmax = -1000.0;
  ws.ymax = -1000.0;
  ws.zmax = -1000.0;
  ws.xWire.resize(nwir);
  ws.xCathode.resize(nwir);
  ws.zWire.resize(nwir);
  for ( Index iwir=0; iwir<nwir; ++iwir ) {
    const WireInfo& dat = data()[iwir];
    if ( dat.x1() < ws.xmin ) ws.xmin = dat.x1();
    if ( dat.y1() < ws.ymin ) ws.ymin = dat.y1();
    if ( dat.z1() < ws.zmin ) ws.zmin = dat.z1();
    if ( dat.x2() > ws.xmax ) ws.xmax = dat.x2();
    if ( dat.y2() > ws.ymax ) ws.ymax = dat.y2();
    if ( dat.z2() > ws.zmax ) ws.zmax = dat.z2();
    ws.xWire[iwir] = dat.x;
    ws.xCathode[iwir] = dat.x + dat.driftMax;
    ws.zWire[iwir] = dat.z;
  }
  return ws;
}

const GeometryCore* WireSelector::geometry ( ) const

inline

Definition at line 106 of file WireSelector.h.

{ return m_pgeo; }

const WireSelector::GeometryCore * WireSelector::geometry

(

)

Definition at line 21 of file WireSelector.cxx.

                                                       {
  if ( m_pgeo == nullptr ) {
    art::ServiceHandle<geo::Geometry> hgeo;
    m_pgeo = hgeo.get();
  }
  return m_pgeo;
}

bool WireSelector::haveData ( ) const

inline

Definition at line 126 of file WireSelector.h.

{ return m_haveData; };

const PlaneIDVector& WireSelector::planeIDs ( ) const

inline

Definition at line 123 of file WireSelector.h.

{ return m_pids; }

void WireSelector::selectCryostats

(

const IndexVector &

a_icrys

)

Definition at line 31 of file WireSelector.cxx.

                                                             {
  m_icrys = a_icrys;
  selectPlanes();
}

void WireSelector::selectDrift	(	double	dmin,
		double	dmax = 1.e20
	)

Definition at line 53 of file WireSelector.cxx.

                                                       {
  m_driftMin = dmin;
  m_driftMax = dmax;
  selectPlanes();
}

void WireSelector::selectPlanes

(

)

Definition at line 75 of file WireSelector.cxx.

                                {
  const double pi = acos(-1.0);
  const double piOver2 = 0.5*pi;
  m_pids.clear();
  clearData();
  for ( Index icry : cryostats() ) {
    geo::CryostatID cid(icry);
    for ( geo::TPCID tid : geometry()->IterateTPCIDs(cid) ) {
      // Check TPC set.
      if ( m_tpcSets.size() ) {
        Index itps = geometry()->TPCtoTPCset(tid).TPCset;
        if ( find(m_tpcSets.begin(), m_tpcSets.end(), itps) == m_tpcSets.end() ) continue;
      }
      // Check drift.
      if ( driftMin() > 0.0 || driftMax() < 1.e20 ) {
        const geo::TPCGeo& gtpc = geometry()->TPC(tid);
        //double driftSize = gtpc.ActiveWidth();
        double driftSize = gtpc.DriftDistance(); // good for the last anode plane
        if ( driftSize < driftMin() ) continue;
        if ( driftSize >= driftMax() ) continue;
      }
      for ( geo::PlaneID pid : geometry()->IteratePlaneIDs(tid) ) {
        const geo::PlaneGeo& gpla = geometry()->Plane(pid);
        // Check view.
        if ( view() != geo::kUnknown && gpla.View() != view() ) continue;
        // Check wire angle.
        if ( wireAngleTolerance() < pi ) {
          double planeWireAngle = gpla.ThetaZ() - piOver2;
          double dang = planeWireAngle - wireAngle();
          while ( dang < -piOver2 ) dang += pi;
          while ( dang >  piOver2 ) dang -= pi;
          if ( fabs(dang) >= wireAngleTolerance() ) continue;
        }
        // Keep this plane.
        m_pids.push_back(pid);
      }
    }
  }
}

void WireSelector::selectTpcSet

(

Index

itps

)

Definition at line 61 of file WireSelector.cxx.

                                          {
  m_tpcSets.clear();
  m_tpcSets.push_back(itps);
  selectPlanes();
}

void WireSelector::selectTpcSets

(

const IndexVector &

itpss

)

Definition at line 68 of file WireSelector.cxx.

                                                         {
  m_tpcSets = itpss;
  selectPlanes();
}

void WireSelector::selectView

(

View

view

)

Definition at line 38 of file WireSelector.cxx.

                                         {
  m_view = a_view;
  selectPlanes();
}

void WireSelector::selectWireAngle	(	double	wireAngle,
		double	tol = 0.001
	)

Definition at line 45 of file WireSelector.cxx.

                                                           {
  m_wireAngle = angle;
  m_wireAngleTolerance = tol;
  selectPlanes();
}

View WireSelector::view ( ) const

inline

Definition at line 112 of file WireSelector.h.

{ return m_view; }

double WireSelector::wireAngle ( ) const

inline

Definition at line 115 of file WireSelector.h.

{ return m_wireAngle; }

double WireSelector::wireAngleTolerance ( ) const

inline

Definition at line 116 of file WireSelector.h.

{ return m_wireAngleTolerance; }

const WireSummary& WireSelector::wireSummary ( ) const

inline

Definition at line 129 of file WireSelector.h.

{ return m_wireSummary; }

Member Data Documentation

WireInfoVector WireSelector::m_data

private

Definition at line 186 of file WireSelector.h.

WireInfoMap WireSelector::m_datamap

private

Definition at line 187 of file WireSelector.h.

double WireSelector::m_driftMax = 1.e20

private

Definition at line 182 of file WireSelector.h.

double WireSelector::m_driftMin = 0.0

private

Definition at line 181 of file WireSelector.h.

bool WireSelector::m_haveData =false

private

Definition at line 185 of file WireSelector.h.

IndexVector WireSelector::m_icrys

private

Definition at line 177 of file WireSelector.h.

const GeometryCore* WireSelector::m_pgeo =nullptr

private

Definition at line 176 of file WireSelector.h.

PlaneIDVector WireSelector::m_pids

private

Definition at line 183 of file WireSelector.h.

IndexVector WireSelector::m_tpcSets

private

Definition at line 184 of file WireSelector.h.

View WireSelector::m_view =geo::kUnknown

private

Definition at line 178 of file WireSelector.h.

double WireSelector::m_wireAngle =999.

private

Definition at line 179 of file WireSelector.h.

double WireSelector::m_wireAngleTolerance =999.

private

Definition at line 180 of file WireSelector.h.

WireSummary WireSelector::m_wireSummary

private

Definition at line 188 of file WireSelector.h.

---

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

• dunecore/dunecore/Geometry/WireSelector.h
• dunecore/dunecore/Geometry/WireSelector.cxx