Circle.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
//-----------------------------------------------------------------------------
// Description : A class to represent a circle
//-----------------------------------------------------------------------------
 
#ifndef TRGCDCCircle_FLAG_
#define TRGCDCCircle_FLAG_
 
#include "trg/trg/Point2D.h"
#include "trg/cdc/TrackBase.h"
 
#ifdef TRGCDC_SHORT_NAMES
#define TCCircle TRGCDCCircle
#endif
 
namespace Belle2 {
  class TRGCDCHoughPlane;
  class TRGCDCCircleFitter;
 
  class TRGCDCCircle : public TRGCDCTrackBase {
 
  public:
 
    explicit TRGCDCCircle(const std::vector<TRGCDCLink*>& links);
 
    TRGCDCCircle(double r,
                 double phi,
                 double charge,
                 const TRGCDCHoughPlane& plane);
 
    virtual ~TRGCDCCircle();
 
  public:
 
    virtual unsigned objectType(void) const override;
 
    const TRGPoint2D& center(void) const;
 
    double pt(void) const override;
 
    double radius(void) const;
 
    virtual void dump(const std::string& message = std::string(""),
                      const std::string& prefix = std::string("")) const override;
 
    void property(double charge,
                  double radius,
                  HepGeom::Point3D<double> center);
//  void property(double charge, double radius, TPoint2D & center);
 
  public: // Utility functions
 
    int approach2D(TRGCDCLink&) const override;
 
  private:
 
    TRGPoint2D _center;
 
    double _radius;
 
    const TRGCDCHoughPlane* _plane;
 
    static const TRGCDCCircleFitter _fitter;
 
  };
 
//-----------------------------------------------------------------------------
 
  inline
  double
  TRGCDCCircle::pt(void) const
  {
    const double ConstantAlpha = 222.376063; // for 1.5T
    return _radius / ConstantAlpha;
  }
 
  inline
  const TRGPoint2D&
  TRGCDCCircle::center(void) const
  {
    return _center;
  }
 
  inline
  double
  TRGCDCCircle::radius(void) const
  {
    return _radius;
  }
 
  inline
  unsigned
  TRGCDCCircle::objectType(void) const
  {
    return TRGCDCCircleType;
  }
 
  inline
  void
  TRGCDCCircle::property(double c, double r, HepGeom::Point3D<double> e)
  {
    charge(c);
    _radius = r;
    _center = e;
  }
 
} // namespace Belle2
 
#endif /* TRGCDCCircle_FLAG_ */