/Users/marc/Developer/EGSnrc/HEN_HOUSE/user_codes/egs_brachy/egs_brachy/phsp.h

Go to the documentation of this file.

/*
################################################################################
#
#  egs_brachy phsp.h
#  Copyright (C) 2016 Rowan Thomson, Dave Rogers, Randle Taylor, and Marc
#  Chamberland
#
#  This file is part of egs_brachy
#
#  egs_brachy 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.
#
#  egs_brachy 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:
#  <http://www.gnu.org/licenses/>.
#
################################################################################
#
#  When egs_brachy is used for publications, please cite our paper:
#  M. J. P. Chamberland, R. E. P. Taylor, D. W. O. Rogers, and R. M. Thomson,
#  egs brachy: a versatile and fast Monte Carlo code for brachytherapy,
#  Phys. Med. Biol. 61, 8214-8231 (2016).
#
################################################################################
#
#  Author:        Randle Taylor, 2016
#
#  Contributors:  Marc Chamberland
#                 Dave Rogers
#                 Rowan Thomson
#
################################################################################
*/
 
#include "egs_functions.h"
#include "egs_input.h"
#include "egs_advanced_application.h"
#include "egs_transformations.h"
#include "pubsub.h"
#include "iaea_phsp.h"
 
class PHSPControl : public Subscriber {
 
    enum ACCESS {
        READ = 1, 
        WRITE = 2, 
        APPEND = 3 
    };
 
    enum PARTICLE_TYPE {
        ALL_TYPES = -1,
        PHOTON = 1,
        ELECTRON = 2,
        POSITRON = 3,
        NEUTRON = 4,
        PROTON = 5
    };
 
    IAEA_I32 mode;  
    string fname;   
    IAEA_I32 id;    
 
    IAEA_I64 num_written; 
 
 
    EGS_AffineTransform *transform; 
    EGS_Float boundary_step; 
    bool kill_after_scoring; 
 
    bool print_header; 
 
    void initSource();
 
    void writeParticle(EGS_Particle *p);
 
 
    short getIAEAParticleType(const EGS_Particle *p);
 
 
 
public:
 
    PHSPControl(EGS_Input *inp, EGS_AffineTransform *trans, EGS_AdvancedApplication *app, Publisher *pub);
 
    void finish(EGS_I64 n_orig_particles);
 
    void destroySource();
 
 
    void outputResults();
 
    void update(EB_Message message, void *particle);
 
};