#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.

       : 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;
   }

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.

Member Function Documentation

◆ channel()

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

inlineoverridevirtual

Implements Topology.

Definition at line 100 of file RectangularPixelTopology.h.

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

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

◆ colsperroc()

int RectangularPixelTopology::colsperroc ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 163 of file RectangularPixelTopology.h.

163 { return m_COLS_PER_ROC; }

References m_COLS_PER_ROC.

◆ containsBigPixelInX()

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 132 of file RectangularPixelTopology.h.

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

References funct::false, and m_upgradeGeometry.

Referenced by PixelCPEGeneric::localError(), PixelCPEBase::setTheClu(), and PixelTemplateSmearerBase::smearHit().

◆ containsBigPixelInY()

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 135 of file RectangularPixelTopology.h.

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

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

Referenced by PixelCPEGeneric::localError(), PixelCPEBase::setTheClu(), and PixelTemplateSmearerBase::smearHit().

◆ isItBigPixelInX()

bool RectangularPixelTopology::isItBigPixelInX ( const int ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 114 of file RectangularPixelTopology.h.

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

References m_upgradeGeometry.

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

◆ isItBigPixelInY()

bool RectangularPixelTopology::isItBigPixelInY ( const int iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 118 of file RectangularPixelTopology.h.

                                                        {
     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;
     }
   }

References m_upgradeGeometry, and UNLIKELY.

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

◆ isItEdgePixel()

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 145 of file RectangularPixelTopology.h.

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

References isItEdgePixelInX(), and isItEdgePixelInY().

◆ isItEdgePixelInX()

bool RectangularPixelTopology::isItEdgePixelInX ( int ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 143 of file RectangularPixelTopology.h.

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

References m_nrows.

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

◆ isItEdgePixelInY()

bool RectangularPixelTopology::isItEdgePixelInY ( int iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 144 of file RectangularPixelTopology.h.

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

References m_ncols.

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

◆ localError()

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

overridevirtual

Implements Topology.

Definition at line 292 of file RectangularPixelTopology.cc.

                                                                                                             {
   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));
 }

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().

◆ localPosition()

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

overridevirtual

Implements Topology.

Definition at line 156 of file RectangularPixelTopology.cc.

                                                                                    {
   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);
 }

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().

◆ localX()

float RectangularPixelTopology::localX ( const float mpX ) const

overridevirtual

Implements PixelTopology.

Definition at line 196 of file RectangularPixelTopology.cc.

                                                             {
   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;
 }

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

Referenced by localPosition().

◆ localY()

float RectangularPixelTopology::localY ( const float mpY ) const

overridevirtual

Implements PixelTopology.

Definition at line 251 of file RectangularPixelTopology.cc.

                                                             {
   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;
 }

References begin, end, dqmMemoryStats::float, createfilelist::int, dqmiolumiharvest::j, LogDebug, cuda_std::lower_bound(), m_COLS_PER_ROC, m_pitchy, m_ROCS_Y, m_upgradeGeometry, m_yoffset, and UNLIKELY.

Referenced by localPosition().

◆ measurementError()

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

overridevirtual

Implements Topology.

Definition at line 308 of file RectangularPixelTopology.cc.

                                                                                                             {
   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));
 }

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().

◆ measurementPosition()

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

inlineoverridevirtual

Implements Topology.

Definition at line 82 of file RectangularPixelTopology.h.

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

References AlCaHLTBitMon_ParallelJobs::p, and pixel().

Referenced by PixelTemplateSmearerBase::smearHit().

◆ ncolumns()

int RectangularPixelTopology::ncolumns ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 155 of file RectangularPixelTopology.h.

155 { return (m_ncols); }

References m_ncols.

Referenced by PixelCPEBase::setTheClu().

◆ nrows()

int RectangularPixelTopology::nrows ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 153 of file RectangularPixelTopology.h.

153 { return (m_nrows); }

References m_nrows.

Referenced by PixelCPEBase::setTheClu().

◆ pitch()

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 151 of file RectangularPixelTopology.h.

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

References m_pitchx, and 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.

                                                                                {
   // 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);
 }