CTPPSPixelSimTopology Class Reference

#include <CTPPSPixelSimTopology.h>

Inheritance diagram for CTPPSPixelSimTopology:

Classes
class	PixelInfo

Public Member Functions
	CTPPSPixelSimTopology ()
PixelInfo	getPixelsInvolved (double x, double y, double sigma, double &hit_pos_x, double &hit_pos_y) const
void	pixelRange (unsigned int arow, unsigned int acol, double &lower_x, double &higher_x, double &lower_y, double &higher_y) const
	~CTPPSPixelSimTopology ()
Public Member Functions inherited from CTPPSPixelTopology
double	activeEdgeSigma () const
	CTPPSPixelTopology ()=default
double	detDeadEdgeWidth () const
double	detPitchSimX () const
double	detPitchSimY () const
unsigned short	detPixelNo () const
unsigned short	detPixelSimXNo () const
unsigned short	detPixelSimYNo () const
double	detThickness () const
double	detXWidth () const
double	detYWidth () const
double	physActiveEdgeDist () const
	~CTPPSPixelTopology ()=default

Private Member Functions
double	activeEdgeFactor (double x, double y) const
unsigned int	col (double y) const
double	distanceFromBottomActiveEdge (double x, double y) const
double	distanceFromLeftActiveEdge (double x, double y) const
double	distanceFromRightActiveEdge (double x, double y) const
double	distanceFromTopActiveEdge (double x, double y) const
void	index2RowCol (unsigned int &arow, unsigned int &acol, unsigned int index) const
unsigned int	row (double x) const
void	rowCol2Index (unsigned int arow, unsigned int acol, unsigned int &index) const

Private Attributes
double	active_edge_x_
double	active_edge_y_

Additional Inherited Members
Static Public Member Functions inherited from CTPPSPixelTopology
static bool	isPixelHit (float xLocalCoordinate, float yLocalCoordinate, bool is3x2=true)
Public Attributes inherited from CTPPSPixelTopology
CTPPSPixelIndices	indices_
Static Public Attributes inherited from CTPPSPixelTopology
static double	active_edge_sigma_ = 0.02
static double	dead_edge_width_ = 200E-3
static unsigned short	no_of_pixels_ = 160 * 156
static unsigned short	no_of_pixels_simX_ = 160
static unsigned short	no_of_pixels_simY_ = 156
static double	phys_active_edge_dist_ = 0.150
static double	pitch_simX_ = 100E-3
static double	pitch_simY_ = 150E-3
static double	simX_width_ = 16.6
static double	simY_width_ = 24.4
static double	thickness_ = 0.23

Detailed Description

Definition at line 7 of file CTPPSPixelSimTopology.h.

Constructor & Destructor Documentation

CTPPSPixelSimTopology::CTPPSPixelSimTopology

(

)

Definition at line 3 of file CTPPSPixelSimTopology.cc.

References active_edge_x_, active_edge_y_, CTPPSPixelTopology::phys_active_edge_dist_, CTPPSPixelTopology::simX_width_, and CTPPSPixelTopology::simY_width_.

                                             {
  active_edge_x_ = simX_width_ * 0.5 - phys_active_edge_dist_;
  active_edge_y_ = simY_width_ * 0.5 - phys_active_edge_dist_;
}

CTPPSPixelSimTopology::~CTPPSPixelSimTopology ( )

inline

Definition at line 57 of file CTPPSPixelSimTopology.h.

References getPixelsInvolved(), x, and y.

{}

Member Function Documentation

double CTPPSPixelSimTopology::activeEdgeFactor ( double  x, double  y  ) const

inlineprivate

Definition at line 133 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::active_edge_sigma_, distanceFromBottomActiveEdge(), distanceFromLeftActiveEdge(), distanceFromRightActiveEdge(), distanceFromTopActiveEdge(), and Exception.

Referenced by getPixelsInvolved().

                                                           {
    const double inv_sigma = 1. / active_edge_sigma_;  // precaching
    const double topEdgeFactor = std::erf(-distanceFromTopActiveEdge(x, y) * inv_sigma) * 0.5 + 0.5;
    const double bottomEdgeFactor = std::erf(-distanceFromBottomActiveEdge(x, y) * inv_sigma) * 0.5 + 0.5;
    const double rightEdgeFactor = std::erf(-distanceFromRightActiveEdge(x, y) * inv_sigma) * 0.5 + 0.5;
    const double leftEdgeFactor = std::erf(-distanceFromLeftActiveEdge(x, y) * inv_sigma) * 0.5 + 0.5;

    const double aEF = topEdgeFactor * bottomEdgeFactor * rightEdgeFactor * leftEdgeFactor;

    if (aEF > 1.)
      throw cms::Exception("CTPPSPixelSimTopology") << " active edge factor > 1";

    return aEF;
  }

unsigned int CTPPSPixelSimTopology::col ( double  y ) const

inlineprivate

Definition at line 183 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::dead_edge_width_, Exception, createfilelist::int, CTPPSPixelTopology::pitch_simY_, and CTPPSPixelTopology::simY_width_.

Referenced by python.rootplot.root2matplotlib.Hist2D::colz(), and getPixelsInvolved().

                                          {
    // y in the G4 simulation system
    unsigned int column;

    // columns (y segmentation)
    // now y in the system centered in the bottom left corner of the sensor (sensor view, rocs behind)
    y = y + simY_width_ / 2.;
    if (y < 0. || y > simY_width_)
      throw cms::Exception("CTPPSPixelSimTopology") << " out of reference frame";

    if (y <= (dead_edge_width_ + pitch_simY_))
      column = 0;
    else if (y <= (dead_edge_width_ + (ROCSizeInY - 1) * pitch_simY_))
      column = int((y - dead_edge_width_ - pitch_simY_) / pitch_simY_) + 1;
    else if (y <= (dead_edge_width_ + (ROCSizeInY + 1) * pitch_simY_))
      column = ROCSizeInY - 1;
    else if (y <= (dead_edge_width_ + (ROCSizeInY + 3) * pitch_simY_))
      column = ROCSizeInY;
    else if (y <= (dead_edge_width_ + (2 * ROCSizeInY + 1) * pitch_simY_))
      column = int((y - dead_edge_width_ - pitch_simY_) / pitch_simY_) - 1;
    else if (y <= (dead_edge_width_ + (2 * ROCSizeInY + 3) * pitch_simY_))
      column = 2 * ROCSizeInY - 1;
    else if (y <= (dead_edge_width_ + (2 * ROCSizeInY + 5) * pitch_simY_))
      column = 2 * ROCSizeInY;
    else if (y <= (dead_edge_width_ + (3 * ROCSizeInY + 3) * pitch_simY_))
      column = int((y - dead_edge_width_ - pitch_simY_) / pitch_simY_) - 3;
    else
      column = (3 * ROCSizeInY - 1);

    return column;
  }

double CTPPSPixelSimTopology::distanceFromBottomActiveEdge ( double  x, double  y  ) const

inlineprivate

Definition at line 149 of file CTPPSPixelSimTopology.h.

Referenced by activeEdgeFactor().

{ return (-y - active_edge_y_); }

double CTPPSPixelSimTopology::distanceFromLeftActiveEdge ( double  x, double  y  ) const

inlineprivate

Definition at line 151 of file CTPPSPixelSimTopology.h.

Referenced by activeEdgeFactor().

{ return (-x - active_edge_x_); }

double CTPPSPixelSimTopology::distanceFromRightActiveEdge ( double  x, double  y  ) const

inlineprivate

Definition at line 150 of file CTPPSPixelSimTopology.h.

Referenced by activeEdgeFactor().

{ return (x - active_edge_x_); }

double CTPPSPixelSimTopology::distanceFromTopActiveEdge ( double  x, double  y  ) const

inlineprivate

Definition at line 148 of file CTPPSPixelSimTopology.h.

Referenced by activeEdgeFactor().

{ return (y - active_edge_y_); }

CTPPSPixelSimTopology::PixelInfo CTPPSPixelSimTopology::getPixelsInvolved	(	double	x,
		double	y,
		double	sigma,
		double &	hit_pos_x,
		double &	hit_pos_y
	)		const

Definition at line 8 of file CTPPSPixelSimTopology.cc.

References activeEdgeFactor(), col(), Exception, pixelRange(), row(), CTPPSPixelTopology::simX_width_, and CTPPSPixelTopology::simY_width_.

Referenced by RPixChargeShare::Share(), and ~CTPPSPixelSimTopology().

                                                                                  {
  //hit position wrt the bottom left corner of the sensor (-8.3, -12.2) in sensor view, rocs behind
  hit_pos_x = x + simX_width_ / 2.;
  hit_pos_y = y + simY_width_ / 2.;
  if (!(hit_pos_x * hit_pos_y > 0))
    throw cms::Exception("CTPPSPixelSimTopology") << "out of reference frame";

  double hit_factor = activeEdgeFactor(x, y);

  unsigned int interested_row = row(x);
  unsigned int interested_col = col(y);
  double low_pixel_range_x, high_pixel_range_x, low_pixel_range_y, high_pixel_range_y;
  pixelRange(
      interested_row, interested_col, low_pixel_range_x, high_pixel_range_x, low_pixel_range_y, high_pixel_range_y);

  return CTPPSPixelSimTopology::PixelInfo(low_pixel_range_x,
                                          high_pixel_range_x,
                                          low_pixel_range_y,
                                          high_pixel_range_y,
                                          hit_factor,
                                          interested_row,
                                          interested_col);
}

void CTPPSPixelSimTopology::index2RowCol ( unsigned int &  arow, unsigned int &  acol, unsigned int  index  ) const

inlineprivate

Definition at line 219 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::no_of_pixels_simX_.

                                                                                             {
    acol = index / no_of_pixels_simX_;
    arow = index % no_of_pixels_simX_;
  }

void CTPPSPixelSimTopology::pixelRange ( unsigned int  arow, unsigned int  acol, double &  lower_x, double &  higher_x, double &  lower_y, double &  higher_y  ) const

inline

Definition at line 61 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::dead_edge_width_, Exception, CTPPSPixelTopology::phys_active_edge_dist_, CTPPSPixelTopology::pitch_simX_, CTPPSPixelTopology::pitch_simY_, CTPPSPixelTopology::simX_width_, and CTPPSPixelTopology::simY_width_.

Referenced by PPSPixelDigiAnalyzer::analyze(), getPixelsInvolved(), RPixClusterToHit::make_hit(), and RPixChargeShare::Share().

                                                 {
    // x and y in the system  of Geant4 SIMULATION
    arow = (2 * ROCSizeInX - 1) - arow;
    if (arow > (2 * ROCSizeInX - 1) || acol > (3 * ROCSizeInY - 1))
      throw cms::Exception("CTPPSPixelSimTopology") << "rows or columns exceeding limits";

    // rows (x segmentation)
    if (arow == 0) {
      lower_x = dead_edge_width_ - phys_active_edge_dist_;  // 50 um
      higher_x = dead_edge_width_ + pitch_simX_;            // 300 um
    } else if (arow <= (ROCSizeInX - 2)) {
      lower_x = dead_edge_width_ + arow * pitch_simX_;
      higher_x = dead_edge_width_ + (arow + 1) * pitch_simX_;
    } else if (arow == (ROCSizeInX - 1)) {
      lower_x = dead_edge_width_ + arow * pitch_simX_;
      higher_x = dead_edge_width_ + (arow + 2) * pitch_simX_;
    } else if (arow == ROCSizeInX) {
      lower_x = dead_edge_width_ + (arow + 1) * pitch_simX_;
      higher_x = dead_edge_width_ + (arow + 3) * pitch_simX_;
    } else if (arow <= (2 * ROCSizeInX - 2)) {
      lower_x = dead_edge_width_ + (arow + 2) * pitch_simX_;
      higher_x = dead_edge_width_ + (arow + 3) * pitch_simX_;
    } else if (arow == (2 * ROCSizeInX - 1)) {
      lower_x = dead_edge_width_ + (arow + 2) * pitch_simX_;
      higher_x = dead_edge_width_ + (arow + 3) * pitch_simX_ + phys_active_edge_dist_;
    }

    // columns (y segmentation)
    if (acol == 0) {
      lower_y = dead_edge_width_ - phys_active_edge_dist_;  // 50 um
      higher_y = dead_edge_width_ + pitch_simY_;            // 350 um
    } else if (acol <= (ROCSizeInY - 2)) {
      lower_y = dead_edge_width_ + acol * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 1) * pitch_simY_;
    } else if (acol == (ROCSizeInY - 1)) {
      lower_y = dead_edge_width_ + acol * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 2) * pitch_simY_;
    } else if (acol == ROCSizeInY) {
      lower_y = dead_edge_width_ + (acol + 1) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 3) * pitch_simY_;
    } else if (acol <= (2 * ROCSizeInY - 2)) {
      lower_y = dead_edge_width_ + (acol + 2) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 3) * pitch_simY_;
    } else if (acol == (2 * ROCSizeInY - 1)) {
      lower_y = dead_edge_width_ + (acol + 2) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 4) * pitch_simY_;
    } else if (acol == (2 * ROCSizeInY)) {
      lower_y = dead_edge_width_ + (acol + 3) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 5) * pitch_simY_;
    } else if (acol <= (3 * ROCSizeInY - 2)) {
      lower_y = dead_edge_width_ + (acol + 4) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 5) * pitch_simY_;
    } else if (acol == (3 * ROCSizeInY - 1)) {
      lower_y = dead_edge_width_ + (acol + 4) * pitch_simY_;
      higher_y = dead_edge_width_ + (acol + 5) * pitch_simY_ + phys_active_edge_dist_;
    }

    lower_x = lower_x - simX_width_ / 2.;
    lower_y = lower_y - simY_width_ / 2.;
    higher_x = higher_x - simX_width_ / 2.;
    higher_y = higher_y - simY_width_ / 2.;
  }

unsigned int CTPPSPixelSimTopology::row ( double  x ) const

inlineprivate

Definition at line 153 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::dead_edge_width_, Exception, createfilelist::int, CTPPSPixelTopology::pitch_simX_, and CTPPSPixelTopology::simX_width_.

Referenced by getPixelsInvolved().

                                          {
    // x in the G4 simulation system
    x = x + simX_width_ / 2.;

    // now x in the system centered in the bottom left corner of the sensor (sensor view, rocs behind)
    if (x < 0. || x > simX_width_)
      throw cms::Exception("CTPPSPixelSimTopology") << "out of reference frame";

    // rows (x segmentation)
    unsigned int arow;
    if (x <= (dead_edge_width_ + pitch_simX_))
      arow = 0;
    else if (x <= (dead_edge_width_ + (ROCSizeInX - 1) * pitch_simX_))
      arow = int((x - dead_edge_width_ - pitch_simX_) / pitch_simX_) + 1;
    else if (x <= (dead_edge_width_ + (ROCSizeInX + 1) * pitch_simX_))
      arow = (ROCSizeInX - 1);
    else if (x <= (dead_edge_width_ + (ROCSizeInX + 3) * pitch_simX_))
      arow = ROCSizeInX;
    else if (x <= (dead_edge_width_ + (2 * ROCSizeInX + 2) * pitch_simX_))
      arow = int((x - dead_edge_width_ - pitch_simX_) / pitch_simX_) - 1;
    else
      arow = (2 * ROCSizeInX - 1);

    arow = (2 * ROCSizeInX - 1) - arow;
    if (arow > (2 * ROCSizeInX - 1))
      throw cms::Exception("CTPPSPixelSimTopology") << "row number exceeding limit";

    return arow;
  }

void CTPPSPixelSimTopology::rowCol2Index ( unsigned int  arow, unsigned int  acol, unsigned int &  index  ) const

inlineprivate

Definition at line 215 of file CTPPSPixelSimTopology.h.

References CTPPSPixelTopology::no_of_pixels_simX_.

                                                                                            {
    index = acol * no_of_pixels_simX_ + arow;
  }

Member Data Documentation

double CTPPSPixelSimTopology::active_edge_x_

private

Definition at line 130 of file CTPPSPixelSimTopology.h.

Referenced by CTPPSPixelSimTopology().

double CTPPSPixelSimTopology::active_edge_y_

private

Definition at line 131 of file CTPPSPixelSimTopology.h.

Referenced by CTPPSPixelSimTopology().