---

Medium.h

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/gismo/gismo/Medium.h,v 1.3 2001/01/23 03:41:32 burnett Exp $
//


#ifndef GISMO_MEDIUM_H
#define GISMO_MEDIUM_H

#ifdef __GNUG__
#pragma interface
#endif

#include "geometry/GeomObject.h"

class Material;
class Detector;
namespace gui {class DisplayRep; }
class Field;
class Shape;
class Ray;
class Vector;
class Point;
class DetectorVisitor;

class Medium  : public GeomObject
{
public:
    
    Medium(Medium* parent =0, float size=100);

    Medium(Medium* parent, Shape* vol, const char* mat= "vacuum", Detector * det= 0);

    Medium(const Medium& old);
    
    virtual ~Medium();
    
    virtual Medium&  addMedium(Medium* nextMedium);
    virtual Medium& removeMedium (Medium* oldMedium);
    
    
    Medium&  setParent(Medium* p)       {_parent = p; return *this;}
    Medium&  setVolume(Shape* vol)      {_volume = vol;return *this;}
    Medium&  setMaterial(Material *mat) {_material = mat;return *this;}
    Medium&  setDetector(Detector *det) {_detector = det;return *this;}
    Medium&  setTitle(const char *newTitle);
    
    virtual Medium& setKECutOff(float keCut);
    virtual Medium& setMaxStep(float mxStep);
    virtual Medium& setField(Field* );
    
    
    
    virtual const char* nameOf() const;
    virtual int  isComposite()const;
    const char* title()  const  {return _title? _title : "no title";}
    
    Shape&      volume()        {return *_volume;}
    const Shape& volume()const   {return *_volume;}
    Field&      field()const    {return *_field;}
    Material&   material()const {return *_material;}
    
    const Medium* getParent() const {return _parent;}
    const Medium* parent()const     {return _parent;}
    
    float       kECutOff()const  {return _keCutOff;}
    float       maxStep()const   {return _maxStep;}
    Detector*   detector()const  {return _detector;}
    
        GeomObject& transform(const CoordTransform& );
    
    //  Methods for tracking
    
    Ray* CreateRay(const Point& position, const Vector& momentum, float charge,float maxdist=0 )const;
    

    virtual double distanceToLeave( const Ray& r,       const Medium*& nextMedium, double maxStep ) const;
    
    virtual const Medium * inside(const Point& r)const;
    
    virtual double distanceToEnter( const Ray& r,               const Medium*& newstuff, double ) const;
    
    static void setLastMedium(const Medium *last) {lastMedium = last;}
    
    
    virtual void clear();
    virtual void generateResponse();
    virtual void accept(DetectorVisitor&);
    virtual void readData(std::istream&);
    virtual void writeData(std::ostream&);
    
    virtual void notify(); 
    
    virtual void printOn( std::ostream& os = std::cout ) const;
    virtual void printResponse(std::ostream& os= std::cout) const;
    
    virtual void createDetectorView(gui::DisplayRep& v);
    virtual void createResponseView(gui::DisplayRep& v);
    
    
protected:
    void set_defaults();
    static unsigned s_count;  // keep informal track of number of Medium's constructed

    // Medium contents:

    Shape*      _volume;        // Container Volume
    Material*   _material;      // Material description
    Detector*   _detector;      // (optional) Detector
    Field*      _field;         // Magnetic Field, null for no field
    
    Medium*     _parent;        // The medium in which this exists, null for root 
    
    // attributes for tracking
    float       _keCutOff;      // Cut-off momentum for swimming
    float       _maxStep;       // Step size limit
    
    char*       _title;         // User setable title
    
    static const Medium *lastMedium; // pointer to the last distanceToLeave Medium
       
};
//
// Operator overloads for <<
inline std::ostream& operator<< (std::ostream& os, const Medium &m){m.printOn(os);return os;}
inline std::ostream& operator<< (std::ostream& os, const Medium *m){m->printOn(os);return os;}

#endif

---