GGEMSRangeCuts.hh

Go to the documentation of this file.

#ifndef GUARD_GGEMS_PHYSICS_GGEMSRANGECUTS_HH
#define GUARD_GGEMS_PHYSICS_GGEMSRANGECUTS_HH
 
// ************************************************************************
// * This file is part of GGEMS.                                          *
// *                                                                      *
// * GGEMS is free software: you can redistribute it and/or modify        *
// * it under the terms of the GNU General Public License as published by *
// * the Free Software Foundation, either version 3 of the License, or    *
// * (at your option) any later version.                                  *
// *                                                                      *
// * GGEMS 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 General Public License for more details.                         *
// *                                                                      *
// * You should have received a copy of the GNU General Public License    *
// * along with GGEMS.  If not, see <https://www.gnu.org/licenses/>.      *
// *                                                                      *
// ************************************************************************
 
#ifdef _MSC_VER
#pragma warning(disable: 4251) // Deleting warning exporting STL members!!!
#endif
 
#include "GGEMS/materials/GGEMSMaterials.hh"
 
class GGEMSMaterials;
class GGEMSLogEnergyTable;
 
typedef std::unordered_map<std::string, GGfloat> EnergyCutUMap; 
class GGEMS_EXPORT GGEMSRangeCuts
{
  public:
    GGEMSRangeCuts(void);
 
    ~GGEMSRangeCuts(void);
 
    GGEMSRangeCuts(GGEMSRangeCuts const& range_cuts) = delete;
 
    GGEMSRangeCuts& operator=(GGEMSRangeCuts const& range_cuts) = delete;
 
    GGEMSRangeCuts(GGEMSRangeCuts const&& range_cuts) = delete;
 
    GGEMSRangeCuts& operator=(GGEMSRangeCuts const&& range_cuts) = delete;
 
    void SetPhotonDistanceCut(GGfloat const& cut);
 
    void SetElectronDistanceCut(GGfloat const& cut);
 
    void SetPositronDistanceCut(GGfloat const& cut);
 
    inline GGfloat GetPhotonDistanceCut(void) const {return distance_cut_photon_;}
 
    inline EnergyCutUMap GetPhotonEnergyCut(void) const {return energy_cuts_photon_;}
 
    inline GGfloat GetElectronDistanceCut(void) const {return distance_cut_electron_;}
 
    inline EnergyCutUMap GetElectronEnergyCut(void) const {return energy_cuts_electron_;}
 
    inline GGfloat GetPositronDistanceCut(void) const {return distance_cut_positron_;}
 
    inline EnergyCutUMap GetPositronEnergyCut(void) const {return energy_cuts_positron_;}
 
    void ConvertCutsFromDistanceToEnergy(GGEMSMaterials* materials);
 
  private:
    GGfloat ConvertToEnergy(GGEMSMaterialTables* material_table, GGushort const& index_mat, std::string const& particle_name);
 
    void BuildElementsLossTable(GGEMSMaterialTables* material_table, GGushort const& index_mat, std::string const& particle_name);
 
    void BuildAbsorptionLengthTable(GGEMSMaterialTables* material_table, GGushort const& index_mat);
 
    void BuildMaterialLossTable(GGEMSMaterialTables* material_table, GGushort const& index_mat);
 
    GGfloat ComputePhotonCrossSection(GGuchar const& atomic_number, GGfloat const& energy) const;
 
    GGfloat ComputeLossElectron(GGuchar const& atomic_number, GGfloat const& energy) const;
 
    GGfloat ComputeLossPositron(GGuchar const& atomic_number, GGfloat const& energy) const;
 
    GGfloat ConvertLengthToEnergyCut(GGfloat const& length_cut) const;
 
  private:
    GGfloat min_energy_; 
    GGfloat max_energy_; 
    GGsize number_of_bins_; 
    // Photon
    GGfloat distance_cut_photon_; 
    EnergyCutUMap energy_cuts_photon_; 
    // Electron
    GGfloat distance_cut_electron_; 
    EnergyCutUMap energy_cuts_electron_; 
    // Positron
    GGfloat distance_cut_positron_; 
    EnergyCutUMap energy_cuts_positron_; 
    GGEMSLogEnergyTable* range_table_material_; 
    GGEMSLogEnergyTable** loss_table_dedx_table_elements_; 
    GGsize number_of_tables_; 
};
 
#endif // GUARD_GGEMS_PHYSICS_GGEMSRANGECUTS_HH