ME0SegFit Class Reference

#include <ME0SegFit.h>

Public Types
typedef std::vector< const ME0RecHit * >	ME0SetOfHits
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
double	chi2 (void) const
AlgebraicSymMatrix	covarianceMatrix (void)
void	fit (void)
bool	fitdone () const
ME0SetOfHits	hits (void) const
LocalPoint	intercept () const
LocalVector	localdir () const
const ME0EtaPartition *	me0etapartition (uint32_t id) const
	ME0SegFit (std::map< uint32_t, const ME0EtaPartition * > me0etapartmap, ME0SetOfHits hits)
int	ndof (void) const
size_t	nhits (void) const
float	Rdev (float x, float y, float z) const
const ME0EtaPartition *	refme0etapart () 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	~ME0SegFit ()

Protected Member Functions
SMatrix12by4	derivativeMatrix (void)
AlgebraicSymMatrix	flipErrors (const SMatrixSym4 &)
void	setOutFromIP (void)
SMatrixSym12	weightMatrix (void)

Protected Attributes
double	chi2_
bool	fitdone_
ME0SetOfHits	hits_
LocalPoint	intercept_
LocalVector	localdir_
std::map< uint32_t, const ME0EtaPartition * >	me0etapartmap_
int	ndof_
uint32_t	refid_
double	scaleXError_
float	uslope_
float	vslope_

Private Member Functions
void	fit2 (void)
void	fitlsq (void)
void	setChi2 (void)

Detailed Description

Definition at line 37 of file ME0SegFit.h.

Member Typedef Documentation

typedef std::vector<const ME0RecHit*> ME0SegFit::ME0SetOfHits

Definition at line 43 of file ME0SegFit.h.

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

Definition at line 49 of file ME0SegFit.h.

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

Definition at line 52 of file ME0SegFit.h.

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

Definition at line 46 of file ME0SegFit.h.

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

Definition at line 56 of file ME0SegFit.h.

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

Definition at line 53 of file ME0SegFit.h.

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

Definition at line 59 of file ME0SegFit.h.

Constructor & Destructor Documentation

ME0SegFit::ME0SegFit ( std::map< uint32_t, const ME0EtaPartition * >  me0etapartmap, ME0SetOfHits  hits  )

inline

Definition at line 65 of file ME0SegFit.h.

References me0etapartmap_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        : 
  me0etapartmap_( me0etapartmap ), hits_( hits ), scaleXError_( 1.0 ), refid_(me0etapartmap_.begin()->first ), fitdone_( false ) 
    {
      // --- LogDebug info about reading of ME0 Eta Partition map ------------------------------------------
      edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::ctor] cached the me0etapartmap";

      // --- LogDebug for ME0 Eta Partition map ------------------------------------------------------------
      std::stringstream gemetapartmapss; gemetapartmapss<<"[ME0SegFit::ctor] :: me0etapartmap :: elements ["<<std::endl;
      for(std::map<uint32_t, const ME0EtaPartition*>::const_iterator mapIt = me0etapartmap_.begin(); mapIt != me0etapartmap_.end(); ++mapIt)
    {
      gemetapartmapss<<"[ME0 DetId "<<mapIt->first<<" ="<<ME0DetId(mapIt->first)<<", ME0 EtaPart "<<mapIt->second<<"],"<<std::endl;
    }
      gemetapartmapss<<"]"<<std::endl;
      std::string gemetapartmapstr = gemetapartmapss.str();
      edm::LogVerbatim("ME0SegFit") << gemetapartmapstr;
      // --- End LogDebug -----------------------------------------------------------------------------------
    }

virtual ME0SegFit::~ME0SegFit ( )

inlinevirtual

Definition at line 83 of file ME0SegFit.h.

{}

Member Function Documentation

double ME0SegFit::chi2 ( void  ) const

inline

Definition at line 107 of file ME0SegFit.h.

References chi2_.

{ return chi2_; }

AlgebraicSymMatrix ME0SegFit::covarianceMatrix

(

void

)

Definition at line 418 of file ME0SegFit.cc.

References HLT_25ns14e33_v1_cff::A, derivativeMatrix(), flipErrors(), convertSQLiteXML::ok, query::result, weightMatrix(), and create_public_pileup_plots::weights.

                                               {
  
  SMatrixSym12 weights = weightMatrix();
  SMatrix12by4 A = derivativeMatrix();
  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::covarianceMatrix] weights matrix W: \n" << weights;      
  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::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);
  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::covarianceMatrix] (AT W A): \n" << result;      
  ok = result.Invert(); // inverts in place
  if ( !ok ) {
    edm::LogVerbatim("ME0Segment|ME0SegFit") << "[ME0SegFit::calculateError] Failed to invert matrix: \n" << result;      
    //    return ok;  //@@ SHOULD PASS THIS BACK TO CALLER?
  }
  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::covarianceMatrix] (AT W A)^-1: \n" << result;      
  
  // reorder components to match TrackingRecHit interface (ME0Segment isa TrackingRecHit)
  // i.e. slopes first, then positions 
  AlgebraicSymMatrix flipped = flipErrors( result );
    
  return flipped;
}

ME0SegFit::SMatrix12by4 ME0SegFit::derivativeMatrix ( void  )

protected

Definition at line 370 of file ME0SegFit.cc.

References ztail::d, hits_, ME0RecHit::localPosition(), makeMuonMisalignmentScenario::matrix, me0etapartition(), TrackingRecHit::rawId(), refme0etapart(), GeomDet::toLocal(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by covarianceMatrix().

                                                  {
  
  SMatrix12by4 matrix; // 12x4, init to 0
  int row = 0;
  
  for( ME0SetOfHits::const_iterator it = hits_.begin(); it != hits_.end(); ++it) {
    
    const ME0RecHit& hit = (**it);
    ME0DetId d = ME0DetId(hit.rawId());
    const ME0EtaPartition* roll = me0etapartition(d);
    GlobalPoint gp = roll->toGlobal(hit.localPosition());
    LocalPoint lp = refme0etapart()->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;
}

void ME0SegFit::fit

(

void

)

Definition at line 15 of file ME0SegFit.cc.

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

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

void ME0SegFit::fit2 ( void  )

private

Definition at line 36 of file ME0SegFit.cc.

References chi2_, fitdone_, hits_, intercept_, ME0DetId::layer(), ME0RecHit::localPosition(), me0etapartition(), ndof_, refme0etapart(), 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().

                         {

  // 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)


  // 1) Check whether hits are on the same layer
  // -------------------------------------------
  ME0SetOfHits::const_iterator ih = hits_.begin();

  ME0DetId d1 = DetId((*ih)->rawId());
  int il1     = d1.layer();
  const ME0RecHit& h1 = (**ih);
  ++ih;    
  ME0DetId d2 = DetId((*ih)->rawId());
  int il2     = d2.layer();
  const ME0RecHit& h2 = (**ih);
    
  // Skip if on same layer, but should be impossible :)
  if (il1 == il2) {
    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit:fit2] - 2 hits on same layer!!";
    return;
  }
    

  // 2) Global Positions of hit 1 and 2 and
  //    Local  Positions of hit 1 and 2 w.r.t. reference ME0 Eta Partition 
  // ---------------------------------------------------------------------
  const ME0EtaPartition* roll1 = me0etapartition(d1);
  GlobalPoint h1glopos = roll1->toGlobal(h1.localPosition());
  const ME0EtaPartition* roll2 = me0etapartition(d2);
  GlobalPoint h2glopos = roll2->toGlobal(h2.localPosition());
    
  // We want hit wrt first me0 eta partition 
  // ( = reference m00 eta partition) 
  // (and local z will be != 0)
  LocalPoint h1pos = refme0etapart()->toLocal(h1glopos);  
  LocalPoint h2pos = refme0etapart()->toLocal(h2glopos);  
    

  // 3) Now make straight line between the two points in local coords
  // ----------------------------------------------------------------
  float dz = h2pos.z()-h1pos.z();
  if(dz != 0.0) {
    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;
}

bool ME0SegFit::fitdone ( ) const

inline

Definition at line 113 of file ME0SegFit.h.

References fitdone_.

Referenced by fit().

{ return fitdone_; }

void ME0SegFit::fitlsq ( void  )

private

Definition at line 101 of file ME0SegFit.cc.

References ztail::d, fitdone_, hits_, ME0EtaPartition::id(), intercept_, ME0RecHit::localPosition(), ME0RecHit::localPositionError(), LogTrace, me0etapartition(), convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), TrackingRecHit::rawId(), refme0etapart(), setChi2(), setOutFromIP(), AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::toLocal(), uslope_, findQualityFiles::v, vslope_, LocalError::xx(), LocalError::xy(), LocalError::yy(), and z.

Referenced by fit().

                           {
  
  // Linear least-squares fit to up to 6 ME0 rechits, one per layer in a ME0 chamber.
  // Comments adapted from Tim Cox' comments in the original  ME0SegAlgoSK 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; 

  ME0SetOfHits::const_iterator ih = hits_.begin();
  
  // LogDebug :: Loop over the TrackingRecHits and print the ME0 Hits  
  for (ih = hits_.begin(); ih != hits_.end(); ++ih) 
    {
      edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::fitlsq] - looping over ME0RecHits";
      const ME0RecHit& hit = (**ih);
      ME0DetId d = ME0DetId(hit.rawId());
      edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::fitlsq] - Tracking RecHit in detid ("<<d.rawId()<<")";
      edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::fitlsq] - ME0DetId ("<<ME0DetId(d.rawId())<<")";
    }

  // Loop over the ME0RecHits and make small (2x2) matrices used to fill the blockdiagonal covariance matrix E^-1
  for (ih = hits_.begin(); ih != hits_.end(); ++ih) {
    const ME0RecHit& hit = (**ih);
    ME0DetId d = DetId(hit.rawId());
    const ME0EtaPartition* roll = me0etapartition(d);
    GlobalPoint gp         = roll->toGlobal(hit.localPosition());
    LocalPoint  lp         = refme0etapart()->toLocal(gp); 
    
    // LogDebug
    #ifdef EDM_ML_DEBUG // have lines below only compiled when in debug mode
    std::stringstream lpss; lpss<<lp; std::string lps = lpss.str();
    std::stringstream gpss; gpss<<gp; std::string gps = gpss.str();
    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::fitlsq] - Tracking RecHit global position "<<std::setw(30)<<gps<<" and local position "<<std::setw(30)<<lps
                     <<" wrt reference ME0 eta partition "<<refme0etapart()->id().rawId()<<" = "<<refme0etapart()->id();
    #endif

    // 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);
    
    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::fit] 2x2 covariance matrix for this ME0RecHit :: [[" << IC(0,0) <<", "<< IC(0,1) <<"]["<< IC(1,0) <<","<<IC(1,1)<<"]]";

    // Invert covariance matrix (and trap if it fails!)
    bool ok = IC.Invert();
    if ( !ok ) {
      edm::LogVerbatim("ME0Segment|ME0SegFit") << "[ME0SegFit::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("ME0Segment|ME0SegFit") << "[ME0SegFit::fit] Failed to invert matrix: \n" << M;
    //    return ok; //@@ SHOULD PASS THIS BACK TO CALLER?
  }
  else {
    p = M * B;
  }

  LogTrace("ME0SegFit") << "[ME0SegFit::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;

}

AlgebraicSymMatrix ME0SegFit::flipErrors ( const SMatrixSym4 &  a )

protected

Definition at line 446 of file ME0SegFit.cc.

References a, i, and j.

Referenced by covarianceMatrix().

                                                               { 
    
  // The ME0Segment needs the error matrix re-arranged to match
  // parameters in order (uz, vz, u0, v0) 
  // where uz, vz = slopes, u0, v0 = intercepts
    
  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::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
    }
  }

  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::flipErrors] after copy:";
  edm::LogVerbatim("ME0SegFit") << "(" << hold(1,1) << "  " << hold(1,2) << "  " << hold(1,3) << "  " << hold(1,4);
  edm::LogVerbatim("ME0SegFit") << " " << hold(2,1) << "  " << hold(2,2) << "  " << hold(2,3) << "  " << hold(2,4);
  edm::LogVerbatim("ME0SegFit") << " " << hold(3,1) << "  " << hold(3,2) << "  " << hold(3,3) << "  " << hold(3,4);
  edm::LogVerbatim("ME0SegFit") << " " << 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)

  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::flipErrors] after flip:";
  edm::LogVerbatim("ME0SegFit") << "(" << hold(1,1) << "  " << hold(1,2) << "  " << hold(1,3) << "  " << hold(1,4);
  edm::LogVerbatim("ME0SegFit") << " " << hold(2,1) << "  " << hold(2,2) << "  " << hold(2,3) << "  " << hold(2,4);
  edm::LogVerbatim("ME0SegFit") << " " << hold(3,1) << "  " << hold(3,2) << "  " << hold(3,3) << "  " << hold(3,4);
  edm::LogVerbatim("ME0SegFit") << " " << hold(4,1) << "  " << hold(4,2) << "  " << hold(4,3) << "  " << hold(4,4) << ")";

  return hold;
}

ME0SetOfHits ME0SegFit::hits ( void  ) const

inline

Definition at line 104 of file ME0SegFit.h.

References hits_.

{ return hits_; }

LocalPoint ME0SegFit::intercept ( ) const

inline

Definition at line 109 of file ME0SegFit.h.

References intercept_.

{ return intercept_;}

LocalVector ME0SegFit::localdir ( ) const

inline

Definition at line 110 of file ME0SegFit.h.

References localdir_.

{ return localdir_;}

const ME0EtaPartition* ME0SegFit::me0etapartition ( uint32_t  id ) const

inline

Definition at line 111 of file ME0SegFit.h.

References me0etapartmap_.

Referenced by derivativeMatrix(), fit2(), fitlsq(), and setChi2().

{ return me0etapartmap_.find(id)->second; }

int ME0SegFit::ndof ( void  ) const

inline

Definition at line 108 of file ME0SegFit.h.

References ndof_.

{ return ndof_; }

size_t ME0SegFit::nhits ( void  ) const

inline

Definition at line 106 of file ME0SegFit.h.

References hits_.

Referenced by fit().

{ return hits_.size(); }

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

Definition at line 513 of file ME0SegFit.cc.

References mathSSE::sqrt(), xdev(), and ydev().

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

const ME0EtaPartition* ME0SegFit::refme0etapart ( ) const

inline

Definition at line 112 of file ME0SegFit.h.

References me0etapartmap_, and refid_.

Referenced by derivativeMatrix(), fit2(), fitlsq(), setChi2(), and setOutFromIP().

{ return me0etapartmap_.find(refid_)->second; }

double ME0SegFit::scaleXError ( void  ) const

inline

Definition at line 105 of file ME0SegFit.h.

References scaleXError_.

Referenced by weightMatrix().

{ return scaleXError_; }

void ME0SegFit::setChi2 ( void  )

private

Definition at line 284 of file ME0SegFit.cc.

References chi2_, ztail::d, hits_, intercept_, ME0RecHit::localPosition(), ME0RecHit::localPositionError(), me0etapartition(), ndof_, convertSQLiteXML::ok, TrackingRecHit::rawId(), refme0etapart(), GeomDet::toLocal(), uslope_, findQualityFiles::v, vslope_, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by fitlsq().

                            {
  
  double chsq = 0.;

  ME0SetOfHits::const_iterator ih;
  for (ih = hits_.begin(); ih != hits_.end(); ++ih) {

    const ME0RecHit& hit = (**ih);
    ME0DetId d = ME0DetId(hit.rawId());
    const ME0EtaPartition* roll = me0etapartition(d);
    GlobalPoint gp         = roll->toGlobal(hit.localPosition());
    LocalPoint lp          = refme0etapart()->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;
    
    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::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);

    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::setChi2] IC before = \n" << IC;

    // Invert covariance matrix
    bool ok = IC.Invert();
    if (!ok ){
      edm::LogVerbatim("ME0Segment|ME0SegFit") << "[ME0SegFit::setChi2] Failed to invert covariance matrix: \n" << IC;
      //      return ok;
    }
    edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::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;

  edm::LogVerbatim("ME0SegFit") << "[ME0SegFit::setChi2] chi2 = " << chi2_ << "/" << ndof_ << " dof";

}

void ME0SegFit::setOutFromIP ( void  )

protected

Definition at line 395 of file ME0SegFit.cc.

References intercept_, localdir_, refme0etapart(), mathSSE::sqrt(), GeomDet::toGlobal(), csvLumiCalc::unit, uslope_, vslope_, and z.

Referenced by fit2(), and fitlsq().

                             {
  // 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    = ( refme0etapart()->toGlobal( intercept_ ) ).z();
  double globalZdir    = ( refme0etapart()->toGlobal( localDir  ) ).z();
  double directionSign = globalZpos * globalZdir;
  localdir_ = (directionSign * localDir ).unit();
}

void ME0SegFit::setScaleXError ( double  factor )

inline

Definition at line 92 of file ME0SegFit.h.

References scaleXError_.

{ scaleXError_ = factor; }

ME0SegFit::SMatrixSym12 ME0SegFit::weightMatrix ( void  )

protected

Definition at line 337 of file ME0SegFit.cc.

References hits_, ME0RecHit::localPositionError(), makeMuonMisalignmentScenario::matrix, convertSQLiteXML::ok, scaleXError(), LocalError::xx(), LocalError::xy(), and LocalError::yy().

Referenced by covarianceMatrix().

                                              {
  
  bool ok = true;

  SMatrixSym12 matrix = ROOT::Math::SMatrixIdentity(); // 12x12, init to 1's on diag

  int row = 0;
  
  for (ME0SetOfHits::const_iterator it = hits_.begin(); it != hits_.end(); ++it) {
    
    const ME0RecHit& 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("ME0SegFit") << "[ME0SegFit::weightMatrix] Failed to invert matrix: \n" << matrix;      
    //    return ok; //@@ SHOULD PASS THIS BACK TO CALLER?
  }
  return matrix;
}

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

Definition at line 505 of file ME0SegFit.cc.

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

Referenced by Rdev().

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

float ME0SegFit::xfit

(

float

z

)

const

Definition at line 495 of file ME0SegFit.cc.

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

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

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

Definition at line 509 of file ME0SegFit.cc.

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

Referenced by Rdev().

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

float ME0SegFit::yfit

(

float

z

)

const

Definition at line 501 of file ME0SegFit.cc.

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

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

Member Data Documentation

double ME0SegFit::chi2_

protected

Definition at line 145 of file ME0SegFit.h.

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

bool ME0SegFit::fitdone_

protected

Definition at line 149 of file ME0SegFit.h.

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

ME0SetOfHits ME0SegFit::hits_

protected

Definition at line 140 of file ME0SegFit.h.

Referenced by derivativeMatrix(), fit2(), fitlsq(), hits(), nhits(), setChi2(), and weightMatrix().

LocalPoint ME0SegFit::intercept_

protected

Definition at line 143 of file ME0SegFit.h.

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

LocalVector ME0SegFit::localdir_

protected

Definition at line 144 of file ME0SegFit.h.

Referenced by localdir(), and setOutFromIP().

std::map<uint32_t, const ME0EtaPartition*> ME0SegFit::me0etapartmap_

protected

Definition at line 138 of file ME0SegFit.h.

Referenced by me0etapartition(), ME0SegFit(), and refme0etapart().

int ME0SegFit::ndof_

protected

Definition at line 146 of file ME0SegFit.h.

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

uint32_t ME0SegFit::refid_

protected

Definition at line 148 of file ME0SegFit.h.

Referenced by refme0etapart().

double ME0SegFit::scaleXError_

protected

Definition at line 147 of file ME0SegFit.h.

Referenced by scaleXError(), and setScaleXError().

float ME0SegFit::uslope_

protected

Definition at line 141 of file ME0SegFit.h.

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

float ME0SegFit::vslope_

protected

Definition at line 142 of file ME0SegFit.h.

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