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 IMdcCalibFunSvc *mdcCalibFunSvc)

Static Public Attributes
static int	m_debug = 0
static TrkContextEv *	_context = NULL
static const IMdcCalibFunSvc *	_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);
}

Referenced by Hough2D().

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 147 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 69 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(), Hough3D::Hough3D(), and HoughStereo::HoughStereo().

getChi2_2D()

double Hough2D::getChi2_2D ( ) const

inline

Definition at line 65 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 Hough3D::Hough3D(), and Hough3D::Hough3D().

getCirX()

double Hough2D::getCirX ( ) const

inline

Definition at line 33 of file Hough2D.h.

{ return _circleX; }

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

getCirY()

double Hough2D::getCirY ( ) const

inline

Definition at line 34 of file Hough2D.h.

{ return _circleY; }

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

getD0()

double Hough2D::getD0 ( ) const

inline

Definition at line 36 of file Hough2D.h.

{ return _d0; }

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

getNfit()

double Hough2D::getNfit ( ) const

inline

Definition at line 64 of file Hough2D.h.

{ return _nfit; }

getOmega()

double Hough2D::getOmega ( ) const

inline

Definition at line 38 of file Hough2D.h.

{ return _omega; }

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

getPhi0()

double Hough2D::getPhi0 ( ) const

inline

Definition at line 37 of file Hough2D.h.

{ return _phi0; }

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

getPt()

double Hough2D::getPt ( ) const

inline

Definition at line 67 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 388 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()

void Hough2D::setCalib ( const IMdcCalibFunSvc * 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()

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 69 of file Hough2D.h.

Referenced by fit(), and setContext().

_mdcCalibFunSvc

const IMdcCalibFunSvc * Hough2D::_mdcCalibFunSvc = NULL

static

Definition at line 70 of file Hough2D.h.

Referenced by fit(), and setCalib().

m_debug

int Hough2D::m_debug = 0

static

Definition at line 68 of file Hough2D.h.

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

m_dropTrkChi2Cut

double Hough2D::m_dropTrkChi2Cut = 99999.

static

Definition at line 74 of file Hough2D.h.

Referenced by fit().

m_dropTrkChi2NdfCut

double Hough2D::m_dropTrkChi2NdfCut = 99999.

static

Definition at line 78 of file Hough2D.h.

Referenced by fit().

m_dropTrkDrCut

double Hough2D::m_dropTrkDrCut = 1

static

Definition at line 75 of file Hough2D.h.

Referenced by fit().

m_dropTrkDzCut

double Hough2D::m_dropTrkDzCut = 10

static

Definition at line 76 of file Hough2D.h.

m_dropTrkNhitCut

double Hough2D::m_dropTrkNhitCut = 5

static

Definition at line 77 of file Hough2D.h.

Referenced by fit().

m_qualityFactor

int Hough2D::m_qualityFactor = 99999

static

Definition at line 73 of file Hough2D.h.

Referenced by fit().

---

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

• Hough2D.h
• Hough2D.cxx