tally_data.h

#ifndef MCXRAYTRANSPORT_TALLY_DATA_H
#define MCXRAYTRANSPORT_TALLY_DATA_H
 
#include "photon.h"
#include "voxel.h"
 
struct TempVoxelData {
    Voxel& voxel;
    double energy_deposited = 0.0;
 
    explicit TempVoxelData(Voxel& voxel) : voxel(voxel) {}
};
 
struct TempSurfaceTallyData {
    double free_path = 0.0;
    bool isInteract = false;
    double entrance_cosine = 0.0;
    bool already_counted = false;
    Photon initial_photon;
};
 
struct TempVolumeTallyData {
    double energy_deposited = 0.0;
    double free_path = 0.0;
    bool isInteract = false;
    bool already_counted = false;
    Photon initial_photon;
    Photon final_photon;
};
 
struct VolumeTallyData {
    Eigen::Vector<double, Eigen::Dynamic> energy_deposition;
    Eigen::Vector<double, Eigen::Dynamic> incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> primary_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> secondary_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> single_coherent_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> single_incoherent_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> multiple_scatter_incident_energy;
    int number_of_particles = 0;
    int number_of_primary_particles = 0;
    int number_of_secondary_particles = 0;
    Eigen::Vector<double, Eigen::Dynamic> entrance_cosines;
    Eigen::Vector<double, Eigen::Dynamic> primary_entrance_cosines;
    int number_of_interactions = 0;
    double area = 0.0;
};
 
struct SurfaceTallyData {
    Eigen::Vector<double, Eigen::Dynamic> energy_deposition;
    Eigen::Vector<double, Eigen::Dynamic> incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> primary_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> secondary_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> single_coherent_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> single_incoherent_incident_energy;
    Eigen::Vector<double, Eigen::Dynamic> multiple_scatter_incident_energy;
    int number_of_particles = 0;
    int number_of_primary_particles = 0;
    int number_of_secondary_particles = 0;
    Eigen::Vector<double, Eigen::Dynamic> entrance_cosines;
    Eigen::Vector<double, Eigen::Dynamic> primary_entrance_cosines;
    int number_of_interactions = 0;
    double area = 0.0;
};
 
#endif //MCXRAYTRANSPORT_TALLY_DATA_H