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 std::string	name (void) const =0
	returns name.

virtual std::string	version (void) const =0
	returns version.

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

virtual int	debugLevel (void) const
	returns debug level.

virtual void	clear (void)=0
	clear internal information.

virtual int	debugLevel (int)
	sets debug level.

virtual bool	doStereo (bool)
	sets flag to reconstruct 3D.

virtual bool	doSalvage (bool)
	sets flag to salvage hits.

virtual int	doit (const AList< TMDCWireHit > &axialHits, const AList< TMDCWireHit > &stereoHits, AList< TTrack > &tracks3D, AList< TTrack > &tracks2D)=0
	finds tracks. 'hits' are used to reconstruct. 'tracks' can be used for both inputs and outputs. Return value = (0, +, -) means (success, warning, fatal error).

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 TFastFinder::doit(), and 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()