ChannelGeo.cxx

Go to the documentation of this file.

// ChannelGeo.cxx

#include "ChannelGeo.h"

#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "larcore/Geometry/Geometry.h"

//**********************************************************************

ChannelGeo::ChannelGeo(Index icha)
: ChannelGeo(icha, &*art::ServiceHandle<geo::Geometry>()) { }

//**********************************************************************

ChannelGeo::ChannelGeo(Index icha, const geo::GeometryCore* pgeo)
: m_icha(icha), m_valid(false) {
  if ( pgeo == nullptr ) return;
  if ( icha >= pgeo->Nchannels() ) return;
  m_valid = true;
  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;
    }
    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);
  }
}

//**********************************************************************