egs_gaussian_shape.cpp

Go to the documentation of this file.

/*
###############################################################################
#
#  EGSnrc egs++ guassian shape
#  Copyright (C) 2015 National Research Council Canada
#
#  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:          Iwan Kawrakow, 2005
#
#  Contributors:    Hannah Gallop
#
###############################################################################
*/
 
 
#include "egs_gaussian_shape.h"
#include "egs_input.h"
#include "egs_functions.h"
 
static bool EGS_GAUSSIAN_SHAPE_LOCAL inputSet = false;
static shared_ptr<EGS_BlockInput> EGS_GAUSSIAN_SHAPE_LOCAL shapeBlockInput = make_shared<EGS_BlockInput>("shape");
 
extern "C" {
 
    static void setInputs() {
        inputSet = true;
 
        shapeBlockInput->addSingleInput("library", true, "The type of shape, loaded by shared library in egs++/dso.", {"egs_gaussian_shape"});
        shapeBlockInput->addSingleInput("sigma", true, "1 or 2 or 3 inputs, the sigma for Gaussian spread along x, y and z, respectively.");
 
        auto shapePtr = shapeBlockInput->addBlockInput("shape");
        setShapeInputs(shapePtr);
    }
 
    EGS_GAUSSIAN_SHAPE_EXPORT string getExample() {
        string example;
        example = {
            R"(
    # Example of egs_gaussian_shape
    #:start shape:
        library = egs_gaussian_shape
        :start shape:
            definition of the shape to be smeared
        :stop shape:
        sigma = 1, 2 or 3 inputs
    :stop shape:
)"};
        return example;
    }
 
    EGS_GAUSSIAN_SHAPE_EXPORT shared_ptr<EGS_BlockInput> getInputs() {
        if(!inputSet) {
            setInputs();
        }
        return shapeBlockInput;
    }
 
    EGS_GAUSSIAN_SHAPE_EXPORT EGS_BaseShape *createShape(EGS_Input *input,
            EGS_ObjectFactory *f) {
        if (!input) {
            egsWarning("createShape(circle): null input?\n");
            return 0;
        }
        vector<EGS_Float> sigma;
        int err = input->getInput("sigma",sigma);
        if (err) {
            egsWarning("createShape(gaussian shape): no 'sigma' input\n");
            return 0;
        }
        EGS_Input *ishape = input->takeInputItem("shape",false);
        EGS_BaseShape *shape = 0;
        if (ishape) {
            shape = EGS_BaseShape::createShape(ishape);
            delete ishape;
        }
        if (!shape) {
            string shape_name;
            int err = input->getInput("shape name",shape_name);
            if (err) {
                egsWarning("createShape(gaussian shape): no inline shape definition"
                           " and no 'shape name' keyword\n");
                return 0;
            }
            shape = EGS_BaseShape::getShape(shape_name);
            if (!shape) {
                egsWarning("createShape(gaussian shape): no shape named %s "
                           "exists\n",shape_name.c_str());
                return 0;
            }
        }
        EGS_GaussianShape *s;
        if (sigma.size() == 1) {
            s = new EGS_GaussianShape(shape,sigma[0],0,0,"",f);
        }
        else if (sigma.size() == 2) {
            s = new EGS_GaussianShape(shape,sigma[0],sigma[1],0,"",f);
        }
        else {
            s = new EGS_GaussianShape(shape,sigma[0],sigma[1],sigma[2],"",f);
        }
        s->setName(input);
        s->setTransformation(input);
        return s;
    }
 
}