ClusterShape Class Reference

#include <ClusterShape.h>

Public Member Functions
	ClusterShape ()
void	determineShape (const PixelGeomDetUnit &pixelDet, const SiPixelRecHit &recHit, ClusterData &data)
void	determineShape (const PixelGeomDetUnit &pixelDet, const SiPixelCluster &cluster, ClusterData &data)
	~ClusterShape ()

Private Member Functions
int	getDirection (int low, int hig, int olow, int ohig)
bool	processColumn (std::pair< int, int > pos, bool inTheLoop)

Private Attributes
int	hig
int	low
int	odir
int	ohig
int	olow
std::vector< SiPixelCluster::Pixel >	pixels_
int	x [2]
int	y [2]

Detailed Description

Definition at line 14 of file ClusterShape.h.

Constructor & Destructor Documentation

ClusterShape::ClusterShape

(

)

Definition at line 15 of file ClusterShape.cc.

{
}

ClusterShape::~ClusterShape

(

)

Definition at line 20 of file ClusterShape.cc.

{
}

Member Function Documentation

void ClusterShape::determineShape	(	const PixelGeomDetUnit &	pixelDet,
		const SiPixelRecHit &	recHit,
		ClusterData &	data
	)

Definition at line 80 of file ClusterShape.cc.

References SiPixelRecHit::cluster().

Referenced by processColumn(), and SiPixelClusterShapeCacheProducer::produce().

{
  determineShape(pixelDet, *(recHit.cluster()), data);
}

void ClusterShape::determineShape	(	const PixelGeomDetUnit &	pixelDet,
		const SiPixelCluster &	cluster,
		ClusterData &	data
	)

Definition at line 87 of file ClusterShape.cc.

References a, create_public_lumi_plots::ax, b, edm::VecArray< T, N >::capacity(), PVValHelper::dx, PVValHelper::dy, ClusterData::hasBigPixelsOnlyInside, mps_fire::i, createfilelist::int, ClusterData::isComplete, PixelTopology::isItBigPixelInX(), PixelTopology::isItBigPixelInY(), PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), ClusterData::isStraight, SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), heppy_hadd::odir, digitizers_cfi::pixel, SiPixelCluster::pixel(), SiPixelCluster::pixelADC(), edm::VecArray< T, N >::push_back_unchecked(), alignCSCRings::s, ClusterData::size, jetUpdater_cfi::sort, PixelGeomDetUnit::specificTopology(), x, and y.

{
  // Topology
  const PixelTopology * theTopology = (&(pixelDet.specificTopology())); 
 
  // Initialize
  data.isStraight = true;
  data.isComplete = true;
 
  x[0] = -1; x[1] = -1;
  y[0] = -1; y[1] = -1;
  olow = -2; ohig = -2; odir = 0;
  low = 0; hig = 0;

  pair<int,int> pos;
 
  // Get sorted pixels
  size_t npixels = cluster.pixelADC().size();
  pixels_.reserve(npixels);
  for(size_t i=0; i<npixels; ++i) {
    pixels_.push_back(cluster.pixel(i));
  }
  sort(pixels_.begin(),pixels_.end(),[](auto const& a, auto const& b){
              // slightly faster by avoiding branches
              return (a.x < b.x) | ((a.x == b.x) & (a.y < b.y));
          });

  // Look at all the pixels
  for(const auto& pixel: pixels_)
  {
    // Position
    pos.first  = (int)pixel.x;
    pos.second = (int)pixel.y;

    // Check if at the edge or big 
    if(theTopology->isItEdgePixelInX(pos.first) ||
       theTopology->isItEdgePixelInY(pos.second))
    { data.isComplete = false; } // break; }
 
    // Check if straight
    if(pos.first > x[1])
    { // column ready
      if(processColumn(pos, true) == false)
      { data.isStraight = false; } // break; }
    }
    else
    { // increasing column
      if(pos.second > hig+1) // at least a pixel is missing
      { data.isStraight = false; } // break; }
 
      hig = pos.second;
    }
  }
  pixels_.clear();

  // Check if straight, process last column
  if(processColumn(pos, false) == false)
    data.isStraight = false;

  // Treat clusters with big pixel(s) inside
  const int minPixelRow = cluster.minPixelRow();
  const int maxPixelRow = cluster.maxPixelRow();
  for(int ix = minPixelRow + 1;
          ix < maxPixelRow; ix++)
    x[1] += theTopology->isItBigPixelInX(ix);
 
  const int minPixelCol = cluster.minPixelCol();
  const int maxPixelCol = cluster.maxPixelCol();
  for(int iy = minPixelCol + 1;
          iy < maxPixelCol; iy++)
    y[1] += theTopology->isItBigPixelInY(iy);

  // Treat clusters with bix pixel(s) outside, FIXME FIXME
  unsigned int px = 0;
  px += theTopology->isItBigPixelInX(minPixelRow);
  px += theTopology->isItBigPixelInX(maxPixelRow);

  unsigned int py = 0;
  py += theTopology->isItBigPixelInY(minPixelCol);
  py += theTopology->isItBigPixelInY(maxPixelCol);

  data.hasBigPixelsOnlyInside = (px <= 0 && py <= 0);

  //if( (px > 0 || py > 0) && odir == 0)
  if( !data.hasBigPixelsOnlyInside && odir == 0)
  {
    // if outside and don't know the direction FIXME?
    data.isComplete = false;
  }
  // else
  { // FIXME do it
    assert((px+1)*(py+1) <= data.size.capacity());
    const int pre_dx = x[1] - x[0];
    const int pre_dy = y[1] - y[0];
    for(unsigned int ax = 0; ax <= px; ax++)
    for(unsigned int ay = 0; ay <= py; ay++)
    {
      int dx = pre_dx + ax;
      int dy = pre_dy + ay;
      if(odir != 0) dy *= odir;
  
      pair<int,int> s(dx,dy);
      data.size.push_back_unchecked(s);
    }
  }
}

int ClusterShape::getDirection ( int  low, int  hig, int  olow, int  ohig  )

private

Definition at line 25 of file ClusterShape.cc.

{
  if(hig == ohig && low == olow) return  0;
  if(hig >= ohig && low >= olow) return  1; 
  if(hig <= ohig && low <= olow) return -1;
   
  return -2;
}

bool ClusterShape::processColumn ( std::pair< int, int >  pos, bool  inTheLoop  )

private

Definition at line 35 of file ClusterShape.cc.

References determineShape(), dir, heppy_hadd::odir, x, and y.

{
  if(x[1] > -1)
  { // Process previous column
    if(low < y[0] || x[1] == x[0]) y[0] = low;
    if(hig > y[1] || x[1] == x[0]) y[1] = hig;

    if(x[1] > x[0])
    { // Determine direction
      int dir = getDirection(low,hig, olow,ohig);

      // no direction
      if(dir == -2) return false;

      if(x[1] > x[0]+1)
      { // Check if direction changes
        if(odir*dir == -1)
        { odir = -2; return false; }
      }

      if(x[1] <= x[0]+1 || odir == 0)
        odir = dir;

    }

    olow = low; ohig = hig;
  }
  else
  { // Very first column, initialize
    x[0] = pos.first;
  }

  // Open new column
  if(inTheLoop)
  {
    x[1] = pos.first;
    low  = pos.second;
    hig  = pos.second;
  }

  return(true);
}

Member Data Documentation

int ClusterShape::hig

private

Definition at line 38 of file ClusterShape.h.

int ClusterShape::low

private

Definition at line 38 of file ClusterShape.h.

int ClusterShape::odir

private

Definition at line 38 of file ClusterShape.h.

int ClusterShape::ohig

private

Definition at line 38 of file ClusterShape.h.

int ClusterShape::olow

private

Definition at line 38 of file ClusterShape.h.

std::vector<SiPixelCluster::Pixel> ClusterShape::pixels_

private

Definition at line 37 of file ClusterShape.h.

int ClusterShape::x[2]

private

Definition at line 38 of file ClusterShape.h.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), Vispa.Gui.WidgetContainer.WidgetContainer::autosize(), Vispa.Gui.VispaWidget.VispaWidget::boundingRect(), geometryXMLparser.Alignable::pos(), Vispa.Gui.ConnectableWidget.ConnectableWidget::positionizeMenuWidget(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

int ClusterShape::y[2]

private

Definition at line 38 of file ClusterShape.h.

Referenced by svgfig.Ellipse::__repr__(), Vispa.Gui.WidgetContainer.WidgetContainer::autosize(), Vispa.Gui.VispaWidget.VispaWidget::boundingRect(), geometryXMLparser.Alignable::pos(), Vispa.Gui.ConnectableWidget.ConnectableWidget::positionizeMenuWidget(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().