---

CompositeMedium.cxx

Go to the documentation of this file.

// $Id: CompositeMedium.cxx,v 1.4 2001/09/04 19:26:01 atwood Exp $
// File: CompositeMedium.cxx
//
// This file is part of Gismo 2 
//
//

#ifdef __GNUG__
#pragma implementation
#endif

#include "gismo/CompositeMedium.h"

#include "gismo/Detector.h"
#include "gismo/Material.h"
#include "gismo/Field.h"
#include "gui/DisplayRep.h"

#include "geometry/Shape.h"
#include "geometry/Volume.h"
#include "geometry/Ray.h"
#include "facilities/error.h"

# include <algorithm> // for find
//              constructors

CompositeMedium::CompositeMedium(Medium * prnt, float size)
: Medium(prnt, size)
{}

CompositeMedium::CompositeMedium(Medium* prnt, Shape* vol, const char* matName, Detector* det)
: Medium(prnt, vol, matName, det)
{}

CompositeMedium::~CompositeMedium()
{
    deleteInnerMedia();
    if(_parent) {
        _parent->removeMedium(this);
       if(&_parent->field() != _field) delete _field;
    }
}
//               set tracking attributes: for this and all inner media
Medium& CompositeMedium::setField(Field* fld)
{
    for(iterator it=begin(); it !=end(); ++it )
       if(&((*it)->field()) == _field) (*it)->setField(fld);
    return Medium::setField(fld);
}

Medium& CompositeMedium::setKECutOff(float keCut)
{
    for(iterator it=begin(); it !=end(); ++it )
       if((*it)->kECutOff() == _keCutOff) (*it)->setKECutOff(keCut);
    return Medium::setKECutOff(keCut);
}

Medium& CompositeMedium::setMaxStep(float mxStep)
{
     for(iterator it=begin(); it !=end(); ++it )
       if( (*it)->maxStep() == maxStep()) (*it)->setMaxStep(mxStep);

    return Medium::setMaxStep(mxStep);
}


//                manage inner media
void CompositeMedium::deleteInnerMedia()
{
    for(reverse_iterator it=rbegin(); it !=rend(); ++it)
        delete *it;
}

Medium& CompositeMedium::addMedium ( Medium* nextMedium)
{
    push_back(nextMedium);
    nextMedium->setParent(this);
    return *nextMedium;
}

Medium& CompositeMedium::removeMedium (Medium* oldMedium)
{
    iterator it = std::find(begin(), end(), oldMedium);
    if( it != end())
         erase(it);
    return *oldMedium;
}

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

int CompositeMedium::isComposite() const {
  return 1;
}

//                       coordinate transformation
GeomObject&
CompositeMedium::transform(const CoordTransform& T)
{
    Medium::transform(T);
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->transform(T);
    return *this;
}

//              Detector data processing methods
void CompositeMedium::clear()
{
    Medium::clear();

    for(iterator it=begin(); it !=end(); ++it)
        (*it)->clear();
}

void CompositeMedium::generateResponse()
{
    Medium::generateResponse();
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->generateResponse();

}
void CompositeMedium::accept(DetectorVisitor& a)
{
    Medium::accept(a);
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->accept(a);
}

void
CompositeMedium::writeData(std::ostream & os)
{
    Medium::writeData(os);
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->writeData(os);
}

void
CompositeMedium::readData(std::istream & is)
{
    Medium::readData(is);
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->readData(is);
}
void CompositeMedium::notify()
{
    Medium::notify();
    for(iterator it=begin(); it !=end(); ++it)
        (*it)->notify();
}


//                    Methods used in Tracking

const Medium *
CompositeMedium::inside(const Point& r)const
{
    // if not inside the container, will just return parent
    const Medium* nextStuff = Medium::inside(r);
    if( nextStuff != this )
          return nextStuff;

    // OK, inside. default will be us
    nextStuff = this;

    for(const_iterator it=begin(); it !=end(); ++it ) {
       nextStuff = (*it)->inside(r) ;
           if(nextStuff != this) break;
    }
    return nextStuff;
}



double
CompositeMedium::distanceToLeave( const Ray& r,
                                const Medium* & newStuff,
                                double maxDist ) const
{

    // first get distance to container wall
    float distance = Medium::distanceToLeave(r, newStuff, maxDist);

    // If within ambig. distance of surface ... just return (WBA 31-8-2001)
    if(distance < Volume::Surface_EPSILON) return distance;

    // then check for inner Media and loop over them in case ray will hit one

    for(const_iterator it=begin(); it !=end(); ++it){
        const Medium* inner = *it;
        
        maxDist = (maxDist > distance) ? distance : maxDist;
        const Medium *nextStuff;

        double testDistance = inner->distanceToEnter( r, nextStuff,  maxDist );
        if(testDistance <= distance) {
            newStuff = nextStuff;
            distance = testDistance;
            if( distance==0 ) // means that it is already inside
            {
              //THB: no, it may be coming in!  newStuff = inner;
              break;
            }
        }
    }
    return distance;
}


// -----------------------------------------------------------------
//   Printing and GUI Methods
//
void CompositeMedium::printOn( std::ostream& os ) const
{
    Medium::printOn(os);
    int i=0;
    for(const_iterator it=begin(); it !=end(); ++it) {

        os << "\n" << title() <<" inner Medium# "<< ++i ;
           (*it)->printOn(os);
    }
}

void
CompositeMedium::printResponse(std::ostream& os)const
{
    Medium::printResponse(os);
    for(const_iterator it=begin(); it !=end(); ++it)
       (*it)->printResponse(os);
}

void CompositeMedium::createDetectorView(gui::DisplayRep& v)
{

    Medium::createDetectorView(v);
    for(iterator it=begin(); it !=end(); ++it )
       (*it)->createDetectorView(v.nested());
}

void CompositeMedium::createResponseView(gui::DisplayRep& v)
{
    Medium::createResponseView(v);
    for(iterator it=begin(); it !=end(); ++it )
       (*it)->createResponseView(v);
}

---