MdcxHit Class Reference

#include <MdcxHit.h>

Public Member Functions
	MdcxHit (const MdcDigi *pdcdatum, float c0=0, float cresol=.0180)
	MdcxHit (const MdcHit *pdchhit, float c0=0, float cresol=.0180)
virtual	~MdcxHit ()
const MdcDigi *	getDigi () const
const MdcHit *	getMdcHit () const
int	SuperLayer () const
int	WireNo () const
int	Layer () const
float	t () const
float	x () const
float	y () const
float	xpos () const
float	ypos () const
float	xneg () const
float	yneg () const
float	wx () const
float	wy () const
float	wz () const
float	pw () const
float	phiMid () const
float	stereo () const
float	v () const
int	type () const
void	SetConstErr (int i)
float	tcor (float zh=0.0, float tof=0.0, float tzero=0.0) const
float	d (MdcxHel &hel) const
float	d (float zh=0.0, float tof=0.0, float tzero=0.0, int wamb=0, float eang=0.0) const
float	e (float dd=0.0) const
float	pull (MdcxHel &hel) const
float	residual (MdcxHel &hel) const
std::vector< float >	derivatives (MdcxHel &hel) const
void	print (std::ostream &o, int i=0) const
void	printAll (std::ostream &o, int i=0) const
void	SetUsedOnHel (int i)
int	GetUsedOnHel () const

Static Public Member Functions
static void	setMdcCalibFunSvc (const MdcCalibFunSvc *calibSvc)
static void	setCountPropTime (bool countPropTime)
static void	setMdcDetector (const MdcDetector *gm)

Protected Member Functions
void	process ()

Protected Attributes
unsigned	_wirenumber
unsigned	_layernumber
unsigned	_superlayer
float	_t
float	_q
unsigned	_iTdc
unsigned	_iAdc
float	_T0Walk
float	_x
float	_y
float	_r
float	_s
float	_L
float	_d
float	_v
int	_type
float	_e
float	_xpos
float	_ypos
float	_xneg
float	_yneg
float	_p
float	_sp
float	_cp
double	_pw
double	_wx
double	_wy
double	_wz
const MdcHit *	_mdcHit
const MdcDigi *	_mdcDigi
int	_consterr
int	usedonhel
float	_c0
float	_cresol

Detailed Description

Definition at line 42 of file MdcxHit.h.

Constructor & Destructor Documentation

MdcxHit() [1/2]

MdcxHit::MdcxHit	(	const MdcDigi *	pdcdatum,
		float	c0 = 0,
		float	cresol = .0180
	)

Definition at line 48 of file MdcxHit.cxx.

    :_mdcHit(0), _mdcDigi(pdcdatum), _c0(c0), _cresol(cresol)
{
  process();
}

MdcxHit() [2/2]

MdcxHit::MdcxHit	(	const MdcHit *	pdchhit,
		float	c0 = 0,
		float	cresol = .0180
	)

Definition at line 54 of file MdcxHit.cxx.

    :_mdcHit(pdchhit),_mdcDigi(pdchhit->digi()), _c0(c0), _cresol(cresol)
{
  process();
}

~MdcxHit()

virtual MdcxHit::~MdcxHit ( )

inlinevirtual

Definition at line 50 of file MdcxHit.h.

{ };

Member Function Documentation

d() [1/2]

float MdcxHit::d	(	float	zh = 0.0,
		float	tof = 0.0,
		float	tzero = 0.0,
		int	wamb = 0,
		float	eang = 0.0
	)		const

Definition at line 135 of file MdcxHit.cxx.

                                                                                    {
 
  //std::cout<<__FILE__<<" "<<__LINE__<<" call entrance "<< entranceAngle<< std::endl;
  //tof = hel.Doca_Tof();//yzhang delete for csmc temp FIXME
  float t_corr = tcor(z,tof,tzero);
  double eAngle = EntranceAngle(entranceAngle);
 
  //yzhang add for 64-bit
  if (fabs(z)>150. || fabs(t_corr)>1500. || fabs(eAngle)>999){
    return 9999.;
  }
  //zhangy
 
  int lrCalib = 2;
  if (wamb==1)lrCalib = 0;
  else if (wamb==-1) lrCalib = 1;
 
  double driftD = 0.1 * m_mdcCalibFunSvc->driftTimeToDist(t_corr,_layernumber,_wirenumber,lrCalib,eAngle);//to cm
  //std::cout<<"MdcxHit ("<<_layernumber<<","<<_wirenumber<<" dd "<<driftD<<" dt "<<t_corr<<" lr "<<lrCalib<<" eAngle "<<eAngle<<std::endl;
 
 
  if (fabs(driftD)<=0.0001) driftD = 0.001;
  return driftD;
} 

d() [2/2]

float MdcxHit::d

(

MdcxHel &

hel

)

const

Definition at line 160 of file MdcxHit.cxx.

                                   {
  float doca=hel.Doca(*this); // changes hel's internal state...
  return d(hel.Doca_Zh(),hel.Doca_Tof(),hel.T0(),
      hel.Doca_Wamb(),hel.Doca_Eang());
}//endof d

Referenced by d(), derivatives(), printAll(), process(), pull(), residual(), and MdcxFindSegs::trial().

derivatives()

std::vector< float > MdcxHit::derivatives

(

MdcxHel &

hel

)

const

Definition at line 192 of file MdcxHit.cxx.

                                                        {
  // compute derivatives for |hel|
  std::vector<float> deriv = hel.derivatives(*this);
  float dtemp = d(hel.Doca_Zh(), hel.Doca_Tof(), hel.T0(),
      hel.Doca_Wamb(), hel.Doca_Eang());
  //float dtemp = d(hel);
  deriv[0] = dtemp - deriv[0];
  //  deriv[0] -= v()*hel.T0();
  float ewire = e(dtemp);
  for (uint32_t i = 0; i < deriv.size(); i++) deriv[i] /= ewire;
  return deriv;
}//endof derivatives

e()

float MdcxHit::e

(

float

dd = 0.0

)

const

Definition at line 166 of file MdcxHit.cxx.

                              {
  //if (0!=_consterr) return _cresol;//yzhang delete 2009-11-03 
  float errTemp = fabs(getSigma(dd));
  //  protect against resolution = 0; set min resol to 1 nm
  float errMin = 1.0e-7;
  return errTemp<errMin?errMin:errTemp;
}

Referenced by derivatives(), printAll(), and pull().

getDigi()

const MdcDigi * MdcxHit::getDigi ( ) const

inline

Definition at line 63 of file MdcxHit.h.

{return _mdcDigi;}

Referenced by print().

getMdcHit()

const MdcHit * MdcxHit::getMdcHit ( ) const

inline

Definition at line 64 of file MdcxHit.h.

{return _mdcHit;}

Referenced by printAll().

GetUsedOnHel()

int MdcxHit::GetUsedOnHel ( ) const

inline

Definition at line 103 of file MdcxHit.h.

{return usedonhel;}

Referenced by MdcxAddHits::GetAssociates().

Layer()

int MdcxHit::Layer ( ) const

inline

Definition at line 67 of file MdcxHit.h.

{return _layernumber;}// layer#

Referenced by MdcxAddHits::GetAssociates(), print(), printAll(), and MdcxFindSegs::process().

phiMid()

float MdcxHit::phiMid ( ) const

inline

Definition at line 79 of file MdcxHit.h.

{return _pw;} //yzhang add

print()

void MdcxHit::print	(	std::ostream &	o,
		int	i = 0
	)		const

Definition at line 205 of file MdcxHit.cxx.

                                           {
  o << "(" << Layer() << "," << WireNo() << ",id "<< getDigi()->getTrackIndex()<<") " ;
}

Referenced by MdcxAddHits::GetAssociates().

printAll()

void MdcxHit::printAll	(	std::ostream &	o,
		int	i = 0
	)		const

Definition at line 209 of file MdcxHit.cxx.

                                              {
  o << " hit" << i << " (" << Layer() << "," << WireNo() << ")  ";
  if (getMdcHit()) o << " phi "<< getMdcHit()->phi();
  o << "dd " << d() << " dde "
    << e() << " rt " << t() << endl;
}

process()

void MdcxHit::process ( )

protected

Definition at line 72 of file MdcxHit.cxx.

                      {
  assert(m_gm);
  assert(m_mdcCalibFunSvc);
  Identifier _id = (*_mdcDigi).identify();
  _layernumber=MdcID::layer(_id);
  _wirenumber=MdcID::wire(_id);
  _superlayer=(_layernumber)/4;
  _iTdc = _mdcDigi->getTimeChannel();
  _iAdc = _mdcDigi->getChargeChannel();
  _t = RawDataUtil::MdcTime(_iTdc);
  _q = RawDataUtil::MdcCharge(_iAdc);
  _T0Walk = m_mdcCalibFunSvc->getT0(_layernumber,_wirenumber)
    + m_mdcCalibFunSvc->getTimeWalk(_layernumber,_iAdc);
  // pointer to layer
  const MdcLayer* layerPtr= m_gm->Layer(_layernumber);
  _L = layerPtr->zLength(); 
  _x = layerPtr->xWire(_wirenumber);  
  _y = layerPtr->yWire(_wirenumber); 
  _r = sqrt(_x*_x + _y*_y);
  _s = layerPtr->stereo(); 
  _p = layerPtr->phiOffset() + _wirenumber*layerPtr->dPhi();
  //double deltaz = m_gm->zOffSet(); //yzhang to BESIII
  //double deltaz = 0.0;
  double tst = _s;
  double pw = atan2(_y, _x);
  _pw = pw; 
  _wx = -tst*sin(pw);
  _wy =  tst*cos(pw); 
  _wz = (tst*tst < 1.0) ? sqrt(1.0-tst*tst) : 0.0;
  _sp = sin(_p);
  _cp = cos(_p);
  _consterr = 0;  //zoujh FIXME
  _e = _cresol;
  _d = d();
  // note _v is a total cludge
  _v = (_t-_c0 > 0.0) ? _d/(_t-_c0) : 0.0013;//yzhang 2009-11-03 
  _xneg = _x + _d*_sp;
  _yneg = _y - _d*_cp;
  _xpos = _x - _d*_sp;
  _ypos = _y + _d*_cp;
  usedonhel = 0;
  _type = Constants::viewOfsLayer[_superlayer];
}

Referenced by MdcxHit().

pull()

float MdcxHit::pull

(

MdcxHel &

hel

)

const

Definition at line 174 of file MdcxHit.cxx.

                                      {
  // compute pulls for |hel|
  float doca=hel.Doca(*this); if(hel.Mode() == 0)doca=fabs(doca);
  return (d(hel.Doca_Zh(),hel.Doca_Tof(),hel.T0(),
    hel.Doca_Wamb(),hel.Doca_Eang())-doca)/e(doca);
  //return residual(hel)/e();
}

Referenced by MdcxAddHits::GetAssociates().

pw()

float MdcxHit::pw ( ) const

inline

Definition at line 78 of file MdcxHit.h.

{return _pw;} 

Referenced by MdcxAddHits::GetAssociates(), and process().

residual()

float MdcxHit::residual

(

MdcxHel &

hel

)

const

Definition at line 182 of file MdcxHit.cxx.

                                          {
  // compute residuals for |hel|
  float doca = hel.Doca(_wx, _wy, _wz, _x, _y);
  if(hel.Mode() == 0) doca = fabs(doca);
  doca += v()*hel.T0();
  return d(hel.Doca_Zh(),hel.Doca_Tof(),hel.T0(),
      hel.Doca_Wamb(),hel.Doca_Eang())-doca;
  //return d(hel) - doca;
}

SetConstErr()

void MdcxHit::SetConstErr ( int  i )

inline

Definition at line 85 of file MdcxHit.h.

{_consterr = i;}

Referenced by MdcxAddHits::GetAssociates().

setCountPropTime()

void MdcxHit::setCountPropTime ( bool  countPropTime )

static

Definition at line 64 of file MdcxHit.cxx.

                                                 {
  m_countPropTime = countPropTime;
}

Referenced by MdcxTrackFinder::initialize().

setMdcCalibFunSvc()

void MdcxHit::setMdcCalibFunSvc ( const MdcCalibFunSvc *  calibSvc )

static

Definition at line 60 of file MdcxHit.cxx.

                                                              {
  m_mdcCalibFunSvc = calibSvc;
}

Referenced by MdcxTrackFinder::initialize().

setMdcDetector()

void MdcxHit::setMdcDetector ( const MdcDetector *  gm )

static

Definition at line 68 of file MdcxHit.cxx.

                                                  {
  m_gm = gm;
}

Referenced by MdcxTrackFinder::beginRun().

SetUsedOnHel()

void MdcxHit::SetUsedOnHel ( int  i )

inline

Definition at line 102 of file MdcxHit.h.

{usedonhel=i;}

Referenced by MdcxAddHits::GetAssociates().

stereo()

float MdcxHit::stereo ( ) const

inline

Definition at line 80 of file MdcxHit.h.

{return _s;} // stereo angle of wire

SuperLayer()

int MdcxHit::SuperLayer ( ) const

inline

Definition at line 65 of file MdcxHit.h.

{return _superlayer;} //SuperLayer#

t()

float MdcxHit::t ( ) const

inline

Definition at line 68 of file MdcxHit.h.

{return _t;} // drift time

Referenced by printAll().

tcor()

float MdcxHit::tcor	(	float	zh = 0.0,
		float	tof = 0.0,
		float	tzero = 0.0
	)		const

Definition at line 116 of file MdcxHit.cxx.

                                                        {
  //drift time 
  double tprop = 0.;
  if (m_countPropTime){ tprop = m_mdcCalibFunSvc->getTprop(_layernumber,z*10); }
  double dt = (_t - _T0Walk -_c0 - tprop - tof - tzero);
/*
  std::cout<<"("<<_layernumber<<","<<_wirenumber<<") dt "<<dt
    <<" rt "<<_t
    <<" t0walk "<< _T0Walk
    <<" z "<<z
    <<" _c0 "<<_c0
    <<" tzero "<<tzero
    <<" tof "<<tof
    <<" tprop "<<tprop
    <<std::endl;
   */ 
  return dt;
}

Referenced by d().

type()

int MdcxHit::type ( ) const

inline

Definition at line 82 of file MdcxHit.h.

{return _type;} // flags |MdcxHit|

v()

float MdcxHit::v ( ) const

inline

Definition at line 81 of file MdcxHit.h.

{return _v;} // drift velocity

Referenced by MdcxHel::derivatives(), and residual().

WireNo()

int MdcxHit::WireNo ( ) const

inline

Definition at line 66 of file MdcxHit.h.

{return _wirenumber;} // Wire#

Referenced by print(), printAll(), and MdcxFindSegs::process().

wx()

float MdcxHit::wx ( ) const

inline

Definition at line 75 of file MdcxHit.h.

{return _wx;} 

Referenced by MdcxHel::Doca().

wy()

float MdcxHit::wy ( ) const

inline

Definition at line 76 of file MdcxHit.h.

{return _wy;} 

Referenced by MdcxHel::Doca().

wz()

float MdcxHit::wz ( ) const

inline

Definition at line 77 of file MdcxHit.h.

{return _wz;} 

Referenced by MdcxHel::Doca().

x()

float MdcxHit::x ( ) const

inline

Definition at line 69 of file MdcxHit.h.

{return _x;} // x of wire

Referenced by MdcxHel::Doca(), and MdcxFindSegs::trial().

xneg()

float MdcxHit::xneg ( ) const

inline

Definition at line 73 of file MdcxHit.h.

{return _xneg;} 

Referenced by MdcxFindSegs::trial().

xpos()

float MdcxHit::xpos ( ) const

inline

Definition at line 71 of file MdcxHit.h.

{return _xpos;}  

Referenced by MdcxFindSegs::trial().

y()

float MdcxHit::y ( ) const

inline

Definition at line 70 of file MdcxHit.h.

{return _y;} // y of wire

Referenced by MdcxHel::Doca(), and MdcxFindSegs::trial().

yneg()

float MdcxHit::yneg ( ) const

inline

Definition at line 74 of file MdcxHit.h.

{return _yneg;} 

Referenced by MdcxFindSegs::trial().

ypos()

float MdcxHit::ypos ( ) const

inline

Definition at line 72 of file MdcxHit.h.

{return _ypos;} 

Referenced by MdcxFindSegs::trial().

Member Data Documentation

_c0

float MdcxHit::_c0

protected

Definition at line 135 of file MdcxHit.h.

Referenced by process(), and tcor().

_consterr

int MdcxHit::_consterr

protected

Definition at line 133 of file MdcxHit.h.

Referenced by process(), and SetConstErr().

_cp

float MdcxHit::_cp

protected

Definition at line 128 of file MdcxHit.h.

Referenced by process().

_cresol

float MdcxHit::_cresol

protected

Definition at line 136 of file MdcxHit.h.

Referenced by process().

_d

float MdcxHit::_d

protected

Definition at line 123 of file MdcxHit.h.

Referenced by process().

_e

float MdcxHit::_e

protected

Definition at line 126 of file MdcxHit.h.

Referenced by process().

_iAdc

unsigned MdcxHit::_iAdc

protected

Definition at line 116 of file MdcxHit.h.

Referenced by process().

_iTdc

unsigned MdcxHit::_iTdc

protected

Definition at line 115 of file MdcxHit.h.

Referenced by process().

_L

float MdcxHit::_L

protected

Definition at line 122 of file MdcxHit.h.

Referenced by process().

_layernumber

unsigned MdcxHit::_layernumber

protected

Definition at line 111 of file MdcxHit.h.

Referenced by d(), Layer(), process(), and tcor().

_mdcDigi

const MdcDigi* MdcxHit::_mdcDigi

protected

Definition at line 132 of file MdcxHit.h.

Referenced by getDigi(), and process().

_mdcHit

const MdcHit* MdcxHit::_mdcHit

protected

Definition at line 131 of file MdcxHit.h.

Referenced by getMdcHit().

_p

float MdcxHit::_p

protected

Definition at line 128 of file MdcxHit.h.

Referenced by process().

_pw

double MdcxHit::_pw

protected

Definition at line 129 of file MdcxHit.h.

Referenced by phiMid(), process(), and pw().

_q

float MdcxHit::_q

protected

Definition at line 114 of file MdcxHit.h.

Referenced by process().

_r

float MdcxHit::_r

protected

Definition at line 120 of file MdcxHit.h.

Referenced by process().

_s

float MdcxHit::_s

protected

Definition at line 121 of file MdcxHit.h.

Referenced by process(), and stereo().

_sp

float MdcxHit::_sp

protected

Definition at line 128 of file MdcxHit.h.

Referenced by process().

_superlayer

unsigned MdcxHit::_superlayer

protected

Definition at line 112 of file MdcxHit.h.

Referenced by process(), and SuperLayer().

_t

float MdcxHit::_t

protected

Definition at line 113 of file MdcxHit.h.

Referenced by process(), t(), and tcor().

_T0Walk

float MdcxHit::_T0Walk

protected

Definition at line 117 of file MdcxHit.h.

Referenced by process(), and tcor().

_type

int MdcxHit::_type

protected

Definition at line 125 of file MdcxHit.h.

Referenced by process(), and type().

_v

float MdcxHit::_v

protected

Definition at line 124 of file MdcxHit.h.

Referenced by process(), and v().

_wirenumber

unsigned MdcxHit::_wirenumber

protected

Definition at line 110 of file MdcxHit.h.

Referenced by d(), process(), and WireNo().

_wx

double MdcxHit::_wx

protected

Definition at line 130 of file MdcxHit.h.

Referenced by process(), residual(), and wx().

_wy

double MdcxHit::_wy

protected

Definition at line 130 of file MdcxHit.h.

Referenced by process(), residual(), and wy().

_wz

double MdcxHit::_wz

protected

Definition at line 130 of file MdcxHit.h.

Referenced by process(), residual(), and wz().

_x

float MdcxHit::_x

protected

Definition at line 118 of file MdcxHit.h.

Referenced by process(), residual(), and x().

_xneg

float MdcxHit::_xneg

protected

Definition at line 127 of file MdcxHit.h.

Referenced by process(), and xneg().

_xpos

float MdcxHit::_xpos

protected

Definition at line 127 of file MdcxHit.h.

Referenced by process(), and xpos().

_y

float MdcxHit::_y

protected

Definition at line 119 of file MdcxHit.h.

Referenced by process(), residual(), and y().

_yneg

float MdcxHit::_yneg

protected

Definition at line 127 of file MdcxHit.h.

Referenced by process(), and yneg().

_ypos

float MdcxHit::_ypos

protected

Definition at line 127 of file MdcxHit.h.

Referenced by process(), and ypos().

usedonhel

int MdcxHit::usedonhel

protected

Definition at line 134 of file MdcxHit.h.

Referenced by GetUsedOnHel(), process(), and SetUsedOnHel().

---

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

• MdcxHit.h
• MdcxHit.cxx