Reconstruction/TrkReco/include/TrkReco/RkFitCylinder.h

Go to the documentation of this file.

//-----------------------------------------------------------------------
// File from RkFit module
//
// Filename : RkFitCylinder.h
//------------------------------------------------------------------------
// Description :
// Cylinder is an Element whose shape is a cylinder.
//------------------------------------------------------------------------
// Modif :
//------------------------------------------------------------------------
 
#ifndef _DEFINE_CYLINDER_H_
#define _DEFINE_CYLINDER_H_
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#include "TrkReco/RkFitMaterial.h"
#include "TrkReco/TRunge.h"
/// Cylinder is an Element whose shape is a cylinder.
class RkFitCylinder {
protected:
  double ro_;
  double ri_;
  double zf_;
  double zb_;
 
public:
  const RkFitMaterial* material_;
  /// Constructor
  RkFitCylinder( const RkFitMaterial* material, double radius, double thick, double length,
                 double z0 ) {
    material_ = material;
    ri_       = radius;
    ro_       = radius + thick;
    zf_       = z0 + length * .5;
    zb_       = z0 - length * .5;
  }
 
  /// Find intersection with Helix
  double intersect( TRunge& track, HepPoint3D& x ) const;
 
  double intersect( TRunge& track, HepPoint3D& x, const HepPoint3D& point ) const;
  void   updateTrack( TRunge& track, double y[6] ) const;
 
  /// Check if the position x is inside the current cylinder
  bool isInside( const HepPoint3D& x ) const;
  bool isInside2( const HepPoint3D& x ) const;
 
  /// Extract radius of the cylinder
  virtual double radius( void ) const { return ro_; }
};
#endif