RectangularMTDTopology Class Reference

#include <RectangularMTDTopology.h>

Inheritance diagram for RectangularMTDTopology:

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
LocalPoint	moduleToPixelLocalPoint (const LocalPoint &mlp) const
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
LocalPoint	pixelToModuleLocalPoint (const LocalPoint &plp, int row, int col) const
LocalPoint	pixelToModuleLocalPoint (const LocalPoint &plp, int channel) const
	RectangularMTDTopology (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
void	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)
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 40 of file RectangularMTDTopology.h.

Constructor & Destructor Documentation

RectangularMTDTopology::RectangularMTDTopology ( 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 45 of file RectangularMTDTopology.h.

References localPosition(), and setOffset().

    : 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 ) {

    setOffset(BIG_PIX_PER_ROC_X, BIG_PIX_PER_ROC_Y, ROWS_PER_ROC, COLS_PER_ROC);

      }

Member Function Documentation

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

inlineoverridevirtual

Implements Topology.

Definition at line 98 of file RectangularMTDTopology.h.

References pixel(), and MTDChannelIdentifier::pixelToChannel().

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

int RectangularMTDTopology::colsperroc ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 198 of file RectangularMTDTopology.h.

References m_COLS_PER_ROC.

                                  {      
    return m_COLS_PER_ROC;   
  }

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 149 of file RectangularMTDTopology.h.

References m_upgradeGeometry.

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

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 152 of file RectangularMTDTopology.h.

References isItBigPixelInY(), and m_upgradeGeometry.

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

bool RectangularMTDTopology::isItBigPixelInX ( const int  ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 132 of file RectangularMTDTopology.h.

References m_upgradeGeometry.

Referenced by localError().

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

bool RectangularMTDTopology::isItBigPixelInY ( const int  iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 136 of file RectangularMTDTopology.h.

References m_upgradeGeometry, and UNLIKELY.

Referenced by containsBigPixelInY(), and localError().

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

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 168 of file RectangularMTDTopology.h.

References isItEdgePixelInX(), and isItEdgePixelInY().

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

bool RectangularMTDTopology::isItEdgePixelInX ( int  ixbin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 162 of file RectangularMTDTopology.h.

References m_nrows.

Referenced by isItEdgePixel().

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

bool RectangularMTDTopology::isItEdgePixelInY ( int  iybin ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 165 of file RectangularMTDTopology.h.

References m_ncols.

Referenced by isItEdgePixel().

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

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

overridevirtual

Implements Topology.

Definition at line 330 of file RectangularMTDTopology.cc.

References objects.autophobj::float, createfilelist::int, isItBigPixelInX(), isItBigPixelInY(), m_pitchx, m_pitchy, MeasurementError::uu(), MeasurementError::vv(), PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

Referenced by measurementPosition().

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

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

overridevirtual

Implements Topology.

Definition at line 187 of file RectangularMTDTopology.cc.

References EPS, objects.autophobj::float, localX(), localY(), LogDebug, m_ncols, m_nrows, PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

Referenced by RectangularMTDTopology().

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

float RectangularMTDTopology::localX ( const float  mpX ) const

overridevirtual

Implements PixelTopology.

Definition at line 233 of file RectangularMTDTopology.cc.

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

Referenced by localPosition(), and pixelToModuleLocalPoint().

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

float RectangularMTDTopology::localY ( const float  mpY ) const

overridevirtual

Implements PixelTopology.

Definition at line 289 of file RectangularMTDTopology.cc.

References begin, constexpr, end, objects.autophobj::float, createfilelist::int, LogDebug, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, m_COLS_PER_ROC, m_pitchy, m_ROCS_Y, m_upgradeGeometry, m_yoffset, and UNLIKELY.

Referenced by localPosition(), and pixelToModuleLocalPoint().

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

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

overridevirtual

Implements Topology.

Definition at line 348 of file RectangularMTDTopology.cc.

References objects.autophobj::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().

Referenced by measurementPosition().

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

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

inlineoverridevirtual

Implements Topology.

Definition at line 77 of file RectangularMTDTopology.h.

References localError(), measurementError(), AlCaHLTBitMon_ParallelJobs::p, and pixel().

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

LocalPoint RectangularMTDTopology::moduleToPixelLocalPoint ( const LocalPoint &  mlp ) const

inline

Definition at line 107 of file RectangularMTDTopology.h.

References createfilelist::int, m_pitchx, m_pitchy, m_xoffset, m_yoffset, pixel(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                  {
    std::pair<float,float> p = pixel( mlp );
    return LocalPoint( mlp.x() - (m_xoffset + (int(p.first)+0.5f)*m_pitchx),
               mlp.y() - (m_yoffset + (int(p.second)+0.5f)*m_pitchy),
               mlp.z());
  }

int RectangularMTDTopology::ncolumns ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 182 of file RectangularMTDTopology.h.

References m_ncols.

Referenced by MTDThresholdClusterizer::setup().

                                {
    return ( m_ncols );
  }

int RectangularMTDTopology::nrows ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 178 of file RectangularMTDTopology.h.

References m_nrows.

Referenced by BtlDigiHitsValidation::analyze(), BtlRecHitsValidation::analyze(), BtlSimHitsValidation::analyze(), BTLTileDeviceSim::getHitsResponse(), BTLBarDeviceSim::getHitsResponse(), and MTDThresholdClusterizer::setup().

                             {
    return ( m_nrows );
  }

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

inlineoverridevirtual

Implements PixelTopology.

Definition at line 174 of file RectangularMTDTopology.h.

References objects.autophobj::float, m_pitchx, and m_pitchy.

Referenced by BTLBarDeviceSim::getHitsResponse(), and MTDTimeCalib::getTimeCalib().

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

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

overridevirtual

Implements PixelTopology.

Definition at line 35 of file RectangularMTDTopology.cc.

References EPSCM, objects.autophobj::float, createfilelist::int, LogDebug, m_COLS_PER_ROC, m_pitchx, m_pitchy, m_upgradeGeometry, m_xoffset, m_yoffset, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by channel(), BTLTileDeviceSim::getHitsResponse(), BTLBarDeviceSim::getHitsResponse(), measurementPosition(), and moduleToPixelLocalPoint().

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

LocalPoint RectangularMTDTopology::pixelToModuleLocalPoint ( const LocalPoint &  plp, int  row, int  col  ) const

inline

Definition at line 113 of file RectangularMTDTopology.h.

References m_pitchx, m_pitchy, m_xoffset, m_yoffset, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by BtlDigiHitsValidation::analyze(), BtlRecHitsValidation::analyze(), BtlSimHitsValidation::analyze(), BTLTileDeviceSim::getHitsResponse(), BTLBarDeviceSim::getHitsResponse(), and pixelToModuleLocalPoint().

                                                                                    {    
    return LocalPoint( plp.x() + (m_xoffset + (row+0.5f)*m_pitchx),
               plp.y() + (m_yoffset + (col+0.5f)*m_pitchy),
               plp.z()); 
  }

LocalPoint RectangularMTDTopology::pixelToModuleLocalPoint ( const LocalPoint &  plp, int  channel  ) const

inline

Definition at line 118 of file RectangularMTDTopology.h.

References MTDChannelIdentifier::channelToPixel(), localX(), localY(), and pixelToModuleLocalPoint().

                                                                               {
    std::pair<int,int> p = MTDChannelIdentifier::channelToPixel(channel);
    return pixelToModuleLocalPoint(plp,p.first,p.second);
  }

int RectangularMTDTopology::rocsX ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 190 of file RectangularMTDTopology.h.

References m_ROCS_X.

                             {   
    return m_ROCS_X;     
  }      

int RectangularMTDTopology::rocsY ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 186 of file RectangularMTDTopology.h.

References m_ROCS_Y.

                             {   
    return m_ROCS_Y;     
  }      

int RectangularMTDTopology::rowsperroc ( ) const

inlineoverridevirtual

Implements PixelTopology.

Definition at line 194 of file RectangularMTDTopology.h.

References m_ROWS_PER_ROC.

                                  {      
    return m_ROWS_PER_ROC;   
  }      

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
	)

Topology for rectangular pixel detector with BIG pixels.

Definition at line 15 of file RectangularMTDTopology.cc.

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.

Referenced by RectangularMTDTopology().

                                                                    {

  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;

}

float RectangularMTDTopology::xoffset ( ) const

inline

Definition at line 201 of file RectangularMTDTopology.h.

References m_xoffset.

                        {
    return m_xoffset;
  }

float RectangularMTDTopology::yoffset ( ) const

inline

Definition at line 204 of file RectangularMTDTopology.h.

References m_yoffset.

                        {
    return m_yoffset;
  }

Member Data Documentation

int RectangularMTDTopology::m_COLS_PER_ROC

private

Definition at line 218 of file RectangularMTDTopology.h.

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

int RectangularMTDTopology::m_ncols

private

Definition at line 216 of file RectangularMTDTopology.h.

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

int RectangularMTDTopology::m_nrows

private

Definition at line 215 of file RectangularMTDTopology.h.

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

float RectangularMTDTopology::m_pitchx

private

Definition at line 211 of file RectangularMTDTopology.h.

Referenced by localError(), localX(), measurementError(), moduleToPixelLocalPoint(), pitch(), pixel(), pixelToModuleLocalPoint(), and setOffset().

float RectangularMTDTopology::m_pitchy

private

Definition at line 212 of file RectangularMTDTopology.h.

Referenced by localError(), localY(), measurementError(), moduleToPixelLocalPoint(), pitch(), pixel(), pixelToModuleLocalPoint(), and setOffset().

int RectangularMTDTopology::m_ROCS_X

private

Definition at line 219 of file RectangularMTDTopology.h.

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

int RectangularMTDTopology::m_ROCS_Y

private

Definition at line 220 of file RectangularMTDTopology.h.

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

int RectangularMTDTopology::m_ROWS_PER_ROC

private

Definition at line 217 of file RectangularMTDTopology.h.

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

bool RectangularMTDTopology::m_upgradeGeometry

private

Definition at line 221 of file RectangularMTDTopology.h.

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

float RectangularMTDTopology::m_xoffset

private

Definition at line 213 of file RectangularMTDTopology.h.

Referenced by localX(), measurementError(), moduleToPixelLocalPoint(), pixel(), pixelToModuleLocalPoint(), setOffset(), and xoffset().

float RectangularMTDTopology::m_yoffset

private

Definition at line 214 of file RectangularMTDTopology.h.

Referenced by localY(), measurementError(), moduleToPixelLocalPoint(), pixel(), pixelToModuleLocalPoint(), setOffset(), and yoffset().