MdcHit Class Reference

#include <MdcHit.h>

Inheritance diagram for MdcHit:

Public Member Functions
	MdcHit (const MdcDigi *&digi, const MdcDetector *&det)
	MdcHit (const MdcHit &)
virtual	~MdcHit ()
MdcHit &	operator= (const MdcHit &)
bool	operator== (const MdcHit &) const
const MdcDigi *	digi () const
const MdcLayer *	layer () const
const MdcSWire *	wire () const
const MdcDetector *	geom () const
const MdcCalibFunSvc *	calibSvc () const
Identifier	mdcId () const
unsigned	layernumber () const
unsigned	wirenumber () const
unsigned	tdcIndex () const
unsigned	adcIndex () const
double	charge () const
double	rawTime () const
double	driftTime (double tof, double z) const
double	driftDist (double, int, double, double, double) const
double	driftDist (double bunchTime, int ambig) const
double	sigma (double, int, double, double, double) const
double	sigma (double driftdist, int ambig=0) const
int	whichView () const
const Trajectory *	hitTraj () const
TrkEnums::TrkViewInfo	whatView () const
double	phi () const
double	x () const
double	y () const
const double	phi (double z) const
double	x (double z) const
double	y (double z) const
double	rMid () const
double	zlen () const
unsigned	status () const
void	setCalibSvc (const MdcCalibFunSvc *calibSvc)
void	setCountPropTime (const bool count)
void	setCosmicFit (const bool cosmicfit)
bool	isCosmicFit () const
void	print (std::ostream &o) const
void	printAll (std::ostream &o) const
Public Member Functions inherited from TrkFundHit
	TrkFundHit ()
virtual	~TrkFundHit ()
int	nUsedHits () const
bool	usedHit (void) const
std::pair< TrkFundHit::hot_iterator, TrkFundHit::hot_iterator >	getUsedHits () const
TrkFundHit::hot_iterator	begin () const
TrkFundHit::hot_iterator	end () const
bool	usedOnTrack (const TrkRecoTrk *t) const
const TrkHitOnTrk *	getHitOnTrack (const TrkRecoTrk *trk) const
const TrkHitOnTrk *	setUsedHit (const TrkHitOnTrk *hit)
const TrkHitOnTrk *	setUnusedHit (const TrkHitOnTrk *hit)

Friends
class	MdcMakeHits

Additional Inherited Members
Public Types inherited from TrkFundHit
typedef TrkHitOnTrkIter< TrkFundHit >	hot_iterator
Protected Types inherited from TrkFundHit
typedef std::vector< constTrkHitOnTrk * >::iterator	iterator_implementation
typedef const TrkHitOnTrk	iterator_value_type
Protected Member Functions inherited from TrkFundHit
TrkFundHit &	operator= (const TrkFundHit &)
Protected Attributes inherited from TrkFundHit
std::vector< const TrkHitOnTrk * >	_hitList

Detailed Description

Definition at line 44 of file MdcHit.h.

Constructor & Destructor Documentation

MdcHit() [1/2]

MdcHit::MdcHit	(	const MdcDigi *&	digi,
		const MdcDetector *&	det )

Definition at line 33 of file MdcHit.cxx.

                                                             :
  TrkFundHit(), _digiPtr(aDigi), _geomPtr(det)
{ 
  _digiPtr  = aDigi;
  _geomPtr  = det;
  assert( _digiPtr != NULL);
  assert( _geomPtr != NULL);
  _id       = aDigi->identify();
  _layerPtr = det->Layer(_id);
  _wirePtr  = det->Wire(_id);
  assert( _layerPtr != NULL);
  assert( _wirePtr  != NULL);
  m_mdcCalibFunSvc = NULL;
  _layer    = MdcID::layer(_id); 
  _wire     = MdcID::wire (_id); 
  _iTdc     = _digiPtr->getTimeChannel();
  _iAdc     = _digiPtr->getChargeChannel();
  _rawTime  = RawDataUtil::MdcTime(_iTdc);
  _charge   = RawDataUtil::MdcCharge(_iAdc);
  _rmid     = _wirePtr->rMid();
  _zlen     = _layerPtr->zLength();   
  _phi      = _wirePtr->phi();
  _cosphi   = cos(_phi);
  _sinphi   = sin(_phi);
  _status   = 0;
}

Referenced by MdcHit(), and MdcHit().

MdcHit() [2/2]

MdcHit::MdcHit

(

const MdcHit &

other

)

Definition at line 60 of file MdcHit.cxx.

                                  :
  TrkFundHit() , _digiPtr(other._digiPtr), _geomPtr(other._geomPtr),
  _layerPtr(other._layerPtr), _wirePtr(other._wirePtr), _id(other._id),
  _layer(other._layer), _wire(other._wire),
  _iTdc(other._iTdc), _iAdc(other._iAdc),
  _rawTime(other._rawTime), _charge(other._charge),
  _rmid(other._rmid), _zlen(other._zlen), 
  _phi(other._phi), _cosphi(other._cosphi), _sinphi(other._sinphi),
  _status(other._status), _T0Walk(other._T0Walk) 
{
  m_mdcCalibFunSvc = other.m_mdcCalibFunSvc;
}

~MdcHit()

MdcHit::~MdcHit ( )

virtual

Definition at line 102 of file MdcHit.cxx.

                {
  // This is ugly and inefficient.  This, along with the rest of 
  //  the hitList mess, should be cleaned up by handling the 
  //  association in an external map
 
  // Not written as a loop because removeHit() modifies TrkFundHit::_hitList
  short count = 0;
  while (nUsedHits() > count) {
    bool removed = _hitList[count]->parentTrack()->hits()->removeHit(this);
    if (!removed) count++;
  }
}

Member Function Documentation

adcIndex()

unsigned MdcHit::adcIndex ( ) const

inline

Definition at line 64 of file MdcHit.h.

{ return _iAdc; }

Referenced by MdcTrack::storeTrack().

calibSvc()

const MdcCalibFunSvc * MdcHit::calibSvc ( ) const

inline

Definition at line 59 of file MdcHit.h.

{ return m_mdcCalibFunSvc; }

Referenced by setCalibSvc().

charge()

double MdcHit::charge ( ) const

inline

Definition at line 65 of file MdcHit.h.

{ return _charge; } 

Referenced by MdcHitOnTrack::charge(), and MdcTrkRecon::fillTrackList().

digi()

const MdcDigi * MdcHit::digi ( ) const

inline

Definition at line 55 of file MdcHit.h.

{ return _digiPtr;  }

Referenced by MdcTrackList::pickHits(), print(), and MdcTrack::storeTrack().

driftDist() [1/2]

double MdcHit::driftDist	(	double	bunchTime,
		int	ambig ) const

Definition at line 178 of file MdcHit.cxx.

                                                    {
  return driftDist(bunchTime+crudeTof(), ambig, 0., 0., 0. );
}

driftDist() [2/2]

double MdcHit::driftDist	(	double	tof,
		int	ambig,
		double	entranceAngle,
		double	,
		double	z ) const

Definition at line 156 of file MdcHit.cxx.

{ 
  double driftD;
  //drift time ns, layer id begin with 0, entrance angle rads,
  //lr ambig: wire ambig 1,-1,0 -> Calib 0,1,2
  int lrCalib=2;
  if (ambig==1) lrCalib = 0;
  else if (ambig==-1) lrCalib = 1;
 
  // tof in s, driftDist in cm, dirftTime in ns
  if (fabs(z)>150. || fabs(driftTime(tof,z))>1500.){
    return 9999.;
  }
  driftD = 0.1 * m_mdcCalibFunSvc->driftTimeToDist(driftTime(tof,z),_layer,_wire,lrCalib,entranceAngle);//to cm
  //std::cout<<"driftDist "<<"("<<_layer <<","<<_wire <<") dd "<<driftD<<" dt "<<driftTime(tof,z) <<" lr "<<lrCalib <<" eAng "<<entranceAngle <<" tof "<<tof*1.e9<<" z "<<z <<" t0walk "<<_T0Walk<<" rawT "<<_rawTime <<" tprop "<< _rawTime - driftTime(tof,z)- _T0Walk-1.e9*tof<<std::endl;
 
  if (abs(driftD)<Constants::epsilon) driftD = 0.00001;
  return driftD;
}

Referenced by MdcSeg::addHits(), driftDist(), MdcxCosmicSewer::execute(), Hough2D::fit(), Hough3D::fit(), MdcTrackList::pickHits(), printAll(), and MdcSegInfoSterO::zPosition().

driftTime()

double MdcHit::driftTime	(	double	tof,
		double	z ) const

Definition at line 142 of file MdcHit.cxx.

                                            {
  // t prop
  double tprop = 0.;
  if (m_countPropTime){ tprop = m_mdcCalibFunSvc->getTprop(_layer,z*10.); }
 
  //tof in s, driftTime in ns, _T0Walk in ns
  double driftT = _rawTime - _T0Walk -1.e9*tof - tprop;
 
  //if(driftT >10000)
  //std::cout<< "driftTime ("<<_layer<<","<<_wire<<") dt "<<driftT<<" rawTime "<<_rawTime<<" z "<<z<<" tprop "<<tprop <<" t0walk "<<_T0Walk<<" tof "<<1.e9*tof<< std::endl;//yzhang debug
  return driftT; 
}  

Referenced by driftDist(), MdcTrackList::dropMultiHotInLayer(), MdcTrkRecon::fillTrackList(), Hough2D::fit(), Hough3D::fit(), printAll(), and MdcTrack::storeTrack().

geom()

const MdcDetector * MdcHit::geom ( ) const

inline

Definition at line 58 of file MdcHit.h.

{ return _geomPtr;  }

hitTraj()

const Trajectory * MdcHit::hitTraj

(

)

const

Definition at line 231 of file MdcHit.cxx.

                      {
  return layer()->makeHitTrajInGlobalCoords(wirenumber(),0.0);
}

isCosmicFit()

bool MdcHit::isCosmicFit ( ) const

inline

Definition at line 88 of file MdcHit.h.

{ return m_cosmicFit; }

layer()

const MdcLayer * MdcHit::layer ( ) const

inline

Definition at line 56 of file MdcHit.h.

{ return _layerPtr; }

Referenced by hitTraj(), MdcHitOnTrack::layer(), MdcHitOnTrack::MdcHitOnTrack(), MdcTrackList::pickHits(), and MdcSegGrouper::transferHits().

layernumber()

unsigned MdcHit::layernumber ( ) const

inline

Definition at line 61 of file MdcHit.h.

{ return _layer; }

Referenced by MdcHitMapGuts::addHit(), MdcSeg::addHits(), MdcSegGrouper::calcParByHits(), MdcTrackList::dropMultiHotInLayer(), MdcTrkRecon::execute(), MdcTrkRecon::fillSegList(), MdcTrkRecon::fillTrackList(), MdcHitOnTrack::layernumber(), MdcTrackList::pickHits(), MdcTrack::storeTrack(), MdcSeg::testCombSeg(), MdcSeg::testCombSegAmbig(), MdcSeg::testCombSegPhi(), MdcSeg::testCombSegPt(), MdcSeg::testCombSegTheta(), and MdcSegInfoSterO::zPosition().

mdcId()

Identifier MdcHit::mdcId ( ) const

inline

Definition at line 60 of file MdcHit.h.

{ return _id; }

operator=()

MdcHit & MdcHit::operator=

(

const MdcHit &

other

)

Definition at line 74 of file MdcHit.cxx.

                                       {
  if(&other != this){
    _digiPtr     = other._digiPtr;
    _geomPtr     = other._geomPtr;
    _layerPtr    = other._layerPtr;
    _wirePtr     = other._wirePtr;
    m_mdcCalibFunSvc = other.m_mdcCalibFunSvc;
    _id          = other._id;
    _layer       = other._layer;
    _wire        = other._wire;
    _iTdc        = other._iTdc;
    _iAdc        = other._iAdc;
    _rawTime     = other._rawTime;
    _charge      = other._charge;
    _digiPtr     = other._digiPtr;
    _rmid        = other._rmid;
    _zlen        = other._zlen;
    _phi         = other._phi;
    _cosphi      = other._cosphi;
    _sinphi      = other._sinphi;
    _status      = other._status;
    _T0Walk      = other._T0Walk;
  }
  return *this;
}

operator==()

bool MdcHit::operator==

(

const MdcHit &

rhs

)

const

Definition at line 116 of file MdcHit.cxx.

                                            {
  return (this == &rhs);
}

phi() [1/2]

double MdcHit::phi ( ) const

inline

Definition at line 75 of file MdcHit.h.

{ return _phi;}       // phi at chamber center

Referenced by MdcSeg::addHits(), MdcHitOnTrack::entranceAngle(), MdcxCosmicSewer::execute(), MdcTrackList::pickHits(), and MdcxHit::printAll().

phi() [2/2]

const double MdcHit::phi ( double z ) const

inline

Definition at line 78 of file MdcHit.h.

{return wire()->phiDC(z);}  // phi at z 

print()

void MdcHit::print

(

std::ostream &

o

)

const

Definition at line 121 of file MdcHit.cxx.

                                {
  o << "(" << _layer << ","<< _wire <<","<<digi()->getTrackIndex()<<")";
}

Referenced by MdcSegInfoSterO::calcStereo(), operator<<(), MdcTrackList::pickHits(), MdcSeg::plotSeg(), MdcSeg::plotSegAll(), MdcSegData::poisonHits(), and MdcSegInfoSterO::zPosition().

printAll()

void MdcHit::printAll ( std::ostream & o ) const

virtual

Reimplemented from TrkFundHit.

Definition at line 126 of file MdcHit.cxx.

                                   {
  o << "Hit:(" <<  _layer << "," << _wire 
    << ") x:" << x() << " y:" << y()
    << " rt:" << rawTime()
    << " dt:" <<driftTime(0.0,0)
    << " dd:" <<driftDist(0.0,1)
    << endl;
}

rawTime()

double MdcHit::rawTime ( ) const

inline

Definition at line 66 of file MdcHit.h.

{ return _rawTime; }

Referenced by Hough2D::fit(), Hough3D::fit(), printAll(), MdcHitOnTrack::rawTime(), and MdcHitOnTrack::setT0().

rMid()

double MdcHit::rMid ( ) const

inline

Definition at line 81 of file MdcHit.h.

{ return _rmid; } // R at chamber center

setCalibSvc()

void MdcHit::setCalibSvc

(

const MdcCalibFunSvc *

calibSvc

)

Definition at line 136 of file MdcHit.cxx.

                                                  { 
  m_mdcCalibFunSvc = calibSvc;  
  _T0Walk = m_mdcCalibFunSvc->getT0(_layer,_wire) +  m_mdcCalibFunSvc->getTimeWalk(_layer, _iAdc);
}

Referenced by MdcMakeHits::execute(), MdcxCosmicSewer::execute(), Hough2D::fit(), Hough3D::fit(), and MdcxCosmicSewer::MdcxHitsToHots().

setCosmicFit()

void MdcHit::setCosmicFit ( const bool cosmicfit )

inline

Definition at line 87 of file MdcHit.h.

{ m_cosmicFit = cosmicfit; }

Referenced by MdcMakeHits::execute(), MdcxCosmicSewer::execute(), and MdcxCosmicSewer::MdcxHitsToHots().

setCountPropTime()

void MdcHit::setCountPropTime ( const bool count )

inline

Definition at line 86 of file MdcHit.h.

{ m_countPropTime = count;}

Referenced by MdcMakeHits::execute(), MdcxCosmicSewer::execute(), Hough2D::fit(), Hough3D::fit(), and MdcxCosmicSewer::MdcxHitsToHots().

sigma() [1/2]

double MdcHit::sigma	(	double	driftdist,
		int	ambig = 0 ) const

Definition at line 226 of file MdcHit.cxx.

                                               {
  return sigma(driftdist, ambig, 0., 0., 0.);// cm
}

sigma() [2/2]

double MdcHit::sigma	(	double	driftdist,
		int	ambig,
		double	entranceAngle,
		double	dipAngle,
		double	z ) const

Definition at line 184 of file MdcHit.cxx.

                                     {
  double sig = 9999.;
 
#ifdef MDCPATREC_RESLAYER
  if (_layer == m_resLayer){
    //give a huge sigma to skip this layer when fit track
    return 9999.;
  }
#endif
  //if(fabs(z)>150. || driftdist==9999. ) {
  //  return 9999.;
  //} 
  if ( m_mdcCalibFunSvc ) {
    //layid begin with 0, entrance angle in rads,
    //distance: cm ->Calib mm //z: cm -> Calib mm
    //Q: MonteCalo eV -> Calib fC FIXME
    //lr ambig: wire ambig 1,-1,0 -> Calib 0,1,2 
    int lrCalib=2;
    if (ambig==1) lrCalib = 0;
    else if (ambig==-1) lrCalib = 1;
 
    //std::cout<<"layer "<<_layer<< " lrCalib "<< lrCalib<< " driftdist "<< driftdist << " eAngle "<<entranceAngle
    //<<" tan(dipAngle) "<<tan(dipAngle)<< " z "<<z <<" iAdc "<<_iAdc<<std::endl;
    sig = 0.1 * m_mdcCalibFunSvc->getSigma(_layer,lrCalib,driftdist*10.,entranceAngle,
    tan(dipAngle),z*10.,_iAdc); //Calib special resolution mm -> cm
 
    if(fabs(sig)<Constants::epsilon){
      sig = 999.;
    }
    //if(sig<=0){
    //std::cout<<__FILE__<<"   "<<__LINE__ <<" sigma "<<sig
    //<<" layer "<<_layer <<" lrCalib "<<lrCalib <<" driftdist "<<driftdist*10
    //<<" eAngle "<<entranceAngle<<" dipAngle "<<(dipAngle) <<" tanl "<<tan(dipAngle)
    //<<" z "<<z <<" iAdc "<<_iAdc <<std::endl;
    //}
  }
  //if (_layer<4) sig /=5.;//yzhang TEST 2011-05-27 
  return sig;
}

Referenced by MdcSeg::addHits(), MdcxCosmicSewer::execute(), MdcTrackList::pickHits(), and MdcSegInfoSterO::zPosition().

status()

unsigned MdcHit::status ( ) const

inline

Definition at line 83 of file MdcHit.h.

{ return _status;}// status

Referenced by MdcRecoHitOnTrack::status().

tdcIndex()

unsigned MdcHit::tdcIndex ( ) const

inline

Definition at line 63 of file MdcHit.h.

{ return _iTdc; }

Referenced by MdcTrack::storeTrack(), and MdcRecoHitOnTrack::tdcIndex().

whatView()

TrkEnums::TrkViewInfo MdcHit::whatView ( ) const

inlinevirtual

Implements TrkFundHit.

Definition at line 74 of file MdcHit.h.

{ return whichView()==0?TrkEnums::xyView:TrkEnums::bothView; }

Referenced by MdcHitOnTrack::whatView().

whichView()

int MdcHit::whichView ( ) const

inline

Definition at line 72 of file MdcHit.h.

{ return _layerPtr->view(); } 

Referenced by whatView(), and MdcHitOnTrack::whichView().

wire()

const MdcSWire * MdcHit::wire ( ) const

inline

Definition at line 57 of file MdcHit.h.

{ return _wirePtr;  }

Referenced by phi(), x(), y(), and MdcSegInfoSterO::zPosition().

wirenumber()

unsigned MdcHit::wirenumber ( ) const

inline

Definition at line 62 of file MdcHit.h.

{ return _wire; }

Referenced by MdcHitMapGuts::addHit(), MdcSeg::addHits(), MdcSegGrouper::calcParByHits(), MdcTrackList::dropMultiHotInLayer(), MdcTrkRecon::execute(), MdcTrkRecon::fillSegList(), MdcTrkRecon::fillTrackList(), hitTraj(), MdcTrackList::pickHits(), MdcTrack::storeTrack(), MdcSeg::testCombSeg(), MdcSeg::testCombSegAmbig(), MdcSeg::testCombSegPhi(), MdcSeg::testCombSegPt(), MdcSeg::testCombSegTheta(), MdcHitOnTrack::wire(), and MdcSegInfoSterO::zPosition().

x() [1/2]

double MdcHit::x ( ) const

inline

Definition at line 76 of file MdcHit.h.

{ return _rmid*_cosphi; } // x at chamber center

Referenced by MdcTrackList::pickHits(), and printAll().

x() [2/2]

double MdcHit::x ( double z ) const

inline

Definition at line 79 of file MdcHit.h.

{return wire()->xWireDC(z);}// x at global z

y() [1/2]

double MdcHit::y ( ) const

inline

Definition at line 77 of file MdcHit.h.

{ return _rmid*_sinphi; } // y at chamber center

Referenced by MdcTrackList::pickHits(), printAll(), and MdcSegGrouper::transferHits().

y() [2/2]

double MdcHit::y ( double z ) const

inline

Definition at line 80 of file MdcHit.h.

{return wire()->yWireDC(z);}// y at global z

zlen()

double MdcHit::zlen ( ) const

inline

Definition at line 82 of file MdcHit.h.

{ return _zlen; } // chamber extent in z

Friends And Related Symbol Documentation

MdcMakeHits

friend class MdcMakeHits

friend

Definition at line 103 of file MdcHit.h.

---

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

• MdcHit.h
• MdcHit.cxx