MdcLine Class Reference

#include <MdcLine.h>

List of all members.

Public Member Functions
	MdcLine (int n)
	~MdcLine ()
int	fit (int nUse)
Public Attributes
double *	x
double *	y
double *	sigma
double *	resid
int	nPoint
double	slope
double	intercept
double	chisq
double	errmat [3]

---

Detailed Description

Definition at line 13 of file MdcLine.h.

---

Constructor & Destructor Documentation

MdcLine::MdcLine

(

int

n

)

[inline]

Definition at line 25 of file MdcLine.h.

References nPoint, resid, sigma, x, and y.

                 { 
    nPoint = n;
    x = new double[n];
    y = new double[n];
    sigma = new double[n];
    resid = new double[n];
  };

MdcLine::~MdcLine

(

)

[inline]

Definition at line 32 of file MdcLine.h.

References resid, sigma, x, and y.

             { 
    delete [] x;
    delete [] y;
    delete [] sigma;
    delete [] resid;
  };

---

Member Function Documentation

int MdcLine::fit

(

int

nUse

)

Definition at line 10 of file MdcLine.cxx.

References chisq, errmat, intercept, nPoint, resid, sigma, slope, weight, x, and y.

Referenced by MdcSegGrouper::calcParByHits(), MdcSegInfoSterO::calcStereo(), and MdcSegFinder::tryPatterns().

                         {

  double s1, sx, sy, sxx, sxy;  // The various sums. 
  double delinv, denom;
  double weight = 0.013;     // weight for hits, calculated from sigmas. 
  int ihit;
  if (nUse > nPoint) nUse = nPoint;
  s1 = sx = sy = sxx = sxy = 0.0;
  chisq = 0.0;
  for (ihit = 0; ihit < nUse; ihit++) {
    
    if (sigma[ihit] < 0.0) continue;
    weight = 1. / (sigma[ihit] * sigma[ihit]);//NEED sigma of MdcHits
    s1 += weight;    
    sx += x[ihit] * weight;
    sy += y[ihit] * weight;
    sxx += x[ihit] * (x[ihit] * weight);
    sxy += y[ihit] * (x[ihit] * weight);

    //std::cout<<__FILE__<<" "<<__LINE__<< " "<<ihit<<" sigma "<<
      //sigma[ihit] <<" weight "
      //<< weight<<" x "
      //<< x[ihit]<<" y "
      //<< y[ihit]<< std::endl;
  }
  
  /* Calculate parameters. */
  denom = s1 * sxx - sx * sx;
  delinv = (denom == 0.0) ? 1.e20 : 1. / denom;
  intercept = (sy * sxx - sx * sxy) * delinv;
  slope = (s1 * sxy - sx * sy) * delinv;
  errmat[0] = sxx * delinv;
  errmat[1] = -sx * delinv;
  errmat[2] = s1 * delinv;
  //std::cout << " After Fit:" << std::endl;//yzhang debug
  /* Calculate residuals. */
  for (ihit = 0; ihit < nUse; ihit++) {
    if (sigma[ihit] < 0.0) continue;
    resid[ihit] = ( y[ihit] - intercept - slope * x[ihit] );
    chisq += resid[ihit] * resid[ihit] / (sigma[ihit] * sigma[ihit]);
    //std::cout << "resid "<<ihit<<" "<<resid[ihit]<<" sigma "<<sigma[ihit]<<" chisq add "<<resid[ihit] * resid[ihit] / (sigma[ihit] * sigma[ihit]) << std::endl;//yzhang debug 
  }

  //std::cout<<__FILE__<<" "<<__LINE__
  //<< " intercept "<<intercept<<" slope "<<slope 
          //<< " chisq "<<chisq<<std::endl;//yzhang debug 
  return 0;
}

---

Member Data Documentation

double MdcLine::chisq

Definition at line 22 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), fit(), MdcSeg::setValues(), and MdcSegFinder::tryPatterns().

double MdcLine::errmat[3]

Definition at line 23 of file MdcLine.h.

Referenced by fit(), and MdcSeg::setValues().

double MdcLine::intercept

Definition at line 21 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), MdcSegInfoSterO::calcStereo(), fit(), MdcSeg::setValues(), and MdcSegFinder::tryPatterns().

int MdcLine::nPoint

Definition at line 19 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), fit(), and MdcLine().

double* MdcLine::resid

Definition at line 18 of file MdcLine.h.

Referenced by fit(), MdcLine(), MdcSegFinder::tryPatterns(), and ~MdcLine().

double* MdcLine::sigma

Definition at line 17 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), MdcSegInfoSterO::calcStereo(), fit(), MdcLine(), MdcSegFinder::tryPatterns(), MdcSegInfoSterO::zPosition(), and ~MdcLine().

double MdcLine::slope

Definition at line 20 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), MdcSegInfoSterO::calcStereo(), fit(), MdcSeg::setValues(), and MdcSegFinder::tryPatterns().

double* MdcLine::x

Definition at line 15 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), MdcSegInfoSterO::calcStereo(), fit(), MdcLine(), MdcSegFinder::tryPatterns(), MdcSegInfoSterO::zPosition(), and ~MdcLine().

double* MdcLine::y

Definition at line 16 of file MdcLine.h.

Referenced by MdcSegGrouper::calcParByHits(), fit(), MdcLine(), MdcSegFinder::tryPatterns(), MdcSegInfoSterO::zPosition(), and ~MdcLine().

---