#include <CSCSegFit.h>

Inheritance diagram for CSCSegFit:

Public Types
typedef std::vector< const CSCRecHit2D * >	CSCSetOfHits

typedef ROOT::Math::SMatrix< double, 12, 4 >	SMatrix12by4

typedef ROOT::Math::SMatrix< double, 4 >	SMatrix4

typedef ROOT::Math::SMatrix< double, 12, 12, ROOT::Math::MatRepSym< double, 12 > >	SMatrixSym12

typedef ROOT::Math::SMatrix< double, 2, 2, ROOT::Math::MatRepSym< double, 2 > >	SMatrixSym2

typedef ROOT::Math::SMatrix< double, 4, 4, ROOT::Math::MatRepSym< double, 4 > >	SMatrixSym4

typedef ROOT::Math::SVector< double, 4 >	SVector4

Public Member Functions
const CSCChamber *	chamber () const

double	chi2 (void) const

AlgebraicSymMatrix	covarianceMatrix (void)

	CSCSegFit (const CSCChamber *csc, CSCSetOfHits hits)

void	fit (void)

bool	fitdone () const

CSCSetOfHits	hits (void) const

LocalPoint	intercept () const

LocalVector	localdir () const

int	ndof (void) const

size_t	nhits (void) const

float	Rdev (float x, float y, float z) const

double	scaleXError (void) const

void	setScaleXError (double factor)

float	xdev (float x, float z) const

float	xfit (float z) const

float	ydev (float y, float z) const

float	yfit (float z) const

virtual	~CSCSegFit ()

Protected Member Functions
SMatrix12by4	derivativeMatrix (void)

AlgebraicSymMatrix	flipErrors (const SMatrixSym4 &)

void	setOutFromIP (void)

SMatrixSym12	weightMatrix (void)

Protected Attributes
const CSCChamber *	chamber_

double	chi2_

bool	fitdone_

CSCSetOfHits	hits_

LocalPoint	intercept_

LocalVector	localdir_

int	ndof_

double	scaleXError_

float	uslope_

float	vslope_

Private Member Functions
void	fit2 (void)

void	fitlsq (void)

void	setChi2 (void)

Detailed Description

Definition at line 30 of file CSCSegFit.h.

Member Typedef Documentation

◆ CSCSetOfHits

typedef std::vector<const CSCRecHit2D*> CSCSegFit::CSCSetOfHits

Definition at line 34 of file CSCSegFit.h.

◆ SMatrix12by4

typedef ROOT::Math::SMatrix<double, 12, 4> CSCSegFit::SMatrix12by4

Definition at line 40 of file CSCSegFit.h.

◆ SMatrix4

typedef ROOT::Math::SMatrix<double, 4> CSCSegFit::SMatrix4

Definition at line 43 of file CSCSegFit.h.

◆ SMatrixSym12

typedef ROOT::Math::SMatrix<double, 12, 12, ROOT::Math::MatRepSym<double, 12> > CSCSegFit::SMatrixSym12

Definition at line 37 of file CSCSegFit.h.

◆ SMatrixSym2

typedef ROOT::Math::SMatrix<double, 2, 2, ROOT::Math::MatRepSym<double, 2> > CSCSegFit::SMatrixSym2

Definition at line 47 of file CSCSegFit.h.

◆ SMatrixSym4

typedef ROOT::Math::SMatrix<double, 4, 4, ROOT::Math::MatRepSym<double, 4> > CSCSegFit::SMatrixSym4

Definition at line 44 of file CSCSegFit.h.

◆ SVector4

typedef ROOT::Math::SVector<double, 4> CSCSegFit::SVector4

Definition at line 50 of file CSCSegFit.h.

Constructor & Destructor Documentation

◆ CSCSegFit()

CSCSegFit::CSCSegFit	(	const CSCChamber *	csc,
		CSCSetOfHits	hits
	)

inline

Definition at line 55 of file CSCSegFit.h.

56 : chamber_(csc), hits_(hits), scaleXError_(1.0), fitdone_(false) {}

◆ ~CSCSegFit()

virtual CSCSegFit::~CSCSegFit ( )

inlinevirtual

Definition at line 58 of file CSCSegFit.h.

58 {}

Member Function Documentation

◆ chamber()

const CSCChamber* CSCSegFit::chamber ( ) const

inline

Definition at line 86 of file CSCSegFit.h.

86 { return chamber_; }

References chamber_.

Referenced by CSCCondSegFit::correctTheCovX(), derivativeMatrix(), CSCCondSegFit::fit(), fit2(), fitlsq(), geometryXMLparser.CSCAlignable::index(), CSCCondSegFit::setChi2(), setChi2(), and setOutFromIP().

◆ chi2()

double CSCSegFit::chi2 ( void ) const

inline

Definition at line 82 of file CSCSegFit.h.

82 { return chi2_; }

References chi2_.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoTC::buildSegments(), CSCSegAlgoSK::buildSegments(), CSCSegAlgoDF::buildSegments(), CSCSegAlgoShowering::compareProtoSegment(), CSCSegAlgoTC::compareProtoSegment(), CSCSegAlgoDF::compareProtoSegment(), CSCSegAlgoSK::compareProtoSegment(), CSCSegAlgoTC::increaseProtoSegment(), CSCSegAlgoSK::increaseProtoSegment(), CSCSegAlgoST::prune_bad_hits(), and CSCSegAlgoShowering::showerSeg().

◆ covarianceMatrix()

AlgebraicSymMatrix CSCSegFit::covarianceMatrix ( void )

Definition at line 352 of file CSCSegFit.cc.

                                                {
   SMatrixSym12 weights = weightMatrix();
   SMatrix12by4 A = derivativeMatrix();
   //  LogTrace("CSCSegFit") << "[CSCSegFit::covarianceMatrix] weights matrix W: \n" << weights;
   //  LogTrace("CSCSegFit") << "[CSCSegFit::covarianceMatrix] derivatives matrix A: \n" << A;
  
   // (AT W A)^-1
   // e.g. See http://www.phys.ufl.edu/~avery/fitting.html, part I
  
   bool ok;
   SMatrixSym4 result = ROOT::Math::SimilarityT(A, weights);
   //  LogTrace("CSCSegFit") << "[CSCSegFit::covarianceMatrix] (AT W A): \n" << result;
   ok = result.Invert();  // inverts in place
   if (!ok) {
     edm::LogVerbatim("CSCSegment|CSCSegFit") << "[CSCSegFit::calculateError] Failed to invert matrix: \n" << result;
     //    return ok;  //@@ SHOULD PASS THIS BACK TO CALLER?
   }
   //  LogTrace("CSCSegFit") << "[CSCSegFit::covarianceMatrix] (AT W A)^-1: \n" << result;
  
   // reorder components to match TrackingRecHit interface (CSCSegment isa TrackingRecHit)
   // i.e. slopes first, then positions
   AlgebraicSymMatrix flipped = flipErrors(result);
  
   return flipped;
 }

References derivativeMatrix(), flipErrors(), convertSQLiteXML::ok, mps_fire::result, weightMatrix(), and HLT_FULL_cff::weights.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoTC::buildSegments(), CSCSegAlgoDF::buildSegments(), CSCSegAlgoSK::buildSegments(), CSCSegAlgoST::prune_bad_hits(), and CSCSegAlgoShowering::showerSeg().

◆ derivativeMatrix()

CSCSegFit::SMatrix12by4 CSCSegFit::derivativeMatrix ( void )

protected

Definition at line 310 of file CSCSegFit.cc.

                                                   {
   SMatrix12by4 matrix;  // 12x4, init to 0
   int row = 0;
  
   for (CSCSetOfHits::const_iterator it = hits_.begin(); it != hits_.end(); ++it) {
     const CSCRecHit2D& hit = (**it);
     const CSCLayer* layer = chamber()->layer(hit.cscDetId().layer());
     GlobalPoint gp = layer->toGlobal(hit.localPosition());
     LocalPoint lp = chamber()->toLocal(gp);
     float z = lp.z();
  
     matrix(row, 0) = 1.;
     matrix(row, 2) = z;
     ++row;
     matrix(row, 1) = 1.;
     matrix(row, 3) = z;
     ++row;
   }
   return matrix;
 }

References chamber(), runTauDisplay::gp, hits_, CSCChamber::layer(), phase1PixelTopology::layer, makeMuonMisalignmentScenario::matrix, GeomDet::toLocal(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by covarianceMatrix().

◆ fit()

void CSCSegFit::fit ( void )

Definition at line 13 of file CSCSegFit.cc.

                         {
   if (fitdone())
     return;  // don't redo fit unnecessarily
   short n = nhits();
   switch (n) {
     case 1:
       edm::LogVerbatim("CSCSegFit") << "[CSCSegFit::fit] - cannot fit just 1 hit!!";
       break;
     case 2:
       fit2();
       break;
     case 3:
     case 4:
     case 5:
     case 6:
       fitlsq();
       break;
     default:
       edm::LogVerbatim("CSCSegFit") << "[CSCSegFit::fit] - cannot fit more than 6 hits!!";
   }
 }

References fit2(), fitdone(), fitlsq(), dqmiodumpmetadata::n, and nhits().

Referenced by CSCSegAlgoShowering::compareProtoSegment(), CSCSegAlgoDF::compareProtoSegment(), trackingPlots.Iteration::modules(), CSCSegAlgoShowering::updateParameters(), CSCSegAlgoDF::updateParameters(), CSCSegAlgoSK::updateParameters(), and CSCSegAlgoTC::updateParameters().

◆ fit2()

void CSCSegFit::fit2 ( void )

private

Definition at line 35 of file CSCSegFit.cc.

                          {
   // Just join the two points
   // Equation of straight line between (x1, y1) and (x2, y2) in xy-plane is
   //       y = mx + c
   // with m = (y2-y1)/(x2-x1)
   // and  c = (y1*x2-x2*y1)/(x2-x1)
  
   CSCSetOfHits::const_iterator ih = hits_.begin();
   int il1 = (*ih)->cscDetId().layer();
   const CSCRecHit2D& h1 = (**ih);
   ++ih;
   int il2 = (*ih)->cscDetId().layer();
   const CSCRecHit2D& h2 = (**ih);
  
   // Skip if on same layer, but should be impossible :)
   if (il1 == il2) {
     edm::LogVerbatim("CSCSegFit") << "[CSCSegFit:fit]2 - 2 hits on same layer!!";
     return;
   }
  
   const CSCLayer* layer1 = chamber()->layer(il1);
   const CSCLayer* layer2 = chamber()->layer(il2);
  
   GlobalPoint h1glopos = layer1->toGlobal(h1.localPosition());
   GlobalPoint h2glopos = layer2->toGlobal(h2.localPosition());
  
   // We want hit wrt chamber (and local z will be != 0)
   LocalPoint h1pos = chamber()->toLocal(h1glopos);
   LocalPoint h2pos = chamber()->toLocal(h2glopos);
  
   float dz = h2pos.z() - h1pos.z();
  
   uslope_ = (h2pos.x() - h1pos.x()) / dz;
   vslope_ = (h2pos.y() - h1pos.y()) / dz;
  
   float uintercept = (h1pos.x() * h2pos.z() - h2pos.x() * h1pos.z()) / dz;
   float vintercept = (h1pos.y() * h2pos.z() - h2pos.y() * h1pos.z()) / dz;
   intercept_ = LocalPoint(uintercept, vintercept, 0.);
  
   setOutFromIP();
  
   //@@ NOT SURE WHAT IS SENSIBLE FOR THESE...
   chi2_ = 0.;
   ndof_ = 0;
  
   fitdone_ = true;
 }

References chamber(), chi2_, CSCRecHit2D::cscDetId(), PVValHelper::dz, fitdone_, hits_, intercept_, CSCChamber::layer(), CSCDetId::layer(), CSCRecHit2D::localPosition(), ndof_, setOutFromIP(), GeomDet::toGlobal(), GeomDet::toLocal(), uslope_, vslope_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by fit().

◆ fitdone()

bool CSCSegFit::fitdone ( ) const

inline

Definition at line 87 of file CSCSegFit.h.

87 { return fitdone_; }

References fitdone_.

Referenced by fit().

◆ fitlsq()

void CSCSegFit::fitlsq ( void )

private

Definition at line 83 of file CSCSegFit.cc.

                            {
   // Linear least-squares fit to up to 6 CSC rechits, one per layer in a CSC.
   // Comments adapted from  mine in original  CSCSegAlgoSK algorithm.
  
   // Fit to the local x, y rechit coordinates in z projection
   // The strip measurement controls the precision of x
   // The wire measurement controls the precision of y.
   // Typical precision: u (strip, sigma~200um), v (wire, sigma~1cm)
  
   // Set up the normal equations for the least-squares fit as a matrix equation
  
   // We have a vector of measurements m, which is a 2n x 1 dim matrix
   // The transpose mT is (u1, v1, u2, v2, ..., un, vn) where
   // ui is the strip-associated measurement and
   // vi is the wire-associated measurement
   // for a given rechit i.
  
   // The fit is to
   // u = u0 + uz * z
   // v = v0 + vz * z
   // where u0, uz, v0, vz are the parameters to be obtained from the fit.
  
   // These are contained in a vector p which is a 4x1 dim matrix, and
   // its transpose pT is (u0, v0, uz, vz). Note the ordering!
  
   // The covariance matrix for each pair of measurements is 2 x 2 and
   // the inverse of this is the error matrix E.
   // The error matrix for the whole set of n measurements is a diagonal
   // matrix with diagonal elements the individual 2 x 2 error matrices
   // (because the inverse of a diagonal matrix is a diagonal matrix
   // with each element the inverse of the original.)
  
   // In function 'weightMatrix()', the variable 'matrix' is filled with this
   // block-diagonal overall covariance matrix. Then 'matrix' is inverted to the
   // block-diagonal error matrix, and returned.
  
   // Define the matrix A as
   //    1   0   z1  0
   //    0   1   0   z1
   //    1   0   z2  0
   //    0   1   0   z2
   //    ..  ..  ..  ..
   //    1   0   zn  0
   //    0   1   0   zn
  
   // This matrix A is set up and returned by function 'derivativeMatrix()'.
  
   // Then the normal equations are described by the matrix equation
   //
   //    (AT E A)p = (AT E)m
   //
   // where AT is the transpose of A.
  
   // Call the combined matrix on the LHS, M, and that on the RHS, B:
   //     M p = B
  
   // We solve this for the parameter vector, p.
   // The elements of M and B then involve sums over the hits
  
   // The covariance matrix of the parameters is obtained by
   // (AT E A)^-1 calculated in 'covarianceMatrix()'.
  
   // NOTE
   // We need local position of a RecHit w.r.t. the CHAMBER
   // and the RecHit itself only knows its local position w.r.t.
   // the LAYER, so we must explicitly transform global position.
  
   SMatrix4 M;  // 4x4, init to 0
   SVector4 B;  // 4x1, init to 0;
  
   CSCSetOfHits::const_iterator ih = hits_.begin();
  
   for (ih = hits_.begin(); ih != hits_.end(); ++ih) {
     const CSCRecHit2D& hit = (**ih);
     const CSCLayer* layer = chamber()->layer(hit.cscDetId().layer());
     GlobalPoint gp = layer->toGlobal(hit.localPosition());
     LocalPoint lp = chamber()->toLocal(gp);
  
     // Local position of hit w.r.t. chamber
     double u = lp.x();
     double v = lp.y();
     double z = lp.z();
  
     // Covariance matrix of local errors
     SMatrixSym2 IC;  // 2x2, init to 0
  
     IC(0, 0) = hit.localPositionError().xx();
     IC(1, 1) = hit.localPositionError().yy();
     //@@ NOT SURE WHICH OFF-DIAGONAL ELEMENT MUST BE DEFINED BUT (1,0) WORKS
     //@@ (and SMatrix enforces symmetry)
     IC(1, 0) = hit.localPositionError().xy();
     // IC(0,1) = IC(1,0);
  
     // Invert covariance matrix (and trap if it fails!)
     bool ok = IC.Invert();
     if (!ok) {
       edm::LogVerbatim("CSCSegment|CSCSegFit") << "[CSCSegFit::fit] Failed to invert covariance matrix: \n" << IC;
       //      return ok;  //@@ SHOULD PASS THIS BACK TO CALLER?
     }
  
     M(0, 0) += IC(0, 0);
     M(0, 1) += IC(0, 1);
     M(0, 2) += IC(0, 0) * z;
     M(0, 3) += IC(0, 1) * z;
     B(0) += u * IC(0, 0) + v * IC(0, 1);
  
     M(1, 0) += IC(1, 0);
     M(1, 1) += IC(1, 1);
     M(1, 2) += IC(1, 0) * z;
     M(1, 3) += IC(1, 1) * z;
     B(1) += u * IC(1, 0) + v * IC(1, 1);
  
     M(2, 0) += IC(0, 0) * z;
     M(2, 1) += IC(0, 1) * z;
     M(2, 2) += IC(0, 0) * z * z;
     M(2, 3) += IC(0, 1) * z * z;
     B(2) += (u * IC(0, 0) + v * IC(0, 1)) * z;
  
     M(3, 0) += IC(1, 0) * z;
     M(3, 1) += IC(1, 1) * z;
     M(3, 2) += IC(1, 0) * z * z;
     M(3, 3) += IC(1, 1) * z * z;
     B(3) += (u * IC(1, 0) + v * IC(1, 1)) * z;
   }
  
   SVector4 p;
   bool ok = M.Invert();
   if (!ok) {
     edm::LogVerbatim("CSCSegment|CSCSegFit") << "[CSCSegFit::fit] Failed to invert matrix: \n" << M;
     //    return ok; //@@ SHOULD PASS THIS BACK TO CALLER?
   } else {
     p = M * B;
   }
  
   //  LogTrace("CSCSegFit") << "[CSCSegFit::fit] p = "
   //        << p(0) << ", " << p(1) << ", " << p(2) << ", " << p(3);
  
   // fill member variables  (note origin has local z = 0)
   //  intercept_
   intercept_ = LocalPoint(p(0), p(1), 0.);
  
   // localdir_ - set so segment points outwards from IP
   uslope_ = p(2);
   vslope_ = p(3);
   setOutFromIP();
  
   // calculate chi2 of fit
   setChi2();
  
   // flag fit has been done
   fitdone_ = true;
 }

References TtFullHadDaughter::B, chamber(), fitdone_, runTauDisplay::gp, hits_, intercept_, CSCChamber::layer(), phase1PixelTopology::layer, convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, setChi2(), setOutFromIP(), GeomDet::toLocal(), uslope_, findQualityFiles::v, vslope_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by fit().

◆ flipErrors()

AlgebraicSymMatrix CSCSegFit::flipErrors ( const SMatrixSym4 & a )

protected

Definition at line 378 of file CSCSegFit.cc.

                                                              {
   // The CSCSegment needs the error matrix re-arranged to match
   // parameters in order (uz, vz, u0, v0)
   // where uz, vz = slopes, u0, v0 = intercepts
  
   //  LogTrace("CSCSegFit") << "[CSCSegFit::flipErrors] input: \n" << a;
  
   AlgebraicSymMatrix hold(4, 0.);
  
   for (short j = 0; j != 4; ++j) {
     for (short i = 0; i != 4; ++i) {
       hold(i + 1, j + 1) = a(i, j);  // SMatrix counts from 0, AlgebraicMatrix from 1
     }
   }
  
   //  LogTrace("CSCSegFit") << "[CSCSegFit::flipErrors] after copy:";
   //  LogTrace("CSCSegFit") << "(" << hold(1,1) << "  " << hold(1,2) << "  " << hold(1,3) << "  " << hold(1,4);
   //  LogTrace("CSCSegFit") << " " << hold(2,1) << "  " << hold(2,2) << "  " << hold(2,3) << "  " << hold(2,4);
   //  LogTrace("CSCSegFit") << " " << hold(3,1) << "  " << hold(3,2) << "  " << hold(3,3) << "  " << hold(3,4);
   //  LogTrace("CSCSegFit") << " " << hold(4,1) << "  " << hold(4,2) << "  " << hold(4,3) << "  " << hold(4,4) << ")";
  
   // errors on slopes into upper left
   hold(1, 1) = a(2, 2);
   hold(1, 2) = a(2, 3);
   hold(2, 1) = a(3, 2);
   hold(2, 2) = a(3, 3);
  
   // errors on positions into lower right
   hold(3, 3) = a(0, 0);
   hold(3, 4) = a(0, 1);
   hold(4, 3) = a(1, 0);
   hold(4, 4) = a(1, 1);
  
   // must also interchange off-diagonal elements of off-diagonal 2x2 submatrices
   hold(4, 1) = a(1, 2);
   hold(3, 2) = a(0, 3);
   hold(2, 3) = a(3, 0);  // = a(0,3)
   hold(1, 4) = a(2, 1);  // = a(1,2)
  
   //  LogTrace("CSCSegFit") << "[CSCSegFit::flipErrors] after flip:";
   //  LogTrace("CSCSegFit") << "(" << hold(1,1) << "  " << hold(1,2) << "  " << hold(1,3) << "  " << hold(1,4);
   //  LogTrace("CSCSegFit") << " " << hold(2,1) << "  " << hold(2,2) << "  " << hold(2,3) << "  " << hold(2,4);
   //  LogTrace("CSCSegFit") << " " << hold(3,1) << "  " << hold(3,2) << "  " << hold(3,3) << "  " << hold(3,4);
   //  LogTrace("CSCSegFit") << " " << hold(4,1) << "  " << hold(4,2) << "  " << hold(4,3) << "  " << hold(4,4) << ")";
  
   return hold;
 }

References a, mps_fire::i, and dqmiolumiharvest::j.

Referenced by covarianceMatrix().

◆ hits()

CSCSetOfHits CSCSegFit::hits ( void ) const

inline

Definition at line 79 of file CSCSegFit.h.

79 { return hits_; }

References hits_.

Referenced by CSCSegAlgoTC::buildSegments(), CSCSegAlgoDF::buildSegments(), CSCSegAlgoSK::buildSegments(), and CSCSegAlgoShowering::showerSeg().

◆ intercept()

LocalPoint CSCSegFit::intercept ( ) const

inline

Definition at line 84 of file CSCSegFit.h.

84 { return intercept_; }

References intercept_.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoTC::buildSegments(), CSCSegAlgoDF::buildSegments(), CSCSegAlgoSK::buildSegments(), CSCSegAlgoSK::phiAtZ(), CSCSegAlgoTC::phiAtZ(), CSCSegAlgoST::prune_bad_hits(), and CSCSegAlgoShowering::showerSeg().

◆ localdir()

LocalVector CSCSegFit::localdir ( ) const

inline

Definition at line 85 of file CSCSegFit.h.

85 { return localdir_; }

References localdir_.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoTC::buildSegments(), CSCSegAlgoDF::buildSegments(), CSCSegAlgoSK::buildSegments(), CSCSegAlgoSK::phiAtZ(), CSCSegAlgoTC::phiAtZ(), CSCSegAlgoST::prune_bad_hits(), and CSCSegAlgoShowering::showerSeg().

◆ ndof()

int CSCSegFit::ndof ( void ) const

inline

Definition at line 83 of file CSCSegFit.h.

83 { return ndof_; }

References ndof_.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoTC::increaseProtoSegment(), CSCSegAlgoSK::increaseProtoSegment(), and CSCSegAlgoST::prune_bad_hits().

◆ nhits()

size_t CSCSegFit::nhits ( void ) const

inline

Definition at line 81 of file CSCSegFit.h.

81 { return hits_.size(); }

References hits_.

Referenced by CSCSegAlgoST::buildSegments(), fit(), and CSCSegAlgoShowering::showerSeg().

◆ Rdev()

float CSCSegFit::Rdev	(	float	x,
		float	y,
		float	z
	)		const

Definition at line 438 of file CSCSegFit.cc.

                                                      {
   return sqrt(xdev(x, z) * xdev(x, z) + ydev(y, z) * ydev(y, z));
 }

References mathSSE::sqrt(), x, xdev(), y, ydev(), and z.

Referenced by CSCSegAlgoShowering::pruneFromResidual(), and CSCSegAlgoDF::pruneFromResidual().

◆ scaleXError()

double CSCSegFit::scaleXError ( void ) const

inline

Definition at line 80 of file CSCSegFit.h.

80 { return scaleXError_; }

References scaleXError_.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoST::prune_bad_hits(), and weightMatrix().

◆ setChi2()

void CSCSegFit::setChi2 ( void )

private

Definition at line 236 of file CSCSegFit.cc.

                             {
   double chsq = 0.;
  
   CSCSetOfHits::const_iterator ih;
   for (ih = hits_.begin(); ih != hits_.end(); ++ih) {
     const CSCRecHit2D& hit = (**ih);
     const CSCLayer* layer = chamber()->layer(hit.cscDetId().layer());
     GlobalPoint gp = layer->toGlobal(hit.localPosition());
     LocalPoint lp = chamber()->toLocal(gp);
  
     double u = lp.x();
     double v = lp.y();
     double z = lp.z();
  
     double du = intercept_.x() + uslope_ * z - u;
     double dv = intercept_.y() + vslope_ * z - v;
  
     //    LogTrace("CSCSegFit") << "[CSCSegFit::setChi2] u, v, z = " << u << ", " << v << ", " << z;
  
     SMatrixSym2 IC;  // 2x2, init to 0
  
     IC(0, 0) = hit.localPositionError().xx();
     //    IC(0,1) = hit.localPositionError().xy();
     IC(1, 0) = hit.localPositionError().xy();
     IC(1, 1) = hit.localPositionError().yy();
     //    IC(1,0) = IC(0,1);
  
     //    LogTrace("CSCSegFit") << "[CSCSegFit::setChi2] IC before = \n" << IC;
  
     // Invert covariance matrix
     bool ok = IC.Invert();
     if (!ok) {
       edm::LogVerbatim("CSCSegment|CSCSegFit") << "[CSCSegFit::setChi2] Failed to invert covariance matrix: \n" << IC;
       //      return ok;
     }
     //    LogTrace("CSCSegFit") << "[CSCSegFit::setChi2] IC after = \n" << IC;
     chsq += du * du * IC(0, 0) + 2. * du * dv * IC(0, 1) + dv * dv * IC(1, 1);
   }
  
   // fill member variables
   chi2_ = chsq;
   ndof_ = 2. * hits_.size() - 4;
  
   //  LogTrace("CSCSegFit") << "[CSCSegFit::setChi2] chi2 = " << chi2_ << "/" << ndof_ << " dof";
 }

References chamber(), chi2_, runTauDisplay::gp, hits_, intercept_, CSCChamber::layer(), phase1PixelTopology::layer, ndof_, convertSQLiteXML::ok, GeomDet::toLocal(), uslope_, findQualityFiles::v, vslope_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by CSCCondSegFit::fit(), and fitlsq().

◆ setOutFromIP()

void CSCSegFit::setOutFromIP ( void )

protected

Definition at line 331 of file CSCSegFit.cc.

                              {
   // Set direction of segment to point from IP outwards
   // (Incorrect for particles not coming from IP, of course.)
  
   double dxdz = uslope_;
   double dydz = vslope_;
   double dz = 1. / sqrt(1. + dxdz * dxdz + dydz * dydz);
   double dx = dz * dxdz;
   double dy = dz * dydz;
   LocalVector localDir(dx, dy, dz);
  
   // localDir sometimes needs a sign flip
   // Examine its direction and origin in global z: to point outward
   // the localDir should always have same sign as global z...
  
   double globalZpos = (chamber()->toGlobal(intercept_)).z();
   double globalZdir = (chamber()->toGlobal(localDir)).z();
   double directionSign = globalZpos * globalZdir;
   localdir_ = (directionSign * localDir).unit();
 }

References chamber(), PVValHelper::dx, BeamSpotPI::dxdz, PVValHelper::dy, BeamSpotPI::dydz, PVValHelper::dz, intercept_, localdir_, mathSSE::sqrt(), GeomDet::toGlobal(), unit(), uslope_, vslope_, and z.

Referenced by CSCCondSegFit::fit(), fit2(), and fitlsq().

◆ setScaleXError()

void CSCSegFit::setScaleXError ( double factor )

inline

Definition at line 67 of file CSCSegFit.h.

67 { scaleXError_ = factor; }

References DQMScaleToClient_cfi::factor, and scaleXError_.

Referenced by CSCSegAlgoST::buildSegments(), CSCCondSegFit::correctTheCovX(), and CSCSegAlgoST::prune_bad_hits().

◆ weightMatrix()

CSCSegFit::SMatrixSym12 CSCSegFit::weightMatrix ( void )

protected

Definition at line 282 of file CSCSegFit.cc.

                                               {
   bool ok = true;
  
   SMatrixSym12 matrix = ROOT::Math::SMatrixIdentity();  // 12x12, init to 1's on diag
  
   int row = 0;
  
   for (CSCSetOfHits::const_iterator it = hits_.begin(); it != hits_.end(); ++it) {
     const CSCRecHit2D& hit = (**it);
  
     // Note scaleXError allows rescaling the x error if necessary
  
     matrix(row, row) = scaleXError() * hit.localPositionError().xx();
     matrix(row, row + 1) = hit.localPositionError().xy();
     ++row;
     matrix(row, row - 1) = hit.localPositionError().xy();
     matrix(row, row) = hit.localPositionError().yy();
     ++row;
   }
  
   ok = matrix.Invert();  // invert in place
   if (!ok) {
     edm::LogVerbatim("CSCSegment|CSCSegFit") << "[CSCSegFit::weightMatrix] Failed to invert matrix: \n" << matrix;
     //    return ok; //@@ SHOULD PASS THIS BACK TO CALLER?
   }
   return matrix;
 }

References hits_, makeMuonMisalignmentScenario::matrix, convertSQLiteXML::ok, and scaleXError().

Referenced by covarianceMatrix().

◆ xdev()

float CSCSegFit::xdev	(	float	x,
		float	z
	)		const

Definition at line 434 of file CSCSegFit.cc.

434 { return intercept_.x() + uslope_ * z - x; }

References intercept_, uslope_, x, PV3DBase< T, PVType, FrameType >::x(), and z.

Referenced by Rdev().

◆ xfit()

float CSCSegFit::xfit ( float z ) const

Definition at line 426 of file CSCSegFit.cc.

                                    {
   //@@ ADD THIS TO EACH ACCESSOR OF FIT RESULTS?
   //  if ( !fitdone() ) fit();
   return intercept_.x() + uslope_ * z;
 }

References intercept_, uslope_, PV3DBase< T, PVType, FrameType >::x(), and z.

Referenced by CSCSegAlgoShowering::isHitNearSegment(), and CSCSegAlgoDF::isHitNearSegment().

◆ ydev()

float CSCSegFit::ydev	(	float	y,
		float	z
	)		const

Definition at line 436 of file CSCSegFit.cc.

436 { return intercept_.y() + vslope_ * z - y; }

References intercept_, vslope_, y, PV3DBase< T, PVType, FrameType >::y(), and z.

Referenced by Rdev().

◆ yfit()

float CSCSegFit::yfit ( float z ) const

Definition at line 432 of file CSCSegFit.cc.

432 { return intercept_.y() + vslope_ * z; }

References intercept_, vslope_, PV3DBase< T, PVType, FrameType >::y(), and z.

Referenced by CSCSegAlgoShowering::isHitNearSegment(), and CSCSegAlgoDF::isHitNearSegment().

Member Data Documentation

◆ chamber_

const CSCChamber* CSCSegFit::chamber_

protected

Definition at line 108 of file CSCSegFit.h.

Referenced by chamber().

◆ chi2_

double CSCSegFit::chi2_

protected

Definition at line 114 of file CSCSegFit.h.

Referenced by chi2(), fit2(), CSCCondSegFit::setChi2(), and setChi2().

◆ fitdone_

bool CSCSegFit::fitdone_

protected

Definition at line 117 of file CSCSegFit.h.

Referenced by fit2(), fitdone(), and fitlsq().

◆ hits_

CSCSetOfHits CSCSegFit::hits_

protected

Definition at line 109 of file CSCSegFit.h.

Referenced by CSCCondSegFit::correctTheCovX(), derivativeMatrix(), CSCCondSegFit::fit(), fit2(), fitlsq(), hits(), nhits(), CSCCondSegFit::setChi2(), setChi2(), and weightMatrix().

◆ intercept_

LocalPoint CSCSegFit::intercept_

protected

Definition at line 112 of file CSCSegFit.h.

Referenced by CSCCondSegFit::fit(), fit2(), fitlsq(), intercept(), CSCCondSegFit::setChi2(), setChi2(), setOutFromIP(), xdev(), xfit(), ydev(), and yfit().

◆ localdir_

LocalVector CSCSegFit::localdir_

protected

Definition at line 113 of file CSCSegFit.h.

Referenced by localdir(), and setOutFromIP().

◆ ndof_

int CSCSegFit::ndof_

protected

Definition at line 115 of file CSCSegFit.h.

Referenced by fit2(), ndof(), CSCCondSegFit::setChi2(), and setChi2().

◆ scaleXError_

double CSCSegFit::scaleXError_

protected

Definition at line 116 of file CSCSegFit.h.

Referenced by scaleXError(), and setScaleXError().

◆ uslope_

float CSCSegFit::uslope_

protected

Definition at line 110 of file CSCSegFit.h.

Referenced by CSCCondSegFit::fit(), fit2(), fitlsq(), CSCCondSegFit::setChi2(), setChi2(), setOutFromIP(), xdev(), and xfit().

◆ vslope_

float CSCSegFit::vslope_

protected

Definition at line 111 of file CSCSegFit.h.

Referenced by CSCCondSegFit::fit(), fit2(), fitlsq(), CSCCondSegFit::setChi2(), setChi2(), setOutFromIP(), ydev(), and yfit().

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

◆ CSCSetOfHits

◆ SMatrix12by4

◆ SMatrix4

◆ SMatrixSym12

◆ SMatrixSym2

◆ SMatrixSym4

◆ SVector4

Constructor & Destructor Documentation

◆ CSCSegFit()

◆ ~CSCSegFit()

Member Function Documentation

◆ chamber()

◆ chi2()

◆ covarianceMatrix()

◆ derivativeMatrix()

◆ fit()

◆ fit2()

◆ fitdone()

◆ fitlsq()

◆ flipErrors()

◆ hits()

◆ intercept()

◆ localdir()

◆ ndof()

◆ nhits()

◆ Rdev()

◆ scaleXError()

◆ setChi2()

◆ setOutFromIP()

◆ setScaleXError()

◆ weightMatrix()

◆ xdev()

◆ xfit()

◆ ydev()

◆ yfit()

Member Data Documentation

◆ chamber_

◆ chi2_

◆ fitdone_

◆ hits_

◆ intercept_

◆ localdir_

◆ ndof_

◆ scaleXError_

◆ uslope_

◆ vslope_