RectangularMTDTopology.cc

Go to the documentation of this file.

#include "Geometry/MTDGeometryBuilder/interface/RectangularMTDTopology.h"

// Modified for the large pixles.
// Danek Kotlinski & Michele Pioppi, 3/06.
// See documentation in the include file.

//--------------------------------------------------------------------
// Calculate the edge of the active sensor with respect to the center,
// that is simply the half-size.
// Take into account large pixels
void
RectangularMTDTopology::setOffset(const int& BIG_PIX_PER_ROC_X, const int& BIG_PIX_PER_ROC_Y,
                  const int& ROWS_PER_ROC, const int& COLS_PER_ROC) {

  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("RectangularMTDTopology")
    << "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 " << m_ROCS_X << ", ROCS_Y " << m_ROCS_Y
    << "\nNROWS " << m_ROWS_PER_ROC * m_ROCS_X << ", NCOL " << m_COLS_PER_ROC * m_ROCS_Y;

}

//--------------------------------------------------------------------
// PixelTopology interface. 
// Transform LocalPoint in cm to measurement in pitch units.
std::pair<float,float>
RectangularMTDTopology::pixel( const LocalPoint& p ) const
{
  // check limits   
  float py = p.y();
  float px = p.x();
  
  LogDebug("RectangularMTDTopology").log( [&](auto & debugstr) {
#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.
      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;
    }   
    }); 
  
  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("RectangularMTDTopology") << " 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("RectangularMTDTopology") << " 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("RectangularMTDTopology") << " 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("RectangularMTDTopology") << " bad pix x " << mpX << "\n"
                     << px << " " << m_xoffset << " " << m_pitchx << " "
                     << newxbin << " " << ixbin << " " << fractionX;
  }
#endif // EDM_ML_DEBUG
  
  return std::pair<float, float>( mpX, mpY );
}

//----------------------------------------------------------------------
// Topology interface, go from Measurement to Local corrdinates
// pixel coordinates (mp) -> cm (LocalPoint)
LocalPoint
RectangularMTDTopology::localPosition( const MeasurementPoint& mp ) const
{
  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("RectangularMTDTopology") << debugstr.str();
#endif // EDM_ML_DEBUG

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

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

//--------------------------------------------------------------------
// 
// measuremet to local transformation for X coordinate
// X coordinate is in the ROC row number direction
float
RectangularMTDTopology::localX( const float mpx ) const
{
  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("RectangularMTDTopology") << " 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("RectangularMTDTopology") << " 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("RectangularMTDTopology") << " bad lp x " << lpX << "\n"
                     << mpx << " " << binoffx << " "
                     << fractionX << " " << local_pitchx << " " << m_xoffset;
  }
#endif // EDM_ML_DEBUG

  return lpX;
} 

// measuremet to local transformation for Y coordinate
// Y is in the ROC column number direction 
float
RectangularMTDTopology::localY( const float mpy ) const
{
  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( "RectangularMTDTopology" ) << " 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( "RectangularMTDTopology" ) << " bad lp y " << lpY << "\n"
                       << mpy << " " << binoffy << " "
                       << fractionY << " " << local_pitchy << " " << m_yoffset;
  }
#endif // EDM_ML_DEBUG

  return lpY;
}

// Get hit errors in LocalPoint coordinates (cm)
LocalError
RectangularMTDTopology::localError( const MeasurementPoint& mp,
                      const MeasurementError& me ) const
{
  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));
}

// Get errors in pixel pitch units.
MeasurementError
RectangularMTDTopology::measurementError( const LocalPoint& lp,
                        const LocalError& le ) const
{
  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));
}