Hough2D Class Reference

#include <Hough2D.h>

Public Member Functions
	Hough2D ()
	Hough2D (const Hough2D &other)
	Hough2D (recHitCol hitCol, double bunchtime)
TrkRecoTrk *	getTrk () const
double	getCirX () const
double	getCirY () const
double	getCirR () const
double	getD0 () const
double	getPhi0 () const
double	getOmega () const
int	getCharge () const
void	setCharge (int charge)
void	setCirX (double cirx)
void	setCirY (double ciry)
void	setCirR (double cirr)
void	setD0 (double d0)
void	setPhi0 (double phi0)
void	setOmega (double omega)
void	setPt (double pt)
void	setBunchT0 (double bunchT0)
void	setRecHit (const recHitCol &rechit)
int	fitLeast ()
virtual int	fit ()
void	print ()
void	printRecHit () const
double	getNfit () const
double	getChi2_2D () const
double	getPt () const

Static Public Member Functions
static void	setContext (TrkContextEv *context)
static void	setCalib (const MdcCalibFunSvc *mdcCalibFunSvc)

Static Public Attributes
static int	m_debug =0
static TrkContextEv *	_context =NULL
static const MdcCalibFunSvc *	_mdcCalibFunSvc =NULL
static int	m_qualityFactor =99999
static double	m_dropTrkChi2Cut =99999.
static double	m_dropTrkDrCut =1
static double	m_dropTrkDzCut =10
static double	m_dropTrkNhitCut =5
static double	m_dropTrkChi2NdfCut =99999.

Detailed Description

Definition at line 24 of file Hough2D.h.

Constructor & Destructor Documentation

Hough2D() [1/3]

Hough2D::Hough2D

(

)

Definition at line 14 of file Hough2D.cxx.

                {
  //_m_gm = MdcDetector::instance(0);
  //_m_gm = MdcDetector::instance(0);
}

Hough2D() [2/3]

Hough2D::Hough2D

(

const Hough2D &

other

)

Definition at line 26 of file Hough2D.cxx.

                                    {
  _circleR=other._circleR;
  _circleX=other._circleX;
  _circleY=other._circleY;
  _d0=other._d0;
  _phi0=other._phi0;
  _omega=other._omega;
  _charge=other._charge;
  _pt=other._pt;
  _nfit=other._nfit;
  _chi2_aver2D=other._chi2_aver2D;
  _bunchT0=other._bunchT0;
  //_m_gm = other._m_gm;
//  _recHitVec=other._recHitVec;
  _recHitVec=other._recHitVec;
 
//  for(int i =0;i<other._MdcHitCol.size();i++){
//  MdcHit* p_hit= new MdcHit( *(other._MdcHitCol[i]));
//  _MdcHitCol.push_back(p_hit);
//  }
}

Hough2D() [3/3]

Hough2D::Hough2D	(	recHitCol	hitCol,
		double	bunchtime
	)

Definition at line 48 of file Hough2D.cxx.

                                                  :_recHitVec(hitCol),_bunchT0(bunchtime){
  //_m_gm = MdcDetector::instance(0);
}

Member Function Documentation

fit()

int Hough2D::fit ( )

virtual

Definition at line 136 of file Hough2D.cxx.

                {
 // replace in houghTrack
 // if(_charge==-1){                                                                      
 //   _d0=_d0;                                                                               
 //   _omega=-1.*fabs(_omega);                                                                   
 // }                                                                                       
 // if(_charge==1){                                                                       
 //   _d0=-_d0;
 //   _omega=1.*fabs(_omega);
 //   _phi0=_phi0-M_PI;                                                                       
 // } 
  TrkExchangePar tt(_d0,_phi0,_omega,0,0);
  if (m_debug>0) cout<<"q d0 phi0 omega: "<<_charge<<" "<<_d0<<" "<<_phi0<<" "<<_omega<<endl;
  TrkCircleMaker circlefactory;
  float chisum =0.;
  newTrack = circlefactory.makeTrack(tt, chisum, *_context, _bunchT0);
  bool permitFlips = true;
  bool lPickHits = true;
  circlefactory.setFlipAndDrop(*newTrack, permitFlips, lPickHits);
  //int nDigi = digiToHots(newTrack);
  TrkHitList* trkHitList = newTrack->hits();
  int digiId=0;
  vector<MdcHit*> t_hitCol;
  std::vector<HoughRecHit>::iterator iter = _recHitVec.begin();
  for (;iter != _recHitVec.end(); iter++, digiId++) {
    if( (*iter).getflag()!=0) continue;
    //if( (*iter).getOuter()==1) continue;
//  const MdcDetector* mgm = Global::m_gm; 
    const MdcDigi* aDigi = (*iter).digi();
    //MdcHit *hit = new MdcHit(aDigi, _m_gm);
    //MdcHit *hit = new MdcHit(aDigi, mgm);
    MdcHit *hit = new MdcHit(aDigi, Global::m_gm);
//  MdcHit mdcHit(aDigi,_m_gm);
//  MdcHit *hit = &mdcHit;
    t_hitCol.push_back(hit);
    Identifier id = aDigi->identify();
    int layer = MdcID::layer(id);
    int wire  = MdcID::wire(id);
    //if(layer>=12) continue;
    //if ((layer>=0&&layer<=7)||(layer>=20&&layer<=35)) {continue;} 
    hit->setCalibSvc(_mdcCalibFunSvc);
    hit->setCountPropTime(true); 
    // new fltLen and ambig
    int newAmbig = 0;
    // calc. position of this point
    MdcRecoHitOnTrack temp(*hit, newAmbig, _bunchT0*1.e9);//m_bunchT0 nano second here
    MdcHitOnTrack* newhot = &temp;
    newhot->setFltLen( (*iter).getRet().second);             //caculate by mc
    double distoTrack= (*iter).getDisToTrack();
 
    //double ddCut=1.0;
    double ddCut=1.0;
    int use_in2d=1;
    if(hit->driftTime(_bunchT0,0)>1000) use_in2d=0;   //   >1000or800?
    if(hit->driftDist(_bunchT0,0)>ddCut)  use_in2d=0; 
    //if(m_debug>0) cout<<"flt : "<<(*iter).getRet().second <<endl;
    if(m_debug>0) std::cout<<" ("<<layer<<","<<wire<<",rt "<<hit->rawTime()<<",dt "<<hit->driftTime(_bunchT0,0)<<",dd "<<hit->driftDist(_bunchT0,0)<<") T0 " << _mdcCalibFunSvc->getT0(layer,wire) <<" timewalk "<<  _mdcCalibFunSvc->getTimeWalk(layer, aDigi->getChargeChannel())<< "dist: "<<hit->driftDist(_bunchT0,0)<<" disToTrk "<< distoTrack<<"  use?: "<<use_in2d<<endl;
    if(use_in2d==0) continue;
    trkHitList->appendHot(newhot);
  }
  //newTrack->printAll(std::cout);
 
  TrkHitList* qhits = newTrack->hits();
  TrkErrCode err=qhits->fit();
  const TrkFit* newFit = newTrack->fitResult();
  int nActiveHit = newTrack->hots()->nActive();
  int fit_stat=false;
  double chi2=-999.;
  if (!newFit || (newFit->nActive()<3)||err.failure()!=0) {
    if(m_debug>0){
      cout << " global 2d fit failed ";
      if(newFit) cout <<" nAct "<<newFit->nActive();
      cout<<"ERR1 failure ="<<err.failure()<<endl;
      fit_stat=0;
    }
  }else{
    TrkExchangePar par = newFit->helix(0.);
    if( abs( 1/(par.omega()) )>0.03) {
      fit_stat = 1;
      chi2=newFit->chisq();
      if(m_debug>0) cout<<"chi2 "<<chi2<<endl;
    }
    else {
      fit_stat = 0;
      chi2=-999;
    }
 
    bool badQuality = false;
    if(fabs(chi2)>m_qualityFactor*m_dropTrkChi2Cut ){
      if(m_debug>0){
        std::cout<<__FILE__<<"   "<<__LINE__<<" drop track by chi2  "<<chi2<<" > "<<m_qualityFactor<<" * "<<m_dropTrkChi2Cut <<std::endl;
      }
      badQuality=1;
    }
    if( fabs(par.d0())>m_qualityFactor*m_dropTrkDrCut) {
      if(m_debug>0){
        std::cout<<__FILE__<<"   "<<__LINE__<<" drop track by d0 "<<par.d0()<<" > "<<m_qualityFactor<<" * "<<m_dropTrkDrCut <<std::endl;
      }
      badQuality=1;
    }
    if( (fabs(chi2)/qhits->nHit()) > m_qualityFactor*m_dropTrkChi2NdfCut) {
      if(m_debug>0){
        std::cout<<__FILE__<<"   "<<__LINE__<<" drop track by chi2/ndf"<<(fabs(chi2)/qhits->nHit()) <<" > "<<m_qualityFactor<<" * "<<m_dropTrkChi2NdfCut<<std::endl;
      }
      badQuality=1;
    }
    if( nActiveHit <= m_dropTrkNhitCut) {
      if(m_debug>0){
        std::cout<<__FILE__<<"   "<<__LINE__<<" drop track by nhit"<<nActiveHit <<" <= "<<m_dropTrkNhitCut<<std::endl;
      }
      badQuality=1;
    }
    if( badQuality) fit_stat=0;
  }
  if( fit_stat==1 ){
    if(m_debug>0){
      cout << " global 2d fit success"<<endl;
      cout<<__FILE__<<__LINE__<<"AFTER fit 2d "<<endl;
      newTrack->print(std::cout);
    }
    TrkExchangePar par = newFit->helix(0.);
    double d0=par.d0();
    double phi0=par.phi0();
    double omega=par.omega();
    double r_temp=-1./par.omega();
    double x_cirtemp = sin(par.phi0()) *(par.d0()+1/par.omega()) * -1.; //z axis verse,x_babar = -x_belle
    double y_cirtemp = -1.*cos(par.phi0()) *(par.d0()+1/par.omega()) * -1;//???
    if(m_debug>0){
      cout<<" circle after globle 2d:  "<<"("<<x_cirtemp<<","<<y_cirtemp<<","<<r_temp<<")"<<endl;
      cout<<"pt_rec:   "<<1./omega/333.567<<endl;
    }
    _pt=1./omega/333.567; 
    _circleX=x_cirtemp;
    _circleY=y_cirtemp;
    _circleR=_charge/omega;
    _d0=par.d0();
    _phi0=par.phi0();
    _omega=par.omega();
 
    int nfit2d=0;
    if(m_debug>1) cout<<" hot list:"<<endl;
    TrkHotList::hot_iterator hotIter= qhits->hotList().begin();
    int lay=((MdcHit*)(hotIter->hit()))->layernumber();
    int wir=((MdcHit*)(hotIter->hit()))->wirenumber();
    int hittemp=lay*1000+wir;
    while (hotIter!=qhits->hotList().end()) {
      if(m_debug>1){ cout <<"(" <<((MdcHit*)(hotIter->hit()))->layernumber()
        <<","<<((MdcHit*)(hotIter->hit()))->wirenumber() 
          <<":"<<hotIter->isActive()<<")  ";
      }
      if (hotIter->isActive()==1) nfit2d++;
      hotIter++;
    }
    _nfit=nfit2d;
  }
  _chi2_aver2D=chi2/_nfit;
 
  for(int i=0;i<t_hitCol.size();i++){
    delete t_hitCol[i];
  }
  delete newTrack;
 
  if(m_debug>0) cout<<" in 2D fit number of Active hits : "<<_nfit<<endl;
  return fit_stat;
}

fitLeast()

int Hough2D::fitLeast

(

)

Definition at line 68 of file Hough2D.cxx.

                     {
  double x_sum=0;
  double y_sum=0;
  double x2_sum=0;
  double y2_sum=0;
  double x3_sum=0;
  double y3_sum=0;
  double x2y_sum=0;
  double xy2_sum=0;
  double xy_sum=0;
  double a1=0;
  double a2=0;
  double b1=0;
  double b2=0;
  double c1=0;
  double c2=0;
  //  double x_aver,y_aver,r2;
  int hitCandiSize=_recHitVec.size();
  if(hitCandiSize>=3){
    for(int i=0;i<hitCandiSize;i++){
      //    if( _recHitVec[i].getflag()!=0 ) continue;
      //cout<<"size of hit in least2d "<<hitCandiSize<<endl;
      //    if(_recHitVec[i]->getSlayerType() !=0 ) continue;
      x_sum=x_sum+( _recHitVec[i].getMidX() );
      y_sum=y_sum+( _recHitVec[i].getMidY() );
      x2_sum=x2_sum+( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidX() );
      y2_sum=y2_sum+( _recHitVec[i].getMidY() )*( _recHitVec[i].getMidY() );
      x3_sum=x3_sum+( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidX() );
      y3_sum=y3_sum+( _recHitVec[i].getMidY() )*( _recHitVec[i].getMidY() )*( _recHitVec[i].getMidY() );
      x2y_sum=x2y_sum+( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidY() );
      xy2_sum=xy2_sum+( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidY() )*( _recHitVec[i].getMidY() );
      xy_sum=xy_sum+( _recHitVec[i].getMidX() )*( _recHitVec[i].getMidY() );
    }
    //hitCandiSize+=1000;   //0 constrain   
    a1=x2_sum-x_sum*x_sum/hitCandiSize;
    a2=xy_sum-x_sum*y_sum/hitCandiSize;
    b1=a2;
    b2=y2_sum-y_sum*y_sum/hitCandiSize;
    c1=(x3_sum+xy2_sum-x_sum*(x2_sum+y2_sum)/hitCandiSize)/2.;
    c2=(y3_sum+x2y_sum-y_sum*(x2_sum+y2_sum)/hitCandiSize)/2.;
    //cout<<"a1 a2 b1 b2 c1 c2 "<<a1<<" "<<a2<<" "<<b1<<" "<<b2<<" "<<c1<<" "<<c2<<endl;
 
    //x_aver=x_sum/hitCandiSize;
    //y_aver=y_sum/hitCandiSize;
 
    _circleX =(b1*c2-b2*c1)/(b1*b1-a1*b2);
    _circleY =(b1*c1-a1*c2)/(b1*b1-a1*b2);
    _circleR=sqrt( (x2_sum+y2_sum-2*_circleX*x_sum-2*_circleY*y_sum)/hitCandiSize+_circleX*_circleX+_circleY*_circleY ); 
    //double r_temp=sqrt(r2);
    // cout<<"circle: ("<<_circleX<<","<<_circleY<<","<<_circleR<<endl;
 
    _d0=sqrt(_circleX*_circleX+_circleY*_circleY)-_circleR;
    _phi0=atan2(_circleY,_circleX)+M_PI/2.;          
    _omega=1/_circleR;               //according with Q
    double  pt_temp=_circleR/333.567;
    if(m_debug>0) cout<<" pt: "<<pt_temp<<endl;
    if(m_debug>0)  cout<<"par: ("<<_d0<<","<<_phi0<<","<<_omega<<")"<<endl;
    _pt=pt_temp;
 
    _circleX_least =(b1*c2-b2*c1)/(b1*b1-a1*b2);
    _circleY_least =(b1*c1-a1*c2)/(b1*b1-a1*b2);
    _circleR_least=sqrt( (x2_sum+y2_sum-2*_circleX*x_sum-2*_circleY*y_sum)/hitCandiSize+_circleX*_circleX+_circleY*_circleY ); 
    _pt_least=pt_temp;
    return 1;
  }
  else {return 0;}
}

getCharge()

int Hough2D::getCharge ( ) const

inline

Definition at line 39 of file Hough2D.h.

{return _charge;}

Referenced by Hough3D::Hough3D(), and HoughStereo::HoughStereo().

getChi2_2D()

double Hough2D::getChi2_2D ( ) const

inline

Definition at line 63 of file Hough2D.h.

{return _chi2_aver2D;}   

getCirR()

double Hough2D::getCirR ( ) const

inline

Definition at line 35 of file Hough2D.h.

{return _circleR;}

Referenced by HoughStereo::cald(), Hough3D::Hough3D(), HoughTrack::trackCharge2D(), and HoughTrack::trackCharge3D().

getCirX()

double Hough2D::getCirX ( ) const

inline

Definition at line 33 of file Hough2D.h.

{return _circleX;}

Referenced by HoughStereo::cald(), Hough3D::Hough3D(), HoughTrack::trackCharge2D(), and HoughTrack::trackCharge3D().

getCirY()

double Hough2D::getCirY ( ) const

inline

Definition at line 34 of file Hough2D.h.

{return _circleY;}

Referenced by HoughStereo::cald(), Hough3D::Hough3D(), HoughTrack::trackCharge2D(), and HoughTrack::trackCharge3D().

getD0()

double Hough2D::getD0 ( ) const

inline

Definition at line 36 of file Hough2D.h.

{return _d0;}

Referenced by Hough3D::Hough3D().

getNfit()

double Hough2D::getNfit ( ) const

inline

Definition at line 62 of file Hough2D.h.

{return _nfit;}

Referenced by HoughTrack::getNfit2D().

getOmega()

double Hough2D::getOmega ( ) const

inline

Definition at line 38 of file Hough2D.h.

{return _omega;}

Referenced by Hough3D::Hough3D().

getPhi0()

double Hough2D::getPhi0 ( ) const

inline

Definition at line 37 of file Hough2D.h.

{return _phi0;}

Referenced by Hough3D::Hough3D().

getPt()

double Hough2D::getPt ( ) const

inline

Definition at line 65 of file Hough2D.h.

{return _pt;}

getTrk()

TrkRecoTrk * Hough2D::getTrk ( ) const

inline

Definition at line 32 of file Hough2D.h.

{return newTrack;}

print()

void Hough2D::print

(

)

Definition at line 341 of file Hough2D.cxx.

                   {
  cout<<" print hough2d "<<endl;
  cout<<"par: "<<_d0<<" "<<_phi0<<" "<<_omega<<endl;
}

printRecHit()

void Hough2D::printRecHit

(

)

const

setBunchT0()

void Hough2D::setBunchT0 ( double  bunchT0 )

inline

Definition at line 52 of file Hough2D.h.

{_bunchT0=bunchT0;}

setCalib()

static void Hough2D::setCalib ( const MdcCalibFunSvc *  mdcCalibFunSvc )

inlinestatic

Definition at line 55 of file Hough2D.h.

{_mdcCalibFunSvc=mdcCalibFunSvc;}

Referenced by MdcHoughFinder::initialize().

setCharge()

void Hough2D::setCharge ( int  charge )

inline

Definition at line 42 of file Hough2D.h.

{_charge=charge;}

setCirR()

void Hough2D::setCirR ( double  cirr )

inline

Definition at line 45 of file Hough2D.h.

{_circleR=cirr;}

setCirX()

void Hough2D::setCirX ( double  cirx )

inline

Definition at line 43 of file Hough2D.h.

{_circleX=cirx;}

setCirY()

void Hough2D::setCirY ( double  ciry )

inline

Definition at line 44 of file Hough2D.h.

{_circleY=ciry;}

setContext()

static void Hough2D::setContext ( TrkContextEv *  context )

inlinestatic

Definition at line 54 of file Hough2D.h.

{_context=context;}

Referenced by MdcHoughFinder::initialize().

setD0()

void Hough2D::setD0 ( double  d0 )

inline

Definition at line 47 of file Hough2D.h.

{_d0=d0;}

setOmega()

void Hough2D::setOmega ( double  omega )

inline

Definition at line 49 of file Hough2D.h.

{_omega=omega;}

setPhi0()

void Hough2D::setPhi0 ( double  phi0 )

inline

Definition at line 48 of file Hough2D.h.

{_phi0=phi0;}

setPt()

void Hough2D::setPt ( double  pt )

inline

Definition at line 50 of file Hough2D.h.

{_pt=pt;}

setRecHit()

void Hough2D::setRecHit ( const recHitCol &  rechit )

inline

Definition at line 53 of file Hough2D.h.

{_recHitVec=rechit;}

Member Data Documentation

_context

TrkContextEv * Hough2D::_context =NULL

static

Definition at line 67 of file Hough2D.h.

Referenced by fit(), and setContext().

_mdcCalibFunSvc

const MdcCalibFunSvc * Hough2D::_mdcCalibFunSvc =NULL

static

Definition at line 68 of file Hough2D.h.

Referenced by fit(), and setCalib().

m_debug

int Hough2D::m_debug =0

static

Definition at line 66 of file Hough2D.h.

Referenced by fit(), fitLeast(), and MdcHoughFinder::initialize().

m_dropTrkChi2Cut

double Hough2D::m_dropTrkChi2Cut =99999.

static

Definition at line 72 of file Hough2D.h.

Referenced by fit().

m_dropTrkChi2NdfCut

double Hough2D::m_dropTrkChi2NdfCut =99999.

static

Definition at line 76 of file Hough2D.h.

Referenced by fit().

m_dropTrkDrCut

double Hough2D::m_dropTrkDrCut =1

static

Definition at line 73 of file Hough2D.h.

Referenced by fit().

m_dropTrkDzCut

double Hough2D::m_dropTrkDzCut =10

static

Definition at line 74 of file Hough2D.h.

m_dropTrkNhitCut

double Hough2D::m_dropTrkNhitCut =5

static

Definition at line 75 of file Hough2D.h.

Referenced by fit().

m_qualityFactor

int Hough2D::m_qualityFactor =99999

static

Definition at line 71 of file Hough2D.h.

Referenced by fit().

---

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

• Hough2D.h
• Hough2D.cxx