#include <IcebergChannelGeo.h>

Public Types
using	Index = unsigned int

using	Point = TVector3

using	EndPoints = std::pair< Point, Point >

using	EndPointsVector = std::vector< EndPoints >

Public Member Functions
	IcebergChannelGeo (Index icha, const geo::GeometryCore *pgeo)

	IcebergChannelGeo (Index icha)

Index	channel () const

bool	isValid () const

Index	nWires () const

const EndPointsVector &	wires () const

const Point &	top () const

const Point &	bottom () const

Private Attributes
Index	m_icha

bool	m_valid

EndPointsVector	m_ends

Point	m_top

Point	m_bot

Detailed Description

Definition at line 23 of file IcebergChannelGeo.h.

Member Typedef Documentation

using IcebergChannelGeo::EndPoints = std::pair<Point, Point>

Definition at line 29 of file IcebergChannelGeo.h.

using IcebergChannelGeo::EndPointsVector = std::vector<EndPoints>

Definition at line 30 of file IcebergChannelGeo.h.

using IcebergChannelGeo::Index = unsigned int

Definition at line 27 of file IcebergChannelGeo.h.

using IcebergChannelGeo::Point = TVector3

Definition at line 28 of file IcebergChannelGeo.h.

Constructor & Destructor Documentation

IcebergChannelGeo::IcebergChannelGeo	(	Index	icha,
		const geo::GeometryCore *	pgeo
	)

Definition at line 21 of file IcebergChannelGeo.cxx.

 : m_icha(icha), m_valid(false) {
   const string myname = "IcebergChannelGeo::ctor: ";
   if ( pgeo == nullptr ) return;
   if ( icha >= pgeo->Nchannels() ) return;
   m_valid = true;
   if ( pgeo->DetectorName() != "iceberg" ) {
     cout << myname << "WARNING: Iceberg mod is applied to non-Iceberg detector channels." << endl;
   }
   double ylo = 173.819 - 95.23;
   for ( geo::WireID wid : pgeo->ChannelToWire(icha) ) {
     EndPoints ends = pgeo->WireEndPoints<TVector3>(wid);
     if ( ends.first.y() > ends.second.y() ) {
       Point firstPoint = ends.first;
       ends.first = ends.second;
       ends.second = firstPoint;
     }
     Point& plo = ends.first;
     Point& phi = ends.second;
     if ( phi.y() < ylo ) continue;
     if ( plo.y() < ylo ) {
       // Set lower limit to ylo.
       double slop = (plo.z() - phi.z())/(plo.y() - phi.y());
       double zlo = phi.z() + slop*(ylo - phi.y());
       plo.SetY(ylo);
       plo.SetZ(zlo);
     }
     if ( m_ends.size() == 0 ) {
        m_bot = ends.first;
        m_top = ends.second;
     } else {
        if ( ends.first.y() < m_bot.y() ) m_bot = ends.first;
        if ( ends.second.y() > m_top.y() ) m_top = ends.second;
     }
     m_ends.push_back(ends);
   }
 }