TMDCUtil.cxx File Reference

#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TMDCWire.h"
#include "TrkReco/TMDCWireHit.h"
#include "TrkReco/TMDCWireHitMC.h"
#include "TrkReco/TMLink.h"

Go to the source code of this file.

Macros
#define	SRTOPI   0.7978846
#define	UPL   340.0
#define	ROOT2I   0.70710678

Functions
float	CathodeSectorId (unsigned id)
	geocdc utilities
void	bitDisplay (unsigned val)
void	bitDisplay (unsigned val, unsigned f, unsigned l)
int	intersection (const HepPoint3D &c1, double r1, const HepPoint3D &c2, double r2, double eps, HepPoint3D &x1, HepPoint3D &x2)
	Circle utilities.
double	chisq2confLevel (int n, double chi2)
	ALPHA = 10000. / 2.99792458 / 15.

Variables
const HepPoint3D	ORIGIN = HepPoint3D( 0., 0., 0. )
	Constants.

Macro Definition Documentation

ROOT2I

#define ROOT2I   0.70710678

Referenced by chisq2confLevel().

SRTOPI

#define SRTOPI   0.7978846

Referenced by chisq2confLevel().

UPL

#define UPL   340.0

Referenced by chisq2confLevel().

Function Documentation

bitDisplay() [1/2]

void bitDisplay

(

unsigned

val

)

Definition at line 82 of file TMDCUtil.cxx.

{ bitDisplay( val, 31, 0 ); }

Referenced by bitDisplay(), TMDCWireHit::dump(), and TTrackMC::dump().

bitDisplay() [2/2]

void bitDisplay	(	unsigned	val,
		unsigned	f,
		unsigned	l )

Definition at line 84 of file TMDCUtil.cxx.

                                                        {
  unsigned i;
  for ( i = 0; i < f - l; i++ )
  {
    if ( ( i % 8 ) == 0 ) std::cout << " ";
    std::cout << ( val >> ( f - i ) ) % 2;
  }
}

CathodeSectorId()

float CathodeSectorId

(

unsigned

id

)

geocdc utilities

Definition at line 59 of file TMDCUtil.cxx.

                                     {
 
  unsigned layer = id / 64;
 
  if ( layer == 0 ) { return int( id / 8 ); }
 
  if ( layer == 1 )
  {
    if ( id >= 127 ) id -= 64;
    if ( ( id - 6 ) % 8 == 0 ) return ( id - 6 ) / 8 + 0.5;
    else return int( ( id + 1 ) / 8 );
  }
 
  if ( layer == 2 )
  {
    if ( id <= 129 ) id += 64;
    return int( ( id - 2 ) / 8 );
  }
 
  return 9999;
}

chisq2confLevel()

double chisq2confLevel	(	int	n,
		double	chi2 )

ALPHA = 10000. / 2.99792458 / 15.

Converter

Definition at line 147 of file TMDCUtil.cxx.

                                             {
#define SRTOPI 0.7978846
#define UPL 340.0
#define ROOT2I 0.70710678
 
  double prob = 0.0;
  double sum, term;
  int    m;
  int    i, k;
  double temp_i, temp_n;
  double srty;
 
  if ( ( n <= 0 ) || ( chi2 < 0.0 ) ) { return prob; }
  if ( n > 60 )
  {
    temp_n = (double)n;
    srty   = sqrt( chi2 ) - sqrt( temp_n - 0.5 );
    if ( srty < 12.0 )
    {
      prob = 0.5 * erfc( srty );
      return prob;
    }
    return prob;
  }
  if ( chi2 > UPL ) { return prob; }
  sum  = exp( -0.5 * chi2 );
  term = sum;
  m    = (int)floor( n / 2. );
 
  if ( 2 * m == n )
  {
    if ( m == 1 )
    {
      prob = sum;
      return prob;
    }
    else
    {
      for ( i = 2; i < m + 1; i++ )
      {
        temp_i = (double)i;
        term   = 0.5 * chi2 * term / ( temp_i - 1.0 );
        sum    = sum + term;
      }
      prob = sum;
      return prob;
    }
  }
  else
  {
    srty = sqrt( chi2 );
    prob = erfc( ROOT2I * srty );
    if ( n == 1 ) { return prob; }
    if ( n == 3 )
    {
      prob = SRTOPI * srty * sum + prob;
      return prob;
    }
    else
    {
      k = m - 1;
      for ( i = 1; i < k + 1; i++ )
      {
        temp_i = (double)i;
        term   = term * chi2 / ( 2.0 * temp_i + 1.0 );
        sum    = sum + term;
      }
      prob = SRTOPI * srty * sum + prob;
      return prob;
    }
  }
}

Referenced by TTrack::confidenceLevel().

intersection()

int intersection	(	const HepPoint3D &	c1,
		double	r1,
		const HepPoint3D &	c2,
		double	r2,
		double	eps,
		HepPoint3D &	x1,
		HepPoint3D &	x2 )

Circle utilities.

Definition at line 93 of file TMDCUtil.cxx.

                                                   {
 
  double c0x  = c2.x() - c1.x();
  double c0y  = c2.y() - c1.y();
  double c0   = sqrt( c0x * c0x + c0y * c0y );
  double rr1  = abs( r1 );
  double rr2  = abs( r2 );
  double Radd = rr1 + rr2;
  double Rsub = abs( rr1 - rr2 );
 
  // no intersections
 
  if ( c0 > Radd + eps || c0 < 0.001 || c0 < Rsub - eps )
  {
    //-- debug
    // std::cout << "Int2Cir return 0 " << std::endl;
    //-- debug end
    return 0;
  }
 
  // single intersection
 
  else
  {
    if ( c0 > Radd - eps )
    {
      x1.setX( c1.x() + rr1 * c0x / c0 );
      x1.setY( c1.y() + rr1 * c0y / c0 );
      x2.setX( 0.0 );
      x2.setY( 0.0 );
      //--debug
      // std::cout << "Int2Cir return 1" << std::endl;
      //--debug end
      return 1;
    }
  }
 
  // two intersections
 
  double chg    = abs( r1 ) / r1;
  double cosPsi = ( c0 * c0 + rr1 * rr1 - rr2 * rr2 ) / ( 2. * c0 * rr1 );
  double sinPsi = -( chg / abs( chg ) ) * sqrt( 1.0 - cosPsi * cosPsi );
  x1.setX( c1.x() + ( rr1 / c0 ) * ( cosPsi * c0x - sinPsi * c0y ) );
  x1.setY( c1.y() + ( rr1 / c0 ) * ( cosPsi * c0y + sinPsi * c0x ) );
  x2.setX( c1.x() + ( rr1 / c0 ) * ( cosPsi * c0x + sinPsi * c0y ) );
  x2.setY( c1.y() + ( rr1 / c0 ) * ( cosPsi * c0y - sinPsi * c0x ) );
  //-- debug
  // std::cout << "Int2Cir return 2" << std::endl;
  //-- debug end
  return 2;
}

Referenced by MucGeoGeneral::FindIntersectGaps(), MucGeoGeneral::FindIntersections(), MucGeoGeneral::FindIntersectStrips(), MucGeoGeneral::FindIntersectStrips(), TTrack::HelCyl(), main(), and RecMucTrack::Project().

Variable Documentation

ORIGIN

const HepPoint3D ORIGIN = HepPoint3D( 0., 0., 0. )

Constants.

Definition at line 47 of file TMDCUtil.cxx.

Referenced by TConformalFinder0::conformalTransformation(), TConformalFinder0::conformalTransformationDriftCircle(), TConformalFinder0::conformalTransformationRphi(), TTrackManager::divideByIp(), TConformalFinder0::doit(), TConformalFinder::doit(), TFastFinder::doit(), TrkReco::execute(), T3DLineFitter::fit(), TRungeFitter::fit(), TTrackManager::maskCurl(), TTrackManager::merge(), TTrack::szPosition(), T3DLine::T3DLine(), TRunge::TRunge(), TTrack::TTrack(), TTrack::TTrack(), and MdcSegInfoSterO::zPosition().