/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/MdcPatRec/MdcGeom/MdcGeom-00-01-17/MdcGeom/MdcSagTraj.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: MdcSagTraj.h,v 1.3 2009/12/17 00:38:40 zhangy Exp $
//
// Description:
//      Class MdcSagTraj.
//      the flight length at the moment is the flight length of a linear 
//      trajectory, the sag is approximated with a parabola (for the
//      moment, maybe at a later time we can replace with a more accurate
//      description, cosh for example)
//      
//
// Environment:
//      Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//      R. Stroili              originator
//      
//
// Copyright Information:
//      Copyright (C) 1998      INFN & Padova University
// 
// History:
//      Migration for BESIII MDC
//
//------------------------------------------------------------------------

#ifndef DCHSAGTRAJ_H
#define DCHSAGTRAJ_H

//-------------
// C Headers --
//-------------
extern "C" {
}

//---------------
// C++ Headers --
//---------------

//----------------------
// Base Class Headers --
//----------------------
#include "MdcGeom/TrkGeomTraj.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#include "CLHEP/Vector/ThreeVector.h"
using CLHEP::Hep3Vector;

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
class Trajectory;

//              ---------------------
//              -- Class Interface --
//              ---------------------

class MdcSagTraj : public TrkGeomTraj {

public:

  // Constructors
  MdcSagTraj(const double sag, const HepPoint3D& point1, const HepPoint3D& point2);
  MdcSagTraj(const MdcSagTraj& traj);
 
  MdcSagTraj* clone() const;
  // Destructor
  virtual ~MdcSagTraj( );

  // Operators
  MdcSagTraj& operator=(const MdcSagTraj&);

  double sag( void ) const { return _sag; }

// needed implementations for intersection with a Surface
  HepPoint3D position( double ) const;
  Hep3Vector direction( double ) const;
  double curvature( double f = 0. ) const;
  Hep3Vector delDirect( double ) const;
  const Hep3Vector& rawDirection( void ) const { return _direction; }
  void  getInfo(double fltLen, HepPoint3D& , Hep3Vector& direction) const;
  void  getInfo(double fltLen, HepPoint3D& , Hep3Vector& direction, 
                        Hep3Vector& delDirect) const;

  virtual double distTo1stError(double s, double tol, int pathDir) const;
  virtual double distTo2ndError(double s, double tol, int pathDir) const;
  double deltaY(double fltLen) const { return (_a*fltLen+_b)*fltLen; }
  // Support Visitor pattern (see TrkGeomTraj.h)
  void accept(TrkGeomTrajVisitor& visitor) const;

//    virtual int operator==( const MdcSagTraj& ) const;
//            int operator!=( const MdcSagTraj& ) const;

private:

  // Data members
  double _sag;
  double _a;
  double _b;
  double _length;
//   double _a_H;
//   double _b_H;
  HepPoint3D _start; // where the trajectory starts
  HepPoint3D _stop; // where the trajectory stops
  Hep3Vector _direction; // direction (unit) vector for 
                         // null sag (straight line)

  Hep3Vector deviation(double) const; // displacement from the line trajectory
                                      // at a given flightlength
  //Note: if your class needs a copy constructor or an assignment operator, 
  //  make one of the following public and implement it.
//   MdcSagTraj( const MdcSagTraj& );       // Copy Constructor
//   MdcSagTraj&       operator= ( const MdcSagTraj& );  // Assignment op

};

#endif // DCHSAGTRAJ_H

---