MdcTrackUtil Class Reference

#include <MdcTrackUtil.h>

Public Member Functions
	MdcTrackUtil ()
	~MdcTrackUtil ()
int	nLayerTrackPassed (const HepVector helix)
int	nLayerTrackPassed (const double helix[5])
HepVector	patRecPar2BesPar (const HepVector &helixPar)
HepSymMatrix	patRecErr2BesErr (const HepSymMatrix &err)
	MdcTrackUtil ()
	~MdcTrackUtil ()
int	nLayerTrackPassed (const HepVector helix)
int	nLayerTrackPassed (const double helix[5])
HepVector	patRecPar2BesPar (const HepVector &helixPar)
HepSymMatrix	patRecErr2BesErr (const HepSymMatrix &err)
	MdcTrackUtil ()
	~MdcTrackUtil ()
int	nLayerTrackPassed (const HepVector helix)
int	nLayerTrackPassed (const double helix[5])
HepVector	patRecPar2BesPar (const HepVector &helixPar)
HepSymMatrix	patRecErr2BesErr (const HepSymMatrix &err)

Static Public Member Functions
static MdcTrackUtil *	instance ()
static MdcTrackUtil *	instance ()
static MdcTrackUtil *	instance ()

Detailed Description

Definition at line 8 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/MdcTrackUtil.h.

Constructor & Destructor Documentation

MdcTrackUtil() [1/3]

MdcTrackUtil::MdcTrackUtil

(

)

Definition at line 35 of file MdcTrackUtil.cxx.

                           {
  // Initalze magnetic flied
  m_pIMF = service( "MagneticFieldSvc" );
  if ( !m_pIMF ) { std::cout << " ERROR Unable to open Magnetic field service " << std::endl; }
  // get Bz for Check TEMP, Bz may be changed by run
  double gaussToTesla = 1000.;
  Bz                  = m_pIMF->getReferField() * gaussToTesla;
 
  // Initialize MdcGeomSvc
  m_mdcGeomSvc = service( "MdcGeomSvc" );
  if ( !m_mdcGeomSvc ) { std::cout << " FATAL Could not load MdcGeomSvc!  " << std::endl; }
}

Referenced by instance(), and patRecErr2BesErr().

~MdcTrackUtil() [1/3]

MdcTrackUtil::~MdcTrackUtil ( )

inline

Definition at line 12 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/MdcTrackUtil.h.

{};

MdcTrackUtil() [2/3]

MdcTrackUtil::MdcTrackUtil

(

)

~MdcTrackUtil() [2/3]

MdcTrackUtil::~MdcTrackUtil ( )

inline

Definition at line 12 of file InstallArea/x86_64-el9-gcc13-opt/include/TrackUtil/MdcTrackUtil.h.

{};

MdcTrackUtil() [3/3]

MdcTrackUtil::MdcTrackUtil

(

)

~MdcTrackUtil() [3/3]

MdcTrackUtil::~MdcTrackUtil ( )

inline

Definition at line 12 of file Reconstruction/TrackUtil/include/TrackUtil/MdcTrackUtil.h.

{};

Member Function Documentation

instance() [1/3]

MdcTrackUtil * MdcTrackUtil::instance ( )

static

Definition at line 30 of file MdcTrackUtil.cxx.

                                     {
  if ( 0 == _myself ) { _myself = new MdcTrackUtil(); }
  return _myself;
}

instance() [2/3]

MdcTrackUtil * MdcTrackUtil::instance ( )

static

instance() [3/3]

MdcTrackUtil * MdcTrackUtil::instance ( )

static

nLayerTrackPassed() [1/6]

int MdcTrackUtil::nLayerTrackPassed

(

const double

helix[5]

)

Definition at line 57 of file MdcTrackUtil.cxx.

                                                           {
  int nLayer = 0;
 
  for ( unsigned iLayer = 0; iLayer < 43; iLayer++ )
  {
    // flightLength is the path length of track in the x-y plane
    // guess flightLength by the radius in middle of the wire.
    double rMidLayer    = m_mdcGeomSvc->Layer( iLayer )->Radius();
    double flightLength = rMidLayer;
 
    HepPoint3D pivot( 0., 0., 0. );
    double     dz    = helix[3];
    double     c     = CLHEP::c_light * 100.; // unit from m/s to cm/s
    double     alpha = 1 / ( c * Bz );        //~333.567
    double     kappa = helix[2];
    double     rc    = ( -1. ) * alpha / kappa;
    // std::cout<<"MdcTrackUtil alpha   "<<alpha<<std::endl;
    double tanl = helix[4];
    double phi0 = helix[1];
    double phi  = flightLength / rc + phi0; // turning angle
    double z    = pivot.z() + dz - ( alpha / kappa ) * tanl * phi;
 
    double layerHalfLength = m_mdcGeomSvc->Layer( iLayer )->Length() / 2.;
 
    // std::cout<<"MdcTrackUtil length  "<<layerHalfLength<<std::endl;
 
    if ( fabs( z ) < fabs( layerHalfLength ) ) ++nLayer;
  }
 
  return nLayer;
}

nLayerTrackPassed() [2/6]

int MdcTrackUtil::nLayerTrackPassed

(

const double

helix[5]

)

nLayerTrackPassed() [3/6]

int MdcTrackUtil::nLayerTrackPassed

(

const double

helix[5]

)

nLayerTrackPassed() [4/6]

int MdcTrackUtil::nLayerTrackPassed

(

const HepVector

helix

)

Definition at line 49 of file MdcTrackUtil.cxx.

                                                           {
  double helixParam[5];
  for ( int i = 0; i < 5; ++i ) helixParam[i] = helix[i];
 
  return nLayerTrackPassed( helixParam );
}

Referenced by nLayerTrackPassed().

nLayerTrackPassed() [5/6]

int MdcTrackUtil::nLayerTrackPassed

(

const HepVector

helix

)

nLayerTrackPassed() [6/6]

int MdcTrackUtil::nLayerTrackPassed

(

const HepVector

helix

)

patRecErr2BesErr() [1/3]

HepSymMatrix MdcTrackUtil::patRecErr2BesErr

(

const HepSymMatrix &

err

)

Definition at line 107 of file MdcTrackUtil.cxx.

                                                                     {
  // error matrix
  // std::cout<<" err1  "<<err<<" "<<err.num_row()<<std::endl;
  // V(Y)=S * V(X) * ST , mS = S , mVy = V(Y) , err() = V(X)
  // int n = err.num_row();
  HepSymMatrix mS( err.num_row(), 0 );
  mS[0][0]         = -1.; // dr0=-d0
  mS[1][1]         = 1.;
  mS[2][2]         = Bz / -333.567; // pxy -> cpa
  mS[3][3]         = 1.;
  mS[4][4]         = 1.;
  HepSymMatrix mVy = err.similarity( mS );
  // std::cout<<" err2  "<<n<<" "<<mVy<<std::endl;
  return mVy;
}

patRecErr2BesErr() [2/3]

HepSymMatrix MdcTrackUtil::patRecErr2BesErr

(

const HepSymMatrix &

err

)

patRecErr2BesErr() [3/3]

HepSymMatrix MdcTrackUtil::patRecErr2BesErr

(

const HepSymMatrix &

err

)

patRecPar2BesPar() [1/3]

HepVector MdcTrackUtil::patRecPar2BesPar

(

const HepVector &

helixPar

)

Definition at line 90 of file MdcTrackUtil.cxx.

                                                                    {
  HepVector helix( 5, 0 );
  double    d0    = -helixPar[0]; // cm
  double    phi0  = helixPar[1] + CLHEP::halfpi;
  double    omega = Bz * helixPar[2] / -333.567;
  double    z0    = helixPar[3]; // cm
  double    tanl  = helixPar[4];
  helix[0]        = d0;
  helix[1]        = phi0;
  helix[2]        = omega;
  helix[3]        = z0;
  helix[4]        = tanl;
  // std::cout<<"helix   "<<helix<<std::endl;
  return helix;
}

patRecPar2BesPar() [2/3]

HepVector MdcTrackUtil::patRecPar2BesPar

(

const HepVector &

helixPar

)

patRecPar2BesPar() [3/3]

HepVector MdcTrackUtil::patRecPar2BesPar

(

const HepVector &

helixPar

)

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TrackUtil/MdcTrackUtil.h
• InstallArea/x86_64-el9-gcc13-opt/include/TrackUtil/MdcTrackUtil.h
• Reconstruction/TrackUtil/include/TrackUtil/MdcTrackUtil.h
• MdcTrackUtil.cxx