Material.h

Go to the documentation of this file.

//
// Class Material
//   Base class for representing materials
//
// Copyright (c) 2019 - 2021, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#ifndef MATERIAL_H
#define MATERIAL_H
 
#include <array>
#include <cmath>
#include <map>
#include <memory>
#include <string>
 
namespace Physics {
    class Material {
    public:
        enum FitCoeffs { A1 = 0, A2, A3, A4, A5, B1, B2, B3, B4, B5 };
 
        Material(
                double atomicNumber, double atomicMass, double massDensity, double radiationLength,
                double meanExcitationEnergy, std::array<double, 10> fitCoefficients)
            : atomicNumber_m(atomicNumber),
              atomicMass_m(atomicMass),
              massDensity_m(massDensity),
              radiationLength_m(radiationLength / massDensity / 100),
              meanExcitationEnergy_m(meanExcitationEnergy),
              stoppingPowerFitCoefficients_m(fitCoefficients) {}
 
        double getAtomicNumber() const;          // [1]
        double getAtomicMass() const;            // [u]
        double getMassDensity() const;           // [g cm^-3]
        double getRadiationLength() const;       // [m]
        double getMeanExcitationEnergy() const;  // [eV]
        double getStoppingPowerFitCoefficients(FitCoeffs n) const;
 
        static std::shared_ptr<Material> getMaterial(const std::string& name);
        static std::shared_ptr<Material> addMaterial(
                const std::string& name, std::shared_ptr<Material> mat_ptr);
 
    private:
        static std::map<std::string, std::shared_ptr<Material> > protoTable_sm;
 
        const double atomicNumber_m;
        const double atomicMass_m;
        const double massDensity_m;
        const double radiationLength_m;
        const double meanExcitationEnergy_m;
        const std::array<double, 10> stoppingPowerFitCoefficients_m;
    };
 
    inline double Material::getAtomicNumber() const { return atomicNumber_m; }
 
    inline double Material::getAtomicMass() const { return atomicMass_m; }
 
    inline double Material::getMassDensity() const { return massDensity_m; }
 
    inline double Material::getRadiationLength() const { return radiationLength_m; }
 
    inline double Material::getMeanExcitationEnergy() const { return meanExcitationEnergy_m; }
 
    inline double Material::getStoppingPowerFitCoefficients(Material::FitCoeffs n) const {
        return stoppingPowerFitCoefficients_m[n];
    }
 
    // inline
    // double Material::getMeanExcitationEnergy() const {
    //     double Z = getAtomicNumber();
    //     if (Z < 13.0)
    //         return 12 * Z + 7;
 
    //     return 9.76 * Z + 58.8 * std::pow(Z, -.19);
    // }
}  // namespace Physics
#endif