RPixClusterToHit Class Reference

#include <RPixClusterToHit.h>

Public Member Functions
void	buildHits (unsigned int detId, const std::vector< CTPPSPixelCluster > &clusters, std::vector< CTPPSPixelRecHit > &hits)
void	make_hit (CTPPSPixelCluster aCluster, std::vector< CTPPSPixelRecHit > &hits)
	RPixClusterToHit (edm::ParameterSet const &conf)
	~RPixClusterToHit ()

Private Attributes
int	verbosity_

Detailed Description

Definition at line 15 of file RPixClusterToHit.h.

Constructor & Destructor Documentation

RPixClusterToHit::RPixClusterToHit

(

edm::ParameterSet const &

conf

)

Definition at line 4 of file RPixClusterToHit.cc.

References edm::ParameterSet::getUntrackedParameter(), and verbosity_.

{
verbosity_ = conf.getUntrackedParameter<int>("RPixVerbosity");
}

RPixClusterToHit::~RPixClusterToHit

(

)

Definition at line 9 of file RPixClusterToHit.cc.

{}

Member Function Documentation

void RPixClusterToHit::buildHits	(	unsigned int	detId,
		const std::vector< CTPPSPixelCluster > &	clusters,
		std::vector< CTPPSPixelRecHit > &	hits
	)

Definition at line 11 of file RPixClusterToHit.cc.

References mps_fire::i, make_hit(), and verbosity_.

Referenced by CTPPSPixelRecHitProducer::run().

{
  
  if(verbosity_) edm::LogInfo("RPixClusterToHit")<<" RPixClusterToHit "<<detId<<" received cluster array of size = "<<clusters.size();
   for(unsigned int i=0; i<clusters.size();i++){
     make_hit(clusters[i],hits);
  }

}

void RPixClusterToHit::make_hit	(	CTPPSPixelCluster	aCluster,
		std::vector< CTPPSPixelRecHit > &	hits
	)

Definition at line 22 of file RPixClusterToHit.cc.

References CTPPSPixelCluster::colSpan(), CTPPSPixelIndices::getDefaultColDetSize(), CTPPSPixelIndices::getDefaultRowDetSize(), CTPPSPixelIndices::getROCId(), mps_fire::i, createfilelist::int, CTPPSPixelCluster::minPixelCol(), CTPPSPixelCluster::minPixelRow(), CTPPSPixelCluster::pixelADC(), CTPPSPixelCluster::pixelCol(), CTPPSPixelSimTopology::pixelRange(), CTPPSPixelCluster::pixelRow(), CTPPSPixelCluster::rowSpan(), CTPPSPixelCluster::size(), CTPPSPixelCluster::sizeCol(), CTPPSPixelCluster::sizeRow(), ecaldqm::topology(), verbosity_, and create_public_lumi_plots::weights.

Referenced by buildHits().

                                                                                              {

// take a cluster, generate a rec hit and push it in the rec hit vector

//call the topology
  CTPPSPixelSimTopology topology;
//call the numbering inside the ROC
  CTPPSPixelIndices pxlInd;
// get information from the cluster 
// get the whole cluster size and row/col size
  unsigned int thisClusterSize = aCluster.size();
  unsigned int thisClusterRowSize = aCluster.sizeRow();
  unsigned int thisClusterColSize = aCluster.sizeCol();

// get the minimum pixel row/col
  unsigned int thisClusterMinRow = aCluster.minPixelRow();
  unsigned int thisClusterMinCol = aCluster.minPixelCol();

// calculate "on edge" flag
  bool anEdgePixel = false;
  if(aCluster.minPixelRow() == 0 || aCluster.minPixelCol() == 0 ||  
     int(aCluster.minPixelRow()+aCluster.rowSpan()) == (pxlInd.getDefaultRowDetSize()-1) ||  
     int(aCluster.minPixelCol()+aCluster.colSpan()) == (pxlInd.getDefaultColDetSize()-1) ) 
    anEdgePixel = true;

// check for bad (ADC=0) pixels in cluster 
  bool aBadPixel = false;
  for(unsigned int i = 0; i < thisClusterSize; i++){
    if(aCluster.pixelADC(i)==0) aBadPixel = true;
  } 

// check for spanning two ROCs
  bool twoRocs = false;
  int currROCId = pxlInd.getROCId(aCluster.pixelCol(0),aCluster.pixelRow(0) ); 

  for(unsigned int i = 1; i < thisClusterSize; i++){
    if(pxlInd.getROCId(aCluster.pixelCol(i),aCluster.pixelRow(i) ) != currROCId){
      twoRocs = true;
      break;
    }
  }

//estimate position and error of the hit
  double avgWLocalX = 0;
  double avgWLocalY = 0;
  double weights = 0;
  double weightedVarianceX = 0.;
  double weightedVarianceY = 0.;

  if(verbosity_)
    edm::LogInfo("RPixClusterToHit") << " hit pixels: ";

  for(unsigned int i = 0; i < thisClusterSize; i++){
    
    if(verbosity_)edm::LogInfo("RPixClusterToHit") << aCluster.pixelRow(i) << " " << aCluster.pixelCol(i)<<" " << aCluster.pixelADC(i);

    double minPxlX = 0;
    double minPxlY = 0;
    double maxPxlX = 0;
    double maxPxlY = 0;
    topology.pixelRange(aCluster.pixelRow(i),aCluster.pixelCol(i),minPxlX,maxPxlX,minPxlY, maxPxlY);
    double halfSizeX = (maxPxlX-minPxlX)/2.;
    double halfSizeY = (maxPxlY-minPxlY)/2.;
    double avgPxlX = minPxlX + halfSizeX;
    double avgPxlY = minPxlY + halfSizeY;
//error propagation
    weightedVarianceX += aCluster.pixelADC(i)*aCluster.pixelADC(i)*halfSizeX*halfSizeX/3.;
    weightedVarianceY += aCluster.pixelADC(i)*aCluster.pixelADC(i)*halfSizeY*halfSizeY/3.;
    
    avgWLocalX += avgPxlX*aCluster.pixelADC(i);
    avgWLocalY += avgPxlY*aCluster.pixelADC(i);
    weights += aCluster.pixelADC(i);

  } 

  if(weights == 0){
    edm::LogError("RPixClusterToHit") << " unexpected weights = 0 for cluster (Row_min, Row_max, Col_min, Col_max) = ("
                      << aCluster.minPixelRow() << ","
                      << aCluster.minPixelRow()+aCluster.rowSpan() << ","
                      << aCluster.minPixelCol() << ","
                      << aCluster.minPixelCol()+aCluster.colSpan()
                      << ")"; 
    return;
  }

  double invWeights = 1./weights;
  double avgLocalX = avgWLocalX*invWeights;
  double avgLocalY = avgWLocalY*invWeights;
  
  double varianceX = weightedVarianceX*invWeights*invWeights;
  double varianceY = weightedVarianceY*invWeights*invWeights;

  LocalPoint lp(avgLocalX,avgLocalY,0);
  LocalError le(varianceX,0,varianceY);
  CTPPSPixelRecHit rh(lp,le,anEdgePixel,aBadPixel,twoRocs,thisClusterMinRow,thisClusterMinCol,thisClusterSize,thisClusterRowSize,thisClusterColSize);
  if(verbosity_)edm::LogInfo("RPixClusterToHit") << lp << " with error " << le ;

  hits.push_back(rh);

  return;

}

Member Data Documentation

int RPixClusterToHit::verbosity_

private

Definition at line 27 of file RPixClusterToHit.h.

Referenced by buildHits(), make_hit(), and RPixClusterToHit().