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

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 MdcSegInfoCsmc.h.

Constructor & Destructor Documentation

MdcSegInfoCsmc()

MdcSegInfoCsmc::MdcSegInfoCsmc ( )

inline

Definition at line 30 of file MdcSegInfoCsmc.h.

{ };

~MdcSegInfoCsmc()

MdcSegInfoCsmc::~MdcSegInfoCsmc ( )

inline

Definition at line 31 of file MdcSegInfoCsmc.h.

{ };

Member Function Documentation

calcStraight() [1/2]

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/2]

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().

d0()

double MdcSegInfoCsmc::d0 ( ) const

inline

Definition at line 33 of file MdcSegInfoCsmc.h.

{return _par0;}

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

parIsAngle()

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
}

phi0()

double MdcSegInfoCsmc::phi0 ( ) const

inline

Definition at line 34 of file MdcSegInfoCsmc.h.

{return _par1;}

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

sigD0()

double MdcSegInfoCsmc::sigD0

(

)

const

Definition at line 117 of file MdcSegInfoCsmc.cxx.

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

Referenced by MdcSegGrouperCsmc::incompWithSeg().

sigPhi0()

double MdcSegInfoCsmc::sigPhi0

(

)

const

Definition at line 111 of file MdcSegInfoCsmc.cxx.

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

Referenced by MdcSegGrouperCsmc::incompWithSeg().

---

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

• MdcSegInfoCsmc.h
• MdcSegInfoCsmc.cxx