A class to find tracks with the conformal method. More...

#include <TConformalFinder0.h>

Inheritance diagram for TConformalFinder0:

Public Member Functions
	TConformalFinder0 (float maxSigma, float fraction, float stereoZ3, float stereoZ4, float stereoChisq3, float stereoChisq4, float stereoMaxSigma, unsigned fittingCorrections, float salvageLevel, bool cosmic)
	Constructor.

virtual	~TConformalFinder0 ()
	Destructor.

std::string	name (void) const
	returns name.

std::string	version (void) const
	returns version.

void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.

const AList< TMLink > &	axialConformalLinks (void) const
	returns a list of conformal links.

const AList< TMLink > &	stereoConformalLinks (void) const
	returns a list of conformal links.

void	clear (void)
	clear internal information.

int	doit (const AList< TMDCWireHit > &axialHits, const AList< TMDCWireHit > &stereoHits, AList< TTrack > &tracks, AList< TTrack > &tracks3D)
	finds tracks.

AList< AList< TSegment0 > >	findSegments (const AList< TMLink > &in) const
	finds segments.

AList< AList< TSegment0 > >	findSegments2 (const AList< TMLink > &in) const

AList< TSegment0 >	findClusters (const THistogram &) const
	finds segments. (obsolete functions)

AList< TSegment0 >	findClusters2 (const THistogram &) const

bool	doStereo (bool)
	sets swtich for stereo reconstruction.

bool	doSalvage (bool)
	sets switch for salvaging.

Public Member Functions inherited from TFinderBase
	TFinderBase ()
	Constructor.

virtual	~TFinderBase ()
	Destructor.

virtual int	debugLevel (void) const
	returns debug level.

virtual int	debugLevel (int)
	sets debug level.

Static Public Member Functions
static void	conformalTransformation (const HepPoint3D &center, const AList< TMDCWireHit > &hits, AList< TMLink > &links)
	transforms hits into a conformal plane. 'center' is a center of the transformation. Z position of 'center' must be zero. Transformed positions are newly created, and are stored in 'links'.

static void	conformalTransformationDriftCircle (const HepPoint3D &center, const AList< TMDCWireHit > &hits, AList< TMLink > &links)
	transforms drift circle of hits into a conformal plane. transformed positions( x0, y0, r) are stored in 'links'.

static void	conformalTransformationRphi (const HepPoint3D &center, const AList< TMDCWireHit > &hits, AList< TMLink > &links)
	transforms hits into a conformal plane. 'center' is a center of the transformation. Z position of 'center' must be zero. Transformed positions are newly created, and are stored in 'links'.

Detailed Description

A class to find tracks with the conformal method.

Definition at line 40 of file TConformalFinder0.h.

Constructor & Destructor Documentation

◆ TConformalFinder0()

TConformalFinder0::TConformalFinder0	(	float	maxSigma,
		float	fraction,
		float	stereoZ3,
		float	stereoZ4,
		float	stereoChisq3,
		float	stereoChisq4,
		float	stereoMaxSigma,
		unsigned	fittingCorrections,
		float	salvageLevel,
		bool	cosmic )

Constructor.

Definition at line 33 of file TConformalFinder0.cxx.

: TFinderBase(),
  _builder(0),
  _doStereo(true),
  _doSalvage(true),  //liuqg
  _fraction(fraction) {
 
    //...Parameters for a track...
    _trackSelector.nLinks(4);
    _trackSelector.nSuperLayers(2);
    _trackSelector.minPt(0.05);
    _trackSelector.maxImpact(100.);
    _trackSelector.maxSigma(maxSigma);
    _trackSelector.nLinksStereo(3);
    _trackSelector.maxDistance(30.);
 
    //...Make a builder...
    if (cosmic) _builder = new TBuilderCosmic("cosmic builder", salvageLevel);
    else        _builder = new TBuilder0("conformal builder",
                     stereoZ3,
                     stereoZ4,
                     stereoChisq3,
                     stereoChisq4,
                     stereoMaxSigma,
                     fittingCorrections,
                     salvageLevel);
 
    //...Set up TBuilder...
    _builder->trackSelector(_trackSelector);
}

◆ ~TConformalFinder0()

TConformalFinder0::~TConformalFinder0 ( )

virtual

Destructor.

Definition at line 73 of file TConformalFinder0.cxx.

                                      {
    delete _builder;
}

Member Function Documentation

◆ axialConformalLinks()

const AList< TMLink > & TConformalFinder0::axialConformalLinks ( void ) const

inline

returns a list of conformal links.

Definition at line 190 of file TConformalFinder0.h.

                                                 {
    return _axialConfLinks;
}

◆ clear()

void TConformalFinder0::clear ( void )

virtual

clear internal information.

Implements TFinderBase.

Definition at line 89 of file TConformalFinder0.cxx.

                             {
    HepAListDeleteAll(_axialConfLinks);
    HepAListDeleteAll(_stereoConfLinks);
    _unusedAxialConfLinks.removeAll();
    _unusedStereoConfLinks.removeAll();
    _goodAxialConfLinks.removeAll();
    HepAListDeleteAll(_circles);
    _tracks.removeAll();
}

◆ conformalTransformation()

void TConformalFinder0::conformalTransformation	(	const HepPoint3D &	center,
		const AList< TMDCWireHit > &	hits,
		AList< TMLink > &	links )

static

transforms hits into a conformal plane. 'center' is a center of the transformation. Z position of 'center' must be zero. Transformed positions are newly created, and are stored in 'links'.

Definition at line 100 of file TConformalFinder0.cxx.

                                             {
 
    unsigned nHits = hits.length();
    if (center == ORIGIN) {
    for (unsigned i = 0; i < nHits; i++) {
        TMDCWireHit * h = hits[i];
        const HepPoint3D & p = h->xyPosition();
        
//zsl       HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2());
        HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
//      std::cout<<" xypo "<<p<<" pos "<<cp<<std::endl;
        links.append(new TMLink(0, h, cp));
    }
    }
    else {
    for (unsigned i = 0; i < nHits; i++) {
        TMDCWireHit * h = hits[i];
        HepPoint3D p(h->xyPosition() - center);
        HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
        links.append(new TMLink(0, h, cp));
    }   
    }
}

◆ conformalTransformationDriftCircle()

void TConformalFinder0::conformalTransformationDriftCircle	(	const HepPoint3D &	center,
		const AList< TMDCWireHit > &	hits,
		AList< TMLink > &	links )

static

transforms drift circle of hits into a conformal plane. transformed positions( x0, y0, r) are stored in 'links'.

Definition at line 127 of file TConformalFinder0.cxx.

                                                                             { //added by Liuqg for Tsf
    unsigned nHits = hits.length();
    if (center == ORIGIN) {
        for (unsigned i = 0; i < nHits; i++) {
            TMDCWireHit * h = hits[i];
            const HepPoint3D & p = h->xyPosition();
                                                                                                                             
            const double r = 0.5*(h->drift(0) + h->drift(1));
            HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
//      HepPoint3D cp2(100 * 2. * p.x() / (p.mag2() - r*r), 100 * 2. * p.y() / (p.mag2() - r*r), 0.);
        HepPoint3D cp2(2. * p.x() / (p.mag2() - r*r), 2. * p.y() / (p.mag2() - r*r), 0.);
//          float cDrift = 100 * 2. * r / (p.mag2() - r*r);
        double cDrift = 2. * r / (p.mag2() - r*r);
                                                                                                                             
        links.append(new TMLink(0, h, cp, cp2, cDrift));
        }
    }
    else {
        for (unsigned i = 0; i < nHits; i++) {
            TMDCWireHit * h = hits[i];
            HepPoint3D p(h->xyPosition() - center);
                                                                                                                             
            const double r = 0.5*(h->drift(0) + h->drift(1));
            HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
        //unit of the following is km-1, origin is cm.
//      HepPoint3D cp2(100 * 2. * p.x() / (p.mag2() - r*r), 100 * 2. * p.y() / (p.mag2() - r*r), 0.);
        HepPoint3D cp2(2. * p.x() / (p.mag2() - r*r), 2. * p.y() / (p.mag2() - r*r), 0.);
//          float cDrift = 100 * 2. * r / (p.mag2() - r*r);
        double cDrift = 2. * r / (p.mag2() - r*r);
                                                                                                                             
        links.append(new TMLink(0, h, cp, cp2, cDrift));
        }
    }
}

Referenced by TConformalFinder::doit(), and TrkReco::execute().

◆ conformalTransformationRphi()

void TConformalFinder0::conformalTransformationRphi	(	const HepPoint3D &	center,
		const AList< TMDCWireHit > &	hits,
		AList< TMLink > &	links )

static

transforms hits into a conformal plane. 'center' is a center of the transformation. Z position of 'center' must be zero. Transformed positions are newly created, and are stored in 'links'.

Definition at line 165 of file TConformalFinder0.cxx.

                                                 {
 
    unsigned nHits = hits.length();
    if (center == ORIGIN) {
    for (unsigned i = 0; i < nHits; i++) {
        TMDCWireHit * h = hits[i];
        const HepPoint3D & p = h->xyPosition();
        HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
        double r = log(cp.mag()) + 4.;
        double phi = atan2(cp.y(), cp.x()) + M_PI;
        HepPoint3D cpt(phi, r, 0.);
        links.append(new TMLink(0, h, cpt));
    }
    }
    else {
    for (unsigned i = 0; i < nHits; i++) {
        TMDCWireHit * h = hits[i];
        HepPoint3D p(h->xyPosition() - center);
        HepPoint3D cp(2. * p.x() / p.mag2(), 2. * p.y() / p.mag2(), 0.);
        double r = log(cp.mag()) + 4.;
        double phi = atan2(cp.y(), cp.x()) + M_PI;
        HepPoint3D cpt(phi, r, 0.);
        links.append(new TMLink(0, h, cpt));
    }   
    }
}

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

◆ doit()

int TConformalFinder0::doit	(	const AList< TMDCWireHit > &	axialHits,
		const AList< TMDCWireHit > &	stereoHits,
		AList< TTrack > &	tracks,
		AList< TTrack > &	tracks3D )

virtual

finds tracks.

Implements TFinderBase.

Definition at line 560 of file TConformalFinder0.cxx.

                             {
 
    //...For debug...
    if (debugLevel()) {
    std::cout << name() << " ... processing" << std::endl;
    std::cout << "    axialHits=" << axialHits.length();
    std::cout << ",stereoHits=" << stereoHits.length();
    std::cout << ",tracks=" << tracks.length();
    std::cout << std::endl;
 
    if (debugLevel() > 1)
        std::cout << name() << " ... conformal transformation0" << std::endl;
    }
 
    //...Conformal transformation with IP constraint...
    conformalTransformationRphi(ORIGIN, axialHits, _axialConfLinks);
    conformalTransformationRphi(ORIGIN, stereoHits, _stereoConfLinks);
    _unusedAxialConfLinks.append(_axialConfLinks);
    _unusedStereoConfLinks.append(_stereoConfLinks);
    AList<TMLink> unusedConfLinks;
    if (_doSalvage) {
    unusedConfLinks.append(_axialConfLinks);
    unusedConfLinks.append(_stereoConfLinks);
    }
 
    //...For debug...
    if (debugLevel() > 1)
    std::cout << name() << " ... selecting good hits" << std::endl;
 
    //...Select good axial hits...
    AList<TMLink> goodHits;
    int nLinks = _axialConfLinks.length();
    for (unsigned i = 0; i < nLinks; i++) {
    TMLink * l = _axialConfLinks[i];
    const TMDCWireHit & h = * l->hit();
//liuqg if ((h.state() & WireHitIsolated) &&
//      (h.state() & WireHitContinuous))
        goodHits.append(l);
    }
    //...Main algorithm...
    standardFinding(goodHits, unusedConfLinks, _fraction);
 
    //...Main algorithm for second trial...
    specialFinding(goodHits, unusedConfLinks, _fraction);
 
    //...For debug...
    if (debugLevel()) {
    std::cout << name() << " ... processed : ";
    std::cout << "good hits=" << goodHits.length();
    std::cout << ",tracks=" << _tracks.length();
    std::cout << std::endl;
    }
 
    tracks.append(_tracks);
    return 0;
}

◆ doSalvage()

bool TConformalFinder0::doSalvage ( bool a )

inlinevirtual

sets switch for salvaging.

Reimplemented from TFinderBase.

Definition at line 208 of file TConformalFinder0.h.

                                   {
    return _doSalvage = a;
}

◆ doStereo()

bool TConformalFinder0::doStereo ( bool a )

inlinevirtual

sets swtich for stereo reconstruction.

Reimplemented from TFinderBase.

Definition at line 202 of file TConformalFinder0.h.

                                  {
    return _doStereo = a;
}

◆ dump()

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

virtual

dumps debug information.

Reimplemented from TFinderBase.

Definition at line 78 of file TConformalFinder0.cxx.

                                                                          {
    std::cout << pre;
    TFinderBase::dump(msg);
    std::cout << pre;
    if (msg.find("state") != std::string::npos) {
    std::cout << "#axialConfPos=" << _axialConfLinks.length();
    std::cout << ",#stereoConfPos=" << _stereoConfLinks.length();
    }
}

◆ findClusters()

AList< TSegment0 > TConformalFinder0::findClusters ( const THistogram & hist ) const

finds segments. (obsolete functions)

Definition at line 195 of file TConformalFinder0.cxx.

                                                             {
 
    //...Obtain raw clusters...
    AList<TSegment0> list = hist.clusters0();
    unsigned n = list.length();
    if (n == 0) return list;
 
#ifdef TRKRECO_DEBUG_DETAIL
    // static TChecker chk0("clusters before splitting");
    // chk0.check(list);
    // chk0.dump("detail", "    ");
#endif
 
    //...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);
 
#ifdef TRKRECO_DEBUG_DETAIL
    // static TChecker chk1("clusters after splitting");
    // chk1.check(list);
    // chk1.dump("detail", "    ");
#endif
 
    return list;
}

Referenced by findSegments().

◆ findClusters2()

AList< TSegment0 > TConformalFinder0::findClusters2 ( const THistogram & hist ) const

Definition at line 243 of file TConformalFinder0.cxx.

                                                              {
 
    //...Obtain raw clusters...
    AList<TSegment0> list = hist.clusters0();
    unsigned n = list.length();
    if (n == 0) return list;
 
#ifdef TRKRECO_DEBUG_DETAIL
    // static TChecker chk0("clusters before splitting (2)");
    // chk0.check(list);
    // chk0.dump("detail", "    ");
#endif
 
    //...Examine each cluster...
    for (unsigned i = 0; i < n; i++) {
    TSegment0 * c = list[i];
    unsigned type = c->clusterType();
 
    if ((type == 1) || (type == 2)) {
        c->dump("hits mc", "    ");
        c->solveDualHits();
    }
    }
 
#ifdef TRKRECO_DEBUG_DETAIL
    // static TChecker chk1("clusters after splitting (2)");
    // chk1.check(list);
    // chk1.dump("detail", "    ");
#endif
 
    return list;
}

◆ findSegments()

AList< AList< TSegment0 > > TConformalFinder0::findSegments ( const AList< TMLink > & in ) const

finds segments.

Definition at line 869 of file TConformalFinder0.cxx.

                                                              {
    AList< AList<TSegment0> > a;
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << name() << " ... finding segments : given hits =" << std::endl;
    Dump(in, "sort");
#endif
 
    //...Create lists of links for each super layer...
    AList<TMLink> links[5];
    unsigned n = in.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * in[i];
//  links[l.wire()->superLayerId() / 2].append(l);
    links[l.wire()->axialStereoLayerId()/4].append(l);
    }
 
    //...Create phi hists and clusters for each super layer...
    THistogram * hist[5];
    hist[0] = new THistogram(76);
    hist[1] = new THistogram(100);
    hist[2] = new THistogram(128);
    hist[3] = new THistogram(256);
    hist[4] = new THistogram(288);
    for (unsigned i = 0; i < 5; i++) {
    hist[i]->fillX(links[i]);
    AList<TSegment0> * b = new AList<TSegment0>();
    a.append(b);
    b->append(findClusters(* hist[i]));
    delete hist[i];
    }
 
    return a;
}

◆ findSegments2()

AList< AList< TSegment0 > > TConformalFinder0::findSegments2 ( const AList< TMLink > & in ) const

◆ name()

std::string TConformalFinder0::name ( void ) const

inlinevirtual

returns name.

Implements TFinderBase.

Definition at line 214 of file TConformalFinder0.h.

                                  {
    return "Conformal Finder 0";
}

Referenced by doit(), and findSegments().

◆ stereoConformalLinks()

const AList< TMLink > & TConformalFinder0::stereoConformalLinks ( void ) const

inline

returns a list of conformal links.

Definition at line 196 of file TConformalFinder0.h.

                                                  {
    return _stereoConfLinks;
}

◆ version()

std::string TConformalFinder0::version ( void ) const

virtual

returns version.

Implements TFinderBase.

Definition at line 29 of file TConformalFinder0.cxx.

                                     {
    return "1.63";
}

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

Public Member Functions

Static Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ TConformalFinder0()

◆ ~TConformalFinder0()

Member Function Documentation

◆ axialConformalLinks()

◆ clear()

◆ conformalTransformation()

◆ conformalTransformationDriftCircle()

◆ conformalTransformationRphi()

◆ doit()

◆ doSalvage()

◆ doStereo()

◆ dump()

◆ findClusters()

◆ findClusters2()

◆ findSegments()

◆ findSegments2()

◆ name()

◆ stereoConformalLinks()

◆ version()