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

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TLine0.h,v 1.5 2010/03/31 09:58:59 liucy Exp $
//-----------------------------------------------------------------------------
// Filename : TLine0.h
// Section  : Tracking
// Owner    : Yoshi Iwasaki
// Email    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// Description : A class to represent a line in tracking.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------

#ifndef TLine0_FLAG_
#define TLine0_FLAG_

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

#include <string>

#define HEP_SHORT_NAMES
#include "TrkReco/TTrackBase.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TLineFitter.h"

class TLine0 : public TTrackBase {

  public:
    TLine0();

    TLine0(const AList<TMLink> &);

    virtual ~TLine0();

  public:// Selectors
    virtual unsigned objectType(void) const;

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

    double a(void) const;

    double b(void) const;

    double chi2(void) const;

    double reducedChi2(void) const;

  public:// Utilities
    double distance(const TMLink &) const;

  public:// Modifiers
    // int fitx(void);

    int fit2();

    int fit2s();

    int fit2p();

    int fit2sp();

    void removeChits();

    void refine(AList<TMLink> & list, float maxSigma);

    void removeSLY(AList<TMLink> & list);

    void appendSLY(AList<TMLink> & list);

    void appendByszdistance(AList<TMLink> & list, unsigned isl, float maxSigma);

    void property(double a, double b, double det);

  private:// Always updated
    mutable bool _fittedUpdated;

  private:// Updated when fitted
    double _a;
    double _b;
    double _det;
    static const TLineFitter _fitter;

  private:// Updated when fitted and accessed
    mutable double _chi2;
    mutable double _reducedChi2;
};

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

#ifdef TLine0_NO_INLINE
#define inline
#else
#undef inline
#define TLine0_INLINE_DEFINE_HERE
#endif

#ifdef TLine0_INLINE_DEFINE_HERE

inline
double
TLine0::a(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::a !!! fit not performed" << std::endl;
#endif
    return _a;
}

inline
double
TLine0::b(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::b !!! fit not performed" << std::endl;
#endif
    return _b;
}

inline
double
TLine0::distance(const TMLink & l) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::distance !!! fit not performed" << std::endl;
#endif
    double dy = fabs(_a * l.position().x() + _b - l.position().y());
    double invCos = sqrt(1. + _a * _a);
    return dy / invCos;
}

inline
void
TLine0::property(double a, double b, double det) {
    _a = a;
    _b = b;
    _det = det;
}

inline
unsigned
TLine0::objectType(void) const {
    return Line;
}

#endif

#undef inline

#endif /* TLine0_FLAG_ */

---