#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.

64{ return _iAdc; }

Referenced by MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ calibSvc()

const MdcCalibFunSvc * MdcHit::calibSvc ( ) const

inline

Definition at line 59 of file MdcHit.h.

59{ return m_mdcCalibFunSvc; }

Referenced by setCalibSvc().

◆ charge()

double MdcHit::charge ( ) const

inline

Definition at line 65 of file MdcHit.h.

65{ return _charge; }

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

◆ digi()

const MdcDigi * MdcHit::digi ( ) const

inline

Definition at line 55 of file MdcHit.h.

55{ return _digiPtr; }

Referenced by MdcTrackList::pickHits(), print(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ 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(), CgemMdcFitAlg::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(), CgemMdcFitAlg::fit(), printAll(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ geom()

const MdcDetector * MdcHit::geom ( ) const

inline

Definition at line 58 of file MdcHit.h.

58{ return _geomPtr; }

◆ hitTraj()

const Trajectory * MdcHit::hitTraj ( ) const

Definition at line 231 of file MdcHit.cxx.

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

Referenced by HoughHit::residual().

◆ isCosmicFit()

bool MdcHit::isCosmicFit ( ) const

inline

Definition at line 88 of file MdcHit.h.

88{ return m_cosmicFit; }

◆ layer()

const MdcLayer * MdcHit::layer ( ) const

inline

Definition at line 56 of file MdcHit.h.

56{ return _layerPtr; }

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

◆ layernumber()

unsigned MdcHit::layernumber ( ) const

inline

◆ mdcId()

Identifier MdcHit::mdcId ( ) const

inline

Definition at line 60 of file MdcHit.h.

60{ 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.

75{ 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.

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

MdcHit::wire

const MdcSWire * wire() const

Definition MdcHit.h:57

MdcSWire::phiDC

double phiDC(double z) const

Definition MdcSWire.h:60

◆ 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.

66{ return _rawTime; }

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

◆ rMid()

double MdcHit::rMid ( ) const

inline

Definition at line 81 of file MdcHit.h.

81{ 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(), CgemMdcFitAlg::fit(), HoughTrack::fitCircle(), HoughTrack::fitHelix(), HoughTrack::fitHelix(), and MdcxCosmicSewer::MdcxHitsToHots().

◆ setCosmicFit()

void MdcHit::setCosmicFit ( const bool cosmicfit )

inline

Definition at line 87 of file MdcHit.h.

87{ 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.

86{ m_countPropTime = count;}

Referenced by MdcMakeHits::execute(), MdcxCosmicSewer::execute(), CgemMdcFitAlg::fit(), HoughTrack::fitCircle(), HoughTrack::fitHelix(), HoughTrack::fitHelix(), 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(), sigma(), and MdcSegInfoSterO::zPosition().

◆ status()

unsigned MdcHit::status ( ) const

inline

Definition at line 83 of file MdcHit.h.

83{ return _status;}// status

Referenced by MdcRecoHitOnTrack::status().

◆ tdcIndex()

unsigned MdcHit::tdcIndex ( ) const

inline

Definition at line 63 of file MdcHit.h.

63{ return _iTdc; }

Referenced by MdcTrack::storeTrack(), HoughFinder::storeTracks(), MdcRecoHitOnTrack::tdcIndex(), and CgemMdcFitAlg::updateTracks().

◆ whatView()

TrkEnums::TrkViewInfo MdcHit::whatView ( ) const

inlinevirtual

Implements TrkFundHit.

Definition at line 74 of file MdcHit.h.

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

MdcHit::whichView

int whichView() const

Definition MdcHit.h:72

TrkEnums::xyView

@ xyView

Definition TrkEnums.h:22

TrkEnums::bothView

@ bothView

Definition TrkEnums.h:22

Referenced by MdcHitOnTrack::whatView().

◆ whichView()

int MdcHit::whichView ( ) const

inline

Definition at line 72 of file MdcHit.h.

72{ return _layerPtr->view(); }

MdcLayer::view

int view(void) const

Definition MdcLayer.h:28

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

◆ wire()

const MdcSWire * MdcHit::wire ( ) const

inline

Definition at line 57 of file MdcHit.h.

57{ return _wirePtr; }

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

◆ wirenumber()

unsigned MdcHit::wirenumber ( ) const

inline

◆ x() [1/2]

double MdcHit::x ( ) const

inline

Definition at line 76 of file MdcHit.h.

76{ 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.

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

MdcSWire::xWireDC

double xWireDC(double z) const

Definition MdcSWire.h:61

◆ y() [1/2]

double MdcHit::y ( ) const

inline

Definition at line 77 of file MdcHit.h.

77{ 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.

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

MdcSWire::yWireDC

double yWireDC(double z) const

Definition MdcSWire.h:62

◆ zlen()

double MdcHit::zlen ( ) const

inline

Definition at line 82 of file MdcHit.h.

82{ 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:

Public Member Functions

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MdcHit() [1/2]

◆ MdcHit() [2/2]

◆ ~MdcHit()

Member Function Documentation

◆ adcIndex()

◆ calibSvc()

◆ charge()

◆ digi()

◆ driftDist() [1/2]

◆ driftDist() [2/2]

◆ driftTime()

◆ geom()

◆ hitTraj()

◆ isCosmicFit()

◆ layer()

◆ layernumber()

◆ mdcId()

◆ operator=()

◆ operator==()

◆ phi() [1/2]

◆ phi() [2/2]

◆ print()

◆ printAll()

◆ rawTime()

◆ rMid()

◆ setCalibSvc()

◆ setCosmicFit()

◆ setCountPropTime()

◆ sigma() [1/2]

◆ sigma() [2/2]

◆ status()

◆ tdcIndex()

◆ whatView()

◆ whichView()

◆ wire()

◆ wirenumber()

◆ x() [1/2]

◆ x() [2/2]

◆ y() [1/2]

◆ y() [2/2]

◆ zlen()

Friends And Related Symbol Documentation

◆ MdcMakeHits