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 18 of file ClusterShape.cc.

{}

Member Function Documentation

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

Definition at line 75 of file ClusterShape.cc.

References SiPixelRecHit::cluster().

Referenced by SiPixelClusterShapeCacheProducer::produce().

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

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

Definition at line 79 of file ClusterShape.cc.

References a, cms::cuda::assert(), b, edm::VecArray< T, N >::capacity(), PVValHelper::dx, PVValHelper::dy, ClusterData::hasBigPixelsOnlyInside, mps_fire::i, 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, 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 21 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 33 of file ClusterShape.cc.

References DeadROC_duringRun::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 35 of file ClusterShape.h.

int ClusterShape::low

private

Definition at line 35 of file ClusterShape.h.

int ClusterShape::odir

private

Definition at line 35 of file ClusterShape.h.

int ClusterShape::ohig

private

Definition at line 35 of file ClusterShape.h.

int ClusterShape::olow

private

Definition at line 35 of file ClusterShape.h.

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

private

Definition at line 34 of file ClusterShape.h.

int ClusterShape::x[2]

private

Definition at line 35 of file ClusterShape.h.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

int ClusterShape::y[2]

private

Definition at line 35 of file ClusterShape.h.

Referenced by svgfig.Ellipse::__repr__(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().