#include <PixelCPEClusterRepair.h>

Inheritance diagram for PixelCPEClusterRepair:

Classes
struct	ClusterParamTemplate

class	Rule

Public Member Functions
	PixelCPEClusterRepair (edm::ParameterSet const &conf, const MagneticField , const TrackerGeometry &, const TrackerTopology &, const SiPixelLorentzAngle , const SiPixelTemplateDBObject , const SiPixel2DTemplateDBObject )

	~PixelCPEClusterRepair () override

Public Member Functions inherited from PixelCPEBase
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 ()

Private Member Functions
void	callTempReco1D (DetParam const &theDetParam, ClusterParamTemplate &theClusterParam, SiPixelTemplateReco::ClusMatrix &clusterPayload, int ID, LocalPoint &lp) const

void	callTempReco2D (DetParam const &theDetParam, ClusterParamTemplate &theClusterParam, SiPixelTemplateReco2D::ClusMatrix &clusterPayload, int ID, LocalPoint &lp) const

void	checkRecommend2D (DetParam const &theDetParam, ClusterParamTemplate &theClusterParam, SiPixelTemplateReco::ClusMatrix &clusterPayload, int ID) const

ClusterParam *	createClusterParam (const SiPixelCluster &cl) const override

void	fill2DTemplIDs ()

LocalError	localError (DetParam const &theDetParam, ClusterParam &theClusterParam) const override

LocalPoint	localPosition (DetParam const &theDetParam, ClusterParam &theClusterParam) const override

Private Attributes
int	barrelTemplateID_

int	forwardTemplateID_

float	maxSizeMismatchInY_

float	minChargeRatio_

std::vector< Rule >	recommend2D_

bool	runDamagedClusters_

int	speed_

const SiPixel2DTemplateDBObject *	templateDBobject2D_

std::string	templateDir_

std::vector< SiPixelTemplateStore2D >	thePixelTemp2D_

std::vector< SiPixelTemplateStore >	thePixelTemp_

bool	UseClusterSplitter_

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

Protected Types inherited from PixelCPEBase
using	DetParams = std::vector< DetParam >

typedef GloballyPositioned< double >	Frame

Protected Member Functions inherited from PixelCPEBase
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

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

Protected Attributes inherited from PixelCPEBase
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_

Detailed Description

Definition at line 35 of file PixelCPEClusterRepair.h.

Constructor & Destructor Documentation

PixelCPEClusterRepair::PixelCPEClusterRepair	(	edm::ParameterSet const &	conf,
		const MagneticField *	mag,
		const TrackerGeometry &	geom,
		const TrackerTopology &	ttopo,
		const SiPixelLorentzAngle *	lorentzAngle,
		const SiPixelTemplateDBObject *	templateDBobject,
		const SiPixel2DTemplateDBObject *	templateDBobject2D
	)

Definition at line 40 of file PixelCPEClusterRepair.cc.

References barrelTemplateID_, Exception, fill2DTemplIDs(), forwardTemplateID_, edm::ParameterSet::getParameter(), PixelCPEBase::LoadTemplatesFromDB_, LogDebug, maxSizeMismatchInY_, minChargeRatio_, SiPixelTemplate2D::pushfile(), SiPixelTemplate::pushfile(), recommend2D_, runDamagedClusters_, speed_, str, templateDBobject2D_, PixelCPEBase::templateDBobject_, templateDir_, thePixelTemp2D_, thePixelTemp_, and UseClusterSplitter_.

 : PixelCPEBase(conf, mag, geom, ttopo, lorentzAngle, nullptr, templateDBobject, nullptr,1)
 {
    LogDebug("PixelCPEClusterRepair::(constructor)") << endl;
 
    //--- Parameter to decide between DB or text file template access
    if ( LoadTemplatesFromDB_ )
    {
       // Initialize template store to the selected ID [Morris, 6/25/08]
       if ( !SiPixelTemplate::pushfile( *templateDBobject_, thePixelTemp_) )
          throw cms::Exception("PixelCPEClusterRepair")
          << "\nERROR: Templates not filled correctly. Check the sqlite file. Using SiPixelTemplateDBObject version "
          << (*templateDBobject_).version() << "\n\n";
 
       // Initialize template store to the selected ID [Morris, 6/25/08]
       if ( !SiPixelTemplate2D::pushfile( *templateDBobject2D , thePixelTemp2D_) )
          throw cms::Exception("PixelCPEClusterRepair")
          << "\nERROR: Templates not filled correctly. Check the sqlite file. Using SiPixelTemplateDBObject2D version "
          << (*templateDBobject2D).version() << "\n\n";
    }
    else
    {
       LogDebug("PixelCPEClusterRepair") << "Loading templates for barrel and forward from ASCII files." << endl;
       //--- (Archaic) Get configurable template IDs.  This code executes only if we loading pixels from ASCII
       //    files, and then they are mandatory.
       barrelTemplateID_  = conf.getParameter<int>( "barrelTemplateID" );
       forwardTemplateID_ = conf.getParameter<int>( "forwardTemplateID" );
       templateDir_       = conf.getParameter<int>( "directoryWithTemplates" );
      
       if ( !SiPixelTemplate::pushfile( barrelTemplateID_  , thePixelTemp_ , templateDir_ ) )
          throw cms::Exception("PixelCPEClusterRepair")
      << "\nERROR: Template ID " << barrelTemplateID_ << " not loaded correctly from text file. Reconstruction will fail.\n\n";
       
       if ( !SiPixelTemplate::pushfile( forwardTemplateID_ , thePixelTemp_ , templateDir_ ) )
          throw cms::Exception("PixelCPEClusterRepair")
      << "\nERROR: Template ID " << forwardTemplateID_ << " not loaded correctly from text file. Reconstruction will fail.\n\n";
    }
    
    templateDBobject2D_ = templateDBobject2D;
    fill2DTemplIDs();
    speed_ = conf.getParameter<int>( "speed");
    LogDebug("PixelCPEClusterRepair::PixelCPEClusterRepair:") <<
    "Template speed = " << speed_ << "\n";
    
    UseClusterSplitter_ = conf.getParameter<bool>("UseClusterSplitter");   
 
 
 
    maxSizeMismatchInY_ = conf.getParameter<double>("MaxSizeMismatchInY");
    minChargeRatio_ = conf.getParameter<double>("MinChargeRatio");
 
    // read sub-detectors and apply rule to recommend 2D 
    // can be: 
    //     XYX (XYZ = PXB, PXE)
    //     XYZ n (XYZ as above, n = layer, wheel or disk = 1 .. 6 ;)
    std::vector<std::string> str_recommend2D = conf.getParameter<std::vector<std::string> >("Recommend2D");
    recommend2D_.reserve(str_recommend2D.size());
    for (auto & str: str_recommend2D){
      recommend2D_.push_back(str);
    }
 
    // run CR on damaged clusters (and not only on edge hits)
    runDamagedClusters_ = conf.getParameter<bool>("RunDamagedClusters");
 }

PixelCPEClusterRepair::~PixelCPEClusterRepair ( )

override

Definition at line 153 of file PixelCPEClusterRepair.cc.

References thePixelTemp2D_, thePixelTemp_, and x.

 {
    for (auto x : thePixelTemp_)   x.destroy();
    for (auto x : thePixelTemp2D_) x.destroy();
 }

Member Function Documentation

void PixelCPEClusterRepair::callTempReco1D	(	DetParam const &	theDetParam,
		ClusterParamTemplate &	theClusterParam,
		SiPixelTemplateReco::ClusMatrix &	clusterPayload,
		int	ID,
		LocalPoint &	lp
	)		const

private

Definition at line 326 of file PixelCPEClusterRepair.cc.

Referenced by localPosition().

 {
    SiPixelTemplate templ(thePixelTemp_);
    
    // Output:
    float nonsense = -99999.9f; // nonsense init value
    theClusterParam.templXrec_ = theClusterParam.templYrec_ = theClusterParam.templSigmaX_ = theClusterParam.templSigmaY_ = nonsense;
    // If the template recontruction fails, we want to return 1.0 for now
    theClusterParam.probabilityX_ = theClusterParam.probabilityY_ = theClusterParam.probabilityQ_ = 1.f;
    theClusterParam.qBin_ = 0;
    // We have a boolean denoting whether the reco failed or not
    theClusterParam.hasFilledProb_ = false;
    
    // ******************************************************************
    //--- Call normal TemplateReco
    //
    float locBz = theDetParam.bz;
    float locBx = theDetParam.bx;
    //
    const bool deadpix = false;
    std::vector<std::pair<int, int> > zeropix;
    int nypix =0, nxpix = 0;
    //
    theClusterParam.ierr =
    PixelTempReco1D( ID, theClusterParam.cotalpha, theClusterParam.cotbeta,
                    locBz, locBx,
                    clusterPayload,
                    templ,
                    theClusterParam.templYrec_, theClusterParam.templSigmaY_, theClusterParam.probabilityY_,
                    theClusterParam.templXrec_, theClusterParam.templSigmaX_, theClusterParam.probabilityX_,
                    theClusterParam.qBin_,
                    speed_, deadpix, zeropix,
                    theClusterParam.probabilityQ_, nypix, nxpix
                    );
    // ******************************************************************
    
    //--- Check exit status
    if UNLIKELY( theClusterParam.ierr != 0 )
    {
       LogDebug("PixelCPEClusterRepair::localPosition") <<
       "reconstruction failed with error " << theClusterParam.ierr << "\n";
 
       theClusterParam.probabilityX_ = theClusterParam.probabilityY_ = theClusterParam.probabilityQ_ = 0.f;
       theClusterParam.qBin_ = 0;
       
       // 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 PixelTempReco1D fails
       float lorentz_drift = -999.9;
       if ( ! GeomDetEnumerators::isEndcap(theDetParam.thePart) )
          lorentz_drift = 60.0f; // in microns
       else
          lorentz_drift = 10.0f; // in microns
       // GG: trk angles needed to correct for bows/kinks
       if ( theClusterParam.with_track_angle )
       {
          theClusterParam.templXrec_ = theDetParam.theTopol->localX( theClusterParam.theCluster->x(), theClusterParam.loc_trk_pred ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
          theClusterParam.templYrec_ = theDetParam.theTopol->localY( theClusterParam.theCluster->y(), theClusterParam.loc_trk_pred );
       }
       else
       {
          edm::LogError("PixelCPEClusterRepair")
          << "@SUB = PixelCPEClusterRepair::localPosition"
          << "Should never be here. PixelCPEClusterRepair should always be called with track angles. This is a bad error !!! ";
          
          theClusterParam.templXrec_ = theDetParam.theTopol->localX( theClusterParam.theCluster->x() ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
          theClusterParam.templYrec_ = theDetParam.theTopol->localY( theClusterParam.theCluster->y() );
       }
    }   
    else 
    {
        //--- Template Reco succeeded.  The probabilities are filled.
        theClusterParam.hasFilledProb_ = true;
 
        
 
            //--- Go from microns to centimeters
       theClusterParam.templXrec_ *= micronsToCm;
       theClusterParam.templYrec_ *= micronsToCm;
       
       //--- Go back to the module coordinate system
       theClusterParam.templXrec_ += lp.x();
       theClusterParam.templYrec_ += lp.y();
             
    }
    return;
 }

void PixelCPEClusterRepair::callTempReco2D	(	DetParam const &	theDetParam,
		ClusterParamTemplate &	theClusterParam,
		SiPixelTemplateReco2D::ClusMatrix &	clusterPayload,
		int	ID,
		LocalPoint &	lp
	)		const

private

Definition at line 424 of file PixelCPEClusterRepair.cc.

Referenced by localPosition().

 {
    SiPixelTemplate2D templ2d(thePixelTemp2D_);
    
    // Output:
    float nonsense = -99999.9f; // nonsense init value
    theClusterParam.templXrec_ = theClusterParam.templYrec_ = theClusterParam.templSigmaX_ = theClusterParam.templSigmaY_ = nonsense;
    // If the template recontruction fails, we want to return 1.0 for now
    theClusterParam.probabilityX_ = theClusterParam.probabilityY_ = theClusterParam.probabilityQ_ = 1.f;
    theClusterParam.qBin_ = 0;
    // We have a boolean denoting whether the reco failed or not
    theClusterParam.hasFilledProb_ = false;
    
    // ******************************************************************
    //--- Call 2D TemplateReco
    //
    float locBz = theDetParam.bz;
    float locBx = theDetParam.bx;
 
    //--- Input:
    //   edgeflagy - (input) flag to indicate the present of edges in y: 
    //           0=none (or interior gap),1=edge at small y, 2=edge at large y, 3=edge at either end
    //
    //   edgeflagx - (input) flag to indicate the present of edges in x: 
    //           0=none, 1=edge at small x, 2=edge at large x
    //
    //   These two variables are calculated in setTheClu() and stored in edgeTypeX_ and edgeTypeY_
    //
    //--- Output:
    //   deltay - (output) template y-length - cluster length [when > 0, possibly missing end]
    //   npixels - ???     &&& Ask Morris
 
 
    float deltay = 0;    // return param
    int npixels = 0;     // return param
 
 
    //For now, turn off edgeX_ flags
    theClusterParam.edgeTypeX_ = 0;
 
    if(clusterPayload.mrow > 4){
        // The cluster is too big, the 2D reco will perform horribly.
        // Better to return immediately in error
        theClusterParam.ierr2 = 8;
 
    }
    else{
        theClusterParam.ierr2 =
        PixelTempReco2D( ID, theClusterParam.cotalpha, theClusterParam.cotbeta,
                        locBz, locBx,
                        theClusterParam.edgeTypeY_ , theClusterParam.edgeTypeX_ ,
                        clusterPayload,
                        templ2d,
                        theClusterParam.templYrec_, theClusterParam.templSigmaY_, 
                        theClusterParam.templXrec_, theClusterParam.templSigmaX_, 
                        theClusterParam.templProbXY_,
                        theClusterParam.probabilityQ_,
                        theClusterParam.qBin_,
                        deltay, npixels
                        );
    }
    // ******************************************************************
    
    //--- Check exit status
    if UNLIKELY( theClusterParam.ierr2 != 0 )
    {
       LogDebug("PixelCPEClusterRepair::localPosition") <<
       "2D reconstruction failed with error " << theClusterParam.ierr2 << "\n";
       
       theClusterParam.probabilityX_ = theClusterParam.probabilityY_ = theClusterParam.probabilityQ_ = 0.f;
       theClusterParam.qBin_ = 0;
       // GG: 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
       float lorentz_drift = -999.9;
       if ( ! GeomDetEnumerators::isEndcap(theDetParam.thePart) )
          lorentz_drift = 60.0f; // in microns  // &&& replace with a constant (globally)
       else
          lorentz_drift = 10.0f; // in microns
       // GG: trk angles needed to correct for bows/kinks
       if ( theClusterParam.with_track_angle )
       {
          theClusterParam.templXrec_ = theDetParam.theTopol->localX( theClusterParam.theCluster->x(), theClusterParam.loc_trk_pred ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
          theClusterParam.templYrec_ = theDetParam.theTopol->localY( theClusterParam.theCluster->y(), theClusterParam.loc_trk_pred );
       }
       else
       {
          edm::LogError("PixelCPEClusterRepair")
          << "@SUB = PixelCPEClusterRepair::localPosition"
          << "Should never be here. PixelCPEClusterRepair should always be called with track angles. This is a bad error !!! ";
          
          theClusterParam.templXrec_ = theDetParam.theTopol->localX( theClusterParam.theCluster->x() ) - lorentz_drift * micronsToCm; // rough Lorentz drift correction
          theClusterParam.templYrec_ = theDetParam.theTopol->localY( theClusterParam.theCluster->y() );
       }
    }   
    else 
    {
       //--- Template Reco succeeded.
       theClusterParam.hasFilledProb_ = true;
 
       //--- Go from microns to centimeters
       theClusterParam.templXrec_ *= micronsToCm;
       theClusterParam.templYrec_ *= micronsToCm;
       
       //--- Go back to the module coordinate system
       theClusterParam.templXrec_ += lp.x();
       theClusterParam.templYrec_ += lp.y();
    }
    return;
 }

void PixelCPEClusterRepair::checkRecommend2D	(	DetParam const &	theDetParam,
		ClusterParamTemplate &	theClusterParam,
		SiPixelTemplateReco::ClusMatrix &	clusterPayload,
		int	ID
	)		const

private

Definition at line 540 of file PixelCPEClusterRepair.cc.

References PixelCPEBase::DetParam::bx, PixelCPEBase::DetParam::bz, SiPixelCluster::charge(), SiPixelTemplate::clsleny(), PixelCPEBase::ClusterParam::cotalpha, PixelCPEBase::ClusterParam::cotbeta, PixelCPEBase::ClusterParam::edgeTypeY_, GeomDet::geographicalId(), PixelCPEBase::ClusterParam::hasBadPixels_, mps_fire::i, SiPixelTemplate::interpolate(), maxSizeMismatchInY_, SiPixelTemplateReco::ClusMatrix::mcol, minChargeRatio_, SiPixelCluster::minPixelCol(), SiPixelTemplate::qavg(), recommend2D_, PixelCPEClusterRepair::ClusterParamTemplate::recommended2D_, runDamagedClusters_, SiPixelCluster::sizeY(), PixelCPEBase::ClusterParam::theCluster, PixelCPEBase::DetParam::theDet, thePixelTemp_, PixelCPEBase::ttopo_, and SiPixelTemplateReco::ClusMatrix::ydouble.

Referenced by localPosition().

 {
     // recommend2D is false by default
     theClusterParam.recommended2D_  = false;
   
     DetId id = (theDetParam.theDet->geographicalId());
 
     bool recommend = false;
     for (auto & rec: recommend2D_){
         recommend = rec.recommend(id, ttopo_);
         if(recommend) break;
     }
 
     // only run on those layers recommended by configuration
     if(!recommend) {
         theClusterParam.recommended2D_ = false;
         return;
     }
 
     if(theClusterParam.edgeTypeY_){
         // For clusters that have edges in Y, run 2d reco.
         // Flags already set by CPEBase
         theClusterParam.recommended2D_ = true;
         return;
     }
     // The 1d pixel template
     SiPixelTemplate templ(thePixelTemp_);
     if(!templ.interpolate(ID, theClusterParam.cotalpha, theClusterParam.cotbeta, theDetParam.bz, theDetParam.bx)){
         //error setting up template, return false
         theClusterParam.recommended2D_ = false;
         return;
     }
 
 
     //length of the cluster taking into account double sized pixels
     float nypix = clusterPayload.mcol;
     for(int i=0; i<clusterPayload.mcol; i++){
         if(clusterPayload.ydouble[i]) nypix+=1.;
     }
 
     // templ.clsleny() is the expected length of the cluster along y axis.
     // templ.qavg() is the expected total charge of the cluster
     // theClusterParam.theCluster->charge() is the total charge of this cluster
     float nydiff = templ.clsleny() - nypix;
     float qratio = theClusterParam.theCluster->charge()/templ.qavg();
 
     if ( nydiff > maxSizeMismatchInY_  && qratio < minChargeRatio_){
         // If the cluster is shorter than expected and has less charge, likely
         // due to truncated cluster, try 2D reco
         
         theClusterParam.recommended2D_ = true;
         theClusterParam.hasBadPixels_ = true;
 
     // if not RunDamagedClusters flag, don't try to fix any clusters
     if(!runDamagedClusters_) {
         theClusterParam.recommended2D_ = false;
     }
 
         // Figure out what edge flags to set for truncated cluster
         // Truncated clusters usually come from dead double columns
         // 
         // If cluster is of even length,  either both of or neither of beginning and ending
         // edge are on a double column, so we cannot figure out the likely edge of
         // truncation, let the 2D algorithm try extending on both sides (option 3)
         if(theClusterParam.theCluster->sizeY() % 2 == 0) theClusterParam.edgeTypeY_ = 3;
         else{
             //If the cluster is of odd length, only one of the edges can end on
             //a double column, this is the likely edge of truncation
             //Double columns always start on even indexes
             int min_col = theClusterParam.theCluster->minPixelCol();
             if(min_col %2 ==0){
                 //begining edge is at a double column (end edge cannot be,
                 //because odd length) so likely truncated at small y (option 1) 
                 theClusterParam.edgeTypeY_ = 1;
             }
             else{ 
                 //end edge is at a double column (beginning edge cannot be,
                 //because odd length) so likely truncated at large y (option 2) 
                 theClusterParam.edgeTypeY_ = 2;
             }
         }
     }
 }

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

overrideprivatevirtual

Implements PixelCPEBase.

Definition at line 159 of file PixelCPEClusterRepair.cc.

 {
    return new ClusterParamTemplate(cl);
 }

void PixelCPEClusterRepair::fill2DTemplIDs ( )

private

Definition at line 114 of file PixelCPEClusterRepair.cc.

References TrackerGeometry::detUnits(), PixelCPEBase::geom_, SiPixel2DTemplateDBObject::getTemplateID(), mps_fire::i, LogDebug, PixelCPEBase::m_DetParams, TrackerGeometry::offsetDU(), AlCaHLTBitMon_ParallelJobs::p, templateDBobject2D_, and GeomDetEnumerators::tkDetEnum.

Referenced by PixelCPEClusterRepair().

 {
    auto const & dus = geom_.detUnits();
    unsigned m_detectors = dus.size();
    for(unsigned int i=1;i<7;++i) {
       LogDebug("PixelCPEClusterRepair:: 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("PixelCPEClusterRepair::fillDetParams():") <<"caching "<<m_detectors<<" pixel detectors"<<endl;
    //Loop through detectors and store 2D template ID
    bool printed_info = false;
    for (unsigned i=0; i!=m_detectors;++i) {
       auto & p=m_DetParams[i];
 
       p.detTemplateId2D = templateDBobject2D_->getTemplateID(p.theDet->geographicalId());
       if(p.detTemplateId != p.detTemplateId2D && !printed_info){
           edm::LogWarning("PixelCPEClusterRepair") << "different template ID between 1D and 2D " << p.detTemplateId << " " << p.detTemplateId2D << endl;
           printed_info = true;
       }
    }
 
 }

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

overrideprivatevirtual

Implements PixelCPEBase.

Definition at line 633 of file PixelCPEClusterRepair.cc.

References PixelCPEBase::ClusterParam::edgeTypeX_, PixelCPEBase::ClusterParam::edgeTypeY_, Exception, f, PixelCPEClusterRepair::ClusterParamTemplate::ierr, PixelCPEClusterRepair::ClusterParamTemplate::ierr2, GeomDetEnumerators::isBarrel(), GeomDetEnumerators::isTrackerPixel(), LogDebug, SiPixelCluster::sizeX(), SiPixelCluster::sizeY(), PixelCPEClusterRepair::ClusterParamTemplate::templSigmaX_, PixelCPEClusterRepair::ClusterParamTemplate::templSigmaY_, PixelCPEBase::ClusterParam::theCluster, PixelCPEBase::DetParam::thePart, PixelCPEBase::theVerboseLevel, and UNLIKELY.

 {
    
    ClusterParamTemplate & theClusterParam = static_cast<ClusterParamTemplate &>(theClusterParamBase);
    
    //--- Default is the maximum error used for edge clusters.
    //--- (never used, in fact: let comment it out, shut up the complains of the static analyzer, and save a few CPU cycles)
    float xerr = 0.0f, yerr = 0.0f;
 
    //--- Check status of both template calls.
    if UNLIKELY ( (theClusterParam.ierr !=0) || (theClusterParam.ierr2 !=0) ) {
      // If reconstruction fails the hit position is calculated from cluster center of gravity
      // corrected in x by average Lorentz drift. Assign huge errors.
      //
      if UNLIKELY (!GeomDetEnumerators::isTrackerPixel(theDetParam.thePart))
        throw cms::Exception("PixelCPEClusterRepair::localPosition :")
      << "A non-pixel detector type in here?";
      
      // Assign better errors based on the residuals for failed template cases
      if ( GeomDetEnumerators::isBarrel(theDetParam.thePart)) {
          xerr = 55.0f * micronsToCm;      // &&& get errors from elsewhere?
      yerr = 36.0f * micronsToCm;
        }
        else {
      xerr = 42.0f * micronsToCm;
      yerr = 39.0f * micronsToCm;
        }
    }
    // Leave commented for now, until we study the interplay of failure modes
    // of 1D template reco and edges.  For edge hits we run 2D reco by default!
    //
    // else if ( theClusterParam.edgeTypeX_ || theClusterParam.edgeTypeY_ )  {
    //   // for edge pixels assign errors according to observed residual RMS
    //   if      ( theClusterParam.edgeTypeX_ && !theClusterParam.edgeTypeY_ ) {
    //     xerr = 23.0f * micronsToCm;
    //     yerr = 39.0f * micronsToCm;
    //   }
    //   else if ( !theClusterParam.edgeTypeX_ && theClusterParam.edgeTypeY_ ) {
    //     xerr = 24.0f * micronsToCm;
    //     yerr = 96.0f * micronsToCm;
    //   }
    //   else if ( theClusterParam.edgeTypeX_ && theClusterParam.edgeTypeY_ ) {
    //     xerr = 31.0f * micronsToCm;
    //     yerr = 90.0f * micronsToCm;
    //   }
    // }
    else {
      xerr = theClusterParam.templSigmaX_ * micronsToCm;
      yerr = theClusterParam.templSigmaY_ * micronsToCm;
      // &&& should also check ierr (saved as class variable) and return
      // &&& nonsense (another class static) if the template fit failed.
    }       
    
    if (theVerboseLevel > 9) {
      LogDebug("PixelCPEClusterRepair") 
        << " Sizex = " << theClusterParam.theCluster->sizeX() 
        << " Sizey = " << theClusterParam.theCluster->sizeY() 
        << " Edgex = " << theClusterParam.edgeTypeX_ 
        << " Edgey = " << theClusterParam.edgeTypeY_ 
        << " ErrX  = " << xerr << " ErrY  = " << yerr;
    }
    
    
    if ( !(xerr > 0.0f) )
       throw cms::Exception("PixelCPEClusterRepair::localError") 
       << "\nERROR: Negative pixel error xerr = " << xerr << "\n\n";
    
    if ( !(yerr > 0.0f) )
       throw cms::Exception("PixelCPEClusterRepair::localError") 
       << "\nERROR: Negative pixel error yerr = " << yerr << "\n\n";
    
    return LocalError(xerr*xerr, 0, yerr*yerr);
 }

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

overrideprivatevirtual

Implements PixelCPEBase.

Definition at line 174 of file PixelCPEClusterRepair.cc.

 {
    
    ClusterParamTemplate & theClusterParam = static_cast<ClusterParamTemplate &>(theClusterParamBase);
    bool filled_from_2d = false;
    
    if(!GeomDetEnumerators::isTrackerPixel(theDetParam.thePart))
       throw cms::Exception("PixelCPEClusterRepair::localPosition :")
       << "A non-pixel detector type in here?";
 
    
    int ID1 = -9999;
    int ID2 = -9999;
    if ( LoadTemplatesFromDB_ ) {
       ID1 = theDetParam.detTemplateId;
       ID2 = theDetParam.detTemplateId2D;
    } else { // from asci file
       if ( ! GeomDetEnumerators::isEndcap(theDetParam.thePart) )
     ID1 = ID2 = barrelTemplateID_  ; // barrel
       else
     ID1 = ID2 = forwardTemplateID_ ; // forward
    }
 
    // &&& PM, note for later: PixelCPEBase calculates minInX,Y, and maxInX,Y
    //     Why can't we simply use that and save time with row_offset, col_offset
    //     and mrow = maxInX-minInX, mcol = maxInY-minInY ... Except that we
    //     also need to take into account cluster_matrix_size_x,y.
 
    
    //--- Preparing to retrieve ADC counts from the SiPixeltheClusterParam.theCluster
    //    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 = theClusterParam.theCluster->minPixelRow();
    int col_offset = theClusterParam.theCluster->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 TemplReco2D
    float tmp_x = float(row_offset) + 0.5f;
    float tmp_y = float(col_offset) + 0.5f;
 
    
    //--- Store these offsets (to be added later) in a LocalPoint after tranforming
    //    them from measurement units (pixel units) to local coordinates (cm)
    LocalPoint lp;
    if ( theClusterParam.with_track_angle )
       //--- Update call with trk angles needed for bow/kink corrections  // Gavril
       lp = theDetParam.theTopol->localPosition( MeasurementPoint(tmp_x, tmp_y), theClusterParam.loc_trk_pred );
    else
    {
       edm::LogError("PixelCPEClusterRepair")
       << "@SUB = PixelCPEClusterRepair::localPosition"
       << "Should never be here. PixelCPEClusterRepair should always be called with track angles. This is a bad error !!! ";
       lp = theDetParam.theTopol->localPosition( MeasurementPoint(tmp_x, tmp_y) );
    }
    
    //--- Compute the size of the matrix which will be passed to TemplateReco.
    //    We'll later make  clustMatrix[ mrow ][ mcol ]
    int mrow=0, mcol=0;
    for (int i=0 ; i!=theClusterParam.theCluster->size(); ++i )
    {
       auto pix = theClusterParam.theCluster->pixel(i);
       int irow = int(pix.x);
       int icol = int(pix.y);
       mrow = std::max(mrow,irow);
       mcol = std::max(mcol,icol);
    }
    mrow -= row_offset; mrow+=1; mrow = std::min(mrow,cluster_matrix_size_x);
    mcol -= col_offset; mcol+=1; mcol = std::min(mcol,cluster_matrix_size_y);
    assert(mrow>0); assert(mcol>0);
 
 
    //--- Make and fill the bool arrays flagging double pixels
    bool xdouble[mrow], ydouble[mcol];
    // x directions (shorter), rows
    for (int irow = 0; irow < mrow; ++irow)
       xdouble[irow] = theDetParam.theRecTopol->isItBigPixelInX( irow+row_offset );
    //
    // y directions (longer), columns
    for (int icol = 0; icol < mcol; ++icol)
       ydouble[icol] = theDetParam.theRecTopol->isItBigPixelInY( icol+col_offset );
 
    //--- C-style matrix.  We'll need it in either case.
    float clustMatrix[mrow][mcol];
    float clustMatrix2[mrow][mcol];
 
    //--- Prepare struct that passes pointers to TemplateReco.  It doesn't own anything.
    SiPixelTemplateReco::ClusMatrix   clusterPayload  { &clustMatrix[0][0], xdouble, ydouble, mrow,mcol};
    SiPixelTemplateReco2D::ClusMatrix clusterPayload2d{ &clustMatrix2[0][0], xdouble, ydouble, mrow,mcol};
 
 
    //--- Copy clust's pixels (calibrated in electrons) into clustMatrix;
    memset( clustMatrix, 0, sizeof(float)*mrow*mcol );   // Wipe it clean.
    for (int i=0 ; i!=theClusterParam.theCluster->size(); ++i )
    {
       auto pix = theClusterParam.theCluster->pixel(i);
       int irow = int(pix.x) - row_offset;
       int icol = int(pix.y) - col_offset;
       // &&& Do we ever get pixels that are out of bounds ???  Need to check.
       if ( (irow<mrow) & (icol<mcol) ) clustMatrix[irow][icol] =  float(pix.adc);
    }
    // &&& Save for later: fillClustMatrix( float * clustMatrix );
 
    //--- Save a copy of clustMatrix into clustMatrix2
    memcpy( clustMatrix2, clustMatrix, sizeof(float)*mrow*mcol);
 
    //--- Set both return statuses, since we may call only one.
    theClusterParam.ierr  = 0;
    theClusterParam.ierr2 = 0;
 
 
 
    //--- Should we run the 2D reco?
    checkRecommend2D(theDetParam, theClusterParam, clusterPayload, ID1);
    if ( theClusterParam.recommended2D_ ) {
      //--- Call the Template Reco 2d with cluster repair
      filled_from_2d = true;
      callTempReco2D( theDetParam, theClusterParam, clusterPayload2d, ID2, lp );
    }
    else {
      //--- Call the vanilla Template Reco
      callTempReco1D( theDetParam, theClusterParam, clusterPayload, ID1, lp );
      filled_from_2d = false;
    }
 
    //--- Make sure cluster repair returns all info about the hit back up to caller
    //--- Necessary because it copied the base class so it does not modify it
    theClusterParamBase.isOnEdge_ = theClusterParam.isOnEdge_;
    theClusterParamBase.hasBadPixels_ = theClusterParam.hasBadPixels_;
    theClusterParamBase.spansTwoROCs_ = theClusterParam.spansTwoROCs_;
    theClusterParamBase.hasFilledProb_ = theClusterParam.hasFilledProb_;
    theClusterParamBase.qBin_ = theClusterParam.qBin_;
    theClusterParamBase.probabilityQ_ = theClusterParam.probabilityQ_;
    theClusterParamBase.filled_from_2d = filled_from_2d;
    if(filled_from_2d){
        theClusterParamBase.probabilityX_ = theClusterParam.templProbXY_;
        theClusterParamBase.probabilityY_ = 0.;
     }
    else{
        theClusterParamBase.probabilityX_ = theClusterParam.probabilityX_;
        theClusterParamBase.probabilityY_ = theClusterParam.probabilityY_;
    }
 
 
    
    return LocalPoint( theClusterParam.templXrec_, theClusterParam.templYrec_ );
 }