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 <map>
#include <array>
#include <string>
#include <memory>
#include <cmath>

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);
    // }
}
#endif