#include <HoughMap.h>

Public Member Functions
	HoughMap (int charge, HoughHitList &houghHitList, int mapHit, int ntheta, int nrho, double rhoMin, double rhoMaxi, int peakWidth, int peakHigh, double hitpro)
	HoughMap ()
	HoughMap (const HoughMap &other)
void	doMap ()
	~HoughMap ()
TH2D *	getMdcHoughFinder () const
const HoughHitList &	getHitList () const
int	getPeakNumber () const
int	getTrackNumber () const
const HoughPeak &	getPeak (int i) const
const HoughTrack &	getTrack (int i) const
const vector< HoughPeak >	getPeakList () const
const vector< HoughTrack >	getTrackList () const
int	getPeakWidth () const
int	getPeakHigh () const
int	getNTheta () const
int	getNRho () const
void	printPeak ()
void	printTrack ()
void	clearMap ()
void	select_slant ()
vector< double >	get_maxlayer_slant () const
double	get_maxlayer () const
vector< double >	get_nomaxlayer_slant () const
vector< int >	get_nomaxlayerid () const
double	exTheta (int, double, double, int)
double	exRho (int, double, double, int)
int	exThetaBin (double)
int	exRhoBin (double)
void	printMapHit ()

Public Attributes
double	Rho
double	Theta
double	Height
int	MAX

Static Public Attributes
static int	m_debug = 0
static int	m_useHalfCir = 0
static int	m_N1
static int	m_N2

Detailed Description

Definition at line 18 of file HoughMap.h.

Constructor & Destructor Documentation

◆ HoughMap() [1/3]

HoughMap::HoughMap	(	int	charge,
		HoughHitList &	houghHitList,
		int	mapHit,
		int	ntheta,
		int	nrho,
		double	rhoMin,
		double	rhoMaxi,
		int	peakWidth,
		int	peakHigh,
		double	hitpro )

Definition at line 17 of file HoughMap.cxx.

                                    {
  _hitList   = houghHitList;
  _mapHit    = mapHit;
  _nTheta    = ntheta;
  _nRho      = nrho;
  _thetaMin  = 0;
  _thetaMax  = M_PI;
  _rhoMin    = rhoMin;
  _rhoMax    = rhoMax;
  _peakWidth = peakWidth;
  _peakHigh  = peakHigh;
  _hitpro    = hitpro;
  _charge    = charge;
 
  //_mapHitList = new vector<const HoughHit*>*[_nTheta];
  //_houghS=  new double*[_nTheta];
  //_houghS2=  new double*[m_N2];
  //  _houghR=  new double*[_nTheta];
  //  _houghRL=  new double*[_nTheta];
  //  _houghNL=  new double*[_nTheta];
 
  //   for ( int i = 0; i < _nTheta; i++ )
  //   {
  //_mapHitList[i] = new vector<const HoughHit*>[_nRho];
  //    _houghS[i] = new double[_nRho];
  //    _houghR[i] = new double[_nRho];
  //    _houghRL[i] = new double[_nRho];
  //    _houghNL[i] = new double[_nRho];
  //    for(int j=0; j<_nRho; j++){
  //      _houghS[i][j]=0;
  //      _houghR[i][j]=0;
  //    }
  //   }
  //  for(int i=0; i<m_N2; i++){
  //    _houghS2[i] = new double[m_N2];
  //    for(int j=0; j<m_N2; j++){
  //      _houghS2[i][j]=0;
  //    }
  //  }
 
  if ( _charge == -1 )
    _houghSpace = new TH2D( "houghspace", "houghspace", _nTheta, _thetaMin, _thetaMax, _nRho,
                            _rhoMin, _rhoMax );
  if ( _charge == 1 )
    _houghSpace = new TH2D( "houghspace2", "houghspace2", _nTheta, _thetaMin, _thetaMax, _nRho,
                            _rhoMin, _rhoMax );
  //  _houghR= new TH2D("houghR","houghR",_nTheta,_thetaMin,_thetaMax,_nRho,_rhoMin,_rhoMax);
 
  doMap();
}

Referenced by HoughMap().

◆ HoughMap() [2/3]

HoughMap::HoughMap ( )

Definition at line 12 of file HoughMap.cxx.

                   {
  //_mapHitList=NULL;
  _houghSpace = NULL;
}

◆ HoughMap() [3/3]

HoughMap::HoughMap ( const HoughMap & other )

Definition at line 70 of file HoughMap.cxx.

                                          {
  _mapHit = other._mapHit, _nTheta = other._nTheta, _nRho = other._nRho,
  _thetaMin = other._thetaMin, _thetaMax = other._thetaMax, _rhoMin = other._rhoMin,
  _rhoMax = other._rhoMax, _hitList = other._hitList, _houghPeakList = other._houghPeakList,
  _houghTrackList = other._houghTrackList, _houghSpace = other._houghSpace,
  _charge = other._charge;
  //_houghR= other._houghR;
 
  //_mapHitList = new vector<const HoughHit*>*[_nTheta];
  //_houghNL= new double*[_nTheta];
  //_houghRL= new double*[_nTheta];
  //_houghS= new double*[_nTheta];
  // for(int i=0; i<_nTheta; i++){
  //_mapHitList[i] = new vector<const HoughHit*>[_nRho];
  //    _houghNL[i] = new double[_nRho];
  //    _houghRL[i] = new double[_nRho];
  //_houghS[i] = new double[_nRho];
  //}
 
  // for(int i=0; i<_nTheta; i++){
  //    for(int j=0; j<_nRho; j++){
  //_mapHitList[i][j]=other._mapHitList[i][j];
  //      _houghNL[i][j]=other._houghNL[i][j];
  //      _houghRL[i][j]=other._houghRL[i][j];
  //      _houghS[i][j]=other._houghS[i][j];
  //    }
  //  }
}

◆ ~HoughMap()

HoughMap::~HoughMap ( )

Definition at line 143 of file HoughMap.cxx.

143{ clearMap(); }

HoughMap::clearMap

void clearMap()

Definition HoughMap.cxx:121

Member Function Documentation

◆ clearMap()

void HoughMap::clearMap ( )

Definition at line 121 of file HoughMap.cxx.

                        {
 
  //  for(int i=0; i<_nTheta; i++){ delete []_mapHitList[i];delete []_houghNL; delete
  //  []_houghRL;delete []_houghS;} for(int i=0; i<_nTheta; i++){ delete []_mapHitList[i]; }
  // delete []_mapHitList;
  delete _houghSpace;
  //  for(int i=0; i<_nTheta; i++){ delete []_houghS[i];}
  //  delete []_houghS;
  //  delete _houghR;
  //  for(int i=0; i<m_N2; i++)   { delete []_houghS2[i];}
  //  delete []_houghS2;
}

Referenced by ~HoughMap().

◆ doMap()

void HoughMap::doMap ( )

Definition at line 99 of file HoughMap.cxx.

                     {
  buildMap();
  // mapDev();
  // cald_layer();
  // select_slant();
  // gravity();
  // findPeaks();
  // if(m_debug>0)  {cout<<" before sort "<<endl;  printPeak();}
  //  {cout<<" before sort "<<endl;  printPeak();}
  sortPeaks();
  hitFinding();
  if ( m_debug > 0 )
  {
    cout << " after sort " << endl;
    printPeak();
  }
  trackFinder();
  //  {cout<<" after sort "<<endl;  printPeak();}
  // candiTrack();
  if ( m_debug > 0 ) printTrack();
}

Referenced by HoughMap().

◆ exRho()

double HoughMap::exRho	(	int	irho,
		double	rhomin,
		double	rhomax,
		int	n )

Definition at line 382 of file HoughMap.cxx.

                                                                      {
  // double rho = _rhoMin+(irho+1/2.)*(_rhoMax-_rhoMin)/_nRho;
  double rho = rhomin + ( irho + 1 / 2. ) * ( rhomax - rhomin ) / n;
  return rho;
}

◆ exRhoBin()

int HoughMap::exRhoBin ( double )

◆ exTheta()

double HoughMap::exTheta	(	int	itheta,
		double	thetamin,
		double	thetamax,
		int	n )

Definition at line 387 of file HoughMap.cxx.

                                                                              {
  double theta = thetamin + ( itheta + 1 / 2. ) * ( thetamax - thetamin ) / n;
  return theta;
}

◆ exThetaBin()

int HoughMap::exThetaBin ( double )

◆ get_maxlayer()

double HoughMap::get_maxlayer ( ) const

inline

Definition at line 55 of file HoughMap.h.

55{ return m_maxlayer; }

◆ get_maxlayer_slant()

vector< double > HoughMap::get_maxlayer_slant ( ) const

inline

Definition at line 54 of file HoughMap.h.

54{ return maxlayer_slant; }

◆ get_nomaxlayer_slant()

vector< double > HoughMap::get_nomaxlayer_slant ( ) const

inline

Definition at line 56 of file HoughMap.h.

56{ return nomaxlayer_slant; }

◆ get_nomaxlayerid()

vector< int > HoughMap::get_nomaxlayerid ( ) const

inline

Definition at line 57 of file HoughMap.h.

57{ return nomaxlayerid; }

◆ getHitList()

const HoughHitList & HoughMap::getHitList ( ) const

inline

Definition at line 29 of file HoughMap.h.

29{ return _hitList; }

◆ getMdcHoughFinder()

TH2D * HoughMap::getMdcHoughFinder ( ) const

inline

Definition at line 27 of file HoughMap.h.

27{ return _houghSpace; }

◆ getNRho()

int HoughMap::getNRho ( ) const

inline

Definition at line 43 of file HoughMap.h.

43{ return _nRho; }

◆ getNTheta()

int HoughMap::getNTheta ( ) const

inline

Definition at line 42 of file HoughMap.h.

42{ return _nTheta; }

◆ getPeak()

const HoughPeak & HoughMap::getPeak ( int i ) const

inline

Definition at line 33 of file HoughMap.h.

33{ return _houghPeakList[i]; }

◆ getPeakHigh()

int HoughMap::getPeakHigh ( ) const

inline

Definition at line 41 of file HoughMap.h.

41{ return _peakHigh; }

◆ getPeakList()

const vector< HoughPeak > HoughMap::getPeakList ( ) const

inline

Definition at line 37 of file HoughMap.h.

37{ return _houghPeakList; }

◆ getPeakNumber()

int HoughMap::getPeakNumber ( ) const

inline

Definition at line 31 of file HoughMap.h.

31{ return _houghPeakList.size(); }

◆ getPeakWidth()

int HoughMap::getPeakWidth ( ) const

inline

Definition at line 40 of file HoughMap.h.

40{ return _peakWidth; }

◆ getTrack()

const HoughTrack & HoughMap::getTrack ( int i ) const

inline

Definition at line 34 of file HoughMap.h.

34{ return _houghTrackList[i]; }

◆ getTrackList()

const vector< HoughTrack > HoughMap::getTrackList ( ) const

inline

Definition at line 38 of file HoughMap.h.

38{ return _houghTrackList; }

Referenced by HoughTrackList::HoughTrackList().

◆ getTrackNumber()

int HoughMap::getTrackNumber ( ) const

inline

Definition at line 32 of file HoughMap.h.

32{ return _houghTrackList.size(); }

◆ printMapHit()

void HoughMap::printMapHit ( )

◆ printPeak()

void HoughMap::printPeak ( )

Definition at line 133 of file HoughMap.cxx.

                         {
  cout << "Print Peak in HoughMap sumPeak: " << _houghPeakList.size() << endl;
  vector<HoughPeak>::iterator iter = _houghPeakList.begin();
  for ( int i = 0; iter != _houghPeakList.end(); iter++, i++ )
  {
    cout << "count of Peak on HoughMap: " << ( *iter ).getHoughHitList().size() << endl;
    ( *iter ).printAllHit();
  }
}

Referenced by doMap().

◆ printTrack()

void HoughMap::printTrack ( )

Definition at line 488 of file HoughMap.cxx.

                          {
  cout << "Print Track in HoughMap: " << endl;
  vector<HoughTrack>::iterator iter = _houghTrackList.begin();
  for ( int i = 0; iter != _houghTrackList.end(); iter++, i++ )
  {
    cout << "Print Track" << i << endl;
    ( *iter ).printRecHit();
  }
}

Referenced by doMap().

◆ select_slant()

void HoughMap::select_slant ( )

Definition at line 614 of file HoughMap.cxx.

                            {
  int binx, biny, binz;
  _houghSpace->GetMaximumBin( binx, biny, binz );
  double                          thetabin = _houghSpace->GetXaxis()->GetBinCenter( binx );
  double                          rhobin   = _houghSpace->GetYaxis()->GetBinCenter( biny );
  vector<int>                     vec_layer;
  std::vector<HoughHit>::iterator iter1 = _hitList.getList().begin();
  for ( ; iter1 != _hitList.getList().end(); iter1++ )
  {
    if ( ( *iter1 ).getSlayerType() != 0 || ( *iter1 ).getStyle() != 0 ||
         ( *iter1 ).getCirList() != 0 )
      continue;
    vec_layer.push_back( ( *iter1 ).getLayerId() );
  }
 
  if ( vec_layer.size() == 0 ) return;
  vector<int>::iterator laymax   = max_element( vec_layer.begin(), vec_layer.end() );
  int                   maxLayer = *laymax;
  m_maxlayer                     = maxLayer;
  // cout<<"max layer "<<maxLayer<<endl;
 
  vector<int>::iterator iterlay = vec_layer.begin();
  for ( ; iterlay != vec_layer.end(); iterlay++ )
  {
    if ( *iterlay == *laymax )
    {
      vec_layer.erase( laymax );
      iterlay--;
    }
  }
  vector<int>::iterator laymax2   = max_element( vec_layer.begin(), vec_layer.end() );
  int                   maxLayer2 = *laymax2;
  //  cout<<"max layer2 "<<maxLayer2<<endl;
 
  std::vector<HoughHit>::iterator iter = _hitList.getList().begin();
  for ( ; iter != _hitList.getList().end(); iter++ )
  {
    HoughHit* hit = &( *iter );
    if ( hit->getSlayerType() != 0 || hit->getStyle() != 0 || hit->getCirList() != 0 )
      continue;
    double slantOfLine = hit->getV() * cos( thetabin ) - hit->getU() * sin( thetabin );
    // cout<<"("<<hit->getLayerId()<<","<<hit->getWireId()<<") "<< slantOfLine<<endl;
    if ( hit->getLayerId() == maxLayer || hit->getLayerId() == maxLayer2 )
    { maxlayer_slant.push_back( slantOfLine ); }
    else
    {
      nomaxlayer_slant.push_back( slantOfLine );
      nomaxlayerid.push_back( hit->getLayerId() );
    }
  }
}

Member Data Documentation

◆ Height

double HoughMap::Height

Definition at line 67 of file HoughMap.h.

◆ m_debug

int HoughMap::m_debug = 0

static

Definition at line 59 of file HoughMap.h.

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

◆ m_N1

int HoughMap::m_N1

static

Definition at line 61 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ m_N2

int HoughMap::m_N2

static

Definition at line 62 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ m_useHalfCir

int HoughMap::m_useHalfCir = 0

static

Definition at line 60 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ MAX

int HoughMap::MAX

Definition at line 69 of file HoughMap.h.

◆ Rho

double HoughMap::Rho

Definition at line 65 of file HoughMap.h.

◆ Theta

double HoughMap::Theta

Definition at line 66 of file HoughMap.h.

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

Public Member Functions

Public Attributes

Static Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HoughMap() [1/3]

◆ HoughMap() [2/3]

◆ HoughMap() [3/3]

◆ ~HoughMap()

Member Function Documentation

◆ clearMap()

◆ doMap()

◆ exRho()

◆ exRhoBin()

◆ exTheta()

◆ exThetaBin()

◆ get_maxlayer()

◆ get_maxlayer_slant()

◆ get_nomaxlayer_slant()

◆ get_nomaxlayerid()

◆ getHitList()

◆ getMdcHoughFinder()

◆ getNRho()

◆ getNTheta()

◆ getPeak()

◆ getPeakHigh()

◆ getPeakList()

◆ getPeakNumber()

◆ getPeakWidth()

◆ getTrack()

◆ getTrackList()

◆ getTrackNumber()

◆ printMapHit()

◆ printPeak()

◆ printTrack()

◆ select_slant()

Member Data Documentation

◆ Height

◆ m_debug

◆ m_N1

◆ m_N2

◆ m_useHalfCir

◆ MAX

◆ Rho

◆ Theta