/data/refman/pasoursint/CMSSW_5_3_3/src/DataFormats/SiPixelCluster/src/SiPixelCluster.cc

Go to the documentation of this file.

#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"

//---------------------------------------------------------------------------
//---------------------------------------------------------------------------


SiPixelCluster::SiPixelCluster( const SiPixelCluster::PixelPos& pix, int adc) :
  detId_(0),     // &&& To be fixed ?? 
  // The center of pixel with row # N is at N+0.5 in the meas. frame!
        //  theSumX( (pix.row()+0.5) * float(adc)), 
        // theSumY( (pix.col()+0.5) * float(adc)),
        //  theCharge( float(adc)),
  theMinPixelRow( pix.row()),
        //  theMaxPixelRow( pix.row()),
  theMinPixelCol( pix.col()),
        //  theMaxPixelCol( pix.col()),
  // ggiurgiu@fnal.gov, 01/05/12
  // Initialize the split cluster errors to un-physical values.
  // The CPE will check these errors and if they are not un-physical, 
  // it will recognize the clusters as split and assign these (increased) 
  // errors to the corresponding rechit. 
  err_x(-99999.9),
  err_y(-99999.9)
{
  // First pixel in this cluster.
  thePixelADC.push_back( adc );
  thePixelOffset.push_back( pix.row() - theMinPixelRow );
  thePixelOffset.push_back( pix.col() - theMinPixelCol );
}

void SiPixelCluster::add( const SiPixelCluster::PixelPos& pix, int adc) {

  // The center of pixel with row # N is at N+0.5 in the meas. frame!
  //theSumX += (pix.row()+0.5) * float(adc); 
  //theSumY += (pix.col()+0.5) * float(adc); 
  //theCharge += float(adc);
        //  thePixels.push_back( Pixel( pix.row(), pix.col(), adc ) );
        
        int minRow = theMinPixelRow;
        int minCol = theMinPixelCol;
        bool recalculate = false;
        
        if (pix.row() < theMinPixelRow) {
                theMinPixelRow = pix.row();
                recalculate = true;
        }
        if (pix.col() < theMinPixelCol) {
                theMinPixelCol = pix.col();
                recalculate = true;
        }
        if (recalculate) {
                int isize = thePixelADC.size();
                for (int i=0; i<isize; ++i) {
                        int xoffset = (thePixelOffset[i*2] ) + minRow - theMinPixelRow;
                        int yoffset = (thePixelOffset[i*2+1] ) + minCol - theMinPixelCol;
                        thePixelOffset[i*2] = xoffset;
                        thePixelOffset[i*2+1] = yoffset;
                }
        }
        
        thePixelADC.push_back( adc );
        thePixelOffset.push_back( (pix.row() - theMinPixelRow) );
        thePixelOffset.push_back( (pix.col() - theMinPixelCol) );
}