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

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// $Id: TCircle.h,v 1.13 2012/05/28 05:16:29 maoh Exp $
//-----------------------------------------------------------------------------
// Filename : TCircle.h
// Section  : Tracking
// Owner    : Yoshi Iwasaki
// Email    : yoshihito.iwasaki@kek.jp
//-----------------------------------------------------------------------------
// Description : A class to represent a circle in tracking.
//               See http://bsunsrv1.kek.jp/~yiwasaki/tracking/
//-----------------------------------------------------------------------------

#ifndef TCircle_FLAG_
#define TCircle_FLAG_

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

#include <string>

#define HEP_SHORT_NAMES

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

#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TTrackBase.h"
#include "TrkReco/TCircleFitter.h"
//#include "lpav/Lpav.h"
//#include "TrkReco/Lpav.h"
#include "TrackUtil/Lpav.h"

class TCircle : public TTrackBase {

  public:
    TCircle(const AList<TMLink> & links);

    virtual ~TCircle();

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

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

    const HepPoint3D & center(void) const;

    double radius(void) const;

    double pt(void) const;

    double impact(void) const;

    double charge(void) const;

    double weight(const TMLink & l) const;

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

    int fitForCurl(int ipConst = 0);

    void property(double charge, double radius, HepPoint3D center);

  private:// Updated when fitted
    Lpav _circle;
    double _charge;
    double _radius;
    HepPoint3D _center;
    static const TCircleFitter _fitter;
};

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

#ifdef TCircle_NO_INLINE
#define inline
#else
#undef inline
#define TCircle_INLINE_DEFINE_HERE
#endif

#ifdef TCircle_INLINE_DEFINE_HERE

inline
const HepPoint3D &
TCircle::center(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::center !!! fit not performed" << std::endl;
#endif
    return _center;
}

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

inline
double
TCircle::pt(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::pt !!! fit not performed" << std::endl;
#endif
   // return fabs(_radius) / Helix::ConstantAlpha;
    return fabs(_radius) / 333.564095;
}

inline
double
TCircle::impact(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::impact !!! fit not performed" << std::endl;
#endif
    return fabs(sqrt(_center.x() * _center.x() +
                     _center.y() * _center.y())
                - fabs(_radius));
}

inline
double
TCircle::charge(void) const {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::charge !!! fit not performed" << std::endl;
#endif
    return _charge;

}

inline
void
TCircle::property(double c, double r, HepPoint3D e) {
    _charge = c;
    _radius = r;
    _center = e;
}

inline
unsigned
TCircle::objectType(void) const {
    return Circle;
}

#endif

#undef inline

#endif /* TCircle_FLAG_ */

---