#include <PixelCPEBase.h>

Inheritance diagram for PixelCPEBase:

Classes
struct	ClusterParam

struct	DetParam

Public Member Functions
ReturnType	getParameters (const SiPixelCluster &cl, const GeomDetUnit &det) const override

ReturnType	getParameters (const SiPixelCluster &cl, const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const override

	PixelCPEBase (edm::ParameterSet const &conf, const MagneticField mag, const TrackerGeometry &geom, const TrackerTopology &ttopo, const SiPixelLorentzAngle lorentzAngle, const SiPixelGenErrorDBObject genErrorDBObject, const SiPixelTemplateDBObject templateDBobject, const SiPixelLorentzAngle *lorentzAngleWidth, int flag=0)

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

virtual ReturnType	getParameters (const SiPixelCluster &cl, const GeomDetUnit &det, const TrajectoryStateOnSurface &tsos) const

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

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

	PixelClusterParameterEstimator ()

virtual	~PixelClusterParameterEstimator ()

Protected Types
using	DetParams = std::vector< DetParam >

typedef GloballyPositioned< double >	Frame

Protected Member Functions
void	computeAnglesFromDetPosition (DetParam const &theDetParam, ClusterParam &theClusterParam) const

void	computeAnglesFromTrajectory (DetParam const &theDetParam, ClusterParam &theClusterParam, const LocalTrajectoryParameters &ltp) const

void	computeLorentzShifts (DetParam &) const

DetParam const &	detParam (const GeomDetUnit &det) const

LocalVector	driftDirection (DetParam &theDetParam, GlobalVector bfield) const

LocalVector	driftDirection (DetParam &theDetParam, LocalVector bfield) const

bool	isFlipped (DetParam const &theDetParam) const

void	setTheClu (DetParam const &, ClusterParam &theClusterParam) const

Protected Attributes
bool	alpha2Order

bool	DoLorentz_

const SiPixelGenErrorDBObject *	genErrorDBObject_

const TrackerGeometry &	geom_

float	lAOffset_

float	lAWidthBPix_

float	lAWidthFPix_

bool	LoadTemplatesFromDB_

const SiPixelLorentzAngle *	lorentzAngle_

const SiPixelLorentzAngle *	lorentzAngleWidth_

DetParams	m_DetParams =DetParams(1440)

const MagneticField *	magfield_

const SiPixelTemplateDBObject *	templateDBobject_

int	theFlag_

int	theVerboseLevel

const TrackerTopology &	ttopo_

bool	useLAOffsetFromConfig_

bool	useLAWidthFromConfig_

bool	useLAWidthFromDB_

Protected Attributes inherited from PixelClusterParameterEstimator
unsigned int	clusterProbComputationFlag_

Private Member Functions
virtual ClusterParam *	createClusterParam (const SiPixelCluster &cl) const =0

void	fillDetParams ()

virtual LocalError	localError (DetParam const &theDetParam, ClusterParam &theClusterParam) const =0

virtual LocalPoint	localPosition (DetParam const &theDetParam, ClusterParam &theClusterParam) const =0

SiPixelRecHitQuality::QualWordType	rawQualityWord (ClusterParam &theClusterParam) const

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

using	ReturnType = std::tuple< LocalPoint, LocalError, SiPixelRecHitQuality::QualWordType >

typedef std::vector< LocalValues >	VLocalValues

Detailed Description

Definition at line 54 of file PixelCPEBase.h.

Member Typedef Documentation

using PixelCPEBase::DetParams = std::vector<DetParam>

protected

Definition at line 282 of file PixelCPEBase.h.

typedef GloballyPositioned<double> PixelCPEBase::Frame

protected

Definition at line 218 of file PixelCPEBase.h.

Constructor & Destructor Documentation

PixelCPEBase::PixelCPEBase	(	edm::ParameterSet const &	conf,
		const MagneticField *	mag,
		const TrackerGeometry &	geom,
		const TrackerTopology &	ttopo,
		const SiPixelLorentzAngle *	lorentzAngle,
		const SiPixelGenErrorDBObject *	genErrorDBObject,
		const SiPixelTemplateDBObject *	templateDBobject,
		const SiPixelLorentzAngle *	lorentzAngleWidth,
		int	flag = `0`
	)

Definition at line 32 of file PixelCPEBase.cc.

References alpha2Order, PixelClusterParameterEstimator::clusterProbComputationFlag_, DoLorentz_, edm::ParameterSet::existsAs(), fillDetParams(), genErrorDBObject_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), createfilelist::int, lAOffset_, lAWidthBPix_, lAWidthFPix_, LoadTemplatesFromDB_, LogDebug, lorentzAngle_, lorentzAngleWidth_, templateDBobject_, theFlag_, theVerboseLevel, useLAOffsetFromConfig_, useLAWidthFromConfig_, and useLAWidthFromDB_.

 : useLAOffsetFromConfig_(false),
 useLAWidthFromConfig_(false), useLAWidthFromDB_(false), theFlag_(flag),
 magfield_(mag), geom_(geom), ttopo_(ttopo)
 {
    
 #ifdef EDM_ML_DEBUG
    nRecHitsTotal_=0;
    nRecHitsUsedEdge_=0,
 #endif
    
    //--- Lorentz angle tangent per Tesla
    lorentzAngle_ = lorentzAngle;
    lorentzAngleWidth_ = lorentzAngleWidth;
    
    //-- GenError Calibration Object (different from SiPixelCPEGenericErrorParm) from DB
    genErrorDBObject_ = genErrorDBObject;
    
    //-- Template Calibration Object from DB
    if(theFlag_!=0) templateDBobject_ = templateDBobject; // flag to check if it is generic or templates
    
    // Configurables
    // For both templates & generic
    
    // Read templates and/or generic errors from DB
    LoadTemplatesFromDB_ = conf.getParameter<bool>("LoadTemplatesFromDB");
    
    //--- Algorithm's verbosity
    theVerboseLevel =
    conf.getUntrackedParameter<int>("VerboseLevel",0);
    
    //-- Switch on/off E.B
    alpha2Order = conf.getParameter<bool>("Alpha2Order");
    
    //--- A flag that could be used to change the behavior of
    //--- clusterProbability() in TSiPixelRecHit (the *transient* one).
    //--- The problem is that the transient hits are made after the CPE runs
    //--- and they don't get the access to the PSet, so we pass it via the
    //--- CPE itself...
    //
    clusterProbComputationFlag_
    = (unsigned int) conf.getParameter<int>("ClusterProbComputationFlag");
    
    // This LA related parameters are only relevant for the Generic algo
    // They still have to be used in Base since the LA computation is in Base
    
    // Use LA-width from DB.
    // If both (this and from config) are false LA-width is calcuated from LA-offset
    useLAWidthFromDB_ = conf.existsAs<bool>("useLAWidthFromDB")?
    conf.getParameter<bool>("useLAWidthFromDB"):false;
    
    // Use Alignment LA-offset in generic
    // (Experimental; leave commented out)
    //useLAAlignmentOffsets_ = conf.existsAs<bool>("useLAAlignmentOffsets")?
    //conf.getParameter<bool>("useLAAlignmentOffsets"):false;
    
    // Used only for testing
    lAOffset_ = conf.existsAs<double>("lAOffset")?  // fixed LA value
    conf.getParameter<double>("lAOffset"):0.0;
    lAWidthBPix_ = conf.existsAs<double>("lAWidthBPix")?   // fixed LA width
    conf.getParameter<double>("lAWidthBPix"):0.0;
    lAWidthFPix_ = conf.existsAs<double>("lAWidthFPix")?   // fixed LA width
    conf.getParameter<double>("lAWidthFPix"):0.0;
    
    // Use LA-offset from config, for testing only
    if (lAOffset_>0.0) 
      useLAOffsetFromConfig_ = true;
    // Use LA-width from config, split into fpix & bpix, for testing only
    if (lAWidthBPix_>0.0 || lAWidthFPix_>0.0) 
      useLAWidthFromConfig_ = true;
    
    
    // For Templates only
    // Compute the Lorentz shifts for this detector element for templates (from Alignment)
    DoLorentz_ = conf.existsAs<bool>("DoLorentz")?conf.getParameter<bool>("DoLorentz"):false;
    
    LogDebug("PixelCPEBase") <<" LA constants - "
                 << lAOffset_ << " " << lAWidthBPix_ << " " <<lAWidthFPix_ << endl; //dk
    
    fillDetParams();
 }

Member Function Documentation

void PixelCPEBase::computeAnglesFromDetPosition	(	DetParam const &	theDetParam,
		ClusterParam &	theClusterParam
	)		const

protected

Definition at line 334 of file PixelCPEBase.cc.

References PixelCPEBase::ClusterParam::cotalpha, PixelCPEBase::ClusterParam::cotbeta, Topology::localPosition(), PixelCPEBase::ClusterParam::theCluster, PixelCPEBase::DetParam::theOrigin, PixelCPEBase::DetParam::theTopol, PixelCPEBase::ClusterParam::with_track_angle, PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::y(), SiPixelCluster::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by computeAnglesFromTrajectory(), and getParameters().

 {
   // &&& PM: darn, we have a lot of junk here.  Need to clean up.
    
    /*
     // get cluster center of gravity (of charge)
     float xcenter = cl.x();
     float ycenter = cl.y();
     
     // get the cluster position in local coordinates (cm)
     
     // ggiurgiu@jhu.edu 12/09/2010 : This function is called without track info, therefore there are no track
     // angles to provide here. Call the default localPosition (without track info)
     LocalPoint lp = theTopol->localPosition( MeasurementPoint(xcenter, ycenter) );
     
     
     // get the cluster position in global coordinates (cm)
     GlobalPoint gp = theDet->surface().toGlobal( lp );
     float gp_mod = sqrt( gp.x()*gp.x() + gp.y()*gp.y() + gp.z()*gp.z() );
     
     // normalize
     float gpx = gp.x()/gp_mod;
     float gpy = gp.y()/gp_mod;
     float gpz = gp.z()/gp_mod;
     
     // make a global vector out of the global point; this vector will point from the
     // origin of the detector to the cluster
     GlobalVector gv(gpx, gpy, gpz);
     
     // make local unit vector along local X axis
     const Local3DVector lvx(1.0, 0.0, 0.0);
     
     // get the unit X vector in global coordinates/
     GlobalVector gvx = theDet->surface().toGlobal( lvx );
     
     // make local unit vector along local Y axis
     const Local3DVector lvy(0.0, 1.0, 0.0);
     
     // get the unit Y vector in global coordinates
     GlobalVector gvy = theDet->surface().toGlobal( lvy );
     
     // make local unit vector along local Z axis
     const Local3DVector lvz(0.0, 0.0, 1.0);
     
     // get the unit Z vector in global coordinates
     GlobalVector gvz = theDet->surface().toGlobal( lvz );
     
     // calculate the components of gv (the unit vector pointing to the cluster)
     // in the local coordinate system given by the basis {gvx, gvy, gvz}
     // note that both gv and the basis {gvx, gvy, gvz} are given in global coordinates
     float gv_dot_gvx = gv.x()*gvx.x() + gv.y()*gvx.y() + gv.z()*gvx.z();
     float gv_dot_gvy = gv.x()*gvy.x() + gv.y()*gvy.y() + gv.z()*gvy.z();
     float gv_dot_gvz = gv.x()*gvz.x() + gv.y()*gvz.y() + gv.z()*gvz.z();
     */
    
    // all the above is equivalent to
    LocalPoint lp = theDetParam.theTopol->localPosition( MeasurementPoint(theClusterParam.theCluster->x(), theClusterParam.theCluster->y()) );
    auto gvx = lp.x()-theDetParam.theOrigin.x();
    auto gvy = lp.y()-theDetParam.theOrigin.y();
    auto gvz = -1.f/theDetParam.theOrigin.z();
    //--- Note that the normalization is not required as only the ratio used  
    
    //theClusterParam.zneg = (gvz < 0); // Not used, AH
    
    // calculate angles
    theClusterParam.cotalpha = gvx*gvz;
    theClusterParam.cotbeta  = gvy*gvz;
    
    theClusterParam.with_track_angle = false;
    
    
    /*
     // used only in dberror param...
     auto alpha = HALF_PI - std::atan(cotalpha_);
     auto beta = HALF_PI - std::atan(cotbeta_);
     if (zneg) { beta -=PI; alpha -=PI;}
     
     auto alpha_ = atan2( gv_dot_gvz, gv_dot_gvx );
     auto beta_  = atan2( gv_dot_gvz, gv_dot_gvy );
     
     assert(std::abs(std::round(alpha*10000.f)-std::round(alpha_*10000.f))<2);
     assert(std::abs(std::round(beta*10000.f)-std::round(beta_*10000.f))<2);
     */
    
 }

void PixelCPEBase::computeAnglesFromTrajectory	(	DetParam const &	theDetParam,
		ClusterParam &	theClusterParam,
		const LocalTrajectoryParameters &	ltp
	)		const

protected

Definition at line 277 of file PixelCPEBase.cc.

References computeAnglesFromDetPosition(), PixelCPEBase::ClusterParam::cotalpha, PixelCPEBase::ClusterParam::cotbeta, LocalTrajectoryParameters::dxdz(), LocalTrajectoryParameters::dydz(), PixelCPEBase::ClusterParam::loc_trk_pred, LocalTrajectoryParameters::position(), PixelCPEBase::ClusterParam::trk_lp_x, PixelCPEBase::ClusterParam::trk_lp_y, PixelCPEBase::ClusterParam::with_track_angle, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by getParameters(), and setTheClu().

 {
   // &&& PM: this comment is a candidate for deletion, but ask around first.
    //cout<<" in PixelCPEBase:computeAnglesFromTrajectory - "<<endl; //dk
    /*
    //theClusterParam.loc_traj_param = ltp;   
    LocalVector localDir = ltp.momentum();
    float locx = localDir.x();
    float locy = localDir.y();
    float locz = localDir.z();
     // Danek's definition
     alpha_ = acos(locx/sqrt(locx*locx+locz*locz));
     if ( isFlipped() )                    // &&& check for FPIX !!!
     alpha_ = PI - alpha_ ;
     beta_ = acos(locy/sqrt(locy*locy+locz*locz));
     */
    
    
    theClusterParam.cotalpha = ltp.dxdz();
    theClusterParam.cotbeta  = ltp.dydz();
    
    
    LocalPoint trk_lp = ltp.position();
    theClusterParam.trk_lp_x = trk_lp.x();
    theClusterParam.trk_lp_y = trk_lp.y();
    
    theClusterParam.with_track_angle = true;
    
    // GG: needed to correct for bows/kinks
    theClusterParam.loc_trk_pred = 
       Topology::LocalTrackPred(theClusterParam.trk_lp_x, theClusterParam.trk_lp_y, 
                                theClusterParam.cotalpha,theClusterParam.cotbeta);
    
 }

void PixelCPEBase::computeLorentzShifts ( DetParam & theDetParam ) const

protected

Definition at line 550 of file PixelCPEBase.cc.

References PixelCPEBase::DetParam::driftDirection, LogDebug, PixelCPEBase::DetParam::lorentzShiftInCmX, PixelCPEBase::DetParam::lorentzShiftInCmY, PixelCPEBase::DetParam::theThickness, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by fillDetParams().

 {
    // Max shift (at the other side of the sensor) in cm
    theDetParam.lorentzShiftInCmX = theDetParam.driftDirection.x()/theDetParam.driftDirection.z() * theDetParam.theThickness;  //
    theDetParam.lorentzShiftInCmY = theDetParam.driftDirection.y()/theDetParam.driftDirection.z() * theDetParam.theThickness;  //
    
    LogDebug("PixelCPEBase::computeLorentzShifts()") << " lorentzShiftsInCmX,Y = "
                             << theDetParam.lorentzShiftInCmX <<" "
                             << theDetParam.lorentzShiftInCmY <<" "
                             << theDetParam.driftDirection;
 }

virtual ClusterParam* PixelCPEBase::createClusterParam ( const SiPixelCluster & cl ) const

privatepure virtual

Implemented in PixelCPEGeneric, PixelCPETemplateReco, and PixelCPEClusterRepair.

Referenced by getParameters().

PixelCPEBase::DetParam const & PixelCPEBase::detParam ( const GeomDetUnit & det ) const

protected

Definition at line 443 of file PixelCPEBase.cc.

References mps_fire::i, GeomDet::index(), m_DetParams, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by getParameters().

 {
    auto i = det.index();
    //cout << "get parameters of detector " << i << endl;
    assert(i<int(m_DetParams.size()));
    //if (i>=int(m_DetParams.size())) m_DetParams.resize(i+1);  // should never happen!
    const DetParam & p = m_DetParams[i];
    return p;
 }

LocalVector PixelCPEBase::driftDirection	(	DetParam &	theDetParam,
		GlobalVector	bfield
	)		const

protected

Definition at line 461 of file PixelCPEBase.cc.

References GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), GeomDet::surface(), and PixelCPEBase::DetParam::theDet.

Referenced by fillDetParams().

 {
    Frame detFrame(theDetParam.theDet->surface().position(), theDetParam.theDet->surface().rotation());
    LocalVector Bfield = detFrame.toLocal(bfield);
    return driftDirection(theDetParam,Bfield);
 }

LocalVector PixelCPEBase::driftDirection	(	DetParam &	theDetParam,
		LocalVector	bfield
	)		const

protected

Definition at line 470 of file PixelCPEBase.cc.

References funct::abs(), alpha2Order, createTree::dd, f, GeomDet::geographicalId(), SiPixelLorentzAngle::getLorentzAngle(), GeomDetEnumerators::isBarrel(), GeomDetEnumerators::isTrackerPixel(), lAOffset_, lAWidthBPix_, lAWidthFPix_, LogDebug, lorentzAngle_, lorentzAngleWidth_, DetId::rawId(), Scenarios_cff::scale, PixelCPEBase::DetParam::theDet, theFlag_, PixelCPEBase::DetParam::thePart, useLAOffsetFromConfig_, useLAWidthFromConfig_, useLAWidthFromDB_, PixelCPEBase::DetParam::widthLAFractionX, PixelCPEBase::DetParam::widthLAFractionY, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

 {
    // Use LA from DB or from config
    float langle = 0.;
    if( !useLAOffsetFromConfig_ ) {  // get it from DB
       if(lorentzAngle_ != nullptr) {  // a real LA object
          langle = lorentzAngle_->getLorentzAngle(theDetParam.theDet->geographicalId().rawId());
          LogDebug("PixelCPEBase::driftDirection()") 
        <<" la "<<langle<<" "<< theDetParam.theDet->geographicalId().rawId() <<endl;
       } else { // no LA, unused
      langle = 0; // set to a fake value
      LogDebug("PixelCPEBase::driftDirection()") <<" LA object is NULL, assume LA = 0"<<endl; //dk
       }
       LogDebug("PixelCPEBase::driftDirection()") << " Will use LA Offset from DB "<<langle<<endl;
    } else {  // from config file
       langle = lAOffset_;
       LogDebug("PixelCPEBase::driftDirection()") << " Will use LA Offset from config "<<langle<<endl;
    }
    
    // Now do the LA width stuff
    theDetParam.widthLAFractionX = 1.; // predefine to 1 (default) if things fail
    theDetParam.widthLAFractionY = 1.;
    
    // Compute the charge width, generic only
    if(theFlag_==0) {
       
       if(useLAWidthFromDB_ && (lorentzAngleWidth_ != nullptr) ) {
          // take it from a separate, special LA DB object (forWidth)
          
          auto langleWidth = lorentzAngleWidth_->getLorentzAngle(theDetParam.theDet->geographicalId().rawId());
          if(langleWidth!=0.0) theDetParam.widthLAFractionX = std::abs(langleWidth/langle);
          // leave the widthLAFractionY=1.
          //cout<<" LAWidth lorentz width "<<theDetParam.widthLAFractionX<<" "<<theDetParam.widthLAFractionY<<endl;
          
       } else if(useLAWidthFromConfig_) { // get from config
          
          double lAWidth=0;
          if( GeomDetEnumerators::isTrackerPixel(theDetParam.thePart) && GeomDetEnumerators::isBarrel(theDetParam.thePart) ) lAWidth = lAWidthBPix_; // barrel
          else lAWidth = lAWidthFPix_;
          
          if(langle!=0.0) theDetParam.widthLAFractionX = std::abs(lAWidth/langle);
          // fix the FractionY at 1
          
          //cout<<" Lorentz width from config "<<theDetParam.widthLAFractionX<<" "<<theDetParam.widthLAFractionY<<endl;
          
       } else { // get if from the offset LA (old method used until 2013)
          // do nothing
          //cout<<" Old default LA width method "<<theDetParam.widthLAFractionX<<" "<<theDetParam.widthLAFractionY<<endl;
          
       }
       
       //cout<<" Final LA fraction  "<<theDetParam.widthLAFractionX<<" "<<theDetParam.widthLAFractionY<<endl;
       
    }  // if flag
    
 
    float alpha2 = alpha2Order ?  langle*langle : 0; //
    
    // **********************************************************************
    // Our convention is the following:
    // +x is defined by the direction of the Lorentz drift!
    // +z is defined by the direction of E field (so electrons always go into -z!)
    // +y is defined by +x and +z, and it turns out to be always opposite to the +B field.
    // **********************************************************************
    
    float dir_x = -( langle * Bfield.y() + alpha2* Bfield.z()* Bfield.x() );
    float dir_y =  ( langle * Bfield.x() - alpha2* Bfield.z()* Bfield.y() );
    float dir_z = -( 1.f + alpha2* Bfield.z()*Bfield.z() );
    auto scale = 1.f/std::abs( dir_z );  // same as 1 + alpha2*Bfield.z()*Bfield.z()
    LocalVector  dd(dir_x*scale, dir_y*scale, -1.f );  // last is -1 !
    
    LogDebug("PixelCPEBase") << " The drift direction in local coordinate is "  << dd   ;
    
    return dd;
 }

void PixelCPEBase::fillDetParams ( )

private

Definition at line 127 of file PixelCPEBase.cc.

References computeLorentzShifts(), TrackerGeometry::detUnits(), DoLorentz_, driftDirection(), genErrorDBObject_, geom_, SiPixelGenErrorDBObject::getGenErrorID(), SiPixelTemplateDBObject::getTemplateID(), mps_fire::i, createfilelist::int, MagneticField::inTesla(), LoadTemplatesFromDB_, LogDebug, m_DetParams, magfield_, TrackerGeometry::offsetDU(), AlCaHLTBitMon_ParallelJobs::p, templateDBobject_, theFlag_, GeomDetEnumerators::tkDetEnum, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by getParameters(), and PixelCPEBase().

 {
   // &&& PM: I have no idea what this code is doing, and what it is doing here!???
   //
    auto const & dus = geom_.detUnits();
    unsigned m_detectors = dus.size();
    for(unsigned int i=1;i<7;++i) {
       LogDebug("PixelCPEBase:: LookingForFirstStrip") << "Subdetector " << i
       << " GeomDetEnumerator " << GeomDetEnumerators::tkDetEnum[i]
       << " offset " << geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])
       << " is it strip? " << (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() ?
                               dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip() : false);
       if(geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() &&
          dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isTrackerStrip()) {
          if(geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) < m_detectors) m_detectors = geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]);
       }
    }
    LogDebug("LookingForFirstStrip") << " Chosen offset: " << m_detectors;
    
    
    m_DetParams.resize(m_detectors);
    LogDebug("PixelCPEBase::fillDetParams():") <<"caching "<<m_detectors<<" pixel detectors"<<endl;
    for (unsigned i=0; i!=m_detectors;++i) {
       auto & p=m_DetParams[i];
       p.theDet = dynamic_cast<const PixelGeomDetUnit*>(dus[i]);
       assert(p.theDet);
       assert(p.theDet->index()==int(i));
       
       p.theOrigin = p.theDet->surface().toLocal(GlobalPoint(0,0,0));
       
       //--- p.theDet->type() returns a GeomDetType, which implements subDetector()
       p.thePart = p.theDet->type().subDetector();
       
       
       //--- The location in of this DetUnit in a cyllindrical coord system (R,Z)
       //--- The call goes via BoundSurface, returned by p.theDet->surface(), but
       //--- position() is implemented in GloballyPositioned<> template
       //--- ( BoundSurface : Surface : GloballyPositioned<float> )
       //p.theDetR = p.theDet->surface().position().perp();  //Not used, AH
       //p.theDetZ = p.theDet->surface().position().z();  //Not used, AH
       //--- Define parameters for chargewidth calculation
       
       //--- bounds() is implemented in BoundSurface itself.
       p.theThickness = p.theDet->surface().bounds().thickness();
       
       // Cache the det id for templates and generic erros
       
       if(theFlag_==0) { // for generic
          if(LoadTemplatesFromDB_ ) // do only if genError requested
             p.detTemplateId = genErrorDBObject_->getGenErrorID(p.theDet->geographicalId().rawId());
       } else {          // for templates
          p.detTemplateId = templateDBobject_->getTemplateID(p.theDet->geographicalId());
       }
       
       // &&& PM: I'm not sure what this does. Ask around.
       // just for testing
       //int i1 = 0;
       //if(theFlag_==0) i1 = genErrorDBObject_->getGenErrorID(p.theDet->geographicalId().rawId());
       //int i2= templateDBobject_->getTemplateID(p.theDet->geographicalId().rawId());
       //int i3= templateDBobject_->getTemplateID(p.theDet->geographicalId());
       //if(i2!=i3) cout<<i2<<" != "<<i3<<endl;
       //cout<<i<<" "<<p.detTemplateId<<" "<<i1<<" "<<i2<<" "<<i3<<endl;
       
       auto topol = &(p.theDet->specificTopology());
       p.theTopol=topol;
       auto const proxyT = dynamic_cast<const ProxyPixelTopology*>(p.theTopol);
       if (proxyT) p.theRecTopol = dynamic_cast<const RectangularPixelTopology*>(&(proxyT->specificTopology()));
       else p.theRecTopol = dynamic_cast<const RectangularPixelTopology*>(p.theTopol);
       assert(p.theRecTopol);
       
       //--- The geometrical description of one module/plaquette
       //p.theNumOfRow = p.theRecTopol->nrows(); // rows in x //Not used, AH. PM: leave commented out.
       //p.theNumOfCol = p.theRecTopol->ncolumns();  // cols in y //Not used, AH. PM: leave commented out.
       std::pair<float,float> pitchxy = p.theRecTopol->pitch();
       p.thePitchX = pitchxy.first;       // pitch along x
       p.thePitchY = pitchxy.second;      // pitch along y
       
       //p.theSign = isFlipped(&p) ? -1 : 1; //Not used, AH.  PM: leave commented out.
       
       LocalVector Bfield = p.theDet->surface().toLocal(magfield_->inTesla(p.theDet->surface().position()));
       p.bz = Bfield.z();
       p.bx = Bfield.x();
       
       
       //---  Compute the Lorentz shifts for this detector element
       if ( (theFlag_==0) || DoLorentz_ ) {  // do always for generic and if(DOLorentz) for templates
          p.driftDirection = driftDirection(p, Bfield );
          computeLorentzShifts(p);
       }
       
       
       LogDebug("PixelCPEBase::fillDetParams()") << "***** PIXEL LAYOUT *****"
       << " thePart = " << p.thePart
       << " theThickness = " << p.theThickness
       << " thePitchX  = " << p.thePitchX
       << " thePitchY  = " << p.thePitchY;
       //                 << " theLShiftX  = " << p.theLShiftX;
       
       
    }
 }

ReturnType PixelCPEBase::getParameters	(	const SiPixelCluster &	cl,
		const GeomDetUnit &	det
	)		const

inlineoverridevirtual

Implements PixelClusterParameterEstimator.

Definition at line 141 of file PixelCPEBase.h.

References computeAnglesFromDetPosition(), createClusterParam(), detParam(), localError(), localPosition(), rawQualityWord(), and setTheClu().

Referenced by cms::SiPixelRecHitConverter::run().

    {
 #ifdef EDM_ML_DEBUG
       nRecHitsTotal_++ ;
       //std::cout<<" in PixelCPEBase:localParameters(all) - "<<nRecHitsTotal_<<std::endl;  //dk
 #endif
       
       DetParam const & theDetParam = detParam(det);
       ClusterParam * theClusterParam = createClusterParam(cl);
       setTheClu( theDetParam, *theClusterParam );
       computeAnglesFromDetPosition(theDetParam, *theClusterParam);
       
       // localPosition( cl, det ) must be called before localError( cl, det ) !!!
       LocalPoint lp = localPosition(theDetParam, *theClusterParam);
       LocalError le = localError(theDetParam, *theClusterParam);
       SiPixelRecHitQuality::QualWordType rqw = rawQualityWord(*theClusterParam);
       auto tuple = std::make_tuple(lp, le , rqw);
       delete theClusterParam;
       
       //std::cout<<" in PixelCPEBase:localParameters(all) - "<<lp.x()<<" "<<lp.y()<<std::endl;  //dk
       return tuple;
    }

ReturnType PixelCPEBase::getParameters	(	const SiPixelCluster &	cl,
		const GeomDetUnit &	det,
		const LocalTrajectoryParameters &	ltp
	)		const

inlineoverridevirtual

Implements PixelClusterParameterEstimator.

Definition at line 168 of file PixelCPEBase.h.

References GetRecoTauVFromDQM_MC_cff::cl, computeAnglesFromTrajectory(), createClusterParam(), detParam(), fillDetParams(), localError(), localPosition(), rawQualityWord(), and setTheClu().

    {
 #ifdef EDM_ML_DEBUG
       nRecHitsTotal_++ ;
       //std::cout<<" in PixelCPEBase:localParameters(on track) - "<<nRecHitsTotal_<<std::endl;  //dk
 #endif
       
       DetParam const & theDetParam = detParam(det);
       ClusterParam *  theClusterParam = createClusterParam(cl);
       setTheClu( theDetParam, *theClusterParam );
       computeAnglesFromTrajectory(theDetParam, *theClusterParam, ltp);
       
       // localPosition( cl, det ) must be called before localError( cl, det ) !!!
       LocalPoint lp = localPosition(theDetParam, *theClusterParam);
       LocalError le = localError(theDetParam, *theClusterParam);
       SiPixelRecHitQuality::QualWordType rqw = rawQualityWord(*theClusterParam);
       auto tuple = std::make_tuple(lp, le , rqw);
       delete theClusterParam;
       
       //std::cout<<" in PixelCPEBase:localParameters(on track) - "<<lp.x()<<" "<<lp.y()<<std::endl;  //dk
       return tuple;
    }

bool PixelCPEBase::isFlipped ( DetParam const & theDetParam ) const

protected

Definition at line 434 of file PixelCPEBase.cc.

References PV3DBase< T, PVType, FrameType >::perp2(), GeomDet::surface(), PixelCPEBase::DetParam::theDet, and Surface::toGlobal().

 {
    // Check the relative position of the local +/- z in global coordinates.
    float tmp1 = theDetParam.theDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp2();
    float tmp2 = theDetParam.theDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp2();
    if ( tmp2 < tmp1 ) return true;
    else return false;
 }

virtual LocalError PixelCPEBase::localError	(	DetParam const &	theDetParam,
		ClusterParam &	theClusterParam
	)		const

privatepure virtual

Implemented in PixelCPEGeneric, PixelCPETemplateReco, and PixelCPEClusterRepair.

Referenced by getParameters().

virtual LocalPoint PixelCPEBase::localPosition	(	DetParam const &	theDetParam,
		ClusterParam &	theClusterParam
	)		const

privatepure virtual

Implemented in PixelCPEGeneric, PixelCPETemplateReco, and PixelCPEClusterRepair.

Referenced by getParameters().

SiPixelRecHitQuality::QualWordType PixelCPEBase::rawQualityWord ( ClusterParam & theClusterParam ) const

private

A convenience method to fill a whole SiPixelRecHitQuality word in one shot. This way, we can keep the details of what is filled within the pixel code and not expose the Transient SiPixelRecHit to it as well. The name of this function is chosen to match the one in SiPixelRecHit.

Definition at line 569 of file PixelCPEBase.cc.

Referenced by getParameters().

 {
    SiPixelRecHitQuality::QualWordType qualWord(0);
    if (theClusterParam.hasFilledProb_) {
      float probabilityXY=0;
      if(theClusterParam.filled_from_2d) probabilityXY = theClusterParam.probabilityX_;
      else if ( theClusterParam.probabilityX_ !=0 && theClusterParam.probabilityY_ !=0 )
        probabilityXY = theClusterParam.probabilityX_ * theClusterParam.probabilityY_ * (1.f - std::log(theClusterParam.probabilityX_ * theClusterParam.probabilityY_) ) ;
      SiPixelRecHitQuality::thePacking.setProbabilityXY ( probabilityXY ,
                                                         qualWord );
    
      SiPixelRecHitQuality::thePacking.setProbabilityQ  ( theClusterParam.probabilityQ_ ,
                                                         qualWord );
    }
    SiPixelRecHitQuality::thePacking.setQBin          ( theClusterParam.qBin_, 
                                                       qualWord );
    
    SiPixelRecHitQuality::thePacking.setIsOnEdge      ( theClusterParam.isOnEdge_,
                                                       qualWord );
    
    SiPixelRecHitQuality::thePacking.setHasBadPixels  ( theClusterParam.hasBadPixels_,
                                                       qualWord );
    
    SiPixelRecHitQuality::thePacking.setSpansTwoROCs  ( theClusterParam.spansTwoROCs_,
                                                       qualWord );
    
    SiPixelRecHitQuality::thePacking.setHasFilledProb ( theClusterParam.hasFilledProb_,
                                                       qualWord );
    
    return qualWord;
 }

void PixelCPEBase::setTheClu	(	DetParam const &	theDetParam,
		ClusterParam &	theClusterParam
	)		const

protected

Definition at line 233 of file PixelCPEBase.cc.

References computeAnglesFromTrajectory(), RectangularPixelTopology::containsBigPixelInX(), RectangularPixelTopology::containsBigPixelInY(), PixelCPEBase::ClusterParam::edgeTypeX_, PixelCPEBase::ClusterParam::edgeTypeY_, RectangularPixelTopology::isItEdgePixelInX(), RectangularPixelTopology::isItEdgePixelInY(), PixelCPEBase::ClusterParam::isOnEdge_, SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), PixelCPEBase::ClusterParam::spansTwoROCs_, PixelCPEBase::ClusterParam::theCluster, and PixelCPEBase::DetParam::theRecTopol.

Referenced by getParameters().

 {
    
    //--- Geometric Quality Information
    int minInX,minInY,maxInX,maxInY=0;
    minInX = theClusterParam.theCluster->minPixelRow();
    minInY = theClusterParam.theCluster->minPixelCol();
    maxInX = theClusterParam.theCluster->maxPixelRow();
    maxInY = theClusterParam.theCluster->maxPixelCol();
 
    if      ( theDetParam.theRecTopol->isItEdgePixelInX(minInX) )
      theClusterParam.edgeTypeX_ = 1;
    else if ( theDetParam.theRecTopol->isItEdgePixelInX(maxInX) )
      theClusterParam.edgeTypeX_ = 2;
    else
      theClusterParam.edgeTypeX_ = 0;
      
    if      ( theDetParam.theRecTopol->isItEdgePixelInY(minInY) )
      theClusterParam.edgeTypeY_ = 1;
    else if ( theDetParam.theRecTopol->isItEdgePixelInY(maxInY) )
      theClusterParam.edgeTypeY_ = 2;
    else
      theClusterParam.edgeTypeY_ = 0;
    
    theClusterParam.isOnEdge_ = ( theClusterParam.edgeTypeX_ || theClusterParam.edgeTypeY_ );
    
    // &&& FOR NOW UNUSED. KEEP IT IN CASE WE WANT TO USE IT IN THE FUTURE
    // Bad Pixels have their charge set to 0 in the clusterizer
    //hasBadPixels_ = false;
    //for(unsigned int i=0; i<theClusterParam.theCluster->pixelADC().size(); ++i) {
    //if(theClusterParam.theCluster->pixelADC()[i] == 0) { hasBadPixels_ = true; break;}
    //}
    
    theClusterParam.spansTwoROCs_ = theDetParam.theRecTopol->containsBigPixelInX(minInX,maxInX) |
    theDetParam.theRecTopol->containsBigPixelInY(minInY,maxInY);
    
 }