#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 MdcCalibFunSvc *calibSvc)

static void	setMdcGeomSvc (MdcGeomSvc *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 19 of file HoughHit.cxx.

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

◆ HoughHit() [2/3]

HoughHit::HoughHit ( const MdcDigi *const digi )

Definition at line 25 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 54 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 161 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 167 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 129 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 55 of file HoughHit.h.

55{return _digiPtr;}

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

◆ driftTime() [1/2]

double HoughHit::driftTime ( ) const

Definition at line 145 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 152 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 56 of file HoughHit.h.

56{return _bunchTime;}

◆ getCharge()

double HoughHit::getCharge ( ) const

inline

Definition at line 67 of file HoughHit.h.

67{return _charge;}

◆ getCir()

CFCir HoughHit::getCir ( int i ) const

inline

Definition at line 48 of file HoughHit.h.

48{return vec_cfcir[i];}

◆ getCirList()

int HoughHit::getCirList ( ) const

inline

Definition at line 105 of file HoughHit.h.

105{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 141 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 135 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 138 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 75 of file HoughHit.h.

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

◆ getDetectorType()

detectorType HoughHit::getDetectorType ( ) const

inline

Definition at line 54 of file HoughHit.h.

54{return _detectorType;}

◆ getDriftDist()

double HoughHit::getDriftDist ( ) const

inline

Definition at line 69 of file HoughHit.h.

69{return _driftDist;}

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

◆ getDriftDistTruth()

double HoughHit::getDriftDistTruth ( ) const

inline

Definition at line 90 of file HoughHit.h.

90{return _truthDrift;}

◆ getDriftTime()

double HoughHit::getDriftTime ( ) const

inline

Definition at line 68 of file HoughHit.h.

68{return _driftTime;}

◆ getEastPoint()

HepPoint3D HoughHit::getEastPoint ( ) const

inline

Definition at line 58 of file HoughHit.h.

58{return _eastPoint;}

Referenced by HoughStereo::cald().

◆ getFltLen()

double HoughHit::getFltLen ( ) const

inline

Definition at line 76 of file HoughHit.h.

76{ return _flightLength; } // R*dtheta

◆ getIdTruth()

int HoughHit::getIdTruth ( ) const

inline

Definition at line 91 of file HoughHit.h.

91{return _truthId;}

◆ getLayerId()

int HoughHit::getLayerId ( ) const

inline

Definition at line 62 of file HoughHit.h.

62{return _layer;}

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

◆ getLrTruth()

int HoughHit::getLrTruth ( ) const

inline

Definition at line 92 of file HoughHit.h.

92{return _truthlr;}

◆ getMidPoint()

HepPoint3D HoughHit::getMidPoint ( ) const

inline

Definition at line 57 of file HoughHit.h.

57{return _midPoint;}

◆ getMidX()

double HoughHit::getMidX ( ) const

inline

Definition at line 60 of file HoughHit.h.

60{return _midPoint.x();}

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

◆ getMidY()

double HoughHit::getMidY ( ) const

inline

Definition at line 61 of file HoughHit.h.

61{return _midPoint.y();}

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

◆ getPointTruth()

HepPoint3D HoughHit::getPointTruth ( ) const

inline

Definition at line 96 of file HoughHit.h.

96{return _truthPoint;}

◆ getR()

double HoughHit::getR ( ) const

inline

Definition at line 74 of file HoughHit.h.

74{return CF_drift;} //drift in CFS

◆ getRTruth()

double HoughHit::getRTruth ( ) const

inline

Definition at line 95 of file HoughHit.h.

95{return _truthR;}

◆ getSlayerType()

int HoughHit::getSlayerType ( ) const

inline

Definition at line 64 of file HoughHit.h.

64{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 106 of file HoughHit.h.

106{return _style;}

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

◆ getU()

double HoughHit::getU ( ) const

inline

Definition at line 72 of file HoughHit.h.

72{return _u;}

Referenced by HoughMap::select_slant().

◆ getUTruth()

double HoughHit::getUTruth ( ) const

inline

Definition at line 93 of file HoughHit.h.

93{return _truthU;}

Referenced by HoughTrack::cald_layer().

◆ getV()

double HoughHit::getV ( ) const

inline

Definition at line 73 of file HoughHit.h.

73{return _v;}

Referenced by HoughMap::select_slant().

◆ getVTruth()

double HoughHit::getVTruth ( ) const

inline

Definition at line 94 of file HoughHit.h.

94{return _truthV;}

Referenced by HoughTrack::cald_layer().

◆ getWestPoint()

HepPoint3D HoughHit::getWestPoint ( ) const

inline

Definition at line 59 of file HoughHit.h.

59{return _westPoint;}

Referenced by HoughStereo::cald().

◆ getWireId()

int HoughHit::getWireId ( ) const

inline

Definition at line 63 of file HoughHit.h.

63{return _wire;}

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

◆ getXTruth()

double HoughHit::getXTruth ( ) const

inline

Definition at line 87 of file HoughHit.h.

87{return _truthPoint.x();}

◆ getYTruth()

double HoughHit::getYTruth ( ) const

inline

Definition at line 88 of file HoughHit.h.

88{return _truthPoint.y();}

◆ getZTruth()

double HoughHit::getZTruth ( ) const

inline

Definition at line 89 of file HoughHit.h.

89{return _truthPoint.z();}

◆ layer()

const MdcLayer * HoughHit::layer ( ) const

inline

Definition at line 65 of file HoughHit.h.

65{ return _layerPtr; }

Referenced by slayerType().

◆ makeCir()

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

Definition at line 210 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 79 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);
}

Referenced by HoughRecHit::operator=().

◆ print()

void HoughHit::print ( ) const

Definition at line 187 of file HoughHit.cxx.

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

◆ printAll()

void HoughHit::printAll ( ) const

Definition at line 191 of file HoughHit.cxx.

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

◆ printTruth()

void HoughHit::printTruth ( ) const

Definition at line 194 of file HoughHit.cxx.

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

◆ setBunchTime()

static 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 107 of file HoughHit.h.

107{_cirlist=cir;}

◆ setDeltaD()

void HoughHit::setDeltaD ( double d )

inline

Definition at line 101 of file HoughHit.h.

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

Referenced by HoughTrack::cald_layer().

◆ setFltLen()

void HoughHit::setFltLen ( double flt )

inline

Definition at line 102 of file HoughHit.h.

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

Referenced by HoughTrack::cald_layer().

◆ setMdcCalibFunSvc()

static void HoughHit::setMdcCalibFunSvc ( const MdcCalibFunSvc * calibSvc )

inlinestatic

Definition at line 35 of file HoughHit.h.

35{_calibPtr = calibSvc;}

Referenced by MdcHoughFinder::initialize().

◆ setMdcGeomSvc()

static void HoughHit::setMdcGeomSvc ( MdcGeomSvc * 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 108 of file HoughHit.h.

108{_style=sty;}

◆ setTruthInfo()

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

Definition at line 104 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 198 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 53 of file HoughHit.h.

53{return _type;}

◆ wire()

const MdcSWire * HoughHit::wire ( ) const

inline

Definition at line 66 of file HoughHit.h.

66{ return _wirePtr; }

Referenced by HoughHit().

Member Data Documentation

◆ _npart

int HoughHit::_npart =360

static

Definition at line 49 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