#include <HoughHit.h>

Inheritance diagram for HoughHit:

Public Member Functions
	HoughHit ()
	HoughHit (const MdcDigi *const digi)
	HoughHit (const HoughHit &other)
HoughHit &	operator= (const HoughHit &other)
void	setTruthInfo (const MdcMcHit *&mcHit)
void	conformalTrans (double x, double y, double r)
double	getConformal_u (double, double, double)
double	getConformal_v (double, double, double)
double	getConformal_r (double, double, double)
void	makeCir (int n, double phi_begin, double phi_last, double r)
CFCir	getCir (int i) const
HoughHitType	type () const
detectorType	getDetectorType () const
const MdcDigi *	digi () const
double	getBunchTime () const
HepPoint3D	getMidPoint () const
HepPoint3D	getEastPoint () const
HepPoint3D	getWestPoint () const
double	getMidX () const
double	getMidY () const
int	getLayerId () const
int	getWireId () const
int	getSlayerType () const
const MdcLayer *	layer () const
const MdcSWire *	wire () const
double	getCharge () const
double	getDriftTime () const
double	getDriftDist () const
double	getU () const
double	getV () const
double	getR () const
double	getDeltaD () const
double	getFltLen () const
double	driftTime () const
double	driftTime (double tof, double z) const
double	calDriftDist (double, int, double, double, double) const
double	calDriftDist (double bunchTime, int ambig) const
int	slayerType (int layer)
double	getXTruth () const
double	getYTruth () const
double	getZTruth () const
double	getDriftDistTruth () const
int	getIdTruth () const
int	getLrTruth () const
double	getUTruth () const
double	getVTruth () const
double	getRTruth () const
HepPoint3D	getPointTruth () const
void	setDeltaD (double d)
void	setFltLen (double flt)
int	getCirList () const
int	getStyle () const
void	setCirList (int cir)
void	setStyle (int sty)
void	print () const
void	printAll () const
void	printTruth () const

Static Public Member Functions
static void	setMdcCalibFunSvc (const IMdcCalibFunSvc *calibSvc)
static void	setMdcGeomSvc (IMdcGeomSvc *geomSvc)
static void	setBunchTime (double t0)

Static Public Attributes
static int	_npart = 360

Detailed Description

Definition at line 24 of file HoughHit.h.

Constructor & Destructor Documentation

◆ HoughHit() [1/3]

HoughHit::HoughHit ( )

Definition at line 18 of file HoughHit.cxx.

                   {
  _det      = Global::m_gm;
  _digiPtr  = NULL;
  _layerPtr = NULL;
  _wirePtr  = NULL;
}

Referenced by HoughHit(), HoughRecHit::HoughRecHit(), HoughRecHit::HoughRecHit(), HoughRecHit::HoughRecHit(), HoughRecHit::HoughRecHit(), and operator=().

◆ HoughHit() [2/3]

HoughHit::HoughHit ( const MdcDigi *const digi )

Definition at line 24 of file HoughHit.cxx.

                                               {
  _det        = Global::m_gm;
  _digiPtr    = aDigi;
  _id         = aDigi->identify();
  _layer      = MdcID::layer( _id );
  _wire       = MdcID::wire( _id );
  _layerPtr   = _det->Layer( _id );
  _wirePtr    = _det->Wire( _id );
  _rawTime    = RawDataUtil::MdcTime( _digiPtr->getTimeChannel() );
  _charge     = _digiPtr->getChargeChannel();
  _driftTime  = driftTime();
  _driftDist  = calDriftDist( _bunchTime, 0 );
  _slayerType = slayerType( _layer );
  _lr         = 0;
  _rmid       = _wirePtr->rMid();
 
  assert( _mdcGeomSvc != NULL );
  const MdcGeoWire* wire = _mdcGeomSvc->Wire( _layer, _wire );
  assert( wire != NULL );
  HepPoint3D eastP = wire->Backward() / 10.;
  HepPoint3D westP = wire->Forward() / 10.;
  _eastPoint       = eastP;
  _westPoint       = westP;
  _midPoint        = ( eastP + westP ) / 2.;
 
  _type = MIDPOINT;
  conformalTrans( _midPoint.x(), _midPoint.y(), _driftDist );
}

◆ HoughHit() [3/3]

HoughHit::HoughHit ( const HoughHit & other )

Definition at line 53 of file HoughHit.cxx.

    : _digiPtr( other._digiPtr )
    , _det( other._det )
    , _layerPtr( other._layerPtr )
    , _wirePtr( other._wirePtr )
    , _id( other._id )
    , _layer( other._layer )
    , _wire( other._wire )
    , _rawTime( other._rawTime )
    , _charge( other._charge )
    , _eastPoint( other._eastPoint )
    , _westPoint( other._westPoint )
    , _midPoint( other._midPoint )
    , _u( other._u )
    , _v( other._v )
    , _type( other._type )
    , _detectorType( other._detectorType )
    , _rmid( other._rmid )
    , _slayerType( other._slayerType )
    , _cirlist( other._cirlist )
    , _style( other._style )
    , _driftTime( other._driftTime )
    , _driftDist( other._driftDist )
    , CF_drift( other.CF_drift )
    , vec_cfcir( other.vec_cfcir )
    ,
    // truth
    _truthU( other._truthU )
    , _truthV( other._truthV )
    , _truthR( other._truthR )
    , _truthlr( other._truthlr )
    , _truthPoint( other._truthPoint )
    , _truthDrift( other._truthDrift )
    , _deltad( other._deltad )
    , _flightLength( other._flightLength ) {}

Member Function Documentation

◆ calDriftDist() [1/2]

double HoughHit::calDriftDist	(	double	bunchTime,
		int	ambig ) const

Definition at line 188 of file HoughHit.cxx.

                                                                 {
 
  //  double crudeTof = 0; //FIXME
  return calDriftDist( bunchTime + crudeTof(), ambig, 0., 0., 0. );
}

◆ calDriftDist() [2/2]

double HoughHit::calDriftDist	(	double	tof,
		int	ambig,
		double	entranceAngle,
		double	,
		double	z ) const

Definition at line 194 of file HoughHit.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 * _calibPtr->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 ) < 0.00001 ) driftD = 0.00001;
  return driftD;
}

Referenced by calDriftDist(), and HoughHit().

◆ conformalTrans()

void HoughHit::conformalTrans	(	double	x,
		double	y,
		double	r )

Definition at line 154 of file HoughHit.cxx.

                                                            {
  _u       = 2 * x / ( x * x + y * y - r * r );
  _v       = 2 * y / ( x * x + y * y - r * r );
  CF_drift = 2 * r / ( x * x + y * y - r * r );
}

Referenced by HoughHit().

◆ digi()

const MdcDigi * HoughHit::digi ( ) const

inline

Definition at line 53 of file HoughHit.h.

53{ return _digiPtr; }

Referenced by HoughTrack::cald_layer(), HoughRecHit::HoughRecHit(), printAll(), and HoughHitList::remove().

◆ driftTime() [1/2]

double HoughHit::driftTime ( ) const

Definition at line 170 of file HoughHit.cxx.

                                 {
  double tprop = _calibPtr->getTprop( _layer, 0 );
  double T0Walk =
      _calibPtr->getT0( _layer, _wire ) + _calibPtr->getTimeWalk( _layer, _charge );
  // tof in ns, driftTime in ns, _T0Walk in ns
  double driftT = _rawTime - T0Walk - 1.e9 * _bunchTime - tprop;
  return driftT;
}

Referenced by HoughHitList::addMdcDigiList(), calDriftDist(), HoughTrack::find_pair_hit(), HoughTrack::find_stereo_hit(), and HoughHit().

◆ driftTime() [2/2]

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

Definition at line 178 of file HoughHit.cxx.

                                                       {
  double tprop = _calibPtr->getTprop( _layer, z * 10. );
  double T0Walk =
      _calibPtr->getT0( _layer, _wire ) + _calibPtr->getTimeWalk( _layer, _charge );
  // tof in ns, driftTime in ns, _T0Walk in ns
  double driftT = _rawTime - T0Walk - 1.e9 * tof - tprop;
 
  return driftT;
}

◆ getBunchTime()

double HoughHit::getBunchTime ( ) const

inline

Definition at line 54 of file HoughHit.h.

54{ return _bunchTime; }

◆ getCharge()

double HoughHit::getCharge ( ) const

inline

Definition at line 65 of file HoughHit.h.

65{ return _charge; }

◆ getCir()

CFCir HoughHit::getCir ( int i ) const

inline

Definition at line 47 of file HoughHit.h.

47{ return vec_cfcir[i]; }

◆ getCirList()

int HoughHit::getCirList ( ) const

inline

Definition at line 102 of file HoughHit.h.

102{ return _cirlist; }

Referenced by HoughTrack::cald_layer(), and HoughMap::select_slant().

◆ getConformal_r()

double HoughHit::getConformal_r	(	double	x,
		double	y,
		double	r )

Definition at line 166 of file HoughHit.cxx.

                                                              {
  return 2 * r / ( x * x + y * y - r * r );
}

Referenced by setTruthInfo().

◆ getConformal_u()

double HoughHit::getConformal_u	(	double	x,
		double	y,
		double	r )

Definition at line 160 of file HoughHit.cxx.

                                                              {
  return 2 * x / ( x * x + y * y - r * r );
}

Referenced by setTruthInfo().

◆ getConformal_v()

double HoughHit::getConformal_v	(	double	x,
		double	y,
		double	r )

Definition at line 163 of file HoughHit.cxx.

                                                              {
  return 2 * x / ( x * x + y * y - r * r );
}

Referenced by setTruthInfo().

◆ getDeltaD()

double HoughHit::getDeltaD ( ) const

inline

Definition at line 73 of file HoughHit.h.

73{ return _deltad; } // l1l2-R+-r

◆ getDetectorType()

detectorType HoughHit::getDetectorType ( ) const

inline

Definition at line 52 of file HoughHit.h.

52{ return _detectorType; }

◆ getDriftDist()

double HoughHit::getDriftDist ( ) const

inline

Definition at line 67 of file HoughHit.h.

67{ return _driftDist; }

Referenced by HoughTrack::cald_layer(), and HoughPeak::collectHits().

◆ getDriftDistTruth()

double HoughHit::getDriftDistTruth ( ) const

inline

Definition at line 87 of file HoughHit.h.

87{ return _truthDrift; }

◆ getDriftTime()

double HoughHit::getDriftTime ( ) const

inline

Definition at line 66 of file HoughHit.h.

66{ return _driftTime; }

◆ getEastPoint()

HepPoint3D HoughHit::getEastPoint ( ) const

inline

Definition at line 56 of file HoughHit.h.

56{ return _eastPoint; }

◆ getFltLen()

double HoughHit::getFltLen ( ) const

inline

Definition at line 74 of file HoughHit.h.

74{ return _flightLength; } // R*dtheta

◆ getIdTruth()

int HoughHit::getIdTruth ( ) const

inline

Definition at line 88 of file HoughHit.h.

88{ return _truthId; }

◆ getLayerId()

int HoughHit::getLayerId ( ) const

inline

Definition at line 60 of file HoughHit.h.

60{ return _layer; }

Referenced by HoughPeak::collectHits(), digi_in_track(), HoughTrack::find_pair_hit(), HoughTrack::find_stereo_hit(), and HoughMap::select_slant().

◆ getLrTruth()

int HoughHit::getLrTruth ( ) const

inline

Definition at line 89 of file HoughHit.h.

89{ return _truthlr; }

◆ getMidPoint()

HepPoint3D HoughHit::getMidPoint ( ) const

inline

Definition at line 55 of file HoughHit.h.

55{ return _midPoint; }

◆ getMidX()

double HoughHit::getMidX ( ) const

inline

Definition at line 58 of file HoughHit.h.

58{ return _midPoint.x(); }

Referenced by HoughTrack::cald_layer(), and HoughPeak::collectHits().

◆ getMidY()

double HoughHit::getMidY ( ) const

inline

Definition at line 59 of file HoughHit.h.

59{ return _midPoint.y(); }

Referenced by HoughTrack::cald_layer(), and HoughPeak::collectHits().

◆ getPointTruth()

HepPoint3D HoughHit::getPointTruth ( ) const

inline

Definition at line 93 of file HoughHit.h.

93{ return _truthPoint; }

◆ getR()

double HoughHit::getR ( ) const

inline

Definition at line 72 of file HoughHit.h.

72{ return CF_drift; } // drift in CFS

◆ getRTruth()

double HoughHit::getRTruth ( ) const

inline

Definition at line 92 of file HoughHit.h.

92{ return _truthR; }

◆ getSlayerType()

int HoughHit::getSlayerType ( ) const

inline

Definition at line 62 of file HoughHit.h.

62{ return _slayerType; }

Referenced by HoughTrack::cald_layer(), HoughPeak::collectHits(), HoughTrack::find_pair_hit(), HoughTrack::find_stereo_hit(), and HoughMap::select_slant().

◆ getStyle()

int HoughHit::getStyle ( ) const

inline

Definition at line 103 of file HoughHit.h.

103{ return _style; }

Referenced by HoughTrack::cald_layer(), and HoughMap::select_slant().

◆ getU()

double HoughHit::getU ( ) const

inline

Definition at line 70 of file HoughHit.h.

70{ return _u; }

Referenced by HoughMap::select_slant().

◆ getUTruth()

double HoughHit::getUTruth ( ) const

inline

Definition at line 90 of file HoughHit.h.

90{ return _truthU; }

Referenced by HoughTrack::cald_layer().

◆ getV()

double HoughHit::getV ( ) const

inline

Definition at line 71 of file HoughHit.h.

71{ return _v; }

Referenced by HoughMap::select_slant().

◆ getVTruth()

double HoughHit::getVTruth ( ) const

inline

Definition at line 91 of file HoughHit.h.

91{ return _truthV; }

Referenced by HoughTrack::cald_layer().

◆ getWestPoint()

HepPoint3D HoughHit::getWestPoint ( ) const

inline

Definition at line 57 of file HoughHit.h.

57{ return _westPoint; }

◆ getWireId()

int HoughHit::getWireId ( ) const

inline

Definition at line 61 of file HoughHit.h.

61{ return _wire; }

Referenced by HoughTrack::find_pair_hit(), HoughTrack::find_stereo_hit(), and small_layer().

◆ getXTruth()

double HoughHit::getXTruth ( ) const

inline

Definition at line 84 of file HoughHit.h.

84{ return _truthPoint.x(); }

◆ getYTruth()

double HoughHit::getYTruth ( ) const

inline

Definition at line 85 of file HoughHit.h.

85{ return _truthPoint.y(); }

◆ getZTruth()

double HoughHit::getZTruth ( ) const

inline

Definition at line 86 of file HoughHit.h.

86{ return _truthPoint.z(); }

◆ layer()

const MdcLayer * HoughHit::layer ( ) const

inline

Definition at line 63 of file HoughHit.h.

63{ return _layerPtr; }

Referenced by slayerType().

◆ makeCir()

void HoughHit::makeCir	(	int	n,
		double	phi_begin,
		double	phi_last,
		double	r )

Definition at line 243 of file HoughHit.cxx.

                                                                           {
  vec_cfcir.clear();
  vector<CFCir>().swap( vec_cfcir );
  for ( int i = 0; i < n; i++ )
  {
    double phi_slice = ( phi_last - phi_begin ) / n;
    double phi       = phi_begin + ( 1 / 2. + i ) * phi_slice;
    double x         = _u + r * cos( phi );
    double y         = _v + r * sin( phi );
    vec_cfcir.push_back( CFCir( x, y, phi, n, _u, _v, r ) );
  }
}

◆ operator=()

HoughHit & HoughHit::operator= ( const HoughHit & other )

Definition at line 89 of file HoughHit.cxx.

                                                   {
  _digiPtr      = ( other._digiPtr );
  _det          = ( other._det );
  _layerPtr     = ( other._layerPtr );
  _wirePtr      = ( other._wirePtr );
  _id           = ( other._id );
  _layer        = ( other._layer );
  _wire         = ( other._wire );
  _rawTime      = ( other._rawTime );
  _charge       = ( other._charge );
  _eastPoint    = ( other._eastPoint );
  _westPoint    = ( other._westPoint );
  _midPoint     = ( other._midPoint );
  _u            = ( other._u );
  _v            = ( other._v );
  _type         = ( other._type );
  _detectorType = ( other._detectorType );
  _rmid         = ( other._rmid );
  _slayerType   = ( other._slayerType );
  _cirlist      = ( other._cirlist );
  _style        = ( other._style );
  _driftTime    = ( other._driftTime );
  _driftDist    = ( other._driftDist );
  CF_drift      = ( other.CF_drift );
  vec_cfcir     = ( other.vec_cfcir );
  // truth
  _truthU       = ( other._truthU );
  _truthV       = ( other._truthV );
  _truthR       = ( other._truthR );
  _truthlr      = ( other._truthlr );
  _truthDrift   = ( other._truthDrift );
  _truthPoint   = ( other._truthPoint );
  _deltad       = ( other._deltad );
  _flightLength = ( other._flightLength );
 
  return *this;
}

Referenced by HoughRecHit::operator=().

◆ print()

void HoughHit::print ( ) const

Definition at line 217 of file HoughHit.cxx.

                           {
  std::cout << "(" << _layer << "," << _wire << ") " << std::endl;
}

◆ printAll()

void HoughHit::printAll ( ) const

Definition at line 221 of file HoughHit.cxx.

                              {
  std::cout << "(" << _layer << "," << _wire << ") id " << this->digi()->getTrackIndex()
            << " xyz " << _midPoint << endl;
}

◆ printTruth()

void HoughHit::printTruth ( ) const

Definition at line 225 of file HoughHit.cxx.

                                {
  std::cout << "(" << _layer << "," << _wire << ") truth " << _truthId << " xyz "
            << _truthPoint << " uv " << _truthU << "," << _truthV << ")  "
            << " cir list " << _cirlist << " ambig " << _truthlr << std::endl;
}

◆ setBunchTime()

void HoughHit::setBunchTime ( double t0 )

inlinestatic

Definition at line 37 of file HoughHit.h.

37{ _bunchTime = t0; }

Referenced by MdcHoughFinder::execute(), and MdcHoughFinder::initialize().

◆ setCirList()

void HoughHit::setCirList ( int cir )

inline

Definition at line 104 of file HoughHit.h.

104{ _cirlist = cir; }

◆ setDeltaD()

void HoughHit::setDeltaD ( double d )

inline

Definition at line 98 of file HoughHit.h.

98{ _deltad = d; } // l1l2-R+-r

Referenced by HoughTrack::cald_layer().

◆ setFltLen()

void HoughHit::setFltLen ( double flt )

inline

Definition at line 99 of file HoughHit.h.

99{ _flightLength = flt; } // R*dtheta

Referenced by HoughTrack::cald_layer().

◆ setMdcCalibFunSvc()

void HoughHit::setMdcCalibFunSvc ( const IMdcCalibFunSvc * calibSvc )

inlinestatic

Definition at line 35 of file HoughHit.h.

35{ _calibPtr = calibSvc; }

Referenced by MdcHoughFinder::initialize().

◆ setMdcGeomSvc()

void HoughHit::setMdcGeomSvc ( IMdcGeomSvc * geomSvc )

inlinestatic

Definition at line 36 of file HoughHit.h.

36{ _mdcGeomSvc = geomSvc; }

Referenced by MdcHoughFinder::initialize().

◆ setStyle()

void HoughHit::setStyle ( int sty )

inline

Definition at line 105 of file HoughHit.h.

105{ _style = sty; }

◆ setTruthInfo()

void HoughHit::setTruthInfo ( const MdcMcHit *& mcHit )

Definition at line 127 of file HoughHit.cxx.

                                                     {
  if ( !aMcHit ) return;
  _truthDrift = aMcHit->getDriftDistance() / 10.;   // mm to cm
  _truthPoint.setX( aMcHit->getPositionX() / 10. ); // mm to cm
  _truthPoint.setY( aMcHit->getPositionY() / 10. );
  _truthPoint.setZ( aMcHit->getPositionZ() / 10. );
 
  int mcLR = aMcHit->getPositionFlag();
  //  if (mcLR == 0) mcLR = -1;//FIXME
  //  for truth reserve from digi
  if ( mcLR == 1 ) mcLR = -1; //
  if ( mcLR == 0 ) mcLR = 1;  //
  _truthId = aMcHit->getTrackIndex();
  _truthlr = mcLR;
 
  // fix
  _truthU = getConformal_u( _truthPoint.x(), _truthPoint.y(), _truthDrift );
  _truthV = getConformal_v( _truthPoint.y(), _truthPoint.x(), _truthDrift );
  _truthR = getConformal_r( _truthPoint.x(), _truthPoint.y(), _truthDrift );
 
  // std::cout<<__FILE__<<"   "<<_layer<<","<<_wire<<" "<<_truthId<<" truth  "<<_truthPoint<<"
  // truthU
  // "<<_truthU<<std::endl;
 
  _type = DIGIWITHTRUTH;
}

◆ slayerType()

int HoughHit::slayerType ( int layer )

Definition at line 231 of file HoughHit.cxx.

                                    {
  // get nominal shift cell of this layer
  double nomShift = _mdcGeomSvc->Layer( layer )->nomShift();
 
  if ( nomShift > 0 ) return 1;
  else if ( nomShift < 0 ) return -1;
  else return 0;
 
  return -999;
}

Referenced by HoughHit().

◆ type()

HoughHitType HoughHit::type ( ) const

inline

Definition at line 51 of file HoughHit.h.

51{ return _type; }

◆ wire()

const MdcSWire * HoughHit::wire ( ) const

inline

Definition at line 64 of file HoughHit.h.

64{ return _wirePtr; }

Referenced by HoughHit().

Member Data Documentation

◆ _npart

int HoughHit::_npart = 360

static

Definition at line 48 of file HoughHit.h.

Referenced by MdcHoughFinder::initialize().

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

Public Member Functions

Static Public Member Functions

Static Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HoughHit() [1/3]

◆ HoughHit() [2/3]

◆ HoughHit() [3/3]

Member Function Documentation

◆ calDriftDist() [1/2]

◆ calDriftDist() [2/2]

◆ conformalTrans()

◆ digi()

◆ driftTime() [1/2]

◆ driftTime() [2/2]

◆ getBunchTime()

◆ getCharge()

◆ getCir()

◆ getCirList()

◆ getConformal_r()

◆ getConformal_u()

◆ getConformal_v()

◆ getDeltaD()

◆ getDetectorType()

◆ getDriftDist()

◆ getDriftDistTruth()

◆ getDriftTime()

◆ getEastPoint()

◆ getFltLen()

◆ getIdTruth()

◆ getLayerId()

◆ getLrTruth()

◆ getMidPoint()

◆ getMidX()

◆ getMidY()

◆ getPointTruth()

◆ getR()

◆ getRTruth()

◆ getSlayerType()

◆ getStyle()

◆ getU()

◆ getUTruth()

◆ getV()

◆ getVTruth()

◆ getWestPoint()

◆ getWireId()

◆ getXTruth()

◆ getYTruth()

◆ getZTruth()

◆ layer()

◆ makeCir()

◆ operator=()

◆ print()

◆ printAll()

◆ printTruth()

◆ setBunchTime()

◆ setCirList()

◆ setDeltaD()

◆ setFltLen()

◆ setMdcCalibFunSvc()

◆ setMdcGeomSvc()

◆ setStyle()

◆ setTruthInfo()

◆ slayerType()

◆ type()

◆ wire()

Member Data Documentation

◆ _npart