TBuilderCosmic Class Reference

A class to build a cosmic track. More...

#include <TBuilderCosmic.h>

Inheritance diagram for TBuilderCosmic:

List of all members.

Public Member Functions
	TBuilderCosmic (const std::string &name, float salvageLevel)
	Constructor.
virtual	~TBuilderCosmic ()
	Destructor.
TTrack *	buildStereo (TTrack &track, const AList< TMLink > &) const
	appends stereo hits to a track.
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.
const TMSelector &	trackSelector (void) const
	returns a track selector.
virtual const TMSelector &	trackSelector (const TMSelector &)
	sets a track selector.
TTrack *	buildRphi (const AList< TMLink > &) const
	builds a r/phi track from TMLinks or from Segments.
TTrack *	buildStereo0 (TTrack &track, const AList< TMLink > &) const
	appends stereo hits to a track. (old version)
void	appendClusters (TTrack &track, const AList< TMLink > &) const
	appends TMLinks in a list.
void	salvage (TTrack &track, AList< TMLink > &list) const
	salvages links in a list. Used links will be removed from a list.
virtual int	fit (TTrackBase &) const
	fits a track using a private fitter.
Protected Attributes
TMSelector	_circleSelector
TMSelector	_trackSelector
TMSelector	_lineSelector
float	_stereoZ3
float	_stereoZ4
float	_stereoChisq3
float	_stereoChisq4
float	_stereoMaxSigma
Private Attributes
TCosmicFitter	_fitter

---

Detailed Description

A class to build a cosmic track.

Definition at line 36 of file TBuilderCosmic.h.

---

Constructor & Destructor Documentation

TBuilderCosmic::TBuilderCosmic	(	const std::string &	name,
		float	salvageLevel
	)

Constructor.

Definition at line 24 of file TBuilderCosmic.cxx.

: TBuilder0(name, salvageLevel), _fitter("TBuilderCosmic Fitter") {
}

TBuilderCosmic::~TBuilderCosmic

(

)

[virtual]

Destructor.

Definition at line 28 of file TBuilderCosmic.cxx.

                                {
}

---

Member Function Documentation

void TBuilder0::appendClusters	(	TTrack &	track,
		const AList< TMLink > &
	)			const [inherited]

appends TMLinks in a list.

Definition at line 1417 of file TBuilder0.cxx.

References TBuilder0::_trackSelector, TTrackBase::appendByApproach(), TBuilder0::fit(), TMSelector::maxSigma(), and TTrackBase::refine().

Referenced by TConformalFinder0::appendClusters2().

                                                           {

    AList<TMLink> tmp = list;

    //...Append them...
    track.appendByApproach(tmp, _trackSelector.maxSigma() * 2. / 3.);

    //...Refine it...
    AList<TMLink> bad;
    fit(track);
    track.refine(bad, _trackSelector.maxSigma());
}

TTrack * TBuilder0::buildRphi

(

const AList< TMLink > &

)

const [inherited]

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

Definition at line 83 of file TBuilder0.cxx.

References TBuilder0::_circleSelector, TBuilder0::_name, TBuilder0::_trackSelector, check_raw_filter::dist, Dump(), showlog::err, TBuilder0::fit(), TTrackBase::fit(), genRecEmupikp::i, ganga-rec::j, TMSelector::maxSigma(), NLayers(), TMSelector::nLinks(), NSuperLayers(), TMSelector::preSelect(), SameSuperLayer(), TMSelector::select(), TBuilder0::selectHits(), SuperLayer(), t(), and WireHitFittingValid.

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

                                                     {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << " ... building a rphi track" << std::endl;
    std::cout << _name << " ... selecting good hits" << std::endl;
//      TTrackBase tmp;
//      tmp.append(list);
//      rphiWindow->append(tmp, leda_red);
#endif

    //...Check # of links...
    if (list.length() < _circleSelector.nLinks()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << " ... rejected by nLinks(";
        std::cout << list.length() << ") < ";
        std::cout << _circleSelector.nLinks() << std::endl;
#endif  
    }

    //...Select core hits...
    AList<TMLink> cores = list;
    selectHits(cores);
    if (cores.length() < 5) cores = list;
    //cout<<"list: "<<list.length()<<"  core: "<<cores.length()<<endl;

    //...Core check...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << " ... checking cores : cores=" << std::endl;
    Dump(cores, "detail");
#endif
    unsigned sLinks = SuperLayer(list);
    unsigned sUsed = SuperLayer(cores);
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    super layer ptn=" << sLinks;
    std::cout << ",super layer used=" << sUsed << std::endl;
#endif
    if (sLinks != sUsed) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "    ... appending hits to cores" << std::endl;
#endif
        unsigned diff = sLinks - sUsed;
        unsigned asl[5] = {2, 3, 4, 9, 10};
        for (unsigned j = 0; j < 5; j++) {
            //if (diff & (1 << ((5 - j) * 2))) {
            if (diff & (1 << asl[5 - j])) {
#ifdef TRKRECO_DEBUG_DETAIL
                std::cout << "    super layer " << (5 - j) * 2 << "searching";
                std::cout << std::endl;
#endif
                AList<TMLink> links = SameSuperLayer(list, asl[5 - j]);
                TMLink * best = NULL;
                double bestD = 999.;
                for (unsigned i = 0; i < links.length(); i++) {
                    double dist = links[i]->hit()->drift();
                    if (dist < 0.02) {
#ifdef TRKRECO_DEBUG_DETAIL
                        std::cout << "    " << links[i]->wire()->name();
                        std::cout << " appended (small dist)" << std::endl;
#endif
                        cores.append(links[i]);
                        continue;
                    }
                    if (dist < bestD) {
                        best = links[i];
                        bestD = dist;
                    }
                }
                if (best) {
                    cores.append(best);
#ifdef TRKRECO_DEBUG_DETAIL
                    std::cout << "    " << best->wire()->name();
                    std::cout << " appended (best)" << std::endl;
#endif
                }
            }
        }
    }

    //...Check cores again...
    unsigned nCores = cores.length();
    AList<TMLink> realCores;
    for (unsigned i = 0; i < nCores; i++) {
        TMLink * l = cores[i];
        if (l->hit()->state() & WireHitFittingValid)
            realCores.append(l);
    }
    if (NSuperLayers(realCores) < 2) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ... rejected by small number of super layers" << std::endl;
#endif  
        return NULL;
    }
    if (NLayers(realCores) < 5) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ... rejected by small number of layers" << std::endl;
#endif  
        return NULL;
    }

    //...Make a circle...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << " ... making a circle : #cores=" << cores.length();
    std::cout << std::endl;
#endif
    AList<TMLink> hits = list;
    hits.remove(cores);
    TCircle c(cores);

    //...Test it...
    if (! _circleSelector.preSelect(c)) return NULL;

    //...Fitting...
    int err = c.fit();
    if (err < 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ... rejected by failure of the 1st fit : ";
        std::cout << "err = " << err << std::endl;
#endif
        return NULL;
    }

    //...Test it...
    if (! _circleSelector.select(c)) return NULL;

    //...Make a track...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << " ... making a track" << std::endl;
#endif
    TTrack * t = new TTrack(c);
    if (! _trackSelector.preSelect(* t)) {
        delete t;
        return NULL;
    }

    //...Fitting...
    AList<TMLink> bad;
    //cout<<"track0: "<<t->nLinks()<<endl;
    err = fit(* t);
    if (err < 0) goto discard;
    t->refine(bad, _trackSelector.maxSigma() * 100.);
    //cout<<"track1: "<<t->nLinks()<<endl;
    err = fit(* t);
    if (err < 0) goto discard;
#ifdef TRKRECO_DEBUG_DETAIL
    t->dump("detail", "    1st> ");
#endif

    //...Test it...
    if (! _trackSelector.select(* t)) goto discard;

    //...Try to append non-core hits...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << " ... appending non-core hits" << std::endl;
#endif
    t->appendByApproach(hits, sqrt(_trackSelector.maxSigma()));
    err = fit(* t);
    if (err < 0) goto discard;
//    cout<<"  2D track1: "<<t->links().length()<<endl;
    t->refine(bad, _trackSelector.maxSigma() * 10.);
    err = fit(* t);
    if (err < 0) goto discard;
//    cout<<"  2D track2: "<<t->links().length()<<endl;
    t->refine(bad, _trackSelector.maxSigma());
    err = fit(* t);
    if (err < 0) goto discard;
//    cout<<"  2D track3: "<<t->links().length()<<endl;
#ifdef TRKRECO_DEBUG_DETAIL
    t->dump("detail", "    2nd> ");
#endif

    //...Test it...
    if (! _trackSelector.select(* t)) goto discard;

    //...OK...
#ifdef TRKRECO_DEBUG_DETAIL
//      rphiWindow->append(* t, leda_blue);
//      rphiWindow->draw();
//      rphiWindow->wait();
//      rphiWindow->remove(tmp);
//      rphiWindow->remove(* t);
#endif

    return t;

    //...Something happened...
discard:
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    ... rejected by fitting failure : ";
    std::cout << " err = " << err << std::endl;
//      rphiWindow->append(* t, leda_blue);
//      rphiWindow->draw();
//      rphiWindow->wait();
//      rphiWindow->remove(tmp);
//      rphiWindow->remove(* t);
#endif
    delete t;
    return NULL;
}

TTrack * TBuilderCosmic::buildStereo	(	TTrack &	track,
		const AList< TMLink > &
	)			const [virtual]

appends stereo hits to a track.

Reimplemented from TBuilder0.

Definition at line 32 of file TBuilderCosmic.cxx.

References _fitter, TTrack::_helix, TBuilder0::_lineSelector, TBuilder0::_trackSelector, TLine0::a(), Helix::a(), TTrackBase::append(), TLine0::b(), TTrack::charge(), showlog::err, TCosmicFitter::fit(), TTrackBase::fit(), T3DLine::helix(), TTrack::helix(), genRecEmupikp::i, TMLink::leftRight(), genRecEmupikp::line, TMLink::link(), TTrackBase::links(), TMSelector::maxSigma(), TBuilder0::name(), TMSelector::nLinks(), TMSelector::nLinksStereo(), TMLink::position(), TTrackBase::refine(), TTrackBase::removeLinks(), TMSelector::select(), SortByWireId(), TTrack::szPosition(), t(), tr, WireHitFittingValid, WireHitLeft, and WireHitRight.

                                                                            {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << "(stereo) ... dump of stereo candidate hits" << std::endl;
    AList<TMLink> tmp = list;
    tmp.sort(SortByWireId);
    std::cout << "    ";
    for (unsigned i = 0; i < tmp.length(); i++) {
        TMLink * l = tmp[i];
        std::cout << l->wire()->layerId() << "-";
        std::cout << l->wire()->localId() << ",";
    }
    std::cout << std::endl;
#endif

    //...Check # of links...
    if (list.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << name() << "(stereo) ... rejected by nLinks(";
        std::cout << list.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        return NULL;
    }

    //...Calculate s and z for every links...
    unsigned n = list.length();
    AList<TMLink> forLine;
    for (unsigned i = 0; i < n; i++) {
        TMLink * l = list[i];

        //... Require Fitting vaildation 
        if(!(l->hit()->state()& WireHitFittingValid)) continue;

        TMLink * t = new TMLink(* l);

        //...Assuming wire position...
        t->leftRight(2);
        int err = track.szPosition(* t);
        if (err) {
            delete t;
            continue;
        }

        //...Store the sz link...
        t->link(l);
        forLine.append(t);
    }

#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << "(stereo) ... dump of sz links" << std::endl;
    std::cout << "    ";
    tmp = forLine;
    tmp.sort(SortByWireId);
    for (unsigned i = 0; i < tmp.length(); i++) {
        TMLink * l = tmp[i];
        std::cout << l->wire()->layerId() << "-";
        std::cout << l->wire()->localId() << ",";
    }
    std::cout << std::endl;
#endif  
    
    //...Check # of sz links...
    if (forLine.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << name() << "(stereo) ... rejected by sz nLinks(";
        std::cout << forLine.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        return NULL;
    }

    //...Make a line...
    unsigned nLine = forLine.length();
    TLine0 line(forLine);
    int err = line.fit();

    //...Linear fit...
    if (err < 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << name() << "(stereo) ... linear fit failure. nLinks(";
        std::cout << forLine.length() << ")" << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        return NULL;
    }

#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << "(stereo) ... dump of left-right" << std::endl;
#endif

    //...Decide Left or Right...
    AList<TMLink> forNewLine;
    for (unsigned i = 0; i < nLine; i++) {
        TMLink * t = forLine[i];
        TMLink * tl = new TMLink(* t);
        TMLink * tr = new TMLink(* t);

        tl->leftRight(WireHitLeft);
        tr->leftRight(WireHitRight);

        int err = track.szPosition(* tl);
        if (err) {
            delete tl;
            tl = NULL;
        }
        err = track.szPosition(* tr);
        if (err) {
            delete tr;
            tr = NULL;
        }
        if ((tl == NULL) && (tr == NULL)) continue;

        TMLink * best;
        if (tl == NULL) best = tr;
        else if (tr == NULL) best = tl;
        else {
            if (line.distance(* tl) < line.distance(* tr)) {
                best = tl;
                delete tr;
            }
            else {
                best = tr;
                delete tl;
            }
        }

#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ";
        std::cout << t->wire()->layerId() << "-";
        std::cout << t->wire()->localId();
        if (tl != NULL)
            std::cout << ",left " << tl->position() << "," << line.distance(* tl);
        if (tr != NULL)
            std::cout << ",right " << tr->position() << "," << line.distance(* tr);
        std::cout << std::endl;
#endif

        best->link(t->link());
        forNewLine.append(best);
    }

    //...Check # of sz links...
    if (forNewLine.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << name() << "(stereo) ... rejected by lr nLinks(";
        std::cout << forNewLine.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Create new line...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << "(stereo) ... creating a new line" << std::endl;
#endif  
    unsigned nNewLine = forNewLine.length();
    TLine0 newLine(forNewLine);

    //...Make a seed track again
    err = newLine.fit();

    //...Linear fit...
    if (err < 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << name() << "(stereo) ... 2nd linear fit failure. nLinks(";
        std::cout << forNewLine.length() << ")" << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Remove bad points...
    AList<TMLink> bad;
//    newLine.refine(bad, 40.);  //Liuqg, meaningless while without magnetic field
//    err = newLine.fit();
//    newLine.refine(bad, 20.);
//    err = newLine.fit();
//    newLine.refine(bad, 10.);
//    err = newLine.fit();

    //...Linear fit again...
    if (err < 0) {
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    appendStereo cut ... new line 2nd linear fit failure. ";
        std::cout << "# of links = " << n << "," << nLine;
        std::cout << "," << nNewLine << std::endl;
#endif  
        return NULL;
    }

    //...3D fit...
    const AList<TMLink> & good = newLine.links();
    unsigned nn = good.length();
    for (unsigned i = 0; i < nn; i++) {
        track.append(* good[i]->link());
    }
    Vector a(5);
    a = track.helix().a();
    a[3] = newLine.b();
    a[4] = track.charge() * newLine.a();
    track._helix->a(a);

#ifdef LINE_COSMIC
    T3DLine * Ltrack = new T3DLine(track);
    
    //...Refine...
    err = _fitter.fit(*Ltrack);
    Ltrack->refine(bad, _trackSelector.maxSigma() * 30.);
    err = _fitter.fit(*Ltrack);
    Ltrack->refine(bad, _trackSelector.maxSigma() * 3.);
    err = _fitter.fit(*Ltrack);
//    Ltrack->refine(bad, _trackSelector.maxSigma() * 0.21);  //liuqg, for prelimilary test
    track.refine(bad, _trackSelector.maxSigma()); 
    err = _fitter.fit(*Ltrack);

    //...Test it...
/*    if (! _trackSelector.select(*Ltrack)) {
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        delete Ltrack;
        return NULL;
    }
*/
    //...Termination...
    HepAListDeleteAll(forLine);
    HepAListDeleteAll(forNewLine);
    track.removeLinks();
    track.append(Ltrack->links());
    Vector a1(5);
    a1 = Ltrack->helix().a();
    track._helix->a(a1);
    delete Ltrack;
    return & track;
#endif

    //...Refine...
    err = _fitter.fit(track);
    track.refine(bad, _trackSelector.maxSigma() * 30.);
    err = _fitter.fit(track);
    track.refine(bad, _trackSelector.maxSigma() * 3.);
    err = _fitter.fit(track);
    track.refine(bad, _trackSelector.maxSigma());  //liuqg, for prelimilary test
//    track.refine(bad, _trackSelector.maxSigma() * 0.21); 
    err = _fitter.fit(track);
    //...Test it...
    if (! _trackSelector.select(track)) {
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Termination...
    HepAListDeleteAll(forLine);
    HepAListDeleteAll(forNewLine);
    return & track;
}

TTrack * TBuilder0::buildStereo0	(	TTrack &	track,
		const AList< TMLink > &
	)			const [inherited]

appends stereo hits to a track. (old version)

Definition at line 307 of file TBuilder0.cxx.

References TTrack::_helix, TBuilder0::_lineSelector, TBuilder0::_name, TBuilder0::_trackSelector, TLine0::a(), Helix::a(), TTrackBase::append(), TTrack::charge(), Helix::curv(), Dump(), showlog::err, TBuilder0::fit(), TTrackBase::fit(), TLine0::fit2p(), TLine0::fit2sp(), TTrack::helix(), genRecEmupikp::i, TMLink::leftRight(), genRecEmupikp::line, TMLink::link(), TTrackBase::links(), TMSelector::maxSigma(), TMSelector::nLinks(), TMSelector::nLinksStereo(), TMLink::position(), rb::R(), TTrackBase::refine(), TLine0::refine(), TLine0::removeChits(), TMSelector::select(), TTrack::szPosition(), t(), tr, WireHitLeft, and WireHitRight.

                                                                        {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << "(stereo) ... dump of stereo candidate hits" << std::endl;
    Dump(list, "sort flag", "    ");
#endif

    //...Check # of links...
    if (list.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "(stereo) ... rejected by nLinks(";
        std::cout << list.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        return NULL;
    }

    //...Calculate s and z for every links...
    unsigned n = list.length();
    AList<TMLink> forLine;
    for (unsigned i = 0; i < n; i++) {
        TMLink * l = list[i];
        TMLink * t = new TMLink(* l);

        //...Assuming wire position...
        t->leftRight(2);
        int err = track.szPosition(* t);
        if (err) {
            delete t;
            continue;
        }

        //...Store the sz link...
        t->link(l);
        forLine.append(t);
    }

#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << "(stereo) ... dump of sz links" << std::endl;
    Dump(forLine, "sort flag", "    ");
#endif    

    //...Check # of sz links...
    if (forLine.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "(stereo) ... rejected by sz nLinks(";
        std::cout << forLine.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        return NULL;
    }

    //...Make a line...
    unsigned nLine = forLine.length();
    TLine0 line(forLine);
    int err = line.fit2sp();
    if (err < 0) {
        err = line.fit2p();
        if (err < 0) err = line.fit();
    }

    //...Linear fit...
    if (err < 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "(stereo) ... linear fit failure. nLinks(";
        std::cout << forLine.length() << ")" << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        return NULL;
    }

#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << "(stereo) ... dump of left-right" << std::endl;
#endif

    //...Decide Left or Right...
    AList<TMLink> forNewLine;
    for (unsigned i = 0; i < nLine; i++) {
        TMLink * t = forLine[i];
        TMLink * tl = new TMLink(* t);
        TMLink * tr = new TMLink(* t);

        tl->leftRight(WireHitLeft);
        tr->leftRight(WireHitRight);

        int err = track.szPosition(* tl);
        if (err) {
            delete tl;
            tl = NULL;
        }
        err = track.szPosition(* tr);
        if (err) {
            delete tr;
            tr = NULL;
        }
        if ((tl == NULL) && (tr == NULL)) continue;

        TMLink * best;
        if (tl == NULL) best = tr;
        else if (tr == NULL) best = tl;
        else {
            if (line.distance(* tl) < line.distance(* tr)) {
                best = tl;
                delete tr;
            }
            else {
                best = tr;
                delete tl;
            }
        }

#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    ";
        std::cout << t->wire()->layerId() << "-";
        std::cout << t->wire()->localId();
        if (tl != NULL)
            std::cout << ",left " << tl->position() << "," << line.distance(* tl);
        if (tr != NULL)
            std::cout << ",right " << tr->position() << "," << line.distance(* tr);
        std::cout << std::endl;
#endif

        best->link(t->link());
        forNewLine.append(best);
    }

    //...Check # of sz links...
    if (forNewLine.length() < _lineSelector.nLinksStereo()) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "(stereo) ... rejected by lr nLinks(";
        std::cout << forNewLine.length() << ") < ";
        std::cout << _lineSelector.nLinks() << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Create new line...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << _name << "(stereo) ... creating a new line" << std::endl;
#endif  
    unsigned nNewLine = forNewLine.length();
    TLine0 newLine(forNewLine);

    //... Remove extremely bad points
    newLine.removeChits();

    //...Make a seed track again
    err = newLine.fit2sp();
    if (err < 0) {
        err = newLine.fit2p();
        if (err < 0) err = newLine.fit();
    }

    //...Linear fit...
    if (err < 0) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << _name << "(stereo) ... 2nd linear fit failure. nLinks(";
        std::cout << forNewLine.length() << ")" << std::endl;
#endif  
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Remove bad points...
    AList<TMLink> bad;
    newLine.refine(bad, 40.);
    err = newLine.fit();
    newLine.refine(bad, 30.);
    err = newLine.fit();
    newLine.refine(bad, 20.);
    err = newLine.fit();
    newLine.refine(bad, 10.);
    err = newLine.fit();
    float R = fabs(track.helix().curv());
    if (R > 80.) {
        newLine.refine(bad, 5.);
        err = newLine.fit();
    }

    //...Linear fit again...
    if (err < 0) {
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    appendStereo cut ... new line 2nd linear fit failure. ";
        std::cout << "# of links = " << n << "," << nLine;
        std::cout << "," << nNewLine << std::endl;
#endif  
        return NULL;
    }

    //...3D fit...
    const AList<TMLink> & good = newLine.links();
    unsigned nn = good.length();
    for (unsigned i = 0; i < nn; i++) {
        track.append(* good[i]->link());
    }
    Vector a(5);
    a = track.helix().a();
    a[4] = track.charge() * newLine.a();
    track._helix->a(a);

    //...Refine...
    err = fit(track);
    track.refine(bad, _trackSelector.maxSigma() * 100.);
    err = fit(track);
    track.refine(bad, _trackSelector.maxSigma() * 10.);
    err = fit(track);
    track.refine(bad, _trackSelector.maxSigma());
    err = fit(track);

    //...Test it...
    if (! _trackSelector.select(track)) {
        HepAListDeleteAll(forLine);
        HepAListDeleteAll(forNewLine);
        return NULL;
    }

    //...Termination...
    HepAListDeleteAll(forLine);
    HepAListDeleteAll(forNewLine);
    return & track;
}

void TBuilder0::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)			const [inherited]

dumps debug information.

Definition at line 79 of file TBuilder0.cxx.

                                                                  {
}

int TBuilder0::fit

(

TTrackBase &

)

const [inline, virtual, inherited]

fits a track using a private fitter.

Definition at line 143 of file TBuilder0.h.

References TBuilder0::_fitter, and TCosmicFitter::fit().

Referenced by TBuilder0::appendClusters(), TBuilder0::buildRphi(), TBuilder0::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::resetHelixFit(), TBuilder0::salvage(), and TBuilder0::salvageNormal().

                                   {
    return _fitter.fit(a);
}

const std::string & TBuilder0::name

(

void

)

const [inline, inherited]

returns name.

Definition at line 137 of file TBuilder0.h.

References TBuilder0::_name.

Referenced by buildStereo(), TBuilder0::salvage(), and TBuilder0::salvageNormal().

                          {
    return _name;
}

void TBuilder0::salvage	(	TTrack &	track,
		AList< TMLink > &	list
	)			const [inherited]

salvages links in a list. Used links will be removed from a list.

Definition at line 1479 of file TBuilder0.cxx.

References TBuilder0::_salvageLevel, TBuilder0::fit(), genRecEmupikp::i, TBuilder0::name(), TBuilder0::salvageNormal(), TMDCWireHit::state(), t(), TrackOldConformalFinder, WireHitConformalFinder, and WireHitUsed.

                                                         {

//    if (t.type() == TrackTypeNormal) {
        salvageNormal(t, hits);   //Pt high enough; fromIP.
        return;
//    }
  //tmply for cosmic ray...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << " ... salvaging" << std::endl;
    std::cout << "    # of given hits=" << hits.length() << std::endl;
#endif

    unsigned nHits = hits.length();
    if (nHits == 0) return;

    //...Hit loop...
    AList<TMLink> candidates;
    for (unsigned i = 0; i < nHits; i++) {
        TMLink & l = * hits[i];
        const TMDCWireHit & h = * l.hit();

        //...Already used?...
        if (h.state() & WireHitUsed) continue;
        candidates.append(l);
    }

    //...Try to append this hit...
//    t.appendByApproach(candidates, 10.);
    t.appendByApproach(candidates, _salvageLevel);
    fit(t);
    hits.remove(candidates);
    t.assign(WireHitConformalFinder);
    t.finder(TrackOldConformalFinder);
    // t.assign(WireHitConformalFinder, TrackOldConformalFinder);
}

const TMSelector & TBuilder0::trackSelector

(

const TMSelector &

)

[virtual, inherited]

sets a track selector.

Definition at line 1400 of file TBuilder0.cxx.

References TBuilder0::_circleSelector, TBuilder0::_lineSelector, TBuilder0::_trackSelector, TMSelector::maxDistance(), TMSelector::maxImpact(), TMSelector::maxSigma(), TMSelector::minPt(), TMSelector::nLinks(), TMSelector::nLinksStereo(), and TMSelector::nSuperLayers().

                                             {
    _circleSelector.nLinks(a.nLinks());
    _circleSelector.nSuperLayers(a.nSuperLayers());
    _circleSelector.minPt(a.minPt());
    _circleSelector.maxImpact(a.maxImpact());

    _trackSelector.nLinks(a.nLinks());
    _trackSelector.nSuperLayers(a.nSuperLayers());
    _trackSelector.maxSigma(a.maxSigma());

    _lineSelector.nLinksStereo(a.nLinksStereo());
    _lineSelector.maxDistance(a.maxDistance());

    return a;
}

const TMSelector & TBuilder0::trackSelector

(

void

)

const [inline, inherited]

returns a track selector.

Definition at line 131 of file TBuilder0.h.

References TBuilder0::_trackSelector.

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

                                   {
    return _trackSelector;
}

---

Member Data Documentation

TMSelector TBuilder0::_circleSelector [protected, inherited]

Definition at line 106 of file TBuilder0.h.

Referenced by TBuilder0::buildRphi(), and TBuilder0::trackSelector().

TCosmicFitter TBuilderCosmic::_fitter [private]

Reimplemented from TBuilder0.

Definition at line 54 of file TBuilderCosmic.h.

Referenced by buildStereo().

TMSelector TBuilder0::_lineSelector [protected, inherited]

Definition at line 110 of file TBuilder0.h.

Referenced by buildStereo(), TBuilder0::buildStereo(), TBuilder0::buildStereo0(), and TBuilder0::trackSelector().

float TBuilder0::_stereoChisq3 [protected, inherited]

Definition at line 113 of file TBuilder0.h.

Referenced by TBuilder0::buildStereo().

float TBuilder0::_stereoChisq4 [protected, inherited]

Definition at line 114 of file TBuilder0.h.

Referenced by TBuilder0::buildStereo().

float TBuilder0::_stereoMaxSigma [protected, inherited]

Definition at line 115 of file TBuilder0.h.

Referenced by TBuilder0::buildStereo().

float TBuilder0::_stereoZ3 [protected, inherited]

Definition at line 111 of file TBuilder0.h.

Referenced by TBuilder0::buildStereo().

float TBuilder0::_stereoZ4 [protected, inherited]

Definition at line 112 of file TBuilder0.h.

Referenced by TBuilder0::buildStereo().

TMSelector TBuilder0::_trackSelector [protected, inherited]

Definition at line 107 of file TBuilder0.h.

Referenced by TBuilder0::appendClusters(), TBuilder0::buildRphi(), buildStereo(), TBuilder0::buildStereo(), TBuilder0::buildStereo0(), and TBuilder0::trackSelector().

---