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 constexpr double	active_edge_sigma_ = 0.02
static constexpr double	dead_edge_width_ = 200E-3
static constexpr unsigned short	no_of_pixels_ = 160 * 156
static constexpr unsigned short	no_of_pixels_simX_ = 160
static constexpr unsigned short	no_of_pixels_simY_ = 156
static constexpr double	phys_active_edge_dist_ = 0.150
static constexpr double	pitch_simX_ = 100E-3
static constexpr double	pitch_simY_ = 150E-3
static constexpr double	simX_width_ = 16.6
static constexpr double	simY_width_ = 24.4
static constexpr double	thickness_ = 0.23

Detailed Description

Definition at line 7 of file CTPPSPixelSimTopology.h.

Constructor & Destructor Documentation

CTPPSPixelSimTopology()

CTPPSPixelSimTopology::CTPPSPixelSimTopology

(

)

Definition at line 3 of file CTPPSPixelSimTopology.cc.

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

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

~CTPPSPixelSimTopology()

CTPPSPixelSimTopology::~CTPPSPixelSimTopology ( )

inline

Definition at line 57 of file CTPPSPixelSimTopology.h.

{}

Member Function Documentation

activeEdgeFactor()

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

inlineprivate

Definition at line 133 of file CTPPSPixelSimTopology.h.

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

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

Referenced by getPixelsInvolved().

col()

unsigned int CTPPSPixelSimTopology::col ( double  y ) const

inlineprivate

Definition at line 183 of file CTPPSPixelSimTopology.h.

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

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

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

distanceFromBottomActiveEdge()

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

inlineprivate

Definition at line 149 of file CTPPSPixelSimTopology.h.

{ return (-y - active_edge_y_); }

References active_edge_y_, and y.

Referenced by activeEdgeFactor().

distanceFromLeftActiveEdge()

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

inlineprivate

Definition at line 151 of file CTPPSPixelSimTopology.h.

{ return (-x - active_edge_x_); }

References active_edge_x_, and x.

Referenced by activeEdgeFactor().

distanceFromRightActiveEdge()

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

inlineprivate

Definition at line 150 of file CTPPSPixelSimTopology.h.

{ return (x - active_edge_x_); }

References active_edge_x_, and x.

Referenced by activeEdgeFactor().

distanceFromTopActiveEdge()

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

inlineprivate

Definition at line 148 of file CTPPSPixelSimTopology.h.

{ return (y - active_edge_y_); }

References active_edge_y_, and y.

Referenced by activeEdgeFactor().

getPixelsInvolved()

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.

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

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

Referenced by RPixChargeShare::Share().

index2RowCol()

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

inlineprivate

Definition at line 219 of file CTPPSPixelSimTopology.h.

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

References CTPPSPixelTopology::no_of_pixels_simX_.

pixelRange()

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.

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

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(), and RPixChargeShare::Share().

row()

unsigned int CTPPSPixelSimTopology::row ( double  x ) const

inlineprivate

Definition at line 153 of file CTPPSPixelSimTopology.h.

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

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

Referenced by getPixelsInvolved().

rowCol2Index()

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

inlineprivate

Definition at line 215 of file CTPPSPixelSimTopology.h.

                                                                                            {
    index = acol * no_of_pixels_simX_ + arow;
  }

References CTPPSPixelTopology::no_of_pixels_simX_.

Member Data Documentation

active_edge_x_

double CTPPSPixelSimTopology::active_edge_x_

private

Definition at line 130 of file CTPPSPixelSimTopology.h.

Referenced by CTPPSPixelSimTopology(), distanceFromLeftActiveEdge(), and distanceFromRightActiveEdge().

active_edge_y_

double CTPPSPixelSimTopology::active_edge_y_

private

Definition at line 131 of file CTPPSPixelSimTopology.h.

Referenced by CTPPSPixelSimTopology(), distanceFromBottomActiveEdge(), and distanceFromTopActiveEdge().