---

LayeredMedium.h

Go to the documentation of this file.

//  $Id: LayeredMedium.h,v 1.4 2001/01/23 03:41:32 burnett Exp $
//

#ifndef LAYEREDMEDIUM_H
#define LAYEREDMEDIUM_H

#include "gismo/CompositeMedium.h"
#include "geometry/Volume.h"

class LayeredMedium : public CompositeMedium {
public:

    virtual ~LayeredMedium();

    virtual double distanceToLeave( const Ray& r, const Medium*& nextMedium, double maxStep ) const;
    // returns distance to leave this medium (perhaps to enter an inner Medium if implemented by
    // a CompositeMedium), setting nextMedium if a boundary is reached. Will be limited to maxStep

    virtual double distanceToEnter( const Ray& r, const Medium*& nextMedium, double maxStep) const;

    const Medium* inside(const Point& x)const;

    virtual void createDetectorView(gui::DisplayRep& v){};


//protected:
    LayeredMedium(Medium* mother, int num_sufaces, double inner_offset=0); 
    // abstract class: ctor only accessible to subclasses 


    Medium& addLayer(Medium * inner_medium, double thickness);
    // add a new layer, enlarging the current structure. inner_medium will be associated
    // with the new layer, and any internal Shapes it has must fit within the layer. Subclass
    // is responsible for shifting it into place.


    virtual const Medium* select(const Point& x)const;
    // return the correct inner medium, given the point x
    // virtual to allow subclass to make efficient calculation

    class Slice; 
    class Volume;

    //-------------------------------------------------------------------
    // Define a special volume that handles the logic of an individual layer
    class Slice : public ::Volume {
    public:
        Slice(LayeredMedium& mother, int index);
        ~Slice();

        const char* nameOf()const{return "LayeredMedium::Slice";}

        double distanceToLeave(const Ray& r, double d)const;
        // override in order to reverse inner surface, flag the mother

        const ::Volume& setMotherVolume()const;
        // set up the mother volume for this slice, return reference to it

        const ::Volume& resetMotherVolume()const;
        // reset the mother volume  to have again initial size 

        void printOn(std::ostream& out)const;

        virtual GeomObject& transform(const CoordTransform& ){return *this;};
        // override since don't own any objects

    private:
        LayeredMedium& m_mother;
        int m_index;    // index of this slice: -1 for the mother
        

    };
    //-------------------------------------------------------------------

    const Surface& referenceSurface()const;
    // access to the inner surface used to define the division

    //-------------------------------------------------------------------
    // special Volume for this Medium 
    class Volume : public ::Volume {
    public:
        Volume(int n);

        void setSurface(int n, Surface* surf){ operator[](n)=surf;}
        Surface& surface(int n){return *operator[](n);}
        // full access to internal surfaces 
        
        void setOffsets(double inner, double outer);
        // adjust the offsets

        const char* nameOf()const{return "LayeredMedium::Volume";}
    private:
    };
    //-------------------------------------------------------------------

    friend class Slice;

    std::vector<double> m_offset;
    // list of offsets for the inner and outer surfaces that define the slices

    Volume* m_vol;
    // separate pointer to the Shape to avoid casts
private:

    int m_select_index;
    int m_last_boundary; 
    // tags: keep track of which layer we are in, communicate which surface left
};

#endif

---