TBuilder Class Reference

A class to build a track. More...

#include <TBuilder.h>

Public Member Functions
	TBuilder (const std::string &name, float maxSigma, float maxSigmaStereo, float salvageLevel, float szSegmentDistance, float szLinkDistance, unsigned fittingFlag)
	Constructor with salvage level.
virtual	~TBuilder ()
	Destructor.
const std::string &	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildRphi (const AList< TSegment > &) const
TTrack *	buildStereoNew (const TTrack &t, AList< TSegment > &goodSegments, AList< TSegment > &badSegments) const
TTrack *	buildStereo (const TTrack &t, AList< TSegment > &) const
TTrack *	buildStereo (TTrack &t, TMLine &l, const AList< TMLink > &links) const
TTrack *	buildStereoOld (TTrack &t, const AList< TMLink > &allLinks, const AList< TMLink > &isolatedLinks) const
TTrack *	build (TTrack &t, const TMLine &l) const
void	salvage (TTrack &t, AList< TMLink > &hits) const
	salvages hits.
void	salvage (TTrack &t, AList< TSegment > &segments) const
TMLine *	initialLine (const TTrack &, AList< TSegment > &) const
	makes a line.
TMLine *	initialLineOld (const TTrack &, AList< TSegment > &) const
TMLine *	initialLine1 (const TTrack &, const AList< TSegment > &, const AList< TMLink > &) const
TMLine *	initialLine2 (const TTrack &, const AList< TMLink > &) const
AList< TMLine >	searchInitialLines (unsigned nSuperLayers) const
AList< TMLine >	searchLines6 (void) const
AList< TMLine >	searchLines5 (void) const
AList< TMLine >	searchLines4 (void) const
AList< TMLine >	searchLines3 (void) const
AList< TMLine >	searchLines2 (void) const
AList< TMLine >	searchLines1 (void) const
TMLine	searchLine (const TMLine &initialLine) const
AList< TSegment >	selectStereoSegment (const TMLine &line, const AList< TSegment > &list, const AList< TMLink > &szList) const
void	removeFarSegment (const TMLine &, AList< TSegment > &, AList< TMLink > &) const

Detailed Description

A class to build a track.

Definition at line 34 of file TBuilder.h.

Constructor & Destructor Documentation

TBuilder()

TBuilder::TBuilder	(	const std::string &	name,
		float	maxSigma,
		float	maxSigmaStereo,
		float	salvageLevel,
		float	szSegmentDistance,
		float	szLinkDistance,
		unsigned	fittingFlag
	)

Constructor with salvage level.

Definition at line 36 of file TBuilder.cxx.

: _name(a),
  _fitter("TBuilder Fitter"),
  _maxSigma(maxSigma),
  _maxSigmaStereo(maxSigmaStereo),
  _salvageLevel(sqrt(salvageLevel)),
//  _minNCores(3), //20060307
  _minNCores(2), //liuqg20070312
  _szSegmentDistance(szSegmentDistance),
  _szLinkDistance(szLinkDistance) {
    if (fittingFlag & 1) _fitter.sag(true);
    if (fittingFlag & 2) _fitter.propagation(true);
    if (fittingFlag & 4) _fitter.tof(true);
    if (fittingFlag & 8) _fitter.freeT0(true);
}

~TBuilder()

TBuilder::~TBuilder ( )

virtual

Destructor.

Definition at line 58 of file TBuilder.cxx.

                    {
}

Member Function Documentation

build()

TTrack * TBuilder::build	(	TTrack &	t,
		const TMLine &	l
	)		const

Definition at line 1255 of file TBuilder.cxx.

                                                  {
    AList<TMLink> links = l.links();
    AList<TMLink> toRemove;
    unsigned n = links.length();
    for (unsigned i = 1; i < n; i++) {
    if (links[i - 1]->link() == links[i]->link()) {
        toRemove.append(links[i]);
        continue;
    }
    if (i < 2)
        continue;
    if (links[i - 2]->link() == links[i]->link())
        toRemove.append(links[i]);
    }
    links.remove(toRemove);
 
    //...Pick up links...
    n = links.length();
    for (unsigned i = 0; i < n; i++)
    t.append(* links[i]->link());
 
    Vector a(5);
    a = t.helix().a();
    a[3] = l.b();
    //a[4] = t.charge() * l.a();//yzhang delete 2012-05-08 
    //    //yzhang add 2012-05-08 
    if(_fitter.getMagneticFieldPointer() == NULL){
      std::cout<<"   "<<__FILE__<<"   "<<__LINE__<<" Could not get MagneticFieldSvc  "<<std::endl;
    }
    const double Bz = -1000 * _fitter.getMagneticFieldPointer()->getReferField();
    if(Bz>0){
      a[4] = t.charge() * l.a();//yzhang add 2012-05-08 ???
    }else{
      a[4] = -1. * t.charge() * l.a();//yzhang add 2012-05-08  ???
    }
    //zhangy
    t._helix->a(a);
 
//#ifdef TRKRECO_WINDOW
//    TTrack * td = new TTrack(t);
//    sz.append(* td, leda_green);
//#endif
#ifdef TRKRECO_DEBUG_DETAIL
    static unsigned nTrk = 0;
    ++nTrk;
//      for (unsigned i = 0; i < t.nCores(); i++)
//      if (t.cores()[i]->wire()->name() == "43=220") {
//          std::cout << "43=220 removed" << std::endl;
//          t.remove(* t.cores()[i]);
//      }
    HepVector3D pp;
    if (links[0]->hit()->mc())
    pp = links[0]->hit()->mc()->hep()->p().vect();
    const HepVector3D p0 = pp;
//      if (nTrk == 1) {
//      Helix tmph(Point3D(0, 0, 0), p0, +1);
//      * t._helix = tmph;
//      }
    t.dump("detail", "before fit");
    std::cout << "Pdif mag=" << (t.p() - p0).mag() << std::endl;
//      t.links().removeAll();
//      for (unsigned i = 0; i < BsCouTab(RECMDC_WIRHIT
//    }
#endif
 
    //...Refine...
    AList<TMLink> bad;
    int err = _fitter.fit(t);
 
#ifdef TRKRECO_DEBUG_DETAIL
    t.dump("detail", "after fit");
    std::cout << "Pdif mag=" << (t.p() - p0).mag() << std::endl;
#endif
 
    t.refine(bad, _maxSigmaStereo * 100.);
    err = _fitter.fit(t);
    t.refine(bad, _maxSigmaStereo * 10.);
    err = _fitter.fit(t);
    t.refine(bad, _maxSigmaStereo);
#ifdef TRKRECO_DEBUG_DETAIL
//    if (nTrk == 2) {
//      Helix tmph(Point3D(0, 0, 0), p0, +1);
//      * t._helix = tmph;
//      std::cout << "initial mom" << p0 << " is set" << std::endl;
//    }
#endif
    err = _fitter.fit(t);
 
#ifdef TRKRECO_DEBUG_DETAIL
    t.dump("detail","         ");
    std::cout << "Pdif mag=" << (t.p() - p0).mag() << std::endl;
#endif
 
    //...Termination...
    return & t;
}

Referenced by buildStereoNew().

buildRphi() [1/2]

TTrack * TBuilder::buildRphi

(

const AList< TMLink > &

list

)

const

builds a r/phi track from TMLinks or from Segments.

Definition at line 107 of file TBuilder.cxx.

                                                    {
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi by links : # of links = ";
    std::cout << list.length() << std::endl;
#endif
 
    //...Classify TMLink's...
    AList<TMLink> cores;
    AList<TMLink> nonCores;
    SeparateCores(list, cores, nonCores);
 
#ifdef TRKRECO_DEBUG
    cout<<"  ... cores..."<<endl;
    for(int ii = 0; ii < cores.length(); ++ii) {
    cout<<"layer: "<<cores[ii]->wire()->layerId()
        <<"  local: "<<cores[ii]->wire()->localId()<<endl;
    }
    cout<<"  ...noncores..."<<endl;
    for(int ii = 0; ii < nonCores.length(); ++ii) {
        cout<<"layer: "<<nonCores[ii]->wire()->layerId()
            <<"  local: "<<nonCores[ii]->wire()->localId()<<endl;
    }
#endif
    //...Check # of links...
    unsigned nCores = cores.length();
    if (nCores < _minNCores) {
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi failure : # of cores(=" << nCores;
    std::cout << ") is less then " << _minNCores << std::endl;
#endif
 
    return NULL;
    }
 
    //...Make a circle...
#ifdef TRKRECO_DEBUG
    std::cout <<"links in list = "<<list.length()<<std::endl;
    std::cout << "... making a circle : # cores =" << cores.length() << std::endl;
#endif
    TCircle c(cores);
    int err = c.fit();
    if (err < 0) {
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi failure : circle fit error = ";
    std::cout << err << std::endl;
#endif
    return NULL;
    }
//  cout<<"radius "<<c.radius()<<" charge "<<c.charge()<<" center  "<<c.center()<<endl;
    //...Make a track...
    TTrack * t = new TTrack(c);
//   t->dump("hits sort flag", "    ");
    AList<TMLink> bad;
    err = _fitter.fit(* t);
    //cout<<"fit 1 done err: "<<err<<endl;
    // t->dump("hits sort flag", "    ");
    if (err < 0) goto discard;
    t->refine(bad, _maxSigma * 100.);
    //cout<<"refine 1 done"<<endl;
    //t->dump("hits sort flag", "    ");
    err = _fitter.fit(* t); 
    //cout<<"fit 2 done err: "<<err<<endl;
    //t->dump("hits sort flag", "    ");
    t->refine(bad, _maxSigma * 10.);
    //cout<<"refine 2 done"<<endl;
    //t->dump("hits sort flag", "    ");
    err = _fitter.fit(* t); 
    //cout<<"fit 3 done err: "<<err<<endl;
    //t->dump("hits sort flag", "    ");
    t->refine(bad, _maxSigma);
    //cout<<"refine 3 done"<<endl;
    //t->dump("hits sort flag", "    ");
    err = _fitter.fit(* t);
    //cout<<"fit 4 done err: "<<err<<endl;
    
    if (err < 0) goto discard;
#ifdef TRKRECO_DEBUG_DETAIL
    c.dump("detail", "    ccl> ");
    t->dump("detail", "    1st> ");
#endif
 
    //...Try to append non-core hits...
#ifdef TRKRECO_DEBUG
    std::cout << "... appending non-core hits : # = " << nonCores.length();
    std::cout << std::endl;
#endif
    t->appendByApproach(nonCores, _salvageLevel);
#ifdef TRKRECO_DEBUG
    t->dump("hits sort flag", "    ");
#endif
//   t->dump("hits sort flag", "    ");
    return t;
    
    //...Something happened...
discard:
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi failure : helix fit error = ";
    std::cout << err << std::endl;
#endif
    delete t;
    return NULL;
}

Referenced by buildRphi().

buildRphi() [2/2]

TTrack * TBuilder::buildRphi

(

const AList< TSegment > &

list

)

const

Definition at line 66 of file TBuilder.cxx.

                                                      {
 
#ifdef TRKRECO_DEBUG
    std::cout << "... building rphi by segments : # of segments = ";
    std::cout << list.length() << std::endl;
    for (unsigned i = 0; i < list.length(); i++)
    list[i]->dump("hits sort flag", "    ");
#endif
    //yuany
    //    for (unsigned i = 0; i < list.length(); i++)
    //    list[i]->dump("", "    ");
 
    //...Pick up links...
    AList<TMLink> links = Links(list);
    //yuany
//    for(unsigned i=0;i<links.length();i++){
//      cout<<"wire name "<<links[i]->wire()->name()<<endl;
//    }
 
    //...Main funtion...
    TTrack * t = buildRphi(links);
 
    //...Check used segments...
    if (t) {
    const AList<TMLink> & usedLinks = t->links();
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
        TSegment & segment = * list[i];
        AList<TMLink> used = Links(segment, * t);
        if (used.length()) {
        t->segments().append(segment);
        segment.tracks().append(t);
        }
    }
    }
 
    return t;
 
}

buildStereo() [1/2]

TTrack * TBuilder::buildStereo	(	const TTrack &	t,
		AList< TSegment > &	segments
	)		const

Definition at line 658 of file TBuilder.cxx.

                                                                         {
 
    TTrack & t = * new TTrack(t0);
 
#ifdef TRKRECO_DEBUG
    std::cout << "... building stereo by segments : # of segments = ";
    std::cout << segments.length() << std::endl;
    for (unsigned i = 0; i < segments.length(); i++)
    segments[i]->dump("hits sort flag", "    ");
#endif
#ifdef TRKRECO_WINDOW
    sz.clear();
    sz.skip(false);
    sz.mode(2);
    sz.appendSz(t, segments, leda_black);
    AList<TSegment> tmps = segments;
    std::string s = "# of segments = " + itostring(int(segments.length()));
#endif
 
    //...Find initial line...
    TMLine * line = initialLineOld(t, segments);
    if (! line) {
#ifdef TRKRECO_DEBUG
    std::cout << "... building stereo failure : no initial line found" << std::endl;
#endif
#ifdef TRKRECO_WINDOW
    s = "no initial line found : " + s;
    sz.text(s);
    sz.wait();
#endif
    return NULL;
    }
#ifdef TRKRECO_WINDOW
    sz.append(* line, leda_red);
    sz.text(s);
    sz.wait();
#endif
 
    //...Pick up links...
    AList<TMLink> links;
    unsigned n = segments.length();
    for (unsigned i = 0; i < n; i++)
    links.append(segments[i]->links());
 
    //...Main funtion...
    TTrack * ts = buildStereo(t, * line, links);
 
    //...Check used segments...
    if (ts) {
    AList<TMLink> usedLinks = StereoHits(t.links());
    for (unsigned i = 0; i < n; i++) {
        TSegment & segment = * segments[i];
        AList<TMLink> used = Links(segment, t);
        if (used.length()) {
        t.segments().append(segment);
        segment.tracks().append(t);
        }
    }
    }
 
    HepAListDeleteAll((AList<TMLink> &) line->links());
    delete line;
    return ts;
}

Referenced by buildStereo().

buildStereo() [2/2]

TTrack * TBuilder::buildStereo	(	TTrack &	t,
		TMLine &	l,
		const AList< TMLink > &	links
	)		const

Definition at line 525 of file TBuilder.cxx.

                                                {
#ifdef TRKRECO_DEBUG
    std::cout << "... building stereo by links : # of links = ";
    std::cout << list.length() << std::endl;
#endif
 
    //...Classify TMLink's...
    AList<TMLink> cores;
    AList<TMLink> nonCores;
    SeparateCores(list, cores, nonCores);
    //...Check # of links...
    unsigned nCores = cores.length();
    if (nCores < _minNCores) {
#ifdef TRKRECO_DEBUG
    std::cout << "... stereo building failure : # of cores(=" << nCores;
    std::cout << ") is less then " << _minNCores << std::endl;
#endif
    return NULL;
    }
#ifdef TRKRECO_WINDOW
    sz.appendSz(track, cores, leda_brown);
#endif
 
    //...Cal. left and right position...
    AList<TMLink> forNewLine;
    for (unsigned i = 0; i < nCores; i++) {
    TMLink & t = * cores[i];
    for (unsigned i = 0; i < 2; i++) {
        TMLink & tt = * new TMLink(t);
        tt.leftRight(i);
        int err = track.szPosition(tt);
        if (err) {
        delete & tt;
        }
        else {
        if (line.distance(tt) < _szLinkDistance) {
            tt.link(& t);
            forNewLine.append(tt);
        }
        else {
            delete & tt;
        }
        }
    }
    }
 
    //...Create new line...
#ifdef TRKRECO_DEBUG
    std::cout << "    ... creating a new line" << std::endl;
#endif  
    unsigned nNewLine = forNewLine.length();
    TMLine newLine(forNewLine);
    static TRobustLineFitter fitter("abc");
    int err = fitter.fit(newLine);
    // int err = newLine.fit();
    if (err < 0) {
    HepAListDeleteAll(forNewLine);
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    ... 2nd linear fit failure. nLinks(";
    std::cout << forNewLine.length() << ")" << std::endl;
#endif  
    return NULL;
    }
 
#ifdef TRKRECO_WINDOW
    sz.append(newLine, leda_green);
    sz.wait();
#endif
 
#ifdef TRKRECO_DEBUG_DETAIL
    Dump(forNewLine, "sort hits stereo", "        ");
#endif
 
    //...Decide left/right...
    AList<TMLink> toRemove;
    TMLink * last = NULL;
    for (unsigned i = 0; i < nNewLine; i++) {
    if (last == NULL) last = forNewLine[i]->link();
    else {
        if (last == forNewLine[i]->link()) {
        if (newLine.distance(* forNewLine[i - 1]) >
            newLine.distance(* forNewLine[i]))
            toRemove.append(forNewLine[i - 1]);
        else
            toRemove.append(forNewLine[i]);
        last = NULL;
        }
        else {
        last = forNewLine[i]->link();
        }
    }
    }
    forNewLine.remove(toRemove);
    nNewLine = forNewLine.length();
 
#ifdef TRKRECO_DEBUG_DETAIL
    Dump(toRemove, "sort hits stereo", "      x ");
    Dump(forNewLine, "sort hits stereo", "        ");
#endif
 
    //...3D fit...
    for (unsigned i = 0; i < nNewLine; i++)
    track.append(* forNewLine[i]->link());
    Vector a(5);
    a = track.helix().a();
    a[3] = newLine.b();
    a[4] = track.charge() * newLine.a();
    track._helix->a(a);
 
    //...Refine...
    AList<TMLink> bad;
    err = _fitter.fit(track);
    track.refine(bad, _maxSigmaStereo * 100.);
    err = _fitter.fit(track);
    track.refine(bad, _maxSigmaStereo * 10.);
    err = _fitter.fit(track);
    track.refine(bad, _maxSigmaStereo);
    err = _fitter.fit(track);
 
#ifdef TRKRECO_WINDOW
    sz.text("stereo finished");
    sz.oneShot(track, leda_blue);
#endif
 
    //...Termination...
    HepAListDeleteAll(toRemove);
    HepAListDeleteAll(forNewLine);
    return & track;
}

buildStereoNew()

TTrack * TBuilder::buildStereoNew	(	const TTrack &	t,
		AList< TSegment > &	goodSegments,
		AList< TSegment > &	badSegments
	)		const

Definition at line 1535 of file TBuilder.cxx.

                                                  {
 
#ifdef TRKRECO_DEBUG
    std::cout << "... building stereo by segments(new) : # of segments = ";
    std::cout << segments.length() << std::endl;
    for (unsigned i = 0; i < segments.length(); i++)
    segments[i]->dump("hits sort flag", "    ");
#endif
#ifdef TRKRECO_WINDOW
    sz.clear();
#endif
    TTrack * bestCandidate = NULL;
    AList< AList<TMLink> > poorSeeds;
    bool ok = initializeForStereo(t, segments, badSegments);
    unsigned nSuperLayers = _nSuperLayers + 1;
//        std::cout<<"nSuperLayers in buildStereoNew:"<<_nSuperLayers<<endl;  
//    cout<<"ok of initializeForStereo: "<<ok<<endl;
 
    if (! ok) goto endOfBuilding;
    //...Main loop...
    while (--nSuperLayers) {
 
    //...Initial line search...
    AList<TMLine> initialLines = searchInitialLines(nSuperLayers);
#ifdef TRKRECO_WINDOW
    sz.clear();
    sz.skip(false);
    sz.mode(2);
    sz.appendSz(t, segments, leda_black);
    AList<TSegment> tmps = segments;
    std::string s = "# of segments = " + itostring(int(segments.length()));
    sz.appendSz(t, _allLinks, leda_black);
    sz.append(_forLine, leda_brown);
    s = "nSprLyr=" + itostring(int(nSuperLayers)) +
        ", # of initial lines = " + itostring(int(initialLines.length()));
    for (unsigned i = 0; i < initialLines.length(); i++)
        sz.append(* initialLines[i], leda_red);
    sz.text(s);
    sz.wait();
#endif
//  cout<<"length of initialLines = "<<initialLines.length()<<endl;
    if (initialLines.length() == 0) continue;
    //...Line loop...
    bool found = false;
    unsigned nInitialLines = initialLines.length();
//        cout<<"nInitialLines:"<<nInitialLines<<endl;
    
    for (unsigned i = 0; i < nInitialLines; i++) {
 
        //...Linear fit...
        const TMLine & line = * initialLines[i];
        TMLine newLine = searchLine(line);
 
        //...Skip if this is a seed of the poor result...
        if (poorSeeds.length()) {
        bool poorCase = false;
        for (unsigned j = 0; j < poorSeeds.length(); j++) {
            if (poorSeeds[j]->length() == newLine.links().length()) {
            AList<TMLink> tmp = * poorSeeds[j];
            tmp.remove(newLine.links());
            if (tmp.length() == 0) {
#ifdef TRKRECO_DEBUG_DETAIL
                std::cout << "    ... This is a poor seed :"
                      << " skipped"
                      << " : # of poor seeds = "
                      << poorSeeds.length() << std::endl;
#endif
                poorCase = true;
                break;
            }
            }
        }
        if (poorCase) continue;
        }
        
        //...Is this a good line...
        if (! stereoQuality(newLine.links()))
        continue;
 
        //...3D fit...
        TTrack * t3d = new TTrack(t);
        t3d = build(* t3d, newLine);
        if (t3d == NULL) continue;
 
        //...Check quality...
        unsigned quality = stereoQuality(t3d->links());
//      cout<<"point3, quality:"<<quality<<endl;
 
        //...Best case...
        if (quality == 2) {
#ifdef TRKRECO_WINDOW
        sz.text("stereo finished");
        sz.oneShot(* t3d, leda_blue);
#endif
        if (bestCandidate) delete bestCandidate;
        bestCandidate = t3d;
        found = true;
        break;
        }
        
        //...Poor case...
        AList<TMLink> * tmpL = new AList<TMLink>();
        tmpL->append(newLine.links());
        poorSeeds.append(tmpL);
        if (quality == 0) {
#ifdef TRKRECO_WINDOW
        sz.text("this candidate discarded");
        sz.oneShot(* t3d, leda_black);
#endif
        delete t3d;
        continue;
        }
        
        //...Not enough...
        if (bestCandidate) {
        if (bestCandidate->cores().length() <
            t3d->cores().length()) {
#ifdef TRKRECO_WINDOW
            sz.text("new candidate");
            sz.append(* bestCandidate, leda_brown);
            sz.oneShot(* t3d, leda_green);
            sz.remove(* bestCandidate);
#endif
            delete bestCandidate;
            bestCandidate = t3d;
        }
        else {
#ifdef TRKRECO_WINDOW
            sz.text("this candidate discarded");
            sz.oneShot(* t3d, leda_black);
#endif
            delete t3d;
        }
        }
        else {
        bestCandidate = t3d;
#ifdef TRKRECO_WINDOW
        sz.text("new candidate");
        sz.oneShot(* bestCandidate, leda_green);
        sz.remove(* t3d);
#endif
        }
    }
    
    //...Termination of a loop...
    HepAListDeleteAll(initialLines);
    if (found) break;
    }
 
endOfBuilding:
    _allLinks.removeAll();
    for (unsigned i = 0; i < 6; i++)
    HepAListDeleteAll(_links[i]);
    HepAListDeleteAll(_forLine);
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    ... # of poor seeds = " << poorSeeds.length() << std::endl;
#endif
#ifdef TRKRECO_WINDOW
    if (bestCandidate == NULL) {
    sz.text("3D failed");
    sz.wait();
    }
#endif
        
    //...Check used segments...
    if (bestCandidate) {
    AList<TMLink> usedLinks = StereoHits(bestCandidate->links());
    for (unsigned i = 0; i < segments.length(); i++) {
        TSegment & segment = * segments[i];
        AList<TMLink> used = Links(segment, * bestCandidate);
        if (used.length()) {
        bestCandidate->segments().append(segment);
        segment.tracks().append(bestCandidate);
        }
    }
    }
    
    if (poorSeeds.length())
    HepAListDeleteAll(poorSeeds);
 
#ifdef TRASA_DEBUG_DETAIL
    if (bestCandidate == NULL)
    std::cout << "... building stereo(new) failed" << std::endl;
    else
    std::cout << "... building stereo(new) ok" << std::endl;
#endif
    return bestCandidate;
}

buildStereoOld()

TTrack * TBuilder::buildStereoOld	(	TTrack &	t,
		const AList< TMLink > &	allLinks,
		const AList< TMLink > &	isolatedLinks
	)		const

dump()

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

dumps debug information.

Definition at line 62 of file TBuilder.cxx.

                                                                 {
}

Referenced by buildRphi(), buildStereo(), buildStereoNew(), and salvage().

initialLine()

TMLine * TBuilder::initialLine	(	const TTrack &	t,
		AList< TSegment > &	list
	)		const

makes a line.

Definition at line 315 of file TBuilder.cxx.

                                                                    {
 
    //...Select good segments in acceptance...
    AList<TSegment> bad;
    AList<TMLink> lList;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * new TMLink();
    int err = t.szPosition(* list[i], l);
    if (err) {
        bad.append(list[i]);
        delete & l;
        continue;
    }
 
    const HepPoint3D & p = l.position();
    if (p.y() > l.wire()->forwardPosition().z() ||
        p.y() < l.wire()->backwardPosition().z()) {
        bad.append(list[i]);
        delete & l;
        continue;
    }
    lList.append(l);
    }
    list.remove(bad);
    n = list.length();
 
    //...Cal. expected # of super layers...
    unsigned nAMaxSL = 0;
    unsigned nAMinSL = 10;
    const AList<TSegment> & segments = t.segments();
    unsigned nA = segments.length();
    for (unsigned i = 0; i < nA; i++) {
    unsigned sl = segments[i]->superLayerId();
    if (sl > nAMaxSL) nAMaxSL = sl;
    if (sl < nAMinSL) nAMinSL = sl;
    }
    unsigned nExpected = (nAMaxSL - nAMinSL) / 2;
#ifdef TRKRECO_DEBUG
    std::cout << "    ... initialLine : axial super layer usage = " << nAMinSL;
    std::cout << " ~ " << nAMaxSL << std::endl;
    std::cout << "                    : expected stereo super layers = ";
    std::cout << nExpected << std::endl;
#endif
 
    TMLine * line = NULL;
    AList<TSegment> list0 = list;
    AList<TMLink> lList0 = lList;
    while (1) {
    bool last = (NSuperLayers(lList) <= nExpected);
 
    TMLine * line = initialLine2(t, lList);
    if (line) {
        AList<TSegment> tmp = selectStereoSegment(* line, list0, lList0);
#ifdef TRKRECO_DEBUG
        std::cout << "    ... initialLine2 : # of selected segments = ";
        std::cout << tmp.length() << std::endl;
#endif
        if ((tmp.length() >= nExpected) ||
        (tmp.length() >= 4)         ||
        last) {
        list = tmp;
        return line;
        }
        delete line;
    }
 
    line = initialLine1(t, list, lList);
    if (line) {
        AList<TSegment> tmp = selectStereoSegment(* line, list0, lList0);
        bool ok = (tmp.length() >= nExpected);
#ifdef TRKRECO_DEBUG
        std::cout << "    ... initialLine2 : # of selected segments = ";
        std::cout << tmp.length() << " : ok, last = ";
        std::cout << ok << ", " << last << std::endl;
#endif
        if (ok) {
        list = tmp;
        return line;
        }
        removeFarSegment(* line, list, lList);
        delete line;
    }
    else {
        return NULL;
    }
    }
}

Referenced by searchLine().

initialLine1()

TMLine * TBuilder::initialLine1	(	const TTrack &	t,
		const AList< TSegment > &	list,
		const AList< TMLink > &	lList
	)		const

Definition at line 247 of file TBuilder.cxx.

                                                  {
 
    //...Check input...
    if (list.length() < 2) return NULL;
 
    //...Check super layer pattern...
    AList<TSegment> sl[5];
    AList<TMLink> tl[5];
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i ++) {
    unsigned j = list[i]->superLayerId() / 2;
    sl[j].append(list[i]);
    tl[j].append(lList[i]);
    }
#ifdef TRKRECO_DEBUG
    std::cout << "        ... initialLine1 : super layer ptn = ";
    for (unsigned i = 0; i < 5; i++) std::cout << sl[i].length();
    std::cout << std::endl;
#endif
 
    //...Count single segment layer...
    unsigned nSingle = 0;
    for (unsigned i = 0; i < 5; i++)
    if (sl[i].length() == 1) ++nSingle;
#ifdef TRKRECO_DEBUG
    std::cout << "        ... # of single segment layer = " << nSingle << std::endl;
#endif
    if (nSingle < 2) return NULL;
 
    //...Happy case...
    AList<TSegment> bestCombination;
    AList<TMLink> forLine;
    for (unsigned i = 0; i < 5; i++) {
    if (sl[i].length() != 1) continue;
    bestCombination.append(sl[i]);
    forLine.append(tl[i]);
    }
    TMLine & line = * new TMLine(forLine);
    int err = line.fit();
    if (err) {
    delete & line;
    return NULL;    
    }
    return & line;
}

Referenced by initialLine(), and initialLineOld().

initialLine2()

TMLine * TBuilder::initialLine2	(	const TTrack &	t,
		const AList< TMLink > &	lList
	)		const

Definition at line 296 of file TBuilder.cxx.

                                                  {
#ifdef TRKRECO_DEBUG
    std::cout << "        ... initlialLine2 : # of links = " << lList.length();
    std::cout << std::endl;
#endif
 
    //...Static object...
    static TRobustLineFitter fitter("can you work?");
    TMLine * line = new TMLine(lList);
    int err = fitter.fit(* line);
    if (err) {
    delete line;
    return NULL;
    }
    return line;
}

Referenced by initialLine(), and initialLineOld().

initialLineOld()

TMLine * TBuilder::initialLineOld	(	const TTrack &	t,
		AList< TSegment > &	list
	)		const

Definition at line 405 of file TBuilder.cxx.

                                                                       {
 
    //...Select good segments in acceptance...
    AList<TSegment> bad;
    AList<TMLink> lList;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * new TMLink();
    int err = t.szPosition(* list[i], l);
    if (err) {
        bad.append(list[i]);
        delete & l;
        continue;
    }
 
    const HepPoint3D & p = l.position();
    if (p.y() > l.wire()->forwardPosition().z() ||
        p.y() < l.wire()->backwardPosition().z()) {
        bad.append(list[i]);
        delete & l;
        continue;
    }
    lList.append(l);
    }
    list.remove(bad);
    n = list.length();
 
    //...Cal. expected # of super layers...
    unsigned nAMaxSL = 0;
    unsigned nAMinSL = 10;
    const AList<TSegment> & segments = t.segments();
    unsigned nA = segments.length();
    for (unsigned i = 0; i < nA; i++) {
    unsigned sl = segments[i]->superLayerId();
    if (sl > nAMaxSL) nAMaxSL = sl;
    if (sl < nAMinSL) nAMinSL = sl;
    }
    unsigned nExpected = (nAMaxSL - nAMinSL) / 2;
#ifdef TRKRECO_DEBUG
    std::cout << "    ... initialLine : axial super layer usage = " << nAMinSL;
    std::cout << " ~ " << nAMaxSL << std::endl;
    std::cout << "                    : expected stereo super layers = ";
    std::cout << nExpected << std::endl;
#endif
 
    TMLine * line = NULL;
    AList<TSegment> list0 = list;
    AList<TMLink> lList0 = lList;
    while (1) {
    TMLine * line = initialLine1(t, list, lList);
    if (line) {
        AList<TSegment> tmp = selectStereoSegment(* line, list0, lList0);
#ifdef TRKRECO_DEBUG
        std::cout << "    ... initialLine1 : # of selected segments = ";
        std::cout << tmp.length() << std::endl;
#endif
        if (tmp.length() >= nExpected) {
        lList0.remove(line->links());
        HepAListDeleteAll(lList0);
        list = tmp;
        return line;
        }
        delete line;
    }
 
    line = initialLine2(t, lList);
    if (line) {
        AList<TSegment> tmp = selectStereoSegment(* line, list0, lList0);
        bool ok = (tmp.length() >= nExpected);
        bool last = (NSuperLayers(lList) <= nExpected);
#ifdef TRKRECO_DEBUG
        std::cout << "    ... initialLine2 : # of selected segments = ";
        std::cout << tmp.length() << " : ok, last = ";
        std::cout << ok << ", " << last << std::endl;
#endif
        if (ok || last) {
        lList0.remove(line->links());
        HepAListDeleteAll(lList0);
        list = tmp;
        return line;
        }
        removeFarSegment(* line, list, lList);
        delete line;
    }
    else {
        HepAListDeleteAll(lList0);
        return NULL;
    }
    }
}

Referenced by buildStereo().

name()

const std::string & TBuilder::name ( void  ) const

inline

returns name.

Definition at line 146 of file TBuilder.h.

                         {
    return _name;
}

removeFarSegment()

void TBuilder::removeFarSegment	(	const TMLine &	line,
		AList< TSegment > &	list,
		AList< TMLink > &	szList
	)		const

Definition at line 510 of file TBuilder.cxx.

                                             {
    unsigned n = list.length();
    float maxDistance = 0.;
    unsigned maxId = 0;
    for (unsigned i = 0; i < n; i++) {
    float distance = line.distance(* szList[i]);
    if (distance > maxDistance) maxId = i;
    }
    list.remove(maxId);
    szList.remove(maxId);
}

Referenced by initialLine(), and initialLineOld().

salvage() [1/2]

void TBuilder::salvage	(	TTrack &	t,
		AList< TMLink > &	hits
	)		const

salvages hits.

Definition at line 211 of file TBuilder.cxx.

                                                        {
 
#ifdef TRKRECO_DEBUG
    std::cout << "... salvaging(TBuilder) : # of given hits=" << hits.length();
    std::cout << ", salvage level=" << _salvageLevel << std::endl;
    Dump(hits, "hits sort flag");
#endif
 
    unsigned nHits = hits.length();
    if (nHits == 0) return;
 
    //...Try to append this hit...
    t.appendByApproach(hits, _salvageLevel);
    _fitter.fit(t);
}

Referenced by salvage().

salvage() [2/2]

void TBuilder::salvage	(	TTrack &	t,
		AList< TSegment > &	segments
	)		const

Definition at line 228 of file TBuilder.cxx.

                                                          {
 
#ifdef TRKRECO_DEBUG
    std::cout << "... salvaging by segments : # of segments = ";
    std::cout << list.length() << std::endl;
    for (unsigned i = 0; i < list.length(); i++)
    list[i]->dump("hits sort flag", "    ");
#endif
 
    //...Pick up links...
    AList<TMLink> links;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++)
    links.append(list[i]->links());
 
    salvage(t, links);
}

searchInitialLines()

AList< TMLine > TBuilder::searchInitialLines

(

unsigned

nSuperLayers

)

const

Definition at line 724 of file TBuilder.cxx.

                                                          {
#ifdef TRKRECO_DEBUG
    std::cout << "    ... searchInitialLines : # of segments in sl : ";
    for (unsigned i = 0; i < 5; i++)
    std::cout << _links[i].length() << ",";
    std::cout << std::endl
     << "                           : max # of super layers="
     << nSuperLayerMax << std::endl;
#endif
 
    AList<TMLine> lines;
    if (nSuperLayerMax > 5)
        lines.append(searchLines6());
    else if (nSuperLayerMax > 4)
    lines.append(searchLines5());
    else if (nSuperLayerMax > 3)
    lines.append(searchLines4());
    else if (nSuperLayerMax > 2)
    lines.append(searchLines3());
    else if (nSuperLayerMax > 1)
    lines.append(searchLines2());
//      else
//      lines.append(searchLines1());
 
    lines.sort(SortByB);
    return lines;
}

Referenced by buildStereoNew().

searchLine()

TMLine TBuilder::searchLine

(

const TMLine &

initialLine

)

const

Definition at line 1176 of file TBuilder.cxx.

                                                     {
#ifdef TRKRECO_WINDOW
    std::string old = sz.text();
#endif
 
    static float _szLinkMaxDistance = 5;
 
    unsigned nLinks = _forLine.length();
    TMLine line0 = initialLine;
    line0.remove(initialLine.links());
#ifdef TRKRECO_DEBUG
    line0.fitted(true);
#endif    
    unsigned nGoodLinks0 = 0;
    AList<TMLink> targets;
    unsigned nLoops = 0;
    while (1) {
    ++nLoops;
 
    unsigned nGoodLinks = 0;
    for (unsigned i = 0; i < nLinks; i++) {
        float distance = line0.distance(* _forLine[i]);
        if (distance < _szLinkDistance)
        targets.append(_forLine[i]);
        if (distance < _szLinkMaxDistance) {
        ++nGoodLinks;
        }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "    ... searchLine : # of close hits(last) = " << nGoodLinks0
         << ", # of close hits = " << nGoodLinks << std::endl;
#endif
 
    //...Check condition...
    if (nGoodLinks0 == nGoodLinks) break;
 
    //...Do fit...
    TMLine newLine(targets);
    static TRobustLineFitter fitter("robust fitter");
    int err = fitter.fit(newLine);
    if (err) {
#ifdef TRKRECO_WINDOW
        std::cout << "    ... searchLine : failed to fit" << std::endl;
        std::string s = "line search failed";
        sz.text(s);
        sz.wait();
#endif
        break;
    }
 
    //...To protect infinite loop...
    if (nLoops > 10) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "        reached to max # of loops(10) : break";
#endif
        break;
    }
 
#ifdef TRKRECO_WINDOW
    std::string ss = ", # of close hits(last) = " + itostring(int(nGoodLinks0))
        + ", # of close hits = " + itostring(int(nGoodLinks));
    sz.append(line0, leda_brown);
    sz.append(newLine, leda_green);
    sz.text(old + ", nLoops = " + itostring(int(nLoops)) + ss);
    sz.wait();
    sz.remove(line0);
    sz.remove(newLine);
#endif
 
    line0 = newLine;
    nGoodLinks0 = nGoodLinks;
    targets.removeAll();
    }
 
    return line0;
}

Referenced by buildStereoNew().

searchLines1()

AList< TMLine > TBuilder::searchLines1

(

void

)

const

Definition at line 1155 of file TBuilder.cxx.

                                 {
    AList<TMLine> lines;
 
    TMLine & line = * new TMLine(_forLine);
    int err = line.fit();
    if (err) {
    delete & line;
    }
    else {
    lines.append(line);
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines1 : # of lines found = "
     << lines.length() << std::endl;
#endif
 
    return lines;
}

searchLines2()

AList< TMLine > TBuilder::searchLines2

(

void

)

const

Definition at line 1069 of file TBuilder.cxx.

                                 {
    AList<TMLine> lines;
 
    const AList<TMLink> * l[2];
    for (unsigned skip = 0; skip < 15; skip++) {
    if (skip == 0) {
        l[0] = & _links[0];
        l[1] = & _links[1];
    }
    else if (skip == 1) {
        l[1] = & _links[2];
    }
    else if (skip == 2) {
        l[1] = & _links[3];
    }
    else if (skip == 3) {
        l[1] = & _links[4];
    }
    else if (skip == 4) {
        l[1] = & _links[5];
    }
    else if (skip == 5) {
        l[0] = & _links[1];
        l[1] = & _links[2];
    }
    else if (skip == 6) {
        l[1] = & _links[3];
    }
    else if (skip == 7) {
        l[1] = & _links[4];
    }
    else if (skip == 8) {
        l[1] = & _links[5];
    }
    else if (skip == 9) {
        l[0] = & _links[2];
        l[1] = & _links[3];
    }
    else if (skip == 10) {
        l[1] = & _links[4];
    }
    else if (skip == 11) {
        l[1] = & _links[5];
    }
    else if (skip == 12) {
        l[0] = & _links[3];
        l[1] = & _links[4];
    }
    else if (skip == 13) {
        l[1] = & _links[5];
    }
    else if (skip == 14) {
        l[0] = & _links[4];
    }
 
    unsigned n[2];
    for (unsigned i = 0; i < 2; i++)
        n[i] = l[i]->length();
 
    for (unsigned i0 = 0; i0 < n[0]; i0++) {
        for (unsigned i1 = 0; i1 < n[1]; i1++) {
        AList<TMLink> forLine;
        forLine.append((* l[0])[i0]);
        forLine.append((* l[1])[i1]);
 
        TMLine & line = * new TMLine(forLine);
        int err = line.fit();
        if (err) {
            delete & line;
            continue;
        }
 
        lines.append(line);
        }
    }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines2 : # of lines found = "
     << lines.length() << std::endl;
#endif
 
    return lines;
}

Referenced by searchInitialLines().

searchLines3()

AList< TMLine > TBuilder::searchLines3

(

void

)

const

Definition at line 964 of file TBuilder.cxx.

                                 {
    AList<TMLine> lines;
 
    const AList<TMLink> * l[3];
    for (unsigned skip = 0; skip < 20; skip++) {
    if (skip == 0) {
        l[0] = & _links[0];
        l[1] = & _links[1];
        l[2] = & _links[2];
    }
    else if (skip == 1) {
        l[2] = & _links[3];
    }
    else if (skip == 2) {
        l[2] = & _links[4];
    }
    else if (skip == 3) {
        l[2] = & _links[5];
    }
    else if (skip == 4) {
        l[1] = & _links[2];
    }
    else if (skip == 5) {
        l[2] = & _links[4];
    }
    else if (skip == 6) {
        l[2] = & _links[3];
    }
    else if (skip == 7) {
        l[1] = & _links[3];
        l[2] = & _links[4];
    }
    else if (skip == 8) {
        l[2] = & _links[5];
    }
    else if (skip == 9) {
        l[1] = & _links[4];
    }
    else if (skip == 10) {
        l[0] = & _links[1];
    }
    else if (skip == 11) {
        l[1] = & _links[3];
    }
    else if (skip == 12) {
        l[2] = & _links[4];
    }
    else if (skip == 13) {
        l[1] = & _links[2];
    }
    else if (skip == 14) {
        l[2] = & _links[3];
    }
    else if (skip == 15) {
        l[2] = & _links[5];
    }
    else if (skip == 16) {
        l[0] = & _links[2];
        l[1] = & _links[3];
    }
    else if (skip == 17) {
        l[2] = & _links[3];
    }
    else if (skip == 18) {
        l[1] = & _links[4];
        l[2] = & _links[5];
    }
    else if (skip == 19) {
        l[0] = & _links[3];
    }
 
    unsigned n[3];
    for (unsigned i = 0; i < 3; i++)
        n[i] = l[i]->length();
 
    for (unsigned i0 = 0; i0 < n[0]; i0++) {
        for (unsigned i1 = 0; i1 < n[1]; i1++) {
        for (unsigned i2 = 0; i2 < n[2]; i2++) {
            AList<TMLink> forLine;
            forLine.append((* l[0])[i0]);
            forLine.append((* l[1])[i1]);
            forLine.append((* l[2])[i2]);
 
            TMLine & line = * new TMLine(forLine);
            int err = line.fit();
            if (err) {
            delete & line;
            continue;
            }
 
            lines.append(line);
        }
        }
    }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines3 : # of lines found = "
     << lines.length() << std::endl;
#endif
 
    return lines;
}

Referenced by searchInitialLines().

searchLines4()

AList< TMLine > TBuilder::searchLines4

(

void

)

const

Definition at line 867 of file TBuilder.cxx.

                                 {
    AList<TMLine> lines;
 
    const AList<TMLink> * l[4];
    for (unsigned skip = 0; skip < 15; skip++) {
    if (skip == 0) {
        l[0] = & _links[0];
        l[1] = & _links[1];
        l[2] = & _links[2];
        l[3] = & _links[3];
    }
    else if (skip == 1) {
        l[3] = & _links[4];
    }
    else if (skip == 2) {
        l[3] = & _links[5];
    }
    else if (skip == 3) {
        l[2] = & _links[3];
        l[3] = & _links[4];
    }
    else if (skip == 4) {
        l[3] = & _links[5];
    }
    else if (skip == 5) {
        l[2] = & _links[4];
    }
    else if (skip == 6) {
        l[1] = & _links[2];
        l[2] = & _links[3];
        l[3] = & _links[4];
    }
    else if (skip == 7) {
        l[3] = & _links[5];
    }
    else if (skip == 8) {
        l[2] = & _links[4];
    }
    else if (skip == 9) {
        l[1] = & _links[3];
    }
    else if (skip == 10) {
        l[0] = & _links[1];
        l[1] = & _links[2];
        l[2] = & _links[3];
        l[3] = & _links[4];
    }
    else if (skip == 11) {
        l[3] = & _links[5];
    }
    else if (skip == 12) {
        l[2] = & _links[4];
    }
    else if (skip == 13) {
        l[1] = & _links[3];
    }
    else if (skip == 14) {
        l[0] = & _links[2];
    }
 
    unsigned n[4];
    for (unsigned i = 0; i < 4; i++)
        n[i] = l[i]->length();
 
    for (unsigned i0 = 0; i0 < n[0]; i0++) {
        for (unsigned i1 = 0; i1 < n[1]; i1++) {
        for (unsigned i2 = 0; i2 < n[2]; i2++) {
            for (unsigned i3 = 0; i3 < n[3]; i3++) {
            AList<TMLink> forLine;
            forLine.append((* l[0])[i0]);
            forLine.append((* l[1])[i1]);
            forLine.append((* l[2])[i2]);
            forLine.append((* l[3])[i3]);
 
            TMLine & line = * new TMLine(forLine);
            int err = line.fit();
            if (err) {
                delete & line;
                continue;
            }
 
            lines.append(line);
            }
        }
        }
    }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines4 : # of lines found = "
     << lines.length() << std::endl;
#endif
 
    return lines;
}

Referenced by searchInitialLines().

searchLines5()

AList< TMLine > TBuilder::searchLines5

(

void

)

const

Definition at line 798 of file TBuilder.cxx.

                                 {
    AList<TMLine> lines;
 
    const AList<TMLink> * l[5];
    for (unsigned skip = 0; skip < 6; skip++) {
    if (skip == 0) {
        l[0] = & _links[1];
        l[1] = & _links[2];
        l[2] = & _links[3];
        l[3] = & _links[4];
        l[4] = & _links[5];
    }
    else if (skip == 1) {
        l[0] = & _links[0];
    }
    else if (skip == 2) {
        l[1] = & _links[1];
    }
    else if (skip == 3) {
        l[2] = & _links[2];
    }
    else if (skip == 4) {
        l[3] = & _links[3];
    }
    else if (skip == 5) {
        l[4] = & _links[4];
    }
 
    unsigned n[5];
    for (unsigned i = 0; i < 5; i++)
    n[i] = l[i]->length();
 
    for (unsigned i0 = 0; i0 < n[0]; i0++) {
    for (unsigned i1 = 0; i1 < n[1]; i1++) {
        for (unsigned i2 = 0; i2 < n[2]; i2++) {
        for (unsigned i3 = 0; i3 < n[3]; i3++) {
            for (unsigned i4 = 0; i4 < n[4]; i4++) {
 
            AList<TMLink> forLine;
            forLine.append((* l[0])[i0]);
            forLine.append((* l[1])[i1]);
            forLine.append((* l[2])[i2]);
            forLine.append((* l[3])[i3]);
            forLine.append((* l[4])[i4]);
 
            TMLine & line = * new TMLine(forLine);
            int err = line.fit();
            if (err) {
                delete & line;
                continue;
            }
 
            lines.append(line);
            }
        }
        }
    }
    }
    }
 
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines5 : # of lines found = "
     << lines.length() << std::endl;
#endif
 
    return lines;
}

Referenced by searchInitialLines().

searchLines6()

AList< TMLine > TBuilder::searchLines6

(

void

)

const

Definition at line 753 of file TBuilder.cxx.

                                 {
    unsigned n[6];
    for (unsigned i = 0; i < 6; i++)
        n[i] = _links[i].length();
                                                                                                                             
    AList<TMLine> lines;
    for (unsigned i0 = 0; i0 < n[0]; i0++) {
        for (unsigned i1 = 0; i1 < n[1]; i1++) {
            for (unsigned i2 = 0; i2 < n[2]; i2++) {
                for (unsigned i3 = 0; i3 < n[3]; i3++) {
                    for (unsigned i4 = 0; i4 < n[4]; i4++) {
            for (unsigned i5 = 0; i5 < n[5]; i5++) {
                                                                                                                             
                        AList<TMLink> forLine;
                        forLine.append(_links[0][i0]);
                        forLine.append(_links[1][i1]);
                        forLine.append(_links[2][i2]);
                        forLine.append(_links[3][i3]);
                        forLine.append(_links[4][i4]);
                        forLine.append(_links[5][i5]);
                                                                                                                             
                        TMLine & line = * new TMLine(forLine);
                        int err = line.fit();
                        if (err) {
                            delete & line;
                            continue;
                        }
                                                                                                                             
                        lines.append(line);
            }
                    }
                }
            }
        }
    }
                                                                                                                             
#ifdef TRKRECO_DEBUG
    std::cout << "        ... searchLines5 : # of lines found = "
         << lines.length() << std::endl;
#endif
                                                                                                                             
    return lines;
}

Referenced by searchInitialLines().

selectStereoSegment()

AList< TSegment > TBuilder::selectStereoSegment	(	const TMLine &	line,
		const AList< TSegment > &	list,
		const AList< TMLink > &	szList
	)		const

Definition at line 497 of file TBuilder.cxx.

                                                      {
    AList<TSegment> outList;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    float distance = line.distance(* szList[i]);
    if (distance < _szSegmentDistance) outList.append(list[i]);
    }
    return outList;
}

Referenced by initialLine(), and initialLineOld().

---

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

• TBuilder.h
• TBuilder.cxx