THistogram Class Reference

A class for a histogram used in tracking. More...

#include <THistogram.h>

Public Member Functions
	THistogram (unsigned nBins)
	Constructor.
virtual	~THistogram ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const AList< TMLink > *const	bin (unsigned i) const
	returns a pointer to i'th AList<TMLink>.
const AList< TMLink > *const	bin (int i) const
	returns a pointer to i'th AList<TMLink>.
unsigned	nBin (unsigned i) const
	returns number in i'th bin.
unsigned	nBins (void) const
	returns # of bins.
AList< TMLink >	contents (unsigned center, unsigned width) const
	returns an AList of TMLinks within 'center' bin +- 'width' bin region.
AList< TMLink >	contents (int start, int end) const
	returns an AList of TMLinks from 'start' bin to 'end' bin region.
const AList< TMLink > &	contents (void) const
	returns an AList<TMLink> of all contents.
AList< TSegment0 >	clusters0 (void) const
	returns an AList<TSegment0> of clusters.
AList< TSegment >	clusters (void) const
	returns an AList<TSegment0> of clusters.
AList< TSegment0 >	segments0 (void) const
	returns an AList<TSegment0> using clusters() function.
AList< TSegment >	segments (void) const
	returns an AList<TSegment0> using clusters() function.
void	fillX (const AList< TMLink > &links)
	fills with hits.
void	fillY (const AList< TMLink > &links)
	fills with hits.
void	fillPhi (const AList< TMLink > &links)
	fills with hits.
void	remove (const AList< TMLink > &links)
	removes links.
void	mask (unsigned binN)
	masks a bin. Masked bin will not be returned by ::findPeak.
void	unmask (void)
	remove masks.

Detailed Description

A class for a histogram used in tracking.

Definition at line 34 of file THistogram.h.

Constructor & Destructor Documentation

THistogram()

THistogram::THistogram

(

unsigned

nBins

)

Constructor.

Definition at line 19 of file THistogram.cxx.

                                     : _nBins(nBins) {
    _binSize = 2. * M_PI / (float) _nBins;
    _bins = (unsigned *) malloc(_nBins * sizeof(unsigned));
    _masks = (bool *) malloc(_nBins * sizeof(bool));
    _links = (AList<TMLink> **) malloc(_nBins * sizeof(AList<TMLink> *));
    for (unsigned i = 0; i < _nBins; i++) {
    _bins[i] = 0;
    _masks[i] = false;
    _links[i] = new AList<TMLink>;
    }
}

~THistogram()

THistogram::~THistogram ( )

virtual

Destructor.

Definition at line 31 of file THistogram.cxx.

                        {
    free(_bins);
    free(_masks);
    for (unsigned i = 0; i < _nBins; i++)
    delete _links[i];
    free(_links);
}

Member Function Documentation

bin() [1/2]

const AList< TMLink > *const THistogram::bin ( int  i ) const

inline

returns a pointer to i'th AList<TMLink>.

Definition at line 129 of file THistogram.h.

                           {
    while (i < 0) i += _nBins;
    return _links[i % _nBins];
}

bin() [2/2]

const AList< TMLink > *const THistogram::bin ( unsigned  i ) const

inline

returns a pointer to i'th AList<TMLink>.

Definition at line 122 of file THistogram.h.

                                {
    if (i < _nBins) return _links[i];
    return 0;
}

clusters()

AList< TSegment > THistogram::clusters

(

void

)

const

returns an AList<TSegment0> of clusters.

Definition at line 229 of file THistogram.cxx.

                               {
    AList<TSegment> list;
//    std::cout<<"enter clusters"<<std::endl;
    //...Serach for empty bin...
    unsigned begin = 0;
    while (_bins[begin] > 0)begin++;
    if (begin == _nBins) return list;
 
    //...Start searching...
    unsigned loop = 0;
    while (loop < _nBins) {
    ++loop;
    unsigned id = (begin + loop) % _nBins;
    if (_bins[id]) {
        unsigned size = 0;
        TSegment * c = new TSegment();
        while (_bins[id]) {
        if (_bins[id]) ++size;
        c->append(* _links[id]);
        ++loop;
        id = (begin + loop) % _nBins;
        if (loop == _nBins) break;
        }
        list.append(c);
    }
    }
    return list;
}

Referenced by segments().

clusters0()

AList< TSegment0 > THistogram::clusters0

(

void

)

const

returns an AList<TSegment0> of clusters.

Definition at line 162 of file THistogram.cxx.

                                {
    AList<TSegment0> list;
 
    //...Serach for empty bin...
    unsigned begin = 0;
    while (_bins[begin] > 0) begin++;
    if (begin == _nBins) return list;
 
    //...Start searching...
    unsigned loop = 0;
    while (loop < _nBins) {
    ++loop;
    unsigned id = (begin + loop) % _nBins;
    if (_bins[id]) {
        unsigned size = 0;
        TSegment0 * c = new TSegment0();
        while (_bins[id]) {
        if (_bins[id]) ++size;
        c->append(* _links[id]);
        ++loop;
        id = (begin + loop) % _nBins;
        if (loop == _nBins) break;
        }
        list.append(c);
    }
    }
    return list;
}

Referenced by TFastFinder::doit(), TConformalFinder0::findClusters(), TConformalFinder0::findClusters2(), and segments0().

contents() [1/3]

AList< TMLink > THistogram::contents	(	int	start,
		int	end
	)		const

returns an AList of TMLinks from 'start' bin to 'end' bin region.

Definition at line 154 of file THistogram.cxx.

                                             {
    AList<TMLink> links;
    for (int i = start; i <= end; i++)
    links.append(* bin(i));
    return links;
}

contents() [2/3]

AList< TMLink > THistogram::contents	(	unsigned	center,
		unsigned	width
	)		const

returns an AList of TMLinks within 'center' bin +- 'width' bin region.

Definition at line 143 of file THistogram.cxx.

                                                          {
    AList<TMLink> links;
    for (int i = - (int) width;
     i <= (int) width;
     i++) {
    links.append(* bin((int) center + i));
    }
    return links;
}

contents() [3/3]

const AList< TMLink > & THistogram::contents ( void  ) const

inline

returns an AList<TMLink> of all contents.

Definition at line 154 of file THistogram.h.

                               {
    return _all;
}

dump()

void THistogram::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)		const

dumps debug information.

Definition at line 40 of file THistogram.cxx.

                                                                   {
    std::cout << pre;
    std::cout << "THistogram dump:#bins=" << _nBins << std::endl;
    unsigned nLoops = _nBins / 15 + 1;
    unsigned n0 = 0;
    unsigned n1 = 0;
    for (unsigned i = 0; i < nLoops; i++) {
    for (unsigned j = 0; j < 15; j++) {
        if (n0 == _nBins) break;
        printf("%4d", n0);
        ++n0;
    }
    std::cout << std::endl;
    for (unsigned j = 0; j < 15; j++) {
        if (n1 == _nBins) break;
        if (! _masks[n1]) printf("%4d", _bins[n1]);
        else              printf("-%3d", _bins[n1]);
        ++n1;
    }   
    std::cout << std::endl;
    }
 
    if (msg.find("detail") != std::string::npos) {
    for (unsigned i = 0; i < _nBins; i++) {
        std::cout << "bin " << i << " : ";
        for (unsigned j = 0; j < _links[i]->length(); j++) {
        std::cout << (* _links[i])[j]->wire()->name() << ",";
        }
        std::cout << std::endl;
    }
    }
 
    return;
}

fillPhi()

void THistogram::fillPhi

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 111 of file THistogram.cxx.

                                               {
    _all = (AList<TMLink> &) links;
    unsigned nLinks = links.length();
    double offset = _binSize / 4.;
    for (unsigned i = 0; i < nLinks; i++) {
    TMLink * l = links[i];
    const HepPoint3D & p = l->position();
    float phi = atan2(p.y(), p.x()) + M_PI;
//  std::cout<<"atan "<<atan2(-1., -1.)<<std::endl;
    unsigned pos = (unsigned) floor((phi + offset) / _binSize);
 
//  std::cout <<"layer "<<l->wire()->layerId()<<" cell "<<l->wire()->localId()
//        <<" x "<<p.x()<<" y "<<p.y()<<" pos "<<pos<<" xypos hit "<<l->hit()->xyPosition()<<" xypos wire "<<l->wire()->xyPosition()<<" drift "<<l->hit()->drift(0)<< "+-" <<l->hit()->dDrift(0)<<std::endl;
//  std::cout<<"binsize "<<_binSize<<" offset "<<offset<<" phi "<<phi<<std::endl;
    //...Why is this needed?...
    pos %= _nBins;
 
    ++_bins[pos];
    _links[pos]->append(l);
    }
}

fillX()

void THistogram::fillX

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 76 of file THistogram.cxx.

                                             {
    _all = (AList<TMLink> &) links;
    unsigned nLinks = links.length();
    double offset = _binSize / 4.;
    for (unsigned i = 0; i < nLinks; i++) {
    TMLink * l = links[i];
    const HepPoint3D & p = l->position();
    unsigned pos = (unsigned) floor((p.x() + offset) / _binSize);
 
    //...Why is this needed?...
    pos %= _nBins;
 
    ++_bins[pos];
    _links[pos]->append(l);
    }
}

Referenced by TFastFinder::doit(), and TConformalFinder0::findSegments().

fillY()

void THistogram::fillY

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 94 of file THistogram.cxx.

                                             {
    _all = (AList<TMLink> &) links;
    unsigned nLinks = links.length();
    for (unsigned i = 0; i < nLinks; i++) {
    TMLink * l = links[i];
    const HepPoint3D & p = l->position();
    unsigned pos = (unsigned) floor(p.y() / _binSize);
 
    //...Why is this needed?...
    pos %= _nBins;
 
    ++_bins[pos];
    _links[pos]->append(l);
    }
}

mask()

void THistogram::mask ( unsigned  binN )

inline

masks a bin. Masked bin will not be returned by ::findPeak.

Definition at line 142 of file THistogram.h.

                           {
    _masks[a] = true;
}

nBin()

unsigned THistogram::nBin ( unsigned  i ) const

inline

returns number in i'th bin.

Definition at line 160 of file THistogram.h.

                                 {
    if (i < _nBins) return _bins[i];
    return 0;
}

nBins()

unsigned THistogram::nBins ( void  ) const

inline

returns # of bins.

Definition at line 136 of file THistogram.h.

                            {
    return _nBins;
}

remove()

void THistogram::remove

(

const AList< TMLink > &

links

)

removes links.

Definition at line 134 of file THistogram.cxx.

                                              {
    for (unsigned i = 0; i < _nBins; i++) {
    _links[i]->remove(links);
    _bins[i] = _links[i]->length();
    }
    _all.remove(links);
}

segments()

AList< TSegment > THistogram::segments

(

void

)

const

returns an AList<TSegment0> using clusters() function.

Definition at line 259 of file THistogram.cxx.

                               {
  //yuany
  /*
  for (unsigned i = 0; i < _nBins; i++) {
      std::cout<<"i "<<i<<" bins[i] "<<_bins[i]<<std::endl;
    }
  */
    //...Obtain raw clusters...
    AList<TSegment> list = clusters();
    unsigned n = list.length();
    if (n == 0) return list;
 
    //...Examine each cluster...
    AList<TSegment> splitted;
    for (unsigned i = 0; i < n; i++) {
    TSegment * c = list[i];
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    base segment : ";
    c->dump("hits");
#endif
 
    AList<TSegment> newClusters = c->split();
    if (newClusters.length() == 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ... Solving dual hits" << std::endl;
#endif
        c->solveDualHits();
        continue;
    }
 
    list.append(newClusters);
    splitted.append(c);
#ifdef TRKRECO_DEBUG_DETAIL
    c->dump("hits", "    ");
    std::cout << "    ... splitted as" << std::endl;
    for (unsigned j = 0; j < newClusters.length(); j++) {
        std::cout << "    " << j << " : ";
        newClusters[j]->dump("hits");
    }
#endif
    }
    list.remove(splitted);
    HepAListDeleteAll(splitted);
    
    //yuany
//    n = list.length();
//    for (unsigned i = 0; i < n; i++) {
//  TSegment * c = list[i];
//  std::cout << "    base segment : ";
//  c->dump("hits");
//    }
 
    return list;
}

segments0()

AList< TSegment0 > THistogram::segments0

(

void

)

const

returns an AList<TSegment0> using clusters() function.

Definition at line 192 of file THistogram.cxx.

                                {
 
    //...Obtain raw clusters...
    AList<TSegment0> list = clusters0();
    unsigned n = list.length();
    if (n == 0) return list;
 
    //...Examine each cluster...
    AList<TSegment0> splitted;
    for (unsigned i = 0; i < n; i++) {
    TSegment0 * c = list[i];
 
    AList<TSegment0> newClusters = c->split();
    if (newClusters.length() == 0) {
        c->solveDualHits();
        continue;
    }
 
    list.append(newClusters);
    splitted.append(c);
#ifdef TRKRECO_DEBUG_DETAIL
    c->dump("hits", "    ");
    std::cout << "    ... splitted as" << std::endl;
    for (unsigned j = 0; j < newClusters.length(); j++) {
        std::cout << "    " << j << " : ";
        newClusters[j]->dump("hits");
    }
#endif
    }
    list.remove(splitted);
    HepAListDeleteAll(splitted);
 
    return list;
 
}

unmask()

void THistogram::unmask ( void  )

inline

remove masks.

Definition at line 148 of file THistogram.h.

                       {
    for (unsigned i = 0; i < _nBins; i++) _masks[i] = false;
}

---

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

• THistogram.h
• THistogram.cxx