MdcSegInfoCsmc Class Reference

#include <MdcSegInfoCsmc.h>

Inheritance diagram for MdcSegInfoCsmc:

Public Member Functions
	MdcSegInfoCsmc ()
	~MdcSegInfoCsmc ()
double	d0 () const
double	phi0 () const
double	sigPhi0 () const
double	sigD0 () const
bool	parIsAngle (int i) const
void	calcStraight (double phi, double slope, double rad, const double *inErr)
void	calcStraight (const MdcSeg *parentSeg)
	MdcSegInfoCsmc ()
	~MdcSegInfoCsmc ()
double	d0 () const
double	phi0 () const
double	sigPhi0 () const
double	sigD0 () const
bool	parIsAngle (int i) const
void	calcStraight (double phi, double slope, double rad, const double *inErr)
void	calcStraight (const MdcSeg *parentSeg)
Public Member Functions inherited from MdcSegInfo
	MdcSegInfo ()
virtual	~MdcSegInfo ()
const double *	errmat () const
const double *	inverr () const
double	par (int i) const
double	arc () const
virtual bool	parIsAngle (int i) const =0
void	plotSegInfo () const
	MdcSegInfo ()
virtual	~MdcSegInfo ()
const double *	errmat () const
const double *	inverr () const
double	par (int i) const
double	arc () const
virtual bool	parIsAngle (int i) const =0
void	plotSegInfo () const

Additional Inherited Members
Protected Attributes inherited from MdcSegInfo
double	_par0
double	_par1
double	_errmat [3]
double	_inverr [3]
double	_arc

Detailed Description

Definition at line 27 of file InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h.

Constructor & Destructor Documentation

MdcSegInfoCsmc() [1/2]

MdcSegInfoCsmc::MdcSegInfoCsmc ( )

inline

Definition at line 30 of file InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h.

{ };

~MdcSegInfoCsmc() [1/2]

MdcSegInfoCsmc::~MdcSegInfoCsmc ( )

inline

Definition at line 31 of file InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h.

{ };

MdcSegInfoCsmc() [2/2]

MdcSegInfoCsmc::MdcSegInfoCsmc ( )

inline

Definition at line 30 of file Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-48/MdcTrkRecon/MdcSegInfoCsmc.h.

{ };

~MdcSegInfoCsmc() [2/2]

MdcSegInfoCsmc::~MdcSegInfoCsmc ( )

inline

Definition at line 31 of file Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-48/MdcTrkRecon/MdcSegInfoCsmc.h.

{ };

Member Function Documentation

calcStraight() [1/4]

void MdcSegInfoCsmc::calcStraight

(

const MdcSeg *

parentSeg

)

Definition at line 99 of file MdcSegInfoCsmc.cxx.

                                                    {
//------------------------------------------------------------------
  double slope = parentSeg->slope();
  double radius = parentSeg->superlayer()->rad0();
  double phi = parentSeg->phi();
  const double *inErr = parentSeg->errmat();
 
  calcStraight(phi, slope, radius, inErr);
}

calcStraight() [2/4]

void MdcSegInfoCsmc::calcStraight

(

const MdcSeg *

parentSeg

)

calcStraight() [3/4]

void MdcSegInfoCsmc::calcStraight	(	double	phi,
		double	slope,
		double	rad,
		const double *	inErr
	)

Definition at line 38 of file MdcSegInfoCsmc.cxx.

                                                   {
//------------------------------------------------------------------
  // Load the segment info obj with phi0 and d0 values, calc. assuming 
  // the segments lie on straight tracks.
  //     
  // d0 defined to be positive for tracks that intercept x axis on positive 
  //   phi (neg x) side of origin.  At least that's what I claim. */
  
  // error matrix conversion: 
  //   Vij(f) = SUM(n,m) ( d(fi)/d(yn) * d(fj)/d(ym) * Vnm(y) )
  //   The derivatives are calculated from:
  //       d0 = r**2 * slope / sqrt( 1 + r**2 * slope**2 )   &
  //       phi0 = phiseg - asin(d0 / r) 
  
  //    d0ovSlope = r**2/sqrt(1 + slope**2 * r**2) = d0 / slope
  //    d0ovrad = d0/r
  //    d0ovradsq = (d0/r)**2
  //    dphi0dphi = d(phi0)/d(seg-phi)
  //    dphi0ds   = d(phi0)/d(seg-slope)
  //    dd0ds   = d(d0)/d(seg-slope)
  //    dd0dphi = 0
  
  double d0ovSlope = -rad * rad / sqrt( 1. + slope*slope * rad*rad );
  double d0 = slope * d0ovSlope;
  _par0 = d0;
 
  double d0ovRad =  d0 / rad ;
  // Don't blow away arcsine:
  d0ovRad = ( d0ovRad > 1.0 || d0ovRad < -1.0) ? 
    (d0ovRad+0.00001)/fabs(d0ovRad) : d0ovRad; 
  BesAngle phi0 = phi - asin( d0ovRad );
  _par1 = phi0.posRad();
 
  // Error matrix (anal calculation).
  //  double dphi0dphi = 1.;
  double dphi0ds = -d0ovSlope * sqrt( 1. - d0ovRad*d0ovRad ) / rad;
  double dd0ds = d0ovSlope * ( 1. - d0ovRad*d0ovRad );
 
  // d0:
  _errmat[0] = inErr[2] * dd0ds * dd0ds;
  if (_errmat[0] < 0.) _errmat[0] = 0.;
 
  // phi0:
  _errmat[2] = inErr[2] * dphi0ds * dphi0ds  + inErr[0]  +  
    inErr[1] * 2. * dphi0ds;
  if (_errmat[2] < 0.) _errmat[2] = 0.;
  // phi0|d0:
  _errmat[1] = inErr[2] * dd0ds * dphi0ds + 
               inErr[1] * dd0ds;
 
  int error = mdcTwoInv(_errmat, _inverr);
  if (error) {
    std::cout << " ErrMsg(warning) " 
      << "Failed to invert matrix -- MdcSegInfo::calcStraight" << endl 
      << _errmat[0] << " " << _errmat[1] << " " << _errmat[2]<< std::endl;
  }
}

Referenced by calcStraight(), and MdcSegGrouperCsmc::fillWithSegs().

calcStraight() [4/4]

void MdcSegInfoCsmc::calcStraight	(	double	phi,
		double	slope,
		double	rad,
		const double *	inErr
	)

d0() [1/2]

double MdcSegInfoCsmc::d0 ( ) const

inline

Definition at line 33 of file InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h.

{return _par0;}

Referenced by calcStraight(), and MdcSegGrouperCsmc::incompWithSeg().

d0() [2/2]

double MdcSegInfoCsmc::d0 ( ) const

inline

Definition at line 33 of file Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-48/MdcTrkRecon/MdcSegInfoCsmc.h.

{return _par0;}

parIsAngle() [1/2]

bool MdcSegInfoCsmc::parIsAngle ( int  i ) const

virtual

Implements MdcSegInfo.

Definition at line 31 of file MdcSegInfoCsmc.cxx.

                                      {
//---------------------------------------------------------------------------
  assert (i >= 0 && i < 2);
  return (0 != i);   //i.e., parameter 1 is an angle
}

parIsAngle() [2/2]

bool MdcSegInfoCsmc::parIsAngle ( int  i ) const

virtual

Implements MdcSegInfo.

phi0() [1/2]

double MdcSegInfoCsmc::phi0 ( ) const

inline

Definition at line 34 of file InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h.

{return _par1;}

Referenced by calcStraight(), MdcSegGrouperCsmc::fillWithSegs(), and MdcSegGrouperCsmc::incompWithSeg().

phi0() [2/2]

double MdcSegInfoCsmc::phi0 ( ) const

inline

Definition at line 34 of file Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-48/MdcTrkRecon/MdcSegInfoCsmc.h.

{return _par1;}

sigD0() [1/2]

double MdcSegInfoCsmc::sigD0

(

)

const

Definition at line 117 of file MdcSegInfoCsmc.cxx.

                            {
//------------------------------------------------------------------
  return sqrt(_errmat[0]);
}

Referenced by MdcSegGrouperCsmc::incompWithSeg().

sigD0() [2/2]

double MdcSegInfoCsmc::sigD0

(

)

const

sigPhi0() [1/2]

double MdcSegInfoCsmc::sigPhi0

(

)

const

Definition at line 111 of file MdcSegInfoCsmc.cxx.

                              {
//------------------------------------------------------------------
  return sqrt(_errmat[2]);
}

Referenced by MdcSegGrouperCsmc::incompWithSeg().

sigPhi0() [2/2]

double MdcSegInfoCsmc::sigPhi0

(

)

const

---

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

• InstallArea/include/MdcTrkRecon/MdcTrkRecon/MdcSegInfoCsmc.h
• Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-48/MdcTrkRecon/MdcSegInfoCsmc.h
• MdcSegInfoCsmc.cxx