#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();
}

◆ 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 62 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 137 of file HoughMap.cxx.

                   {
  clearMap();
}

Member Function Documentation

◆ clearMap()

void HoughMap::clearMap ( )

Definition at line 116 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 98 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 340 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 345 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 54 of file HoughMap.h.

54{return m_maxlayer;}

◆ get_maxlayer_slant()

vector< double > HoughMap::get_maxlayer_slant ( ) const

inline

Definition at line 53 of file HoughMap.h.

53{return maxlayer_slant;}

◆ get_nomaxlayer_slant()

vector< double > HoughMap::get_nomaxlayer_slant ( ) const

inline

Definition at line 55 of file HoughMap.h.

55{return nomaxlayer_slant;}

◆ get_nomaxlayerid()

vector< int > HoughMap::get_nomaxlayerid ( ) const

inline

Definition at line 56 of file HoughMap.h.

56{return nomaxlayerid;}

◆ getHitList()

const HoughHitList & HoughMap::getHitList ( ) const

inline

Definition at line 28 of file HoughMap.h.

28{return _hitList;}

◆ getMdcHoughFinder()

TH2D * HoughMap::getMdcHoughFinder ( ) const

inline

Definition at line 26 of file HoughMap.h.

26{return _houghSpace;}

◆ getNRho()

int HoughMap::getNRho ( ) const

inline

Definition at line 42 of file HoughMap.h.

42{return _nRho;}

◆ getNTheta()

int HoughMap::getNTheta ( ) const

inline

Definition at line 41 of file HoughMap.h.

41{return _nTheta;}

◆ getPeak()

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

inline

Definition at line 32 of file HoughMap.h.

32{return _houghPeakList[i];}

◆ getPeakHigh()

int HoughMap::getPeakHigh ( ) const

inline

Definition at line 40 of file HoughMap.h.

40{return _peakHigh;}

◆ getPeakList()

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

inline

Definition at line 36 of file HoughMap.h.

36{return _houghPeakList;}

◆ getPeakNumber()

int HoughMap::getPeakNumber ( ) const

inline

Definition at line 30 of file HoughMap.h.

30{return _houghPeakList.size();}

◆ getPeakWidth()

int HoughMap::getPeakWidth ( ) const

inline

Definition at line 39 of file HoughMap.h.

39{return _peakWidth;}

◆ getTrack()

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

inline

Definition at line 33 of file HoughMap.h.

33{return _houghTrackList[i];}

◆ getTrackList()

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

inline

Definition at line 37 of file HoughMap.h.

37{return _houghTrackList;}

Referenced by HoughTrackList::HoughTrackList().

◆ getTrackNumber()

int HoughMap::getTrackNumber ( ) const

inline

Definition at line 31 of file HoughMap.h.

31{return _houghTrackList.size();}

◆ printMapHit()

void HoughMap::printMapHit ( )

◆ printPeak()

void HoughMap::printPeak ( )

Definition at line 128 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 433 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 554 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 66 of file HoughMap.h.

◆ m_debug

int HoughMap::m_debug =0

static

Definition at line 58 of file HoughMap.h.

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

◆ m_N1

int HoughMap::m_N1

static

Definition at line 60 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ m_N2

int HoughMap::m_N2

static

Definition at line 61 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ m_useHalfCir

int HoughMap::m_useHalfCir =0

static

Definition at line 59 of file HoughMap.h.

Referenced by MdcHoughFinder::initialize().

◆ MAX

int HoughMap::MAX

Definition at line 68 of file HoughMap.h.

◆ Rho

double HoughMap::Rho

Definition at line 64 of file HoughMap.h.

◆ Theta

double HoughMap::Theta

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