---

LayeredTube.h

Go to the documentation of this file.

// $Id: LayeredTube.h,v 1.3 2001/01/25 22:24:32 burnett Exp $
//
// This file is part of Gismo2
//
// Define a LayeredTube, a medium that is a tube with cylindrical layers
// The layers are added by the function member
//
//     addLayer(float thickness, const char* material, Detector detector);
//
// The layers are members of a private class LayeredTube::Layer, itself a Medium
// subclass.
//
//
#ifndef LAYEREDTUBE_H
#define LAYEREDTUBE_H



#include "gismo/CompositeMedium.h"

#include "geometry/Point.h"
#include "geometry/Tube.h"


class LayeredTube : public CompositeMedium
{
  public:

   LayeredTube(Medium* parent, float dz, float rzero);
   // create a new LayeredTube
   // +r direction

   Medium& addLayer(float thickness, const char* material="vacuum",
                         Detector* det=0);
   // add a new layer after the previous one. The enclosing Box will be increased
   // in the z-direction by thickness.
   // returns reference to the new medium for further modification

   const Medium* inside(const Point& p)const;
   // return the parent, if outside, or the appropriate layer if inside 


   const Tube& tube()const{ return (const Tube&)volume();}
   Tube& tube(){ return (Tube&)volume();}
   const Point& origin()const{ return m_origin;}
   const Vector& axis()const{ return tube().axis();}

   GeomObject& transform(const CoordTransform& T);
   void createDetectorView(gui::DisplayRep& view);
   const char* nameOf()const{return "LayeredTube";}


protected:
   double distanceToEnter(const Ray&, const Medium*&, double)const;
   double distanceToLeave(const Ray&, const Medium*& ,double) const;

private:
   float m_dz, m_r_inner, m_r_outer;  // dimensions

   class Layer : public Medium
   {
      friend class LayeredTube;
      Layer(Medium*, const char*, Detector*);

      const Medium* inside(const Point&)const{return 0;}
      const char* nameOf()const{return "LayeredTube::Layer";}
      double distanceToLeave(const Ray& r, const Medium*& , double )const;

      GeomObject& transform(const CoordTransform& T);
      void createDetectorView(gui::DisplayRep& v);

      void printOn(std::ostream&)const;

      const LayeredTube* parent()const;
      // override inherited function to get right type of parent

      Layer *previous, *next;
      float min_r, max_r;  // range of r

      const Point& origin()const;
      const Vector& axis()const;
      const Tube& parentTube()const;
      // access to stuff of the parent

      double localRadius(const Point& p)const;
      // calcualtes the radius with respect to axis in local coords
   };
    // internal class used to store linked list of layers
   Layer * first, *last;

   Point m_origin;  // keep track of origin for convenience

   const Layer * select(const Point& x)const;
   friend class Layer;

};

//                        inlines
inline const LayeredTube*
LayeredTube::Layer::parent()const{  return (const LayeredTube*)_parent;}

inline const Point&
LayeredTube::Layer::origin()const{  return parent()->origin();}

inline const Vector&
LayeredTube::Layer::axis()const{ return parent()->axis();}

inline const Tube&
LayeredTube::Layer::parentTube()const { return parent()->tube();}

inline double
LayeredTube::Layer::localRadius(const Point& x) const
{
    Vector r = x - origin();
    return sqrt( r.mag2() - sqr(r*axis()) );
}

#endif

---