RPixChargeShare Class Reference

#include <RPixChargeShare.h>

Public Member Functions
	RPixChargeShare (const edm::ParameterSet &params, uint32_t det_id)
std::map< unsigned short, double >	Share (const std::vector< RPixSignalPoint > &charge_map)

Private Attributes
double	chargeMap2E_ [4][60][40]
std::string	ChargeMapFile2E_ [4]
uint32_t	det_id_
int	no_of_pixels_
const int	pxlColSize_ = pxlInd.getDefaultColDetSize()
CTPPSPixelIndices	pxlInd
const int	pxlRowSize_ = pxlInd.getDefaultRowDetSize()
std::vector< double >	signalCoupling_
CTPPSPixelSimTopology	theRPixDetTopology_
int	verbosity_
const int	xBinMax_ [4] = {29, 59, 29, 59}
const int	yBinMax_ [4] = {19, 19, 39, 39}

Detailed Description

Definition at line 9 of file RPixChargeShare.h.

Constructor & Destructor Documentation

RPixChargeShare()

RPixChargeShare::RPixChargeShare	(	const edm::ParameterSet &	params,
		uint32_t	det_id
	)

Definition at line 5 of file RPixChargeShare.cc.

    : det_id_(det_id), theRPixDetTopology_() {
  verbosity_ = params.getParameter<int>("RPixVerbosity");
  signalCoupling_.clear();
  ChargeMapFile2E_[0] = params.getParameter<std::string>("ChargeMapFile2E");
  ChargeMapFile2E_[1] = params.getParameter<std::string>("ChargeMapFile2E_2X");
  ChargeMapFile2E_[2] = params.getParameter<std::string>("ChargeMapFile2E_2Y");
  ChargeMapFile2E_[3] = params.getParameter<std::string>("ChargeMapFile2E_2X2Y");
 
  double coupling_constant_ = params.getParameter<double>("RPixCoupling");
 
  signalCoupling_.push_back(coupling_constant_);
  signalCoupling_.push_back((1.0 - coupling_constant_) / 2);
 
  no_of_pixels_ = theRPixDetTopology_.detPixelNo();
 
  double xMap, yMap;
  double chargeprobcollect;
  int xUpper[] = {75, 150, 75, 150};
  int yUpper[] = {50, 50, 100, 100};
  int ix, iy;
  for (int i = 0; i < 4; i++) {
    std::ifstream fChargeMap(ChargeMapFile2E_[i]);
    if (fChargeMap.is_open()) {
      while (fChargeMap >> xMap >> yMap >> chargeprobcollect) {
        ix = int((xMap + xUpper[i]) / 5);
        iy = int((yMap + yUpper[i]) / 5);
        chargeMap2E_[i][ix][iy] = chargeprobcollect;
      }
      fChargeMap.close();
    } else
      throw cms::Exception("RPixChargeShare") << "Charge map file not found";
  }
}

References chargeMap2E_, ChargeMapFile2E_, CTPPSPixelTopology::detPixelNo(), Exception, mps_fire::i, createfilelist::int, no_of_pixels_, CalibrationSummaryClient_cfi::params, signalCoupling_, AlCaHLTBitMon_QueryRunRegistry::string, theRPixDetTopology_, verbosity_, and edmStreamStallGrapher::xUpper.

Member Function Documentation

Share()

std::map< unsigned short, double > RPixChargeShare::Share

(

const std::vector< RPixSignalPoint > &

charge_map

)

Definition at line 40 of file RPixChargeShare.cc.

                                                                                                  {
  std::map<unsigned short, double> thePixelChargeMap;
  if (verbosity_ > 1)
    edm::LogInfo("RPixChargeShare") << det_id_ << " : Clouds to be induced= " << charge_map.size();
 
  double CH = 0;
 
  for (std::vector<RPixSignalPoint>::const_iterator i = charge_map.begin(); i != charge_map.end(); ++i) {
    double hit_pos_x, hit_pos_y;
    // Used to avoid the abort due to hits out of detector
    if (((*i).Position().x() + 16.6 / 2) < 0 || ((*i).Position().x() + 16.6 / 2) > 16.6) {
      edm::LogInfo("RPixChargeShare")
          << "**** Attention ((*i).Position().x()+simX_width_/2.)<0||((*i).Position().x()+simX_width_/2.)>simX_width  ";
      edm::LogInfo("RPixChargeShare") << "(*i).Position().x() = " << (*i).Position().x();
      continue;
    }
    if (((*i).Position().y() + 24.4 / 2.) < 0 || ((*i).Position().y() + 24.4 / 2.) > 24.4) {
      edm::LogInfo("RPixChargeShare")
          << "**** Attention ((*i).Position().y()+simY_width_/2.)<0||((*i).Position().y()+simY_width_/2.)>simY_width  ";
      edm::LogInfo("RPixChargeShare") << "(*i).Position().y() = " << (*i).Position().y();
      continue;
    }
 
    CTPPSPixelSimTopology::PixelInfo relevant_pixels = theRPixDetTopology_.getPixelsInvolved(
        (*i).Position().x(), (*i).Position().y(), (*i).Sigma(), hit_pos_x, hit_pos_y);
    double effic = relevant_pixels.effFactor();
 
    unsigned short pixel_no = relevant_pixels.pixelIndex();
 
    double charge_in_pixel = (*i).Charge() * effic;
 
    CH += charge_in_pixel;
 
    if (verbosity_ > 1)
      edm::LogInfo("RPixChargeShare") << "Efficiency in detector " << det_id_ << " and pixel no " << pixel_no << "  : "
                                      << effic << "  ch: " << charge_in_pixel << "   CHtot: " << CH;
 
    //        QUI SI POTREBBE INTRODURRE IL CHARGE SHARING TRA I PIXELS ..................................
 
    if (signalCoupling_[0] == 0.) {
      thePixelChargeMap[pixel_no] += charge_in_pixel;
    } else {
      int pixel_row = relevant_pixels.pixelRowNo();
      int pixel_col = relevant_pixels.pixelColNo();
      double pixel_lower_x = 0;
      double pixel_lower_y = 0;
      double pixel_upper_x = 0;
      double pixel_upper_y = 0;
      int psize = 0;
      theRPixDetTopology_.pixelRange(pixel_row, pixel_col, pixel_lower_x, pixel_upper_x, pixel_lower_y, pixel_upper_y);
      double pixel_width_x = pixel_upper_x - pixel_lower_x;
      double pixel_width_y = pixel_upper_y - pixel_lower_y;
      if (pixel_row == 0 || pixel_row == pxlRowSize_ - 1)
        pixel_width_x = 0.1;  // Correct edge pixel width
      if (pixel_col == 0 || pixel_col == pxlColSize_ - 1)
        pixel_width_y = 0.15;  //
      double pixel_center_x = pixel_lower_x + (pixel_width_x) / 2.;
      double pixel_center_y = pixel_lower_y + (pixel_width_y) / 2.;
      // xbin and ybin are coordinates (um) ??nside the pixel as in the test beam, swapped wrt plane coodinates.
      int xbin = int((((*i).Position().y() - pixel_center_y) + pixel_width_y / 2.) * 1.e3 / 5.);
      int ybin = int((((*i).Position().x() - pixel_center_x) + pixel_width_x / 2.) * 1.e3 / 5.);
      if (pixel_width_x < 0.11 && pixel_width_y < 0.151) {  // pixel 100x150 um^2
        psize = 0;
        if (xbin > xBinMax_[psize] || ybin > yBinMax_[psize])
          continue;
      }
      if (pixel_width_x > 0.11 && pixel_width_y < 0.151) {  // pixel 200x150 um^2
        psize = 2;
        if (xbin > xBinMax_[psize] || ybin > yBinMax_[psize])
          continue;
      }
      if (pixel_width_x < 0.11 && pixel_width_y > 0.151) {  // pixel 100x300 um^2
        psize = 1;
        if (xbin > xBinMax_[psize] || ybin > yBinMax_[psize])
          continue;
      }
      if (pixel_width_x > 0.11 && pixel_width_y > 0.151) {  // pixel 200x300 um^2
        psize = 3;
        if (xbin > xBinMax_[psize] || ybin > yBinMax_[psize])
          continue;
      }
      double hit2neighbour[8];
      double collect_prob = chargeMap2E_[psize][xbin][ybin];
      thePixelChargeMap[pixel_no] += charge_in_pixel * collect_prob;
      unsigned short neighbour_no[8];
      unsigned short m = 0;
      double closer_neighbour = 0;
      unsigned short closer_no = 0;
      //      Considering the 8 neighbours to share charge
      for (int k = pixel_row - 1; k <= pixel_row + 1; k++) {
        for (int l = pixel_col - 1; l <= pixel_col + 1; l++) {
          if ((k < 0) || k > pxlRowSize_ - 1 || l < 0 || l > pxlColSize_ - 1)
            continue;
          if ((k == pixel_row) && (l == pixel_col))
            continue;
          double neighbour_pixel_lower_x = 0;
          double neighbour_pixel_lower_y = 0;
          double neighbour_pixel_upper_x = 0;
          double neighbour_pixel_upper_y = 0;
          double neighbour_pixel_center_x = 0;
          double neighbour_pixel_center_y = 0;
          //      Check the hit approach to the neighbours
          theRPixDetTopology_.pixelRange(
              k, l, neighbour_pixel_lower_x, neighbour_pixel_upper_x, neighbour_pixel_lower_y, neighbour_pixel_upper_y);
          neighbour_pixel_center_x = neighbour_pixel_lower_x + (neighbour_pixel_upper_x - neighbour_pixel_lower_x) / 2.;
          neighbour_pixel_center_y = neighbour_pixel_lower_y + (neighbour_pixel_upper_y - neighbour_pixel_lower_y) / 2.;
          hit2neighbour[m] = sqrt(pow((*i).Position().x() - neighbour_pixel_center_x, 2.) +
                                  pow((*i).Position().y() - neighbour_pixel_center_y, -2.));
          neighbour_no[m] = l * pxlRowSize_ + k;
          if (hit2neighbour[m] > closer_neighbour) {
            closer_neighbour = hit2neighbour[m];
            closer_no = neighbour_no[m];
          }
          m++;
        }
      }
      double chargetransfereff = (1 - collect_prob) * signalCoupling_[0];
      thePixelChargeMap[closer_no] += charge_in_pixel * chargetransfereff;
    }
  }
 
  return thePixelChargeMap;
}

References L1DTConfigFromDB_cfi::CH, chargeMap2E_, det_id_, CTPPSPixelSimTopology::PixelInfo::effFactor(), CTPPSPixelSimTopology::getPixelsInvolved(), mps_fire::i, createfilelist::int, dqmdumpme::k, cmsLHEtoEOSManager::l, visualization-live-secondInstance_cfg::m, CTPPSPixelSimTopology::PixelInfo::pixelColNo(), CTPPSPixelSimTopology::PixelInfo::pixelIndex(), CTPPSPixelSimTopology::pixelRange(), CTPPSPixelSimTopology::PixelInfo::pixelRowNo(), funct::pow(), pxlColSize_, pxlRowSize_, signalCoupling_, mathSSE::sqrt(), theRPixDetTopology_, verbosity_, xBinMax_, and yBinMax_.

Member Data Documentation

chargeMap2E_

double RPixChargeShare::chargeMap2E_[4][60][40]

private

Definition at line 26 of file RPixChargeShare.h.

Referenced by RPixChargeShare(), and Share().

ChargeMapFile2E_

std::string RPixChargeShare::ChargeMapFile2E_[4]

private

Definition at line 25 of file RPixChargeShare.h.

Referenced by RPixChargeShare().

det_id_

uint32_t RPixChargeShare::det_id_

private

Definition at line 15 of file RPixChargeShare.h.

Referenced by Share().

no_of_pixels_

int RPixChargeShare::no_of_pixels_

private

Definition at line 21 of file RPixChargeShare.h.

Referenced by RPixChargeShare().

pxlColSize_

const int RPixChargeShare::pxlColSize_ = pxlInd.getDefaultColDetSize()

private

Definition at line 20 of file RPixChargeShare.h.

Referenced by Share().

pxlInd

CTPPSPixelIndices RPixChargeShare::pxlInd

private

Definition at line 18 of file RPixChargeShare.h.

pxlRowSize_

const int RPixChargeShare::pxlRowSize_ = pxlInd.getDefaultRowDetSize()

private

Definition at line 19 of file RPixChargeShare.h.

Referenced by Share().

signalCoupling_

std::vector<double> RPixChargeShare::signalCoupling_

private

Definition at line 16 of file RPixChargeShare.h.

Referenced by RPixChargeShare(), and Share().

theRPixDetTopology_

CTPPSPixelSimTopology RPixChargeShare::theRPixDetTopology_

private

Definition at line 17 of file RPixChargeShare.h.

Referenced by RPixChargeShare(), and Share().

verbosity_

int RPixChargeShare::verbosity_

private

Definition at line 23 of file RPixChargeShare.h.

Referenced by RPixChargeShare(), and Share().

xBinMax_

const int RPixChargeShare::xBinMax_[4] = {29, 59, 29, 59}

private

Definition at line 27 of file RPixChargeShare.h.

Referenced by Share().

yBinMax_

const int RPixChargeShare::yBinMax_[4] = {19, 19, 39, 39}

private

Definition at line 28 of file RPixChargeShare.h.

Referenced by Share().