egs_rz.h

Go to the documentation of this file.

/*
###############################################################################
#
#  EGSnrc egs++ rz geometry headers
#  Copyright (C) 2016 Randle E. P. Taylor, Rowan M. Thomson,
#  Marc J. P. Chamberland, Dave W. O. Rogers
#
#  This file is part of EGSnrc.
#
#  EGSnrc is free software: you can redistribute it and/or modify it under
#  the terms of the GNU Affero General Public License as published by the
#  Free Software Foundation, either version 3 of the License, or (at your
#  option) any later version.
#
#  EGSnrc is distributed in the hope that it will be useful, but WITHOUT ANY
#  WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
#  FOR A PARTICULAR PURPOSE.  See the GNU Affero General Public License for
#  more details.
#
#  You should have received a copy of the GNU Affero General Public License
#  along with EGSnrc. If not, see <http://www.gnu.org/licenses/>.
#
###############################################################################
#
#  Author:          Randle Taylor, 2016
#
#  Contributors:    Marc Chamberland
#                   Rowan Thomson
#                   Dave Rogers
#
###############################################################################
#
#  egs_rz was developed for the Carleton Laboratory for Radiotherapy Physics.
#
###############################################################################
*/


#include <string>
#include <algorithm>
#include <fstream>
#include "../egs_nd_geometry/egs_nd_geometry.h"

#ifndef EGS_RZ_
#define EGS_RZ_

#ifdef WIN32

    #define EGS_RZ_EXPORT __declspec(dllexport)
    #define EGS_RZ_LOCAL

#else

    #ifdef HAVE_VISIBILITY
        #define EGS_RZ_EXPORT __attribute__ ((visibility ("default")))
        #define EGS_RZ_LOCAL  __attribute__ ((visibility ("hidden")))
    #else
        #define EGS_RZ_EXPORT
        #define EGS_RZ_LOCAL
    #endif

#endif

class EGS_RZ_EXPORT EGS_RZGeometry: public EGS_NDGeometry {

private:

    vector<EGS_Float> radii; 
    vector<EGS_Float> zbounds;
    vector<EGS_Float> reg_vol; /* calculated vol of each region */
    static string RZType;

public:

    enum DIRS {ZDIR, RDIR};

    EGS_RZGeometry(vector<EGS_BaseGeometry *> geoms,
                   vector<EGS_Float> rads, vector<EGS_Float> zbs, const string &name = "");

    ~EGS_RZGeometry() {};

    const string &getType() const {
        return RZType;
    };

    EGS_Float getBound(int idir, int ind);


    int getNRegDir(int dir);

    EGS_Float getVolume(int ireg);

};


#endif