gmaterial.h

Go to the documentation of this file.

#pragma once
 
// gemc
#include "gbase.h"
 
// c++
#include <string>
#include <vector>
#include <numeric>
 
class GMaterial : public GBase<GMaterial>
{
public:
    GMaterial(const std::string& system, std::vector<std::string> pars, const std::shared_ptr<GLogger>& logger);
 
    [[nodiscard]] virtual std::unique_ptr<GMaterial> clone() const {
        return std::make_unique<GMaterial>(*this);
    }
 
    virtual ~GMaterial() = default;
 
private:
    std::string system;      
    std::string name;        
    std::string description; 
 
    double                   density{};    
    std::vector<std::string> components;   
    std::vector<double>      amounts;      
 
    std::vector<double> photonEnergy;       
    std::vector<double> indexOfRefraction;  
    std::vector<double> absorptionLength;   
    std::vector<double> reflectivity;       
    std::vector<double> efficiency;         
 
    std::vector<double> fastcomponent;      
    std::vector<double> slowcomponent;      
 
    double scintillationyield{};            
    double resolutionscale{};               
    double fasttimeconstant{};              
    double slowtimeconstant{};              
    double yieldratio{};                    
    double birksConstant{};                 
 
    std::vector<double> rayleigh;           
 
    void setComponentsFromString(const std::string& composition);
 
    void getMaterialPropertyFromString(const std::string& parameter, const std::string& propertyName);
 
    friend std::ostream& operator<<(std::ostream& stream, const GMaterial&);
 
    bool assign_if_set(const std::vector<std::string>& pars, size_t& i, double& out);
 
public:
    [[nodiscard]] std::string              getName() const { return name; }
    [[nodiscard]] std::string              getDescription() const { return description; }
    [[nodiscard]] double                   getDensity() const { return density; }
    [[nodiscard]] std::vector<std::string> getComponents() const { return components; }
    [[nodiscard]] std::vector<double>      getAmounts() const { return amounts; }
 
    [[nodiscard]] std::vector<double> getPhotonEnergy() const { return photonEnergy; }
    [[nodiscard]] std::vector<double> getIndexOfRefraction() const { return indexOfRefraction; }
    [[nodiscard]] std::vector<double> getAbsorptionLength() const { return absorptionLength; }
    [[nodiscard]] std::vector<double> getReflectivity() const { return reflectivity; }
    [[nodiscard]] std::vector<double> getEfficiency() const { return efficiency; }
 
    [[nodiscard]] std::vector<double> getFastcomponent() const { return fastcomponent; }
    [[nodiscard]] std::vector<double> getSlowcomponent() const { return slowcomponent; }
    [[nodiscard]] double              getScintillationyield() const { return scintillationyield; }
    [[nodiscard]] double              getResolutionscale() const { return resolutionscale; }
    [[nodiscard]] double              getFasttimeconstant() const { return fasttimeconstant; }
    [[nodiscard]] double              getSlowtimeconstant() const { return slowtimeconstant; }
    [[nodiscard]] double              getYieldratio() const { return yieldratio; }
    [[nodiscard]] double              getBirksConstant() const { return birksConstant; }
 
    [[nodiscard]] std::vector<double> getRayleigh() const { return rayleigh; }
 
    [[nodiscard]] bool isChemicalFormula() const {
        return std::accumulate(amounts.begin(), amounts.end(), 0.0) > 1;
    }
};