---

LayeredBox.cxx

Go to the documentation of this file.

// $Id: LayeredBox.cxx,v 1.3 2001/01/23 03:41:32 burnett Exp $
//  Author: Toby Burnett

#include "gismo/LayeredBox.h"
#include "geometry/Plane.h"
#include "gismo/Detector.h"

#include "gui/DisplayRep.h"
#include <cassert>
#include <cstdlib>
#include <stdio.h>

LayeredBox::LayeredBox(Medium* mother, double dx, double dy)
        :LayeredMedium(mother, 6)
        , m_depth(dy)
        , m_width(dx)
        , m_height(0)
    {

        const Point& p = m_vol->center();
        int n=0;
        m_vol->setSurface(n++, new Plane(p, Vector(0,0,-1),0) );
        m_vol->setSurface(n++, new Plane(p, Vector(0,0, 1),0) );
        m_vol->setSurface(n++, new Plane(p, Vector(-dx/2,0,0)) );
        m_vol->setSurface(n++, new Plane(p, Vector( dx/2,0,0)) );
        m_vol->setSurface(n++, new Plane(p, Vector(0,-dy/2,0)) );
        m_vol->setSurface(n++, new Plane(p, Vector(0, dy/2,0)) );

    }

Medium&
LayeredBox::addLayer(double thickness, const char* material, Detector* det)
{
    // create medium without a volume
    Medium* med = new Medium(this, 0, material, det);


    // now have LayeredMedium add it to its list
    return addLayer(med, thickness);
}
static inline char* itoa(int i)  //little local thing to format integers
{   static char buf[20];
    sprintf(buf,"%d",i);
    return buf;
}

Medium&
LayeredBox::addLayer(Medium* inner_medium, double thickness, const char* mat)
{
    assert( inner_medium->parent() == this); // must have been added

    Medium* added_medium = inner_medium;

    if( inner_medium->isComposite() ){
        // have to isolate: create new Composite to go between
        pop_back(); // unlink from list
        added_medium = new CompositeMedium(this, (Shape*)0, mat);
        added_medium->addMedium(inner_medium);
    }

    // set default title to inner medium
    added_medium->setTitle(
        (std::string("BoxLayer#")+itoa(innerMediaCount())).c_str() );

    // follow a procedure here so that the enlarged medium winds up centered

    // first translate the new guy to the edge
    added_medium->moveZ(m_height/2+thickness/2);

    // translate the medium back, including new one
    moveZ(-thickness/2.0);

    m_height+=thickness;

    LayeredMedium::addLayer(added_medium, thickness);
    return *inner_medium;
}

void
LayeredBox::createDetectorView(gui::DisplayRep& view)
{
    Vector vx = m_vol->surface(3).direction() * (m_width/2),
           vy = m_vol->surface(5).direction() * (m_depth/2),
           vz = m_vol->surface(1).direction() * (m_height);
    const Point& c = m_vol->center();
    Point vertex[] =
    {   c+vx+vy,
        c-vx+vy,
        c+vx-vy,
        c-vx-vy,
        c+vx+vy+vz,
        c-vx+vy+vz,
        c+vx-vy+vz,
        c-vx-vy+vz};
    view.move_to(vertex[0]);
    view.line_to(vertex[2]);
    view.line_to(vertex[3]);
    view.line_to(vertex[1]);
    view.line_to(vertex[0]);

    view.line_to(vertex[4]);
    view.line_to(vertex[6]);
    view.line_to(vertex[7]);
    view.line_to(vertex[5]);
    view.line_to(vertex[4]);

    view.move_to(vertex[3]);
    view.line_to(vertex[7]);

    view.move_to(vertex[1]);
    view.line_to(vertex[5]);

    view.move_to(vertex[2]);
    view.line_to(vertex[6]);
    //-------------------------------------------------------
    // draw inner boundaries, as a nested view
    gui::DisplayRep& nview = view.nested();
    for( unsigned int i=1; i< size(); ++i){
        Vector vz = referenceSurface().direction()* (-m_offset[i]);
        nview.move_to(vertex[0]+vz);
        nview.line_to(vertex[2]+vz);
        nview.line_to(vertex[3]+vz);
        nview.line_to(vertex[1]+vz);
        nview.line_to(vertex[0]+vz);
    }

    // draw inner media, nested again
    for(iterator it=begin(); it !=end(); ++it )
       (*it)->createDetectorView(nview.nested());
}

---