PixelCPEBase Class Reference

#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 ()

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

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
void	setTheClu (DetParam const &, ClusterParam &theClusterParam) const

Protected Attributes
bool	alpha2Order
bool	doLorentzFromAlignment_
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	useLAFromDB_
bool	useLAOffsetFromConfig_
bool	useLAWidthFromConfig_
bool	useLAWidthFromDB_
Protected Attributes inherited from PixelClusterParameterEstimator
unsigned int	clusterProbComputationFlag_

Static Protected Attributes
static constexpr float	bothEdgeXError_ = 31.0f
static constexpr float	bothEdgeYError_ = 90.0f
static constexpr float	clusterSplitMaxError_ = 7777.7f
static constexpr float	xEdgeXError_ = 23.0f
static constexpr float	xEdgeYError_ = 39.0f
static constexpr float	yEdgeXError_ = 24.0f
static constexpr float	yEdgeYError_ = 96.0f

Private Member Functions
virtual std::unique_ptr< 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 43 of file PixelCPEBase.h.

Member Typedef Documentation

DetParams

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

protected

Definition at line 279 of file PixelCPEBase.h.

Frame

typedef GloballyPositioned<double> PixelCPEBase::Frame

protected

Definition at line 203 of file PixelCPEBase.h.

Constructor & Destructor Documentation

PixelCPEBase()

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 30 of file PixelCPEBase.cc.

References alpha2Order, PixelClusterParameterEstimator::clusterProbComputationFlag_, doLorentzFromAlignment_, fillDetParams(), genErrorDBObject_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), createfilelist::int, lAOffset_, lAWidthBPix_, lAWidthFPix_, LoadTemplatesFromDB_, LogDebug, lorentzAngle_, lorentzAngleWidth_, templateDBobject_, theFlag_, theVerboseLevel, useLAFromDB_, 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 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.getParameter<bool>("useLAWidthFromDB");

  // 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.getParameter<double>("lAOffset");
  lAWidthBPix_ = conf.getParameter<double>("lAWidthBPix");
  lAWidthFPix_ = conf.getParameter<double>("lAWidthFPix");

  // 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)
  doLorentzFromAlignment_ = conf.getParameter<bool>("doLorentzFromAlignment");
  useLAFromDB_ = conf.getParameter<bool>("useLAFromDB");

  LogDebug("PixelCPEBase") << " LA constants - " << lAOffset_ << " " << lAWidthBPix_ << " " << lAWidthFPix_
                           << endl;  //dk

  fillDetParams();
}

Member Function Documentation

computeAnglesFromDetPosition()

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

protected

Definition at line 291 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 getParameters().

                                                                                                                {
  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);
    */
}

computeAnglesFromTrajectory()

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

protected

Definition at line 257 of file PixelCPEBase.cc.

References 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().

                                                                                           {
  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);
}

computeLorentzShifts()

void PixelCPEBase::computeLorentzShifts ( DetParam &  theDetParam ) const

protected

Definition at line 423 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;
}

createClusterParam()

virtual std::unique_ptr<ClusterParam> PixelCPEBase::createClusterParam ( const SiPixelCluster &  cl ) const

privatepure virtual

Implemented in PixelCPEClusterRepair, PixelCPETemplateReco, and PixelCPEGenericBase.

Referenced by getParameters().

detParam()

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

protected

Definition at line 322 of file PixelCPEBase.cc.

References cms::cuda::assert(), 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;
}

driftDirection() [1/2]

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

protected

Definition at line 338 of file PixelCPEBase.cc.

References mkfit::Config::Bfield, 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);
}

driftDirection() [2/2]

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

protected

Definition at line 344 of file PixelCPEBase.cc.

References funct::abs(), alpha2Order, mkfit::Config::Bfield, createTree::dd, f, GeomDet::geographicalId(), SiPixelLorentzAngle::getLorentzAngle(), GeomDetEnumerators::isBarrel(), GeomDetEnumerators::isTrackerPixel(), lAOffset_, lAWidthBPix_, lAWidthFPix_, LogDebug, lorentzAngle_, lorentzAngleWidth_, DetId::rawId(), L1EGammaClusterEmuProducer_cfi::scale, PixelCPEBase::DetParam::theDet, theFlag_, PixelCPEBase::DetParam::thePart, useLAOffsetFromConfig_, useLAWidthFromConfig_, useLAWidthFromDB_, PixelCPEBase::DetParam::widthLAFractionX, and PixelCPEBase::DetParam::widthLAFractionY.

                                                                                        {
  // 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;
}

fillDetParams()

void PixelCPEBase::fillDetParams ( )

private

Definition at line 121 of file PixelCPEBase.cc.

References cms::cuda::assert(), mkfit::Config::Bfield, computeLorentzShifts(), TrackerGeometry::detUnits(), doLorentzFromAlignment_, 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, and useLAFromDB_.

Referenced by PixelCPEBase().

                                 {
  // MM: this code finds the last Pixel (Inner Tracker) DetUnit to loop upon when filling the det params, by looking the first detId which is from
  // the Outer Tracker (Strips in case of the Phase-0/1 detector).

  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().isOuterTracker()
                : false);

    if (geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i]) != dus.size() &&
        dus[geom_.offsetDU(GeomDetEnumerators::tkDetEnum[i])]->type().isOuterTracker()) {
      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());
    }

    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

    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) || useLAFromDB_ ||
        doLorentzFromAlignment_) {  // do always for generic and if using LA from DB or alignment 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;
  }
}

fillPSetDescription()

void PixelCPEBase::fillPSetDescription ( edm::ParameterSetDescription &  desc )

static

Definition at line 467 of file PixelCPEBase.cc.

References submitPVResolutionJobs::desc.

Referenced by PixelCPEFastESProducer::fillDescriptions(), PixelCPETemplateRecoESProducer::fillDescriptions(), PixelCPEClusterRepairESProducer::fillDescriptions(), and PixelCPEGenericESProducer::fillDescriptions().

                                                                       {
  desc.add<bool>("LoadTemplatesFromDB", true);
  desc.add<bool>("Alpha2Order", true);
  desc.add<int>("ClusterProbComputationFlag", 0);
  desc.add<bool>("useLAWidthFromDB", true);
  desc.add<double>("lAOffset", 0.0);
  desc.add<double>("lAWidthBPix", 0.0);
  desc.add<double>("lAWidthFPix", 0.0);
  desc.add<bool>("doLorentzFromAlignment", false);
  desc.add<bool>("useLAFromDB", true);
}

getParameters() [1/2]

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

inlineoverridevirtual

Implements PixelClusterParameterEstimator.

Definition at line 133 of file PixelCPEBase.h.

References haddnano::cl, computeAnglesFromDetPosition(), createClusterParam(), detParam(), localError(), localPosition(), LogDebug, rawQualityWord(), setTheClu(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

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

                                                                                                   {
#ifdef EDM_ML_DEBUG
    nRecHitsTotal_++;
    LogDebug("PixelCPEBase") << " in PixelCPEBase:localParameters(all) - " << nRecHitsTotal_;
#endif

    DetParam const& theDetParam = detParam(det);
    std::unique_ptr<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);

    LogDebug("PixelCPEBase") << " in PixelCPEBase:localParameters(all) - " << lp.x() << " " << lp.y();
    return tuple;
  }

getParameters() [2/2]

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

inlineoverridevirtual

Implements PixelClusterParameterEstimator.

Definition at line 157 of file PixelCPEBase.h.

References haddnano::cl, computeAnglesFromTrajectory(), createClusterParam(), detParam(), localError(), localPosition(), LogDebug, rawQualityWord(), setTheClu(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                       {
#ifdef EDM_ML_DEBUG
    nRecHitsTotal_++;
    LogDebug("PixelCPEBase") << " in PixelCPEBase:localParameters(on track) - " << nRecHitsTotal_;
#endif

    DetParam const& theDetParam = detParam(det);
    std::unique_ptr<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);

    LogDebug("PixelCPEBase") << " in PixelCPEBase:localParameters(on track) - " << lp.x() << " " << lp.y();
    return tuple;
  }

localError()

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

privatepure virtual

Implemented in PixelCPETemplateReco, PixelCPEClusterRepair, PixelCPEGeneric, and PixelCPEFast.

Referenced by getParameters().

localPosition()

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

privatepure virtual

Implemented in PixelCPETemplateReco, PixelCPEClusterRepair, PixelCPEGeneric, PixelCPEGenericForBricked, and PixelCPEFast.

Referenced by getParameters().

rawQualityWord()

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 441 of file PixelCPEBase.cc.

References PixelCPEBase::ClusterParam::filled_from_2d, PixelCPEBase::ClusterParam::hasBadPixels_, PixelCPEBase::ClusterParam::hasFilledProb_, PixelCPEBase::ClusterParam::isOnEdge_, dqm-mbProfile::log, PixelCPEBase::ClusterParam::probabilityQ_, PixelCPEBase::ClusterParam::probabilityX_, PixelCPEBase::ClusterParam::probabilityY_, PixelCPEBase::ClusterParam::qBin_, SiPixelRecHitQuality::Packing::setHasBadPixels(), SiPixelRecHitQuality::Packing::setHasFilledProb(), SiPixelRecHitQuality::Packing::setIsOnEdge(), SiPixelRecHitQuality::Packing::setProbabilityQ(), SiPixelRecHitQuality::Packing::setProbabilityXY(), SiPixelRecHitQuality::Packing::setQBin(), SiPixelRecHitQuality::Packing::setSpansTwoROCs(), PixelCPEBase::ClusterParam::spansTwoROCs_, and SiPixelRecHitQuality::thePacking.

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;
}

setTheClu()

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

protected

Definition at line 214 of file PixelCPEBase.cc.

References RectangularPixelTopology::containsBigPixelInX(), RectangularPixelTopology::containsBigPixelInY(), PixelCPEBase::ClusterParam::edgeTypeX_, PixelCPEBase::ClusterParam::edgeTypeY_, PixelCPEBase::ClusterParam::isOnEdge_, SiPixelCluster::maxPixelCol(), SiPixelCluster::maxPixelRow(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), RectangularPixelTopology::ncolumns(), RectangularPixelTopology::nrows(), 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();

  int min_row(0), min_col(0);
  int max_row = theDetParam.theRecTopol->nrows() - 1;
  int max_col = theDetParam.theRecTopol->ncolumns() - 1;

  if (minInX == min_row)
    theClusterParam.edgeTypeX_ = 1;
  else if (maxInX == max_row)
    theClusterParam.edgeTypeX_ = 2;
  else
    theClusterParam.edgeTypeX_ = 0;

  if (minInY == min_col)
    theClusterParam.edgeTypeY_ = 1;
  else if (maxInY == max_col)
    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);
}

Member Data Documentation

alpha2Order

bool PixelCPEBase::alpha2Order

protected

Definition at line 237 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

bothEdgeXError_

constexpr float PixelCPEBase::bothEdgeXError_ = 31.0f

staticprotected

Definition at line 252 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

bothEdgeYError_

constexpr float PixelCPEBase::bothEdgeYError_ = 90.0f

staticprotected

Definition at line 253 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

clusterSplitMaxError_

constexpr float PixelCPEBase::clusterSplitMaxError_ = 7777.7f

staticprotected

Definition at line 255 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

doLorentzFromAlignment_

bool PixelCPEBase::doLorentzFromAlignment_

protected

Definition at line 241 of file PixelCPEBase.h.

Referenced by fillDetParams(), PixelCPETemplateReco::localPosition(), and PixelCPEBase().

genErrorDBObject_

const SiPixelGenErrorDBObject* PixelCPEBase::genErrorDBObject_

protected

Definition at line 235 of file PixelCPEBase.h.

Referenced by fillDetParams(), PixelCPEBase(), PixelCPEFast::PixelCPEFast(), and PixelCPEGeneric::PixelCPEGeneric().

geom_

const TrackerGeometry& PixelCPEBase::geom_

protected

Definition at line 229 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::fill2DTemplIDs(), and fillDetParams().

lAOffset_

float PixelCPEBase::lAOffset_

protected

Definition at line 216 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

lAWidthBPix_

float PixelCPEBase::lAWidthBPix_

protected

Definition at line 217 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

lAWidthFPix_

float PixelCPEBase::lAWidthFPix_

protected

Definition at line 218 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

LoadTemplatesFromDB_

bool PixelCPEBase::LoadTemplatesFromDB_

protected

Definition at line 242 of file PixelCPEBase.h.

Referenced by fillDetParams(), PixelCPEClusterRepair::localPosition(), PixelCPETemplateReco::localPosition(), PixelCPEBase(), PixelCPEClusterRepair::PixelCPEClusterRepair(), PixelCPEGeneric::PixelCPEGeneric(), PixelCPEGenericForBricked::PixelCPEGenericForBricked(), and PixelCPETemplateReco::PixelCPETemplateReco().

lorentzAngle_

const SiPixelLorentzAngle* PixelCPEBase::lorentzAngle_

protected

Definition at line 232 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

lorentzAngleWidth_

const SiPixelLorentzAngle* PixelCPEBase::lorentzAngleWidth_

protected

Definition at line 233 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

m_DetParams

DetParams PixelCPEBase::m_DetParams = DetParams(1440)

protected

Definition at line 281 of file PixelCPEBase.h.

Referenced by detParam(), PixelCPEClusterRepair::fill2DTemplIDs(), fillDetParams(), and PixelCPEFast::fillParamsForGpu().

magfield_

const MagneticField* PixelCPEBase::magfield_

protected

Definition at line 228 of file PixelCPEBase.h.

Referenced by fillDetParams().

templateDBobject_

const SiPixelTemplateDBObject* PixelCPEBase::templateDBobject_

protected

Definition at line 236 of file PixelCPEBase.h.

Referenced by fillDetParams(), PixelCPETemplateReco::localPosition(), PixelCPEBase(), PixelCPEClusterRepair::PixelCPEClusterRepair(), and PixelCPETemplateReco::PixelCPETemplateReco().

theFlag_

int PixelCPEBase::theFlag_

protected

Definition at line 226 of file PixelCPEBase.h.

Referenced by driftDirection(), fillDetParams(), and PixelCPEBase().

theVerboseLevel

int PixelCPEBase::theVerboseLevel

protected

Definition at line 225 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), PixelCPETemplateReco::localError(), PixelCPEGenericForBricked::localPosition(), PixelCPEBase(), PixelCPEGeneric::PixelCPEGeneric(), and PixelCPEGenericForBricked::PixelCPEGenericForBricked().

ttopo_

const TrackerTopology& PixelCPEBase::ttopo_

protected

Definition at line 230 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::checkRecommend2D(), PixelCPEFast::fillParamsForGpu(), and PixelCPEGenericBase::setXYErrors().

useLAFromDB_

bool PixelCPEBase::useLAFromDB_

protected

Definition at line 239 of file PixelCPEBase.h.

Referenced by fillDetParams(), and PixelCPEBase().

useLAOffsetFromConfig_

bool PixelCPEBase::useLAOffsetFromConfig_

protected

Definition at line 220 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

useLAWidthFromConfig_

bool PixelCPEBase::useLAWidthFromConfig_

protected

Definition at line 221 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

useLAWidthFromDB_

bool PixelCPEBase::useLAWidthFromDB_

protected

Definition at line 222 of file PixelCPEBase.h.

Referenced by driftDirection(), and PixelCPEBase().

xEdgeXError_

constexpr float PixelCPEBase::xEdgeXError_ = 23.0f

staticprotected

Definition at line 246 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

xEdgeYError_

constexpr float PixelCPEBase::xEdgeYError_ = 39.0f

staticprotected

Definition at line 247 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

yEdgeXError_

constexpr float PixelCPEBase::yEdgeXError_ = 24.0f

staticprotected

Definition at line 249 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().

yEdgeYError_

constexpr float PixelCPEBase::yEdgeYError_ = 96.0f

staticprotected

Definition at line 250 of file PixelCPEBase.h.

Referenced by PixelCPEClusterRepair::localError(), and PixelCPETemplateReco::localError().