#include <PixelCPETemplateReco.h>

Inheritance diagram for PixelCPETemplateReco:

Public Member Functions
LocalError	localError (const SiPixelCluster &cl, const GeomDetUnit &det) const

LocalPoint	localPosition (const SiPixelCluster &cluster, const GeomDetUnit &det) const

MeasurementPoint	measurementPosition (const SiPixelCluster &, const GeomDetUnit &det) const

	PixelCPETemplateReco (edm::ParameterSet const &conf, const MagneticField , const SiPixelLorentzAngle , const SiPixelTemplateDBObject *)

	~PixelCPETemplateReco ()

Public Member Functions inherited from PixelCPEBase
unsigned int	clusterProbComputationFlag () const

bool	hasBadPixels () const

bool	hasFilledProb () const

bool	isOnEdge () const

LocalValues	localParameters (const SiPixelCluster &cl, const GeomDetUnit &det) const

LocalValues	localParameters (const SiPixelCluster &cl, const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const

LocalValues	localParameters (const SiPixelCluster &cl, const GeomDetUnit &det, float alpha, float beta) const

	PixelCPEBase (edm::ParameterSet const &conf, const MagneticField mag=0, const SiPixelLorentzAngle lorentzAngle=0, const SiPixelCPEGenericErrorParm genErrorParm=0, const SiPixelTemplateDBObject templateDBobject=0)

float	probabilityQ () const

float	probabilityX () const

float	probabilityXY () const

float	probabilityY () const

float	qBin () const

SiPixelRecHitQuality::QualWordType	rawQualityWord () const

void	setMagField (const MagneticField *mag) const

bool	spansTwoRocks () const

Public Member Functions inherited from PixelClusterParameterEstimator
unsigned int	clusterProbComputationFlag () const

	PixelClusterParameterEstimator ()

Public Member Functions inherited from ClusterParameterEstimator< SiPixelCluster >
virtual void	clearParameters () const

virtual void	enterLocalParameters (unsigned int id, std::pair< int, int > &row_col, LocalValues pos_err_info) const

virtual void	enterLocalParameters (uint32_t id, uint16_t firstStrip, LocalValues pos_err_info) const

virtual LocalValues	localParameters (const SiPixelCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual VLocalValues	localParametersV (const SiPixelCluster &cluster, const GeomDetUnit &gd) const

virtual VLocalValues	localParametersV (const SiPixelCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual	~ClusterParameterEstimator ()

Protected Member Functions
float	xpos (const SiPixelCluster &) const

float	ypos (const SiPixelCluster &) const

Protected Member Functions inherited from PixelCPEBase
void	computeAnglesFromDetPosition (const SiPixelCluster &cl, const GeomDetUnit &det) const

void	computeAnglesFromTrajectory (const SiPixelCluster &cl, const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const

void	computeLorentzShifts () const

LocalVector	driftDirection (GlobalVector bfield) const

LocalVector	driftDirectionCorrect (GlobalVector bfield) const

bool	isFlipped () const

virtual float	lorentzShiftX () const

virtual float	lorentzShiftY () const

MeasurementError	measurementError (const SiPixelCluster &, const GeomDetUnit &det) const

MeasurementPoint	measurementPosition (const SiPixelCluster &, const GeomDetUnit &det) const

void	setTheDet (const GeomDetUnit &det, const SiPixelCluster &cluster) const

void	xCharge (const std::vector< SiPixelCluster::Pixel > &, const int &, const int &, float &q1, float &q2) const

void	yCharge (const std::vector< SiPixelCluster::Pixel > &, const int &, const int &, float &q1, float &q2) const

Private Attributes
bool	DoCosmics_

int	ierr

bool	LoadTemplatesFromDB_

int	speed_

SiPixelTemplate	templ_

int	templID_

float	templProbQ_

float	templProbX_

float	templProbY_

int	templQbin_

float	templSigmaX_

float	templSigmaY_

float	templXrec_

float	templYrec_

bool	UseClusterSplitter_

Additional Inherited Members
Public Types inherited from ClusterParameterEstimator< SiPixelCluster >
typedef std::pair< LocalPoint, LocalError >	LocalValues

typedef std::vector< LocalValues >	VLocalValues

Protected Types inherited from PixelCPEBase
typedef GloballyPositioned < double >	Frame

Protected Attributes inherited from PixelCPEBase
bool	alpha2Order

float	alpha_

float	beta_

unsigned int	clusterProbComputationFlag_

float	cotalpha_

float	cotbeta_

LocalVector	driftDirection_

const SiPixelCPEGenericErrorParm *	genErrorParm_

bool	hasBadPixels_

bool	hasFilledProb_

bool	isOnEdge_

LocalTrajectoryParameters	loc_traj_param_

Topology::LocalTrackPred	loc_trk_pred_

const SiPixelLorentzAngle *	lorentzAngle_

double	lorentzShiftInCmX_

double	lorentzShiftInCmY_

double	lorentzShiftX_

double	lorentzShiftY_

const MagneticField *	magfield_

int	nRecHitsTotal_

int	nRecHitsUsedEdge_

float	probabilityQ_

float	probabilityX_

float	probabilityY_

float	qBin_

bool	spansTwoROCs_

const SiPixelTemplateDBObject *	templateDBobject_

const PixelGeomDetUnit *	theDet

float	theDetR

float	theDetZ

float	theLShiftX

float	theLShiftY

float	theNumOfCol

float	theNumOfRow

GeomDetType::SubDetector	thePart

float	thePitchX

float	thePitchY

float	theSign

float	theThickness

const PixelTopology *	theTopol

int	theVerboseLevel

float	trk_lp_x

float	trk_lp_y

bool	with_track_angle

Protected Attributes inherited from PixelClusterParameterEstimator
unsigned int	clusterProbComputationFlag_

Detailed Description

Definition at line 34 of file PixelCPETemplateReco.h.

Constructor & Destructor Documentation

PixelCPETemplateReco::PixelCPETemplateReco	(	edm::ParameterSet const &	conf,
		const MagneticField *	mag,
		const SiPixelLorentzAngle *	lorentzAngle,
		const SiPixelTemplateDBObject *	templateDBobject
	)

Definition at line 46 of file PixelCPETemplateReco.cc.

References DoCosmics_, edm::ParameterSet::getParameter(), LoadTemplatesFromDB_, LogDebug, SiPixelTemplate::pushfile(), speed_, templ_, PixelCPEBase::templateDBobject_, and UseClusterSplitter_.

   : PixelCPEBase(conf, mag, lorentzAngle, 0, templateDBobject)
 {
   //cout << endl;
   //cout << "Constructing PixelCPETemplateReco::PixelCPETemplateReco(...)................................................." << endl;
   //cout << endl;
 
   // &&& initialize the templates, etc.
   
   //-- Use Magnetic field at (0,0,0) to select a template ID [Morris, 6/25/08] (temporary until we implement DB access)
   
   DoCosmics_ = conf.getParameter<bool>("DoCosmics");
   LoadTemplatesFromDB_ = conf.getParameter<bool>("LoadTemplatesFromDB");
   //cout << " PixelCPETemplateReco : (int)LoadTemplatesFromDB_ = " << (int)LoadTemplatesFromDB_ << endl;
   //cout << "field_magnitude = " << field_magnitude << endl;
   
   // ggiurgiu@fnal.gov, 12/17/2008: use configuration parameter to decide between DB or text file template access
   if ( LoadTemplatesFromDB_ )
     {
       //cout << "PixelCPETemplateReco: Loading templates from database (DB) ------------------------------- " << endl;
       
       // Initialize template store to the selected ID [Morris, 6/25/08]  
       if ( !templ_.pushfile( *templateDBobject_) )
     throw cms::Exception("PixelCPETemplateReco") 
       << "\nERROR: Templates not filled correctly. Check the sqlite file. Using SiPixelTemplateDBObject version " 
       << (*templateDBobject_).version() << "\n\n";
     }
   else 
     {
       //cout << "PixelCPETemplateReco : Loading templates 40 and 41 from ASCII files ------------------------" << endl;
 
       if ( !templ_.pushfile( 40 ) )
     throw cms::Exception("PixelCPETemplateReco") 
       << "\nERROR: Templates 40 not loaded correctly from text file. Reconstruction will fail.\n\n";
     
       if ( !templ_.pushfile( 41 ) )
     throw cms::Exception("PixelCPETemplateReco") 
       << "\nERROR: Templates 41 not loaded correctly from text file. Reconstruction will fail.\n\n";
     }
   
   speed_ = conf.getParameter<int>( "speed");
   LogDebug("PixelCPETemplateReco::PixelCPETemplateReco:") <<
     "Template speed = " << speed_ << "\n";
   
   UseClusterSplitter_ = conf.getParameter<bool>("UseClusterSplitter");
   
 }

PixelCPETemplateReco::~PixelCPETemplateReco ( )

Definition at line 99 of file PixelCPETemplateReco.cc.

 {
   // &&& delete template store?
 }

Member Function Documentation

LocalError PixelCPETemplateReco::localError	(	const SiPixelCluster &	cl,
		const GeomDetUnit &	det
	)		const

virtual

Implements PixelCPEBase.

Definition at line 495 of file PixelCPETemplateReco.cc.

References edm::hlt::Exception, SiPixelCluster::getSplitClusterErrorX(), SiPixelCluster::getSplitClusterErrorY(), ierr, PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), LogDebug, SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), micronsToCm, SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, PixelCPEBase::setTheDet(), SiPixelCluster::sizeX(), SiPixelCluster::sizeY(), mathSSE::sqrt(), templSigmaX_, templSigmaY_, PixelCPEBase::thePart, PixelCPEBase::thePitchX, PixelCPEBase::thePitchY, PixelCPEBase::theTopol, and PixelCPEBase::theVerboseLevel.

 {
   //cout << endl;
   //cout << "Set PixelCPETemplate errors .............................................." << endl;
 
   //cout << "CPETemplate : " << endl;
 
   //--- Default is the maximum error used for edge clusters.
   float xerr = thePitchX / sqrt(12.0);
   float yerr = thePitchY / sqrt(12.0);
   
   // Check if the errors were already set at the clusters splitting level
   if ( cluster.getSplitClusterErrorX() > 0.0 && cluster.getSplitClusterErrorX() < 7777.7 && 
        cluster.getSplitClusterErrorY() > 0.0 && cluster.getSplitClusterErrorY() < 7777.7 )
     {
       xerr = cluster.getSplitClusterErrorX() * micronsToCm;
       yerr = cluster.getSplitClusterErrorY() * micronsToCm;
 
       //cout << "Errors set at cluster splitting level : " << endl;
       //cout << "xerr = " << xerr << endl;
       //cout << "yerr = " << yerr << endl;
     }
   else
     {
       // If errors are not split at the cluster splitting level, set the errors here
 
       //cout  << "Errors are not split at the cluster splitting level, set the errors here : " << endl; 
   
       setTheDet( det, cluster );
       
       int maxPixelCol = cluster.maxPixelCol();
       int maxPixelRow = cluster.maxPixelRow();
       int minPixelCol = cluster.minPixelCol();
       int minPixelRow = cluster.minPixelRow();
       
       //--- Are we near either of the edges?
       bool edgex = ( theTopol->isItEdgePixelInX( minPixelRow ) || theTopol->isItEdgePixelInX( maxPixelRow ) );
       bool edgey = ( theTopol->isItEdgePixelInY( minPixelCol ) || theTopol->isItEdgePixelInY( maxPixelCol ) );
       
       if ( ierr !=0 ) 
     {
       // If reconstruction fails the hit position is calculated from cluster center of gravity 
       // corrected in x by average Lorentz drift. Assign huge errors.
       //xerr = 10.0 * (float)cluster.sizeX() * xerr;
       //yerr = 10.0 * (float)cluster.sizeX() * yerr;
       
       // Assign better errors based on the residuals for failed template cases
       if ( thePart == GeomDetEnumerators::PixelBarrel )
         {
           xerr = 55.0 * micronsToCm;
           yerr = 36.0 * micronsToCm;
         }
       else if ( thePart == GeomDetEnumerators::PixelEndcap )
         {
           xerr = 42.0 * micronsToCm;
           yerr = 39.0 * micronsToCm;
         }
       else 
         throw cms::Exception("PixelCPETemplateReco::localError :") << "A non-pixel detector type in here?" ;
 
       //cout << "xerr = " << xerr << endl;
       //cout << "yerr = " << yerr << endl;
       
       //return LocalError(xerr*xerr, 0, yerr*yerr);
     }
       else if ( edgex || edgey )
     {
       // for edge pixels assign errors according to observed residual RMS 
       if      ( edgex && !edgey )
         {
           xerr = 23.0 * micronsToCm;
           yerr = 39.0 * micronsToCm;
         }
       else if ( !edgex && edgey )
         {
           xerr = 24.0 * micronsToCm;
           yerr = 96.0 * micronsToCm;
         }
       else if ( edgex && edgey )
         {
           xerr = 31.0 * micronsToCm;
           yerr = 90.0 * micronsToCm;
         }
       else
         {
           throw cms::Exception(" PixelCPETemplateReco::localError: Something wrong with pixel edge flag !!!");
         }
 
       //cout << "xerr = " << xerr << endl;
       //cout << "yerr = " << yerr << endl;
     }
       else 
     {
       // &&& need a class const
       //const float micronsToCm = 1.0e-4;
       
       xerr = templSigmaX_ * micronsToCm;
       yerr = templSigmaY_ * micronsToCm;
       
       //cout << "xerr = " << xerr << endl;
       //cout << "yerr = " << yerr << endl;
 
       // &&& should also check ierr (saved as class variable) and return
       // &&& nonsense (another class static) if the template fit failed.
     }       
       
       if (theVerboseLevel > 9) 
     {
       LogDebug("PixelCPETemplateReco") <<
         " Sizex = " << cluster.sizeX() << " Sizey = " << cluster.sizeY() << " Edgex = " << edgex << " Edgey = " << edgey << 
         " ErrX  = " << xerr            << " ErrY  = " << yerr;
     }
 
     } // else
   
   if ( !(xerr > 0.0) )
     throw cms::Exception("PixelCPETemplateReco::localError") 
       << "\nERROR: Negative pixel error xerr = " << xerr << "\n\n";
   
   if ( !(yerr > 0.0) )
     throw cms::Exception("PixelCPETemplateReco::localError") 
       << "\nERROR: Negative pixel error yerr = " << yerr << "\n\n";
 
   //cout << "Final errors set to: " << endl;
   //cout << "xerr = " << xerr << endl;
   //cout << "yerr = " << yerr << endl;
   //cout << "Out of PixelCPETemplateREco..........................................................................." << endl;
   //cout << endl;
 
   return LocalError(xerr*xerr, 0, yerr*yerr);
 }

LocalPoint PixelCPETemplateReco::localPosition	(	const SiPixelCluster &	cluster,
		const GeomDetUnit &	det
	)		const

virtual

Reimplemented from PixelCPEBase.

Definition at line 134 of file PixelCPETemplateReco.cc.

References cluster_matrix_size_x, cluster_matrix_size_y, PixelCPEBase::cotalpha_, PixelCPEBase::cotbeta_, gather_cfg::cout, edm::hlt::Exception, GeomDet::geographicalId(), SiPixelTemplateDBObject::getTemplateID(), PixelCPEBase::hasFilledProb_, ierr, MagneticField::inTesla(), PixelTopology::isItBigPixelInX(), PixelTopology::isItBigPixelInY(), LoadTemplatesFromDB_, PixelCPEBase::loc_trk_pred_, Topology::localPosition(), PixelTopology::localX(), PixelTopology::localY(), LogDebug, PixelCPEBase::magfield_, micronsToCm, SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, SiPixelCluster::pixels(), SiPixelTemplateReco::PixelTempReco2D(), SiPixelTemplateSplit::PixelTempSplit(), GloballyPositioned< T >::position(), PixelCPEBase::probabilityQ_, PixelCPEBase::probabilityX_, PixelCPEBase::probabilityY_, SiPixelTemplate2D::pushfile(), PixelCPEBase::qBin_, GloballyPositioned< T >::rotation(), PixelCPEBase::setTheDet(), speed_, mathSSE::sqrt(), GeomDet::surface(), templ_, PixelCPEBase::templateDBobject_, templProbQ_, templProbX_, templProbY_, templQbin_, templSigmaX_, templSigmaY_, templXrec_, templYrec_, PixelCPEBase::theDet, PixelCPEBase::thePart, PixelCPEBase::theTopol, create_public_lumi_plots::tmp_x, create_public_lumi_plots::tmp_y, PixelCPEBase::trk_lp_x, PixelCPEBase::trk_lp_y, UseClusterSplitter_, PixelCPEBase::with_track_angle, SiPixelCluster::x(), SiPixelCluster::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by measurementPosition().

 {
   setTheDet( det, cluster );
   
   bool fpix;  //  barrel(false) or forward(true)
   if ( thePart == GeomDetEnumerators::PixelBarrel )   
     fpix = false;    // no, it's not forward -- it's barrel
   else                                              
     fpix = true;     // yes, it's forward
   
   int ID = -9999;
 
   if ( LoadTemplatesFromDB_ )
     {
       ID = templateDBobject_->getTemplateID(theDet->geographicalId());
     }
   else
     {
       if ( !fpix )
     ID = 40; // barrel
       else 
     ID = 41; // endcap
     }
   
   //cout << "PixelCPETemplateReco : ID = " << ID << endl;
   
 
   // Make from cluster (a SiPixelCluster) a boost multi_array_2d called 
   // clust_array_2d.
   boost::multi_array<float, 2> clust_array_2d(boost::extents[cluster_matrix_size_x][cluster_matrix_size_y]);
   
   // Preparing to retrieve ADC counts from the SiPixelCluster.  In the cluster,
   // we have the following:
   //   int minPixelRow(); // Minimum pixel index in the x direction (low edge).
   //   int maxPixelRow(); // Maximum pixel index in the x direction (top edge).
   //   int minPixelCol(); // Minimum pixel index in the y direction (left edge).
   //   int maxPixelCol(); // Maximum pixel index in the y direction (right edge).
   // So the pixels from minPixelRow() will go into clust_array_2d[0][*],
   // and the pixels from minPixelCol() will go into clust_array_2d[*][0].
   int row_offset = cluster.minPixelRow();
   int col_offset = cluster.minPixelCol();
   
   // Store the coordinates of the center of the (0,0) pixel of the array that 
   // gets passed to PixelTempReco2D
   // Will add these values to the output of  PixelTempReco2D
   float tmp_x = float(cluster.minPixelRow()) + 0.5;
   float tmp_y = float(cluster.minPixelCol()) + 0.5;
   
   // Store these offsets (to be added later) in a LocalPoint after tranforming 
   // them from measurement units (pixel units) to local coordinates (cm)
 
   // ggiurgiu@jhu.edu 12/09/2010 : update call with trk angles needed for bow/kink corrections
   LocalPoint lp;
   
   if ( with_track_angle )
     lp = theTopol->localPosition( MeasurementPoint(tmp_x, tmp_y), loc_trk_pred_ );
   else
     {
       edm::LogError("PixelCPETemplateReco") 
     << "@SUB = PixelCPETemplateReco::localPosition" 
     << "Should never be here. PixelCPETemplateReco should always be called with track angles. This is a bad error !!! ";
       
       lp = theTopol->localPosition( MeasurementPoint(tmp_x, tmp_y) );
     }
     
   const std::vector<SiPixelCluster::Pixel> & pixVec = cluster.pixels();
   std::vector<SiPixelCluster::Pixel>::const_iterator 
     pixIter = pixVec.begin(), pixEnd = pixVec.end();
   
   // Visualize large clusters ---------------------------------------------------------
   // From Petar: maybe this should be moved into a method in the base class?
   /*
     char cluster_matrix[100][100];
     for (int i=0; i<100; i++)
     for (int j=0; j<100; j++)
     cluster_matrix[i][j] = '.';
     
     if ( cluster.sizeX()>cluster_matrix_size_x || cluster.sizeY()>cluster_matrix_size_y )
     {       
     cout << "cluster.size()  = " << cluster.size()  << endl;
     cout << "cluster.sizeX() = " << cluster.sizeX() << endl;
     cout << "cluster.sizeY() = " << cluster.sizeY() << endl;
     
     for ( std::vector<SiPixelCluster::Pixel>::const_iterator pix = pixVec.begin(); pix != pixVec.end(); ++pix )
     {
     int i = (int)(pix->x) - row_offset;
     int j = (int)(pix->y) - col_offset;
     cluster_matrix[i][j] = '*';
     }
     
     for (int i=0; i<(int)cluster.sizeX()+2; i++)
     {
     for (int j=0; j<(int)cluster.sizeY()+2; j++)
     cout << cluster_matrix[i][j];
     cout << endl;
     }
     } // if ( cluster.sizeX()>cluster_matrix_size_x || cluster.sizeY()>cluster_matrix_size_y )
   */
   // End Visualize clusters ---------------------------------------------------------
   
 
   // Copy clust's pixels (calibrated in electrons) into clust_array_2d;
   for ( ; pixIter != pixEnd; ++pixIter ) 
     {
       // *pixIter dereferences to Pixel struct, with public vars x, y, adc (all float)
       // 02/13/2008 ggiurgiu@fnal.gov: type of x, y and adc has been changed to unsigned char, unsigned short, unsigned short
       // in DataFormats/SiPixelCluster/interface/SiPixelCluster.h so the type cast to int is redundant. Leave it there, it 
       // won't hurt. 
       int irow = int(pixIter->x) - row_offset;   // &&& do we need +0.5 ???
       int icol = int(pixIter->y) - col_offset;   // &&& do we need +0.5 ???
       
       // Gavril : what do we do here if the row/column is larger than cluster_matrix_size_x/cluster_matrix_size_y = 7/21 ?
       // Ignore them for the moment...
       if ( irow<cluster_matrix_size_x && icol<cluster_matrix_size_y )
     // 02/13/2008 ggiurgiu@fnal.gov typecast pixIter->adc to float
     clust_array_2d[irow][icol] = (float)pixIter->adc;
     }
   
   // Make and fill the bool arrays flagging double pixels
   std::vector<bool> ydouble(cluster_matrix_size_y), xdouble(cluster_matrix_size_x);
   // x directions (shorter), rows
   for (int irow = 0; irow < cluster_matrix_size_x; ++irow)
     {
       xdouble[irow] = theTopol->isItBigPixelInX( irow+row_offset );
     }
       
   // y directions (longer), columns
   for (int icol = 0; icol < cluster_matrix_size_y; ++icol) 
     {
       ydouble[icol] = theTopol->isItBigPixelInY( icol+col_offset );
     }
 
   // Output:
   float nonsense = -99999.9; // nonsense init value
   templXrec_ = templYrec_ = templSigmaX_ = templSigmaY_ = nonsense;
     // If the template recontruction fails, we want to return 1.0 for now
     templProbY_ = templProbX_ = templProbQ_ = 1.0;
     templQbin_ = 0;
     // We have a boolean denoting whether the reco failed or not
     hasFilledProb_ = false;
     
   float templYrec1_ = nonsense;
   float templXrec1_ = nonsense;
   float templYrec2_ = nonsense;
   float templXrec2_ = nonsense;
 
   // ******************************************************************
   // Do it! Use cotalpha_ and cotbeta_ calculated in PixelCPEBase
 
   GlobalVector bfield = magfield_->inTesla( theDet->surface().position() ); 
   
   Frame detFrame( theDet->surface().position(), theDet->surface().rotation() );
   LocalVector Bfield = detFrame.toLocal( bfield );
   float locBz = Bfield.z();
     
   ierr =
     PixelTempReco2D( ID, cotalpha_, cotbeta_,
              locBz, 
              clust_array_2d, ydouble, xdouble,
              templ_,
              templYrec_, templSigmaY_, templProbY_,
              templXrec_, templSigmaX_, templProbX_, 
              templQbin_, 
              speed_,
              templProbQ_
              );
 
   // ******************************************************************
 
   // Check exit status
   if ( ierr != 0 ) 
     {
       LogDebug("PixelCPETemplateReco::localPosition") <<
     "reconstruction failed with error " << ierr << "\n";
 
       // Gavril: what do we do in this case ? For now, just return the cluster center of gravity in microns
       // In the x case, apply a rough Lorentz drift average correction
       // To do: call PixelCPEGeneric whenever PixelTempReco2D fails
       double lorentz_drift = -999.9;
       if ( thePart == GeomDetEnumerators::PixelBarrel )
     lorentz_drift = 60.0; // in microns
       else if ( thePart == GeomDetEnumerators::PixelEndcap )
     lorentz_drift = 10.0; // in microns
       else 
     throw cms::Exception("PixelCPETemplateReco::localPosition :") 
       << "A non-pixel detector type in here?" << "\n";
       // ggiurgiu@jhu.edu, 21/09/2010 : trk angles needed to correct for bows/kinks
       if ( with_track_angle )
     {
       templXrec_ = theTopol->localX( cluster.x(), loc_trk_pred_ ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
       templYrec_ = theTopol->localY( cluster.y(), loc_trk_pred_ ); 
     }
       else
     {
       edm::LogError("PixelCPETemplateReco") 
         << "@SUB = PixelCPETemplateReco::localPosition" 
         << "Should never be here. PixelCPETemplateReco should always be called with track angles. This is a bad error !!! ";
       
       templXrec_ = theTopol->localX( cluster.x() ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
       templYrec_ = theTopol->localY( cluster.y() );
     }
     }
   else if ( UseClusterSplitter_ && templQbin_ == 0 )
     {
       cout << " PixelCPETemplateReco : We should never be here !!!!!!!!!!!!!!!!!!!!!!" << endl;
       cout << "                 (int)UseClusterSplitter_ = " << (int)UseClusterSplitter_ << endl;
 
       //ierr = 
       //PixelTempSplit( ID, fpix, cotalpha_, cotbeta_, 
       //        clust_array_2d, ydouble, xdouble, 
       //        templ_, 
       //        templYrec1_, templYrec2_, templSigmaY_, templProbY_,
       //        templXrec1_, templXrec2_, templSigmaX_, templProbX_, 
       //        templQbin_ );
 
 
       float dchisq;
       float templProbQ_;
       SiPixelTemplate2D templ2D_;
       templ2D_.pushfile(ID);
         
       ierr =
     SiPixelTemplateSplit::PixelTempSplit( ID, cotalpha_, cotbeta_,
                           clust_array_2d, 
                           ydouble, xdouble,
                           templ_,
                           templYrec1_, templYrec2_, templSigmaY_, templProbY_,
                           templXrec1_, templXrec2_, templSigmaX_, templProbX_,
                           templQbin_, 
                           templProbQ_,
                           true, 
                           dchisq, 
                           templ2D_ );
       
 
       if ( ierr != 0 )
     {
       LogDebug("PixelCPETemplateReco::localPosition") <<
         "reconstruction failed with error " << ierr << "\n";
       
       // Gavril: what do we do in this case ? For now, just return the cluster center of gravity in microns
       // In the x case, apply a rough Lorentz drift average correction
       // To do: call PixelCPEGeneric whenever PixelTempReco2D fails
       double lorentz_drift = -999.9;
       if ( thePart == GeomDetEnumerators::PixelBarrel )
         lorentz_drift = 60.0; // in microns
       else if ( thePart == GeomDetEnumerators::PixelEndcap )
         lorentz_drift = 10.0; // in microns
       else 
         throw cms::Exception("PixelCPETemplateReco::localPosition :") 
           << "A non-pixel detector type in here?" << "\n";
 
       // ggiurgiu@jhu.edu, 12/09/2010 : trk angles needed to correct for bows/kinks
       if ( with_track_angle )
         {
           templXrec_ = theTopol->localX( cluster.x(),loc_trk_pred_ ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
           templYrec_ = theTopol->localY( cluster.y(),loc_trk_pred_ );
         }
       else
         {
           edm::LogError("PixelCPETemplateReco") 
         << "@SUB = PixelCPETemplateReco::localPosition" 
         << "Should never be here. PixelCPETemplateReco should always be called with track angles. This is a bad error !!! ";
           
           templXrec_ = theTopol->localX( cluster.x() ) - lorentz_drift * micronsToCm; // very rough Lorentz drift correction
           templYrec_ = theTopol->localY( cluster.y() );
           
         }
     }
       else
     {
       // go from micrometer to centimeter      
       templXrec1_ *= micronsToCm;
       templYrec1_ *= micronsToCm;     
       templXrec2_ *= micronsToCm;
       templYrec2_ *= micronsToCm;
       
       // go back to the module coordinate system   
       templXrec1_ += lp.x();
       templYrec1_ += lp.y();
       templXrec2_ += lp.x();
       templYrec2_ += lp.y();
       
       // calculate distance from each hit to the track and choose the 
       // hit closest to the track
       float distance11 = sqrt( (templXrec1_ - trk_lp_x)*(templXrec1_ - trk_lp_x) + 
                    (templYrec1_ - trk_lp_y)*(templYrec1_ - trk_lp_y) );
       
       float distance12 = sqrt( (templXrec1_ - trk_lp_x)*(templXrec1_ - trk_lp_x) + 
                    (templYrec2_ - trk_lp_y)*(templYrec2_ - trk_lp_y) );
       
       float distance21 = sqrt( (templXrec2_ - trk_lp_x)*(templXrec2_ - trk_lp_x) + 
                    (templYrec1_ - trk_lp_y)*(templYrec1_ - trk_lp_y) );
       
       float distance22 = sqrt( (templXrec2_ - trk_lp_x)*(templXrec2_ - trk_lp_x) + 
                    (templYrec2_ - trk_lp_y)*(templYrec2_ - trk_lp_y) );
       
       float min_templXrec_ = -999.9;
       float min_templYrec_ = -999.9;
       float distance_min = 9999999999.9;
       if ( distance11 < distance_min )
         {
           distance_min = distance11;
           min_templXrec_ = templXrec1_;
           min_templYrec_ = templYrec1_;
         }
       if ( distance12 < distance_min )
         {
           distance_min = distance12;
           min_templXrec_ = templXrec1_;
           min_templYrec_ = templYrec2_;
         }
       if ( distance21 < distance_min )
         {
           distance_min = distance21;
           min_templXrec_ = templXrec2_;
           min_templYrec_ = templYrec1_;
         }
       if ( distance22 < distance_min )
         {
           distance_min = distance22;
           min_templXrec_ = templXrec2_;
           min_templYrec_ = templYrec2_;
         }
           
       templXrec_ = min_templXrec_;
       templYrec_ = min_templYrec_;
     }
     } // else if ( UseClusterSplitter_ && templQbin_ == 0 )
   else 
     {
       // go from micrometer to centimeter      
       templXrec_ *= micronsToCm;
       templYrec_ *= micronsToCm;
       
       // go back to the module coordinate system 
       templXrec_ += lp.x();
       templYrec_ += lp.y();
     }
     
   // Save probabilities and qBin in the quantities given to us by the base class
   // (for which there are also inline getters).  &&& templProbX_ etc. should be retired...
   probabilityX_  = templProbX_;
   probabilityY_  = templProbY_;
   probabilityQ_  = templProbQ_;
   qBin_          = templQbin_;
   
   if ( ierr == 0 ) 
     hasFilledProb_ = true;
   
   LocalPoint template_lp = LocalPoint( nonsense, nonsense );
   template_lp = LocalPoint( templXrec_, templYrec_ );      
   
   return template_lp;
   
 }

MeasurementPoint PixelCPETemplateReco::measurementPosition	(	const SiPixelCluster &	cluster,
		const GeomDetUnit &	det
	)		const

Definition at line 106 of file PixelCPETemplateReco.cc.

References LocalTrajectoryParameters::dxdz(), LocalTrajectoryParameters::dydz(), PixelCPEBase::loc_traj_param_, localPosition(), Topology::measurementPosition(), PixelCPEBase::theTopol, and PixelCPEBase::with_track_angle.

 {
   LocalPoint lp = localPosition(cluster,det);
 
   // ggiurgiu@jhu.edu 12/09/2010 : trk angles needed to correct for bows/kinks
   if ( with_track_angle )
     return theTopol->measurementPosition(lp, Topology::LocalTrackAngles( loc_traj_param_.dxdz(), loc_traj_param_.dydz() ) );
   else
     {
       edm::LogError("PixelCPETemplateReco") 
     << "@SUB = PixelCPETemplateReco::measurementPosition" 
     << "Should never be here. PixelCPETemplateReco should always be called with track angles. This is a bad error !!! ";
       
       return theTopol->measurementPosition( lp ); 
     }
 
 }