FastVertexFit Class Reference

#include <FastVertexFit.h>

Public Member Functions
	~FastVertexFit ()
void	init ()
void	addTrack (const int number, const HepVector helix, const HepSymMatrix err)
bool	Fit ()
double	chisq () const
VertexParameter	vtx () const
HepVector	Vx () const
HepSymMatrix	Evx () const
HepVector	Pull () const
	~FastVertexFit ()
void	init ()
void	addTrack (const int number, const HepVector helix, const HepSymMatrix err)
bool	Fit ()
double	chisq () const
VertexParameter	vtx () const
HepVector	Vx () const
HepSymMatrix	Evx () const
HepVector	Pull () const

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

Detailed Description

Definition at line 18 of file Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h.

Constructor & Destructor Documentation

~FastVertexFit() [1/2]

FastVertexFit::~FastVertexFit

(

)

Definition at line 15 of file FastVertexFit.cxx.

{;}

~FastVertexFit() [2/2]

FastVertexFit::~FastVertexFit

(

)

Member Function Documentation

addTrack() [1/2]

void FastVertexFit::addTrack	(	const int	number,
		const HepVector	helix,
		const HepSymMatrix	err
	)

Definition at line 29 of file FastVertexFit.cxx.

                                                                                            {
  int ifail;
  int ierr;
  //HepSymMatrix err0(5, 0); //inverse of err;
  //err0 = err.inverse(ifail);
  HepSymMatrix Wc(2, 0);
  Wc[0][0] = err[0][0];
  Wc[0][1] = Wc[1][0] = err[0][3];
  Wc[1][1] = err[3][3];
  Wc = Wc.inverse(ifail);
 
  HepVector p0(3, 0), v0(3, 0);
  double pxy = 1./fabs(helix[2]);
  p0[0] = 0 - pxy*sin(helix[1]); //px
  p0[1] = pxy*cos(helix[1]);     //py
  p0[2] = pxy * helix[4];        //pz
  v0[0] = helix[0]*cos(helix[1]);//x
  v0[1] = helix[0]*sin(helix[1]);//y
  v0[2] = helix[3];              //z  
 
  HepPoint3D vv0(v0[0], v0[1], v0[2]);
  double bField = VertexFitBField::instance()->getBFieldZRef();
 
  int charge = (helix[2]>0 ? +1 :-1);
  double a = alpha * bField * charge;
  double T0 = sqrt((p0[0]+a*p0[1])*(p0[0]+a*p0[1])+(p0[1]-a*p0[0])*(p0[1]-a*p0[0]));
 
  HepMatrix Dc(2, 3, 0);
  Dc[0][0] = (p0[1] - a*v0[0])/T0;
  Dc[0][1] = 0 - (p0[0] + a*v0[1])/T0;
  Dc[1][0] = 0 - (p0[2]/T0) * (p0[0] + a*v0[1])/T0;
  Dc[1][1] = 0 - (p0[2]/T0) * (p0[1] - a*v0[0])/T0;
  Dc[1][2] = 1;
 
  m_D.push_back(Dc);
  m_W.push_back(Wc);
 
  HepSymMatrix DTWD(3, 0);
  DTWD = Wc.similarity(Dc.T());
  HepVector qTrk(2, 0);
  qTrk[0] = helix[0];
  qTrk[1] = helix[3];
  //HepVector xp(3, 0);
  //xp = Dc.inverse(ierr) * qTrk;
  m_DTWD.push_back(DTWD);
  m_xp.push_back(v0);
}

Referenced by DQASelHadron::execute().

addTrack() [2/2]

void FastVertexFit::addTrack	(	const int	number,
		const HepVector	helix,
		const HepSymMatrix	err
	)

chisq() [1/2]

double FastVertexFit::chisq ( ) const

inline

Definition at line 38 of file Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h.

{return m_chisq;}

Referenced by Fit().

chisq() [2/2]

double FastVertexFit::chisq ( ) const

inline

Definition at line 38 of file InstallArea/include/VertexFit/VertexFit/FastVertexFit.h.

{return m_chisq;}

Evx() [1/2]

HepSymMatrix FastVertexFit::Evx ( ) const

inline

Definition at line 41 of file Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h.

{return m_Evx;}   

Evx() [2/2]

HepSymMatrix FastVertexFit::Evx ( ) const

inline

Definition at line 41 of file InstallArea/include/VertexFit/VertexFit/FastVertexFit.h.

{return m_Evx;}   

Fit() [1/2]

bool FastVertexFit::Fit

(

)

Definition at line 77 of file FastVertexFit.cxx.

                        {
  bool fitResult = false;
 
  int ifail;
  HepSymMatrix total_DTWD(3, 0);
  HepVector total_xp(3, 0);
 
  for(int i = 0; i < m_DTWD.size(); i++) {
    total_DTWD += m_DTWD[i];  
    total_xp += m_DTWD[i]*m_xp[i]; 
  }
  m_Vx = total_DTWD.inverse(ifail) * total_xp;
  m_Evx = total_DTWD.inverse(ifail);
 
  double chisq = 0;
  for(int i = 0; i < m_xp.size(); i++) {
    double chi2 = (m_DTWD[i].similarity((m_xp[i] - m_Vx).T()))[0][0];  
    chisq += chi2;
  }
  m_chisq = chisq;
 
  fitResult = true;
  return fitResult;
}

Referenced by DQASelHadron::execute().

Fit() [2/2]

bool FastVertexFit::Fit

(

)

init() [1/2]

void FastVertexFit::init

(

)

Definition at line 19 of file FastVertexFit.cxx.

                         {
  m_D.clear();
  m_W.clear();
  m_DTWD.clear();
  m_xp.clear();
  m_chisq = 999;
  m_Vx = HepVector(3, 0);
  m_Evx = HepSymMatrix(3,0);
}

Referenced by DQASelHadron::execute().

init() [2/2]

void FastVertexFit::init

(

)

instance() [1/2]

FastVertexFit * FastVertexFit::instance ( )

static

Definition at line 9 of file FastVertexFit.cxx.

                                        {
  if(m_pointer) return m_pointer;
  m_pointer = new FastVertexFit();
  return m_pointer;
}

Referenced by DQASelHadron::execute().

instance() [2/2]

static FastVertexFit * FastVertexFit::instance ( )

static

Pull() [1/2]

HepVector FastVertexFit::Pull

(

)

const

Definition at line 102 of file FastVertexFit.cxx.

                                    {
  HepVector pull(3, 0);
  pull[0] = m_Vx[0]/sqrt(m_Evx[0][0]);
  pull[1] = m_Vx[1]/sqrt(m_Evx[1][1]);
  pull[2] = m_Vx[2]/sqrt(m_Evx[2][2]);
  return pull;
}

Pull() [2/2]

HepVector FastVertexFit::Pull

(

)

const

vtx() [1/2]

VertexParameter FastVertexFit::vtx ( ) const

inline

Definition at line 39 of file Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h.

{return m_vtx_infit;}

vtx() [2/2]

VertexParameter FastVertexFit::vtx ( ) const

inline

Definition at line 39 of file InstallArea/include/VertexFit/VertexFit/FastVertexFit.h.

{return m_vtx_infit;}

Vx() [1/2]

HepVector FastVertexFit::Vx ( ) const

inline

Definition at line 40 of file Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h.

{return m_Vx;}

Referenced by DQASelHadron::execute().

Vx() [2/2]

HepVector FastVertexFit::Vx ( ) const

inline

Definition at line 40 of file InstallArea/include/VertexFit/VertexFit/FastVertexFit.h.

{return m_Vx;}

---

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

• Analysis/VertexFit/VertexFit-00-02-82/VertexFit/FastVertexFit.h
• InstallArea/include/VertexFit/VertexFit/FastVertexFit.h
• FastVertexFit.cxx