RectangularPixelTopology Class Reference

#include <RectangularPixelTopology.h>

Inheritance diagram for RectangularPixelTopology:

Public Member Functions
int	channel (const LocalPoint &lp) const override
int	colsperroc () const override
bool	containsBigPixelInX (int ixmin, int ixmax) const override
bool	containsBigPixelInY (int iymin, int iymax) const override
bool	isItBigPixelInX (const int ixbin) const override
bool	isItBigPixelInY (const int iybin) const override
bool	isItEdgePixel (int ixbin, int iybin) const override
bool	isItEdgePixelInX (int ixbin) const override
bool	isItEdgePixelInY (int iybin) const override
LocalError	localError (const MeasurementPoint &, const MeasurementError &) const override
LocalPoint	localPosition (const MeasurementPoint &mp) const override
float	localX (const float mpX) const override
float	localY (const float mpY) const override
MeasurementError	measurementError (const LocalPoint &, const LocalError &) const override
MeasurementPoint	measurementPosition (const LocalPoint &lp) const override
int	ncolumns () const override
int	nrows () const override
std::pair< float, float >	pitch () const override
std::pair< float, float >	pixel (const LocalPoint &p) const override
	RectangularPixelTopology (int nrows, int ncols, float pitchx, float pitchy, bool upgradeGeometry, int ROWS_PER_ROC, int COLS_PER_ROC, int BIG_PIX_PER_ROC_X, int BIG_PIX_PER_ROC_Y, int ROCS_X, int ROCS_Y)
int	rocsX () const override
int	rocsY () const override
int	rowsperroc () const override
float	xoffset () const
float	yoffset () const
Public Member Functions inherited from PixelTopology
virtual float	localX (const float mpX, const Topology::LocalTrackPred &) const
virtual float	localY (const float mpY, const Topology::LocalTrackPred &) const
virtual std::pair< float, float >	pixel (const LocalPoint &p, const Topology::LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
	~PixelTopology () override
Public Member Functions inherited from Topology
virtual int	channel (const LocalPoint &lp, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual LocalError	localError (const MeasurementPoint &mp, const MeasurementError &me, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual LocalPoint	localPosition (const MeasurementPoint &mp, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual MeasurementError	measurementError (const LocalPoint &lp, const LocalError &le, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual MeasurementPoint	measurementPosition (const LocalPoint &lp, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual	~Topology ()

Private Attributes
int	m_COLS_PER_ROC
int	m_ncols
int	m_nrows
float	m_pitchx
float	m_pitchy
int	m_ROCS_X
int	m_ROCS_Y
int	m_ROWS_PER_ROC
bool	m_upgradeGeometry
float	m_xoffset
float	m_yoffset

Additional Inherited Members
Public Types inherited from Topology
typedef Point2DBase< double, LocalTag >	Local2DPoint
typedef Vector2D::MathVector	MathVector2D
typedef Basic2DVector< double >	Vector2D

Detailed Description

Topology for rectangular pixel detector with BIG pixels.

Definition at line 39 of file RectangularPixelTopology.h.

Constructor & Destructor Documentation

RectangularPixelTopology()

RectangularPixelTopology::RectangularPixelTopology ( int  nrows, int  ncols, float  pitchx, float  pitchy, bool  upgradeGeometry, int  ROWS_PER_ROC, int  COLS_PER_ROC, int  BIG_PIX_PER_ROC_X, int  BIG_PIX_PER_ROC_Y, int  ROCS_X, int  ROCS_Y  )

inline

Definition at line 42 of file RectangularPixelTopology.h.

References LogDebug, m_COLS_PER_ROC, m_ncols, m_nrows, m_pitchx, m_pitchy, m_ROCS_X, m_ROCS_Y, m_ROWS_PER_ROC, m_xoffset, and m_yoffset.

      : m_pitchx(pitchx),
        m_pitchy(pitchy),
        m_nrows(nrows),
        m_ncols(ncols),
        m_ROWS_PER_ROC(ROWS_PER_ROC),  // Num of Rows per ROC
        m_COLS_PER_ROC(COLS_PER_ROC),  // Num of Cols per ROC
        m_ROCS_X(ROCS_X),              // 2 for SLHC
        m_ROCS_Y(ROCS_Y),              // 8 for SLHC
        m_upgradeGeometry(upgradeGeometry) {
    // Calculate the edge of the active sensor with respect to the center,
    // that is simply the half-size.
    // Take into account large pixels
    m_xoffset = -(m_nrows + BIG_PIX_PER_ROC_X * m_nrows / ROWS_PER_ROC) / 2. * m_pitchx;
    m_yoffset = -(m_ncols + BIG_PIX_PER_ROC_Y * m_ncols / COLS_PER_ROC) / 2. * m_pitchy;

    LogDebug("RectangularPixelTopology") << "nrows " << m_nrows << ", ncols " << m_ncols << ", pitchx " << m_pitchx
                                         << ", pitchy " << m_pitchy << ", xoffset " << m_xoffset << ", yoffset "
                                         << m_yoffset << ", BIG_PIX_PER_ROC_X " << BIG_PIX_PER_ROC_X
                                         << ", BIG_PIX_PER_ROC_Y " << BIG_PIX_PER_ROC_Y << ", ROWS_PER_ROC "
                                         << ROWS_PER_ROC << ", COLS_PER_ROC " << COLS_PER_ROC << ", ROCS_X " << ROCS_X
                                         << ", ROCS_Y " << ROCS_Y << "\nNROWS " << m_ROWS_PER_ROC * m_ROCS_X
                                         << ", NCOL " << m_COLS_PER_ROC * m_ROCS_Y;
  }

Member Function Documentation

channel()

int RectangularPixelTopology::channel ( const LocalPoint &  lp ) const

inlineoverridevirtual

Implements Topology.

Definition at line 100 of file RectangularPixelTopology.h.

References AlCaHLTBitMon_ParallelJobs::p, pixel(), and PixelChannelIdentifier::pixelToChannel().

                                                   {
    std::pair<float, float> p = pixel(lp);
    return PixelChannelIdentifier::pixelToChannel(int(p.first), int(p.second));
  }

colsperroc()

int RectangularPixelTopology::colsperroc ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 163 of file RectangularPixelTopology.h.

References m_COLS_PER_ROC.

{ return m_COLS_PER_ROC; }

containsBigPixelInX()

bool RectangularPixelTopology::containsBigPixelInX ( int  ixmin, int  ixmax  ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 132 of file RectangularPixelTopology.h.

References funct::false, and m_upgradeGeometry.

Referenced by PixelCPEGenericBase::initializeLocalErrorVariables(), PixelCPEBase::setTheClu(), and PixelTemplateSmearerBase::smearHit().

                                                                {
    return m_upgradeGeometry ? false : ((ixmin <= 80) & (ixmax >= 79));
  }

containsBigPixelInY()

bool RectangularPixelTopology::containsBigPixelInY ( int  iymin, int  iymax  ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 135 of file RectangularPixelTopology.h.

References funct::false, isItBigPixelInY(), and m_upgradeGeometry.

Referenced by PixelCPEGenericBase::initializeLocalErrorVariables(), PixelCPEBase::setTheClu(), and PixelTemplateSmearerBase::smearHit().

                                                                {
    return m_upgradeGeometry ? false
                             : (isItBigPixelInY(iymin) || isItBigPixelInY(iymax) || (iymin / 52) != (iymax / 52));
  }

isItBigPixelInX()

bool RectangularPixelTopology::isItBigPixelInX ( const int  ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 114 of file RectangularPixelTopology.h.

References m_upgradeGeometry.

Referenced by PixelCPEGeneric::localError(), localError(), PixelCPEGeneric::localPosition(), PixelCPEClusterRepair::localPosition(), PixelCPETemplateReco::localPosition(), and PixelTemplateSmearerBase::smearHit().

                                                       {
    return ((m_upgradeGeometry) ? (false) : ((ixbin == 79) | (ixbin == 80)));
  }

isItBigPixelInY()

bool RectangularPixelTopology::isItBigPixelInY ( const int  iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 118 of file RectangularPixelTopology.h.

References m_upgradeGeometry, and UNLIKELY.

Referenced by containsBigPixelInY(), PixelCPEGeneric::localError(), localError(), PixelCPEGeneric::localPosition(), PixelCPEClusterRepair::localPosition(), PixelCPETemplateReco::localPosition(), and PixelTemplateSmearerBase::smearHit().

                                                       {
    if UNLIKELY (m_upgradeGeometry)
      return false;
    else {
      int iybin0 = iybin % 52;
      return ((iybin0 == 0) | (iybin0 == 51));
      // constexpr int bigYIndeces[]{0,51,52,103,104,155,156,207,208,259,260,311,312,363,364,415,416,511};
      // return *std::lower_bound(std::begin(bigYIndeces),std::end(bigYIndeces),iybin) == iybin;
    }
  }

isItEdgePixel()

bool RectangularPixelTopology::isItEdgePixel ( int  ixbin, int  iybin  ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 145 of file RectangularPixelTopology.h.

References isItEdgePixelInX(), and isItEdgePixelInY().

                                                          {
    return (isItEdgePixelInX(ixbin) || isItEdgePixelInY(iybin));
  }

isItEdgePixelInX()

bool RectangularPixelTopology::isItEdgePixelInX ( int  ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 143 of file RectangularPixelTopology.h.

References m_nrows.

Referenced by PixelCPEGenericBase::initializeLocalErrorVariables(), isItEdgePixel(), PixelCPETemplateReco::localError(), and PixelTemplateSmearerBase::smearHit().

{ return ((ixbin == 0) | (ixbin == (m_nrows - 1))); }

isItEdgePixelInY()

bool RectangularPixelTopology::isItEdgePixelInY ( int  iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 144 of file RectangularPixelTopology.h.

References m_ncols.

Referenced by PixelCPEGenericBase::initializeLocalErrorVariables(), isItEdgePixel(), PixelCPETemplateReco::localError(), and PixelTemplateSmearerBase::smearHit().

{ return ((iybin == 0) | (iybin == (m_ncols - 1))); }

localError()

LocalError RectangularPixelTopology::localError ( const MeasurementPoint &  mp, const MeasurementError &  me  ) const

overridevirtual

Implements Topology.

Definition at line 288 of file RectangularPixelTopology.cc.

References dqmMemoryStats::float, createfilelist::int, isItBigPixelInX(), isItBigPixelInY(), m_pitchx, m_pitchy, hlt_dqm_clientPB-live_cfg::me, PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

                                                                                                            {
  float pitchy = m_pitchy;
  int binoffy = int(mp.y());
  if (isItBigPixelInY(binoffy))
    pitchy = 2. * m_pitchy;

  float pitchx = m_pitchx;
  int binoffx = int(mp.x());
  if (isItBigPixelInX(binoffx))
    pitchx = 2. * m_pitchx;

  return LocalError(me.uu() * float(pitchx * pitchx), 0, me.vv() * float(pitchy * pitchy));
}

localPosition()

LocalPoint RectangularPixelTopology::localPosition ( const MeasurementPoint &  mp ) const

overridevirtual

Implements Topology.

Definition at line 156 of file RectangularPixelTopology.cc.

References hcal_dqm_sourceclient-live_cfg::debugstr, EPS, dqmMemoryStats::float, localX(), localY(), LogDebug, m_ncols, m_nrows, PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

Referenced by PixelTemplateSmearerBase::smearHit().

                                                                                   {
  float mpy = mp.y();  // measurements
  float mpx = mp.x();

#ifdef EDM_ML_DEBUG
#define EPS 0
  // check limits
  std::ostringstream debugstr;

  if (mpy < 0.) {
    debugstr << " wrong mp y, fix " << mpy << " " << 0 << "\n";
    mpy = 0.;
  }
  if (mpy >= m_ncols) {
    debugstr << " wrong mp y, fix " << mpy << " " << m_ncols << "\n";
    mpy = float(m_ncols) - EPS;  // EPS is a small number
  }
  if (mpx < 0.) {
    debugstr << " wrong mp x, fix " << mpx << " " << 0 << "\n";
    mpx = 0.;
  }
  if (mpx >= m_nrows) {
    debugstr << " wrong mp x, fix " << mpx << " " << m_nrows << "\n";
    mpx = float(m_nrows) - EPS;  // EPS is a small number
  }
  if (!debugstr.str().empty())
    LogDebug("RectangularPixelTopology") << debugstr.str();
#endif  // EDM_ML_DEBUG

  float lpY = localY(mpy);
  float lpX = localX(mpx);

  // Return it as a LocalPoint
  return LocalPoint(lpX, lpY);
}

localX()

float RectangularPixelTopology::localX ( const float  mpX ) const

overridevirtual

Implements PixelTopology.

Definition at line 196 of file RectangularPixelTopology.cc.

References dqmMemoryStats::float, createfilelist::int, LogDebug, m_pitchx, m_ROCS_X, m_ROWS_PER_ROC, m_upgradeGeometry, m_xoffset, and UNLIKELY.

Referenced by localPosition().

                                                            {
  int binoffx = int(mpx);                  // truncate to int
  float fractionX = mpx - float(binoffx);  // find the fraction
  float local_pitchx = m_pitchx;           // defaultpitch

  if UNLIKELY (m_upgradeGeometry) {
#ifdef EDM_ML_DEBUG
    if (binoffx > m_ROWS_PER_ROC * m_ROCS_X)  // too large
    {
      LogDebug("RectangularPixelTopology")
          << " very bad, binx " << binoffx << "\n"
          << mpx << " " << binoffx << " " << fractionX << " " << local_pitchx << " " << m_xoffset << "\n";
    }
#endif
  } else {
    if (binoffx > 80) {  // ROC 1 - handles x on edge cluster
      binoffx = binoffx + 2;
    } else if (binoffx == 80) {  // ROC 1
      binoffx = binoffx + 1;
      local_pitchx *= 2;
    } else if (binoffx == 79) {  // ROC 0
      binoffx = binoffx + 0;
      local_pitchx *= 2;
    }
    // else if (binoffx>=0) {       // ROC 0
    //  binoffx=binoffx+0;
    // }

#ifdef EDM_ML_DEBUG
    if (binoffx < 0)  // too small
      LogDebug("RectangularPixelTopology")
          << " very bad, binx " << binoffx << "\n"
          << mpx << " " << binoffx << " " << fractionX << " " << local_pitchx << " " << m_xoffset;
#endif
  }

  // The final position in local coordinates
  float lpX = float(binoffx * m_pitchx) + fractionX * local_pitchx + m_xoffset;

#ifdef EDM_ML_DEBUG

  if (lpX < m_xoffset || lpX > (-m_xoffset)) {
    LogDebug("RectangularPixelTopology") << " bad lp x " << lpX << "\n"
                                         << mpx << " " << binoffx << " " << fractionX << " " << local_pitchx << " "
                                         << m_xoffset;
  }
#endif  // EDM_ML_DEBUG

  return lpX;
}

localY()

float RectangularPixelTopology::localY ( const float  mpY ) const

overridevirtual

Implements PixelTopology.

Definition at line 249 of file RectangularPixelTopology.cc.

References dqmMemoryStats::float, createfilelist::int, dqmiolumiharvest::j, LogDebug, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, m_COLS_PER_ROC, m_pitchy, m_ROCS_Y, m_upgradeGeometry, m_yoffset, and UNLIKELY.

Referenced by localPosition().

                                                            {
  int binoffy = int(mpy);                  // truncate to int
  float fractionY = mpy - float(binoffy);  // find the fraction
  float local_pitchy = m_pitchy;           // defaultpitch

  if UNLIKELY (m_upgradeGeometry) {
#ifdef EDM_ML_DEBUG
    if (binoffy > m_ROCS_Y * m_COLS_PER_ROC)  // too large
    {
      LogDebug("RectangularPixelTopology")
          << " very bad, biny " << binoffy << "\n"
          << mpy << " " << binoffy << " " << fractionY << " " << local_pitchy << " " << m_yoffset;
    }
#endif
  } else {  // 415 is last big pixel, 416 and above do not exists!
    constexpr int bigYIndeces[]{0, 51, 52, 103, 104, 155, 156, 207, 208, 259, 260, 311, 312, 363, 364, 415, 416, 511};
    auto const j = std::lower_bound(std::begin(bigYIndeces), std::end(bigYIndeces), binoffy);
    if (*j == binoffy)
      local_pitchy *= 2;
    binoffy += (j - bigYIndeces);
  }

  // The final position in local coordinates
  float lpY = float(binoffy * m_pitchy) + fractionY * local_pitchy + m_yoffset;

#ifdef EDM_ML_DEBUG

  if (lpY < m_yoffset || lpY > (-m_yoffset)) {
    LogDebug("RectangularPixelTopology") << " bad lp y " << lpY << "\n"
                                         << mpy << " " << binoffy << " " << fractionY << " " << local_pitchy << " "
                                         << m_yoffset;
  }
#endif  // EDM_ML_DEBUG

  return lpY;
}

measurementError()

MeasurementError RectangularPixelTopology::measurementError ( const LocalPoint &  lp, const LocalError &  le  ) const

overridevirtual

Implements Topology.

Definition at line 304 of file RectangularPixelTopology.cc.

References dqmMemoryStats::float, createfilelist::int, LIKELY, m_pitchx, m_pitchy, m_upgradeGeometry, m_xoffset, m_yoffset, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

                                                                                                            {
  float pitchy = m_pitchy;
  float pitchx = m_pitchx;

  if LIKELY (!m_upgradeGeometry) {
    int iybin = int((lp.y() - m_yoffset) / m_pitchy);  //get bin for equal picth
    int iybin0 = iybin % 54;                           //This is just to avoid many ifs by using the periodicy
    //quasi bins 0,1,52,53 fall into larger pixels
    if ((iybin0 <= 1) | (iybin0 >= 52))
      pitchy = 2.f * m_pitchy;

    int ixbin = int((lp.x() - m_xoffset) / m_pitchx);  //get bin for equal pitch
    //quasi bins 79,80,81,82 fall into the 2 larger pixels
    if ((ixbin >= 79) & (ixbin <= 82))
      pitchx = 2.f * m_pitchx;
  }

  return MeasurementError(le.xx() / float(pitchx * pitchx), 0, le.yy() / float(pitchy * pitchy));
}

measurementPosition()

MeasurementPoint RectangularPixelTopology::measurementPosition ( const LocalPoint &  lp ) const

inlineoverridevirtual

Implements Topology.

Definition at line 82 of file RectangularPixelTopology.h.

References AlCaHLTBitMon_ParallelJobs::p, and pixel().

Referenced by PixelTemplateSmearerBase::smearHit().

                                                                            {
    std::pair<float, float> p = pixel(lp);
    return MeasurementPoint(p.first, p.second);
  }

ncolumns()

int RectangularPixelTopology::ncolumns ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 155 of file RectangularPixelTopology.h.

References m_ncols.

Referenced by PixelCPEBase::setTheClu().

{ return (m_ncols); }

nrows()

int RectangularPixelTopology::nrows ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 153 of file RectangularPixelTopology.h.

References m_nrows.

Referenced by PixelCPEBase::setTheClu().

{ return (m_nrows); }

pitch()

std::pair<float, float> RectangularPixelTopology::pitch ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 151 of file RectangularPixelTopology.h.

References dqmMemoryStats::float, m_pitchx, and m_pitchy.

{ return std::pair<float, float>(float(m_pitchx), float(m_pitchy)); }

pixel()

std::pair< float, float > RectangularPixelTopology::pixel ( const LocalPoint &  p ) const

overridevirtual

Topology for rectangular pixel detector with BIG pixels.

Implements PixelTopology.

Definition at line 13 of file RectangularPixelTopology.cc.

References hcal_dqm_sourceclient-live_cfg::debugstr, dqmMemoryStats::float, createfilelist::int, LogDebug, m_COLS_PER_ROC, m_pitchx, m_pitchy, m_upgradeGeometry, m_xoffset, m_yoffset, AlCaHLTBitMon_ParallelJobs::p, multPhiCorr_741_25nsDY_cfi::px, and multPhiCorr_741_25nsDY_cfi::py.

Referenced by channel(), and measurementPosition().

                                                                               {
  // check limits
  float py = p.y();
  float px = p.x();

#ifdef EDM_ML_DEBUG
#define EPSCM 0
#define EPS 0
  // This will catch points which are outside the active sensor area.
  // In the digitizer during the early induce_signal phase non valid
  // location are passed here. They are cleaned later.

  std::ostringstream debugstr;
  debugstr << "py = " << py << ", m_yoffset = " << m_yoffset << "px = " << px << ", m_xoffset = " << m_xoffset << "\n";

  if (py < m_yoffset)  // m_yoffset is negative
  {
    debugstr << " wrong lp y " << py << " " << m_yoffset << "\n";
    py = m_yoffset + EPSCM;  // make sure it is in, add an EPS in cm
  }
  if (py > -m_yoffset) {
    debugstr << " wrong lp y " << py << " " << -m_yoffset << "\n";
    py = -m_yoffset - EPSCM;
  }
  if (px < m_xoffset)  // m_xoffset is negative
  {
    debugstr << " wrong lp x " << px << " " << m_xoffset << "\n";
    px = m_xoffset + EPSCM;
  }
  if (px > -m_xoffset) {
    debugstr << " wrong lp x " << px << " " << -m_xoffset << "\n";
    px = -m_xoffset - EPSCM;
  }

  if (!debugstr.str().empty())
    LogDebug("RectangularPixelTopology") << debugstr.str();
#endif  // EDM_ML_DEBUG

  float newybin = (py - m_yoffset) / m_pitchy;
  int iybin = int(newybin);
  float fractionY = newybin - iybin;

  // Normalize it all to 1 ROC
  int iybin0 = 0;
  int numROC = 0;
  float mpY = 0.;

  if (m_upgradeGeometry) {
    iybin0 = (iybin % m_COLS_PER_ROC);  // 0-51
    numROC = iybin / m_COLS_PER_ROC;    // 0-7
    mpY = float(numROC * m_COLS_PER_ROC + iybin0) + fractionY;

#ifdef EDM_ML_DEBUG

    if (iybin0 > m_COLS_PER_ROC) {
      LogDebug("RectangularPixelTopology") << " very bad, newbiny " << iybin0 << "\n"
                                           << py << " " << m_yoffset << " " << m_pitchy << " " << newybin << " "
                                           << iybin << " " << fractionY << " " << iybin0 << " " << numROC;
    }
#endif  // EDM_ML_DEBUG

  } else {
    iybin0 = (iybin % 54);  // 0-53
    numROC = iybin / 54;    // 0-7

    if (iybin0 == 53) {  // inside big pixel
      iybin0 = 51;
      fractionY = (fractionY + 1.) / 2.;
    } else if (iybin0 == 52) {  // inside big pixel
      iybin0 = 51;
      fractionY = fractionY / 2.;
    } else if (iybin0 > 1) {  // inside normal pixel
      iybin0 = iybin0 - 1;
    } else if (iybin0 == 1) {  // inside big pixel
      iybin0 = 0;
      fractionY = (fractionY + 1.) / 2.;
    } else if (iybin0 == 0) {  // inside big pixel
      iybin0 = 0;
      fractionY = fractionY / 2.;
    }

    mpY = float(numROC * 52. + iybin0) + fractionY;
  }

#ifdef EDM_ML_DEBUG

  if (mpY < 0. || mpY >= 416.) {
    LogDebug("RectangularPixelTopology") << " bad pix y " << mpY << "\n"
                                         << py << " " << m_yoffset << " " << m_pitchy << " " << newybin << " " << iybin
                                         << " " << fractionY << " " << iybin0 << " " << numROC;
  }
#endif  // EDM_ML_DEBUG

  // In X
  float newxbin = (px - m_xoffset) / m_pitchx;
  int ixbin = int(newxbin);
  float fractionX = newxbin - ixbin;

#ifdef EDM_ML_DEBUG

  if (ixbin > 161 || ixbin < 0)  //  ixbin < 0 outside range
  {
    LogDebug("RectangularPixelTopology") << " very bad, newbinx " << ixbin << "\n"
                                         << px << " " << m_xoffset << " " << m_pitchx << " " << newxbin << " " << ixbin
                                         << " " << fractionX;
  }
#endif  // EDM_ML_DEBUG

  if (!m_upgradeGeometry) {
    if (ixbin > 82) {  // inside normal pixel, ROC 1
      ixbin = ixbin - 2;
    } else if (ixbin == 82) {  // inside bin pixel
      ixbin = 80;
      fractionX = (fractionX + 1.) / 2.;
    } else if (ixbin == 81) {  // inside big pixel
      ixbin = 80;
      fractionX = fractionX / 2.;
    } else if (ixbin == 80) {  // inside bin pixel, ROC 0
      ixbin = 79;
      fractionX = (fractionX + 1.) / 2.;
    } else if (ixbin == 79) {  // inside big pixel
      ixbin = 79;
      fractionX = fractionX / 2.;
    }
  }

  float mpX = float(ixbin) + fractionX;

#ifdef EDM_ML_DEBUG

  if (mpX < 0. || mpX >= 160.) {
    LogDebug("RectangularPixelTopology") << " bad pix x " << mpX << "\n"
                                         << px << " " << m_xoffset << " " << m_pitchx << " " << newxbin << " " << ixbin
                                         << " " << fractionX;
  }
#endif  // EDM_ML_DEBUG

  return std::pair<float, float>(mpX, mpY);
}

rocsX()

int RectangularPixelTopology::rocsX ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 159 of file RectangularPixelTopology.h.

References m_ROCS_X.

{ return m_ROCS_X; }

rocsY()

int RectangularPixelTopology::rocsY ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 157 of file RectangularPixelTopology.h.

References m_ROCS_Y.

{ return m_ROCS_Y; }

rowsperroc()

int RectangularPixelTopology::rowsperroc ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 161 of file RectangularPixelTopology.h.

References m_ROWS_PER_ROC.

{ return m_ROWS_PER_ROC; }

xoffset()

float RectangularPixelTopology::xoffset ( ) const

inline

Definition at line 164 of file RectangularPixelTopology.h.

References m_xoffset.

{ return m_xoffset; }

yoffset()

float RectangularPixelTopology::yoffset ( ) const

inline

Definition at line 165 of file RectangularPixelTopology.h.

References m_yoffset.

{ return m_yoffset; }

Member Data Documentation

m_COLS_PER_ROC

int RectangularPixelTopology::m_COLS_PER_ROC

private

Definition at line 175 of file RectangularPixelTopology.h.

Referenced by colsperroc(), localY(), pixel(), and RectangularPixelTopology().

m_ncols

int RectangularPixelTopology::m_ncols

private

Definition at line 173 of file RectangularPixelTopology.h.

Referenced by isItEdgePixelInY(), localPosition(), ncolumns(), and RectangularPixelTopology().

m_nrows

int RectangularPixelTopology::m_nrows

private

Definition at line 172 of file RectangularPixelTopology.h.

Referenced by isItEdgePixelInX(), localPosition(), nrows(), and RectangularPixelTopology().

m_pitchx

float RectangularPixelTopology::m_pitchx

private

Definition at line 168 of file RectangularPixelTopology.h.

Referenced by localError(), localX(), measurementError(), pitch(), pixel(), and RectangularPixelTopology().

m_pitchy

float RectangularPixelTopology::m_pitchy

private

Definition at line 169 of file RectangularPixelTopology.h.

Referenced by localError(), localY(), measurementError(), pitch(), pixel(), and RectangularPixelTopology().

m_ROCS_X

int RectangularPixelTopology::m_ROCS_X

private

Definition at line 176 of file RectangularPixelTopology.h.

Referenced by localX(), RectangularPixelTopology(), and rocsX().

m_ROCS_Y

int RectangularPixelTopology::m_ROCS_Y

private

Definition at line 177 of file RectangularPixelTopology.h.

Referenced by localY(), RectangularPixelTopology(), and rocsY().

m_ROWS_PER_ROC

int RectangularPixelTopology::m_ROWS_PER_ROC

private

Definition at line 174 of file RectangularPixelTopology.h.

Referenced by localX(), RectangularPixelTopology(), and rowsperroc().

m_upgradeGeometry

bool RectangularPixelTopology::m_upgradeGeometry

private

Definition at line 178 of file RectangularPixelTopology.h.

Referenced by containsBigPixelInX(), containsBigPixelInY(), isItBigPixelInX(), isItBigPixelInY(), localX(), localY(), measurementError(), and pixel().

m_xoffset

float RectangularPixelTopology::m_xoffset

private

Definition at line 170 of file RectangularPixelTopology.h.

Referenced by localX(), measurementError(), pixel(), RectangularPixelTopology(), and xoffset().

m_yoffset

float RectangularPixelTopology::m_yoffset

private

Definition at line 171 of file RectangularPixelTopology.h.

Referenced by localY(), measurementError(), pixel(), RectangularPixelTopology(), and yoffset().