MucRecLineFit Class Reference

#include <MucRecLineFit.h>

Public Member Functions
	MucRecLineFit ()
	Constructor.
	~MucRecLineFit ()
	Destructor.
int	LineFit (float x[], float y[], float w[], int n, float *a, float *b, float *chisq, float *siga, float *sigb)
int	LineFit (float x[], float y[], float w[], int part, int seg, int orient, int n, float *a, float *b, float *chisq, float *siga, float *sigb)
	MucRecLineFit ()
	Constructor.
	~MucRecLineFit ()
	Destructor.
int	LineFit (float x[], float y[], float w[], int n, float *a, float *b, float *chisq, float *siga, float *sigb)
int	LineFit (float x[], float y[], float w[], int part, int seg, int orient, int n, float *a, float *b, float *chisq, float *siga, float *sigb)
	MucRecLineFit ()
	Constructor.
	~MucRecLineFit ()
	Destructor.
int	LineFit (float x[], float y[], float w[], int n, float *a, float *b, float *chisq, float *siga, float *sigb)
int	LineFit (float x[], float y[], float w[], int part, int seg, int orient, int n, float *a, float *b, float *chisq, float *siga, float *sigb)

Detailed Description

Definition at line 15 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/MucRecLineFit.h.

Constructor & Destructor Documentation

MucRecLineFit() [1/3]

MucRecLineFit::MucRecLineFit

(

)

Constructor.

Definition at line 19 of file MucRecLineFit.cxx.

                             {
  // Constructor.
}

~MucRecLineFit() [1/3]

MucRecLineFit::~MucRecLineFit

(

)

Destructor.

Definition at line 23 of file MucRecLineFit.cxx.

                              {
  // Destructor.
}

MucRecLineFit() [2/3]

MucRecLineFit::MucRecLineFit

(

)

Constructor.

~MucRecLineFit() [2/3]

MucRecLineFit::~MucRecLineFit

(

)

Destructor.

MucRecLineFit() [3/3]

MucRecLineFit::MucRecLineFit

(

)

Constructor.

~MucRecLineFit() [3/3]

MucRecLineFit::~MucRecLineFit

(

)

Destructor.

Member Function Documentation

LineFit() [1/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

Definition at line 167 of file MucRecLineFit.cxx.

{
  double sum, sx, sy, sxx, sxy, syy, det;
  float  chi;
 
  /* The variable "i" is declared 8 bytes long to avoid triggering an
     apparent alignment bug in optimized code produced by gcc-2.95.3.
     The bug doesn't seem to be present in gcc-3.2. */
  long i;
 
  chi = 99999999.0;
 
  /* N must be >= 2 for this guy to work */
  if ( n < 2 )
  {
    cout << "utiLineFit-W: Too few points for line fit \n" << endl;
    return -1;
  }
 
  /* Initialization */
  sum = 0.0;
  sx  = 0.0;
  sy  = 0.0;
  sxx = 0.0;
  sxy = 0.0;
  syy = 0.0;
 
  /* Find sum , sumx ,sumy, sumxx, sumxy */
  for ( i = 0; i < n; i++ )
  {
    // cout<<"x: "<<x[i]<<"  y: "<<y[i]<<"  w: "<<w[i]<<endl;
    sum = sum + w[i];
    sx  = sx + w[i] * x[i];
    sy  = sy + w[i] * y[i];
    sxx = sxx + w[i] * x[i] * x[i];
    sxy = sxy + w[i] * x[i] * y[i];
    syy = syy + w[i] * y[i] * y[i];
  }
 
  det = sum * sxx - sx * sx;
  if ( fabs( det ) < 1.0e-20 )
  {
    *a     = 1.0e20;
    *b     = x[0];
    *chisq = 0.0;
    *siga  = 0.0;
    *sigb  = 0.0;
 
    return 0;
  }
 
  /* compute the best fitted parameters A,B */
  *a = ( sum * sxy - sx * sy ) / det;
  *b = ( sy * sxx - sxy * sx ) / det;
 
  // cout<<"a: "<<*a<<"  b: "<<*b<<endl;
  /* calculate chi-square */
  chi = 0.0;
  for ( i = 0; i < n; i++ )
  {
    chi =
        chi + ( w[i] ) * ( ( y[i] ) - *a * ( x[i] ) - *b ) * ( ( y[i] ) - *a * ( x[i] ) - *b );
  }
 
  *siga = sqrt( fabs( sum / det ) );
  *sigb = sqrt( fabs( sxx / det ) );
 
  *chisq = chi;
  return 0;
}

Referenced by MucRec2DRoad::Fit(), LineFit(), MucRec2DRoad::SimpleFit(), and MucRec2DRoad::SimpleFitNoCurrentGap().

LineFit() [2/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

LineFit() [3/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

LineFit() [4/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	part,
		int	seg,
		int	orient,
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

Definition at line 63 of file MucRecLineFit.cxx.

                                          {
  if ( part > 2 || part < 0 )
  {
    cout << "BAD part id from MucRecLineFit" << endl;
    return -1;
  }
 
  int   status = -1;
  float c, d, sigc, sigd;
 
  if ( part != 1 )
  {
    for ( int i = 0; i < n; i++ )
    {
      float gap0z = ( MucGeoGeneral::Instance()->GetGap( part, seg, 0 )->GetCenter() ).z();
      // cout<<"in MucRecLineFit change x from "<<x[i]<<" to ";
      x[i] -= gap0z;
      // cout<<x[i]<<endl;
    }
    status = LineFit( x, y, w, n, a, b, chisq, siga, sigb );
    // cout<<"in MucRecLineFit sigb = "<<*sigb<<endl;
  }
  else
  {
    if ( orient == 0 )
    { // need change x and y because x has error now.
 
      for ( int i = 0; i < n; i++ )
      {
        float gap0r = ( MucGeoGeneral::Instance()->GetGap( part, 0, 0 )->GetCenter() ).x();
        y[i] -= gap0r;
      }
      status = LineFit( y, x, w, n, &c, &d, chisq, &sigc, &sigd );
      // need to recalculate a,b,siga,sigb now
      //  y=ax+b;  x=cy+d
      //  a = 1/c; b = -d/c
      //  siga^2 = a^2/c^2 * sigc^2
      //  sigb^2 = b^2/c^2 * sigc^2 + b^2/d^2 * sigd^2;
 
      *a    = 1 / c;
      *b    = -1 * d / c;
      *siga = 1 / c / c * sigc;
      //*sigb = sqrt(fabs(1/c/c*sigd*sigd + d*d/c/c/c/c*sigc*sigc));
      *sigb = sigd; ////////////
    }
    else
    { // barrel... seg0,4 need not change; seg2,6 change x<->y; other seg need more complicate
      // operation!
      if ( seg == 0 || seg == 4 )
      {
        float gap0x = ( MucGeoGeneral::Instance()->GetGap( part, seg, 0 )->GetCenter() ).x();
        for ( int i = 0; i < n; i++ ) { x[i] -= gap0x; }
        status = LineFit( x, y, w, n, a, b, chisq, siga, sigb );
      }
      else if ( seg == 2 || seg == 6 )
      {
        float gap0y = ( MucGeoGeneral::Instance()->GetGap( part, seg, 0 )->GetCenter() ).y();
        for ( int i = 0; i < n; i++ ) { y[i] -= gap0y; }
        status = LineFit( y, x, w, n, &c, &d, chisq, &sigc, &sigd );
        // need to recalculate a,b,siga,sigb now
        *a    = 1 / c;
        *b    = -1 * d / c;
        *siga = 1 / c / c * sigc;
        //*sigb = sqrt(fabs(1/c/c*sigd*sigd + d*d/c/c/c/c*sigc*sigc));
        *sigb = sigd;
      }
      else
      {
        for ( int i = 0; i < n; i++ )
        {
          float temp = ( y[i] + x[i] ) / sqrt( 2.0 );
          y[i]       = ( y[i] - x[i] ) / sqrt( 2.0 );
          x[i]       = temp;
        }
 
        if ( seg == 1 || seg == 5 )
        {
          for ( int i = 0; i < n; i++ )
          {
            float gap0x =
                ( MucGeoGeneral::Instance()->GetGap( part, seg - 1, 0 )->GetCenter() ).x();
            x[i] -= gap0x;
          }
          status = LineFit( x, y, w, n, a, b, chisq, siga, sigb );
        }
        else if ( seg == 3 || seg == 7 )
        {
          for ( int i = 0; i < n; i++ )
          {
            float gap0y =
                ( MucGeoGeneral::Instance()->GetGap( part, seg - 1, 0 )->GetCenter() ).y();
            y[i] -= gap0y;
          }
          status = LineFit( y, x, w, n, a, b, chisq, siga, sigb );
        }
      }
    }
  }
 
  return status;
}

LineFit() [5/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	part,
		int	seg,
		int	orient,
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

LineFit() [6/6]

int MucRecLineFit::LineFit	(	float	x[],
		float	y[],
		float	w[],
		int	part,
		int	seg,
		int	orient,
		int	n,
		float *	a,
		float *	b,
		float *	chisq,
		float *	siga,
		float *	sigb )

---

The documentation for this class was generated from the following files:

• InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/MucRecLineFit.h
• InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/MucRecLineFit.h
• Muc/MucRecEvent/include/MucRecEvent/MucRecLineFit.h
• MucRecLineFit.cxx