---

SlicedCone.cxx

Go to the documentation of this file.

//     $Id: SlicedCone.cxx,v 1.3 2001/01/23 03:41:32 burnett Exp $
//  Author: Toby Burnett
#include "gismo/SlicedCone.h"
#include "geometry/Plane.h"
#include "geometry/Conic.h"
#include "gismo/Detector.h"

#include "gui/DisplayRep.h"
#include "geomrep/ArcRep.h"
#include "geomrep/ConeRep.h"
#include <cassert>
#include <strstream>
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
SlicedCone::SlicedCone(Medium* mother, const Cone& cone)
//------------------------------------------------------
    :LayeredMedium(mother, cone.surfaceCount() )
    , m_length(0)
{
    const Point& p = m_vol->center();
    static Vector z_axis(0,0,1);
    int n=0;
    // front and back planes first
    m_vol->setSurface(n++, new Plane(p, Vector(0,0,-1),0) );
    m_vol->setSurface(n++, new Plane(p, Vector(0,0, 1),0) );

    // now transverse surfaces
    double len = cone.length();
    m_outer_radius = cone.outerRadius(-len/2);
    m_inner_radius = cone.innerRadius(-len/2);
    m_outer_slope = (cone.outerRadius(1-len/2)-m_outer_radius);
    m_inner_slope = (cone.innerRadius(1-len/2)-m_inner_radius);
        
    m_vol->setSurface(n++,
             new Conic(p, z_axis,  m_outer_radius, m_outer_slope) );
    if( cone.surfaceCount()>3) {
        m_vol->setSurface(n++,
             new Conic(p, z_axis, -m_inner_radius,-m_inner_slope) );

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

    // set default title to inner medium
    std::ostrstream title;
    title << "SlicedCone Layer#" << innerMediaCount() << '\0';
    med->setTitle(title.str());

    // follow a procedure here so that the enlarged medium winds up centered
    if( m_length>0) moveZ(-m_length/2);

    // now have LayeredMedium add it to its list
    Medium & t = LayeredMedium::addLayer(med, thickness);
    moveZ(m_length/2-thickness/2);
    m_length += thickness;
    return t;
}


//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
void
SlicedCone::createDetectorView(gui::DisplayRep& v)
//-------------------------------------------------
{
    Vector a  = axis();
    Point c = center();


    //create a temporary Cone just to use its drawing
    double h = m_length/2;
    Cone cone(m_length,
        m_inner_radius, m_outer_radius,
        m_inner_radius+2*h*m_inner_slope, m_outer_radius+2*h*m_outer_slope);

    Vector z_axis(0,0,1);

    if( a.z() != 1.0 ) {
        // has been rotated: rotate this cone
        Vector r= z_axis.cross(a);
        if( r.mag()<1e-6) r = Vector(0,1,0); // if 2pi
        cone.rotate( acos(a.z()), r.unit());
    }
    cone.move(c+0.5*m_length*a);
    v.append(ConeRep(cone));
    //-------------------------------------------------------
    // draw inner boundaries, as a nested view
    Vector xhat (1,0,0);
    xhat = xhat.cross(a).unit();
    Vector yhat=xhat.cross(a).unit();

    gui::DisplayRep& nv = v.nested();

    for( std::vector<double>::const_iterator it=m_offset.begin();
                it < m_offset.end(); ++it){
        double d = *it;
        double  radius = m_outer_radius + (d)*m_outer_slope;
        ArcRep arc(c+a*d, radius, xhat, yhat, 2*M_PI);
        arc.update();
        nv.append(arc);
        if( m_inner_radius ==0) continue;
        radius = m_inner_radius + (d)*m_inner_slope;
        ArcRep inner_arc(c+a*d, radius, xhat, yhat, 2*M_PI);
        inner_arc.update();
        nv.append(inner_arc);
    }

}

---