In This Package:

SolidTubs.h

Go to the documentation of this file.

// $Id: SolidTubs.h,v 1.20 2008/10/28 12:04:37 cattanem Exp $
// ===========================================================================
#ifndef     DETDESC_SOLIDTUBS_H
#define     DETDESC_SOLIDTUBS_H 1  
// STD and STL includes 
#include <cmath>
#include <iostream>
// Units
#include "GaudiKernel/SystemOfUnits.h" 
// Geometry Definitions includes 
#include "GaudiKernel/Point3DTypes.h" 
#include "GaudiKernel/Vector3DTypes.h"
// DetDesc 
#include "DetDesc/SolidBase.h" 
// forward declarations
template <class TYPE>
class SolidFactory;
class MsgStream;
//

class SolidTubs: public virtual SolidBase
{
  friend class SolidFactory<SolidTubs>;
  
public:
  
  SolidTubs( const std::string& name                             ,
             const double       ZHalfLength                      , 
             const double       OuterRadius                      ,
             const double       InnerRadius   =   0.0            , 
             const double       StartPhiAngle =   0.0 * Gaudi::Units::degree, 
             const double       DeltaPhiAngle = 360.0 * Gaudi::Units::degree,
             const int          CoverModel    =   0              );     

  virtual ~SolidTubs();

  inline std::string typeName () const { return "SolidTubs"; };

  bool isInside ( const Gaudi::XYZPoint& point ) const; 
  bool isInside ( const Gaudi::Polar3DPoint& point ) const ;
  bool isInside ( const Gaudi::RhoZPhiPoint& point ) const ;
  
  const ISolid* cover () const;
  
  virtual std::ostream& printOut ( std::ostream& os = std::cout ) const;
  
  virtual MsgStream&    printOut ( MsgStream&    os             ) const;
  
  virtual unsigned int intersectionTicks( const Gaudi::XYZPoint & Point,
                                          const Gaudi::XYZVector& Vector,
                                          ISolid::Ticks& ticks  ) const ;

  virtual unsigned int intersectionTicks( const Gaudi::Polar3DPoint  & Point,
                                          const Gaudi::Polar3DVector & Vector,
                                          ISolid::Ticks     & ticks) const ; 

  virtual unsigned int intersectionTicks( const Gaudi::RhoZPhiPoint  & Point,
                                          const Gaudi::RhoZPhiVector & Vector,
                                          ISolid::Ticks     & ticks) const ;   

  virtual unsigned int intersectionTicks( const Gaudi::XYZPoint& Point,
                                          const Gaudi::XYZVector & Vector,
                                          const ISolid::Tick& tickMin,
                                          const ISolid::Tick& tickMax,
                                          ISolid::Ticks& ticks   ) const ; 

  virtual unsigned int intersectionTicks( const Gaudi::Polar3DPoint& Point,
                                          const Gaudi::Polar3DVector & Vector,
                                          const ISolid::Tick& tickMin,
                                          const ISolid::Tick& tickMax,
                                          ISolid::Ticks& ticks   ) const ; 

  virtual unsigned int intersectionTicks( const Gaudi::RhoZPhiPoint& Point,
                                          const Gaudi::RhoZPhiVector & Vector,
                                          const ISolid::Tick& tickMin,
                                          const ISolid::Tick& tickMax,
                                          ISolid::Ticks& ticks   ) const ; 
    

  inline double innerRadius     () const 
  { return m_tubs_innerRadius   ; };
  inline double outerRadius     () const 
  { return m_tubs_outerRadius   ; };
  inline double zHalfLength     () const 
  { return m_tubs_zHalfLength   ; };
  inline double startPhiAngle   () const 
  { return m_tubs_startPhiAngle ; }; 
  inline double deltaPhiAngle   () const 
  { return m_tubs_deltaPhiAngle ; }; 
  inline double innerDiameter   () const 
  { return m_tubs_innerRadius   * 2 ; };
  inline double outerDiameter   () const 
  { return m_tubs_outerRadius   * 2 ; };
  inline double zLength         () const 
  { return m_tubs_zHalfLength   * 2 ; };
  inline double endPhiAngle     () const 
  { return m_tubs_startPhiAngle + m_tubs_deltaPhiAngle ; }; 

  Ticks::size_type maxNumberOfTicks() const { return 4 ; }

protected:
  
  template <class aPoint>
  inline bool insideRho ( const aPoint& point ) const ;
  
  inline bool insidePhi ( const double phi /* [-pi,pi] */ ) const ;

  template <class aPoint>
  inline bool insidePhi ( const aPoint& point ) const ;  
  
  inline bool noPhiGap() const { return m_noPhiGap ; }
  
protected:
  
  SolidTubs( const std::string& name = "Anonymous Tubs");

  void setBP();
  
private:
  
  SolidTubs           ( const SolidTubs & );  
  SolidTubs& operator=( const SolidTubs & );  

  template <class aPoint>
  bool isInsideImpl(const aPoint& point) const;  

  template<class aPoint, class aVector>
  unsigned int intersectionTicksImpl( const aPoint  & Point,
                                      const aVector & Vector,
                                      const ISolid::Tick& tickMin,
                                      const ISolid::Tick& tickMax,
                                      ISolid::Ticks&  ticks) const;

  template<class aPoint, class aVector>
  unsigned int intersectionTicksImpl( const aPoint  & Point,
                                      const aVector & Vector,
                                      ISolid::Ticks& ticks ) const;

private:

  // parameters of tube segment 
  double                  m_tubs_zHalfLength   ;
  double                  m_tubs_outerRadius   ;
  double                  m_tubs_innerRadius   ;
  double                  m_tubs_startPhiAngle ;
  double                  m_tubs_deltaPhiAngle ;
  // model for constructing the cover solid 
  int                     m_tubs_coverModel    ;  
  //
  bool                    m_noPhiGap           ;
};

  
// ===========================================================================
// ===========================================================================
template<class aPoint>
inline bool SolidTubs::insideRho ( const aPoint& point ) const 
{
  const double rho2 = point.perp2();
  if( rho2 > outerRadius() * outerRadius() ) { return false ; }
  if( 0    < innerRadius() && 
      rho2 < innerRadius() * innerRadius() ) { return false ; }
  // inside! 
  return true ;
};
// ===========================================================================


// ===========================================================================
// ===========================================================================
inline bool SolidTubs::insidePhi ( const double phi /* [-pi,pi] */ ) const
{
  return 
    noPhiGap() ||
    ( startPhiAngle ()                   <= phi &&
      startPhiAngle () + deltaPhiAngle() >= phi     ) ||
    ( startPhiAngle ()                   <= phi + 2*M_PI &&
      startPhiAngle () + deltaPhiAngle() >= phi + 2*M_PI ) ;
};

// ===========================================================================
// ===========================================================================
template< class aPoint>
inline bool SolidTubs::insidePhi ( const aPoint& point ) const 
{
  return noPhiGap() || insidePhi( point.phi() ) ;
};
// ===========================================================================
// The END 
// ===========================================================================
#endif 
// ===========================================================================

---