/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/TrkReco/TConformalFinder0.h

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TConformalFinder0.h,v 1.8 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TConformalFinder0.h
// Section  : Tracking
// Owner    : Yoshi Iwasaki
// Email    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// Description : A class to find tracks with the conformal method.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------

#ifndef TConformalFinder0_FLAG_
#define TConformalFinder0_FLAG_

#ifdef TRKRECO_DEBUG_DETAIL
#ifndef TRKRECO_DEBUG
#define TRKRECO_DEBUG
#endif
#endif
#define HEP_SHORT_NAMES

#include "TrkReco/TFinderBase.h"
#include "TrkReco/TBuilder0.h"

#ifndef CLHEP_POINT3D_H
#include "CLHEP/Geometry/Point3D.h"
#endif
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
    typedef HepGeom::Point3D<double> HepPoint3D;
#endif

//class HepPoint3D;
class TMLink;
class THistogram;
class TCircle;
class TSegment0;

class TConformalFinder0 : public TFinderBase {

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

    virtual ~TConformalFinder0();

  public:// Selectors
    std::string name(void) const;

    std::string version(void) const;

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

    const AList<TMLink> & axialConformalLinks(void) const;

    const AList<TMLink> & stereoConformalLinks(void) const;

  public:// Modifiers
    void clear(void);

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

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

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

    static void conformalTransformationRphi(const HepPoint3D & center,
                                            const AList<TMDCWireHit> & hits,
                                            AList<TMLink> & links);
    
    AList< AList<TSegment0> > findSegments(const AList<TMLink> & in) const;
    AList< AList<TSegment0> > findSegments2(const AList<TMLink> & in) const;

    AList<TSegment0> findClusters(const THistogram &) const;
    AList<TSegment0> findClusters2(const THistogram &) const;


    bool doStereo(bool);

    bool doSalvage(bool);

  private:
    AList<TMLink> findCloseHits(const AList<TMLink> & in,
                               const TTrack & track) const;

    TSegment0 * findBestLink(const TSegment0 &,
                            const AList<TSegment0> & in) const;

    TSegment0 * appendCluster(TTrack & track, AList<TSegment0> & in) const;

    AList<TSegment0> findClusterLink(TSegment0 &,
                                    const AList<TSegment0> * const) const;

    TTrack * makeTrack(const AList<TSegment0> &) const;

    void appendClusters2(TTrack & track, AList<TSegment0> &) const;

    AList<TSegment0> findCloseClusters(const TTrack &,
                                      const AList<TSegment0> &,
                                      double maxDistance) const;

    void standardFinding(AList<TMLink> & seeds,
                         AList<TMLink> & unusedLinks,
                         double fraction);

    void specialFinding(AList<TMLink> & seeds,
                        AList<TMLink> & unusedLinks,
                        double fraction);

  private:
    TBuilder0 * _builder;

    TMSelector _circleSelector;
    TMSelector _trackSelector;
    float _fraction;

    AList<TMLink> _axialConfLinks;
    AList<TMLink> _stereoConfLinks;
    AList<TMLink> _unusedAxialConfLinks;
    AList<TMLink> _unusedStereoConfLinks;
    AList<TMLink> _goodAxialConfLinks;

    AList<TCircle> _circles;
    AList<TTrack> _tracks;

    bool _doStereo;
    bool _doSalvage;
};

//-----------------------------------------------------------------------------

#ifdef TConformalFinder0_NO_INLINE
#define inline
#else
#undef inline
#define TConformalFinder0_INLINE_DEFINE_HERE
#endif

#ifdef TConformalFinder0_INLINE_DEFINE_HERE

inline
const AList<TMLink> &
TConformalFinder0::axialConformalLinks(void) const {
    return _axialConfLinks;
}

inline
const AList<TMLink> &
TConformalFinder0::stereoConformalLinks(void) const {
    return _stereoConfLinks;
}

inline
bool
TConformalFinder0::doStereo(bool a) {
    return _doStereo = a;
}

inline
bool
TConformalFinder0::doSalvage(bool a) {
    return _doSalvage = a;
}

inline
std::string
TConformalFinder0::name(void) const {
    return "Conformal Finder 0";
}

#endif

#undef inline

#endif /* TConformalFinder0_FLAG_ */

---