THistogram Class Reference

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

#include <THistogram.h>

List of all members.

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.
Private Attributes
unsigned	_nBins
unsigned *	_bins
bool *	_masks
AList< TMLink > **	_links
AList< TMLink >	_all
float	_binSize

---

Detailed Description

A class for a histogram used in tracking.

Definition at line 34 of file THistogram.h.

---

Constructor & Destructor Documentation

THistogram::THistogram

(

unsigned

nBins

)

Constructor.

Definition at line 19 of file THistogram.cxx.

References _bins, _binSize, _links, _masks, _nBins, genRecEmupikp::i, and M_PI.

                                     : _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

(

)

[virtual]

Destructor.

Definition at line 31 of file THistogram.cxx.

References _bins, _links, _masks, _nBins, and genRecEmupikp::i.

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

---

Member Function Documentation

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.

References _links, and _nBins.

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

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.

References _links, and _nBins.

Referenced by contents().

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

AList< TSegment > THistogram::clusters

(

void

)

const

returns an AList<TSegment0> of clusters.

Definition at line 229 of file THistogram.cxx.

References _bins, _links, _nBins, TTrackBase::append(), and delete_small_size::size.

Referenced by segments().

                               {
    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;
}

AList< TSegment0 > THistogram::clusters0

(

void

)

const

returns an AList<TSegment0> of clusters.

Definition at line 162 of file THistogram.cxx.

References _bins, _links, _nBins, TTrackBase::append(), and delete_small_size::size.

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

                                {
    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;
}

const AList< TMLink > & THistogram::contents

(

void

)

const [inline]

returns an AList<TMLink> of all contents.

Definition at line 154 of file THistogram.h.

References _all.

                               {
    return _all;
}

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.

References bin(), and genRecEmupikp::i.

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

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.

References bin(), and genRecEmupikp::i.

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

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.

References _bins, _links, _masks, _nBins, genRecEmupikp::i, ganga-rec::j, n1, and RealDBUtil::npos.

                                                                   {
    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;
}

void THistogram::fillPhi

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 111 of file THistogram.cxx.

References _all, _bins, _binSize, _links, _nBins, genRecEmupikp::i, M_PI, and boss::pos.

Referenced by TConformalFinder::findSegments().

                                               {
    _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);
    }
}

void THistogram::fillX

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 76 of file THistogram.cxx.

References _all, _bins, _binSize, _links, _nBins, genRecEmupikp::i, and boss::pos.

Referenced by TFastFinder::doit().

                                             {
    _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);
    }
}

void THistogram::fillY

(

const AList< TMLink > &

links

)

fills with hits.

Definition at line 94 of file THistogram.cxx.

References _all, _bins, _binSize, _links, _nBins, genRecEmupikp::i, and boss::pos.

                                             {
    _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);
    }
}

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.

References _masks.

                           {
    _masks[a] = true;
}

unsigned THistogram::nBin

(

unsigned

i

)

const [inline]

returns number in i'th bin.

Definition at line 160 of file THistogram.h.

References _bins, and _nBins.

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

unsigned THistogram::nBins

(

void

)

const [inline]

returns # of bins.

Definition at line 136 of file THistogram.h.

References _nBins.

                            {
    return _nBins;
}

void THistogram::remove

(

const AList< TMLink > &

links

)

removes links.

Definition at line 134 of file THistogram.cxx.

References _all, _bins, _links, _nBins, and genRecEmupikp::i.

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

AList< TSegment > THistogram::segments

(

void

)

const

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

Definition at line 259 of file THistogram.cxx.

References clusters(), genRecEmupikp::i, and ganga-rec::j.

Referenced by TConformalFinder::findSegments().

                               {
  //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;
}

AList< TSegment0 > THistogram::segments0

(

void

)

const

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

Definition at line 192 of file THistogram.cxx.

References clusters0(), genRecEmupikp::i, and ganga-rec::j.

                                {

    //...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;

}

void THistogram::unmask

(

void

)

[inline]

remove masks.

Definition at line 148 of file THistogram.h.

References _masks, _nBins, and genRecEmupikp::i.

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

---

Member Data Documentation

AList<TMLink> THistogram::_all [private]

Definition at line 105 of file THistogram.h.

Referenced by contents(), fillPhi(), fillX(), fillY(), and remove().

unsigned* THistogram::_bins [private]

Definition at line 102 of file THistogram.h.

Referenced by clusters(), clusters0(), dump(), fillPhi(), fillX(), fillY(), nBin(), remove(), THistogram(), and ~THistogram().

float THistogram::_binSize [private]

Definition at line 106 of file THistogram.h.

Referenced by fillPhi(), fillX(), fillY(), and THistogram().

AList<TMLink>** THistogram::_links [private]

Definition at line 104 of file THistogram.h.

Referenced by bin(), clusters(), clusters0(), dump(), fillPhi(), fillX(), fillY(), remove(), THistogram(), and ~THistogram().

bool* THistogram::_masks [private]

Definition at line 103 of file THistogram.h.

Referenced by dump(), mask(), THistogram(), unmask(), and ~THistogram().

unsigned THistogram::_nBins [private]

Definition at line 101 of file THistogram.h.

Referenced by bin(), clusters(), clusters0(), dump(), fillPhi(), fillX(), fillY(), nBin(), nBins(), remove(), THistogram(), unmask(), and ~THistogram().

---