MatterDef.h

Go to the documentation of this file.

#ifndef MATTER_DEF_H
#define MATTER_DEF_H
 
#include "wcpplib/matter/AtomDef.h"
 
namespace Heed {
 
/// Definition of matter (material or any media).
/// Only the basic information: the name, the notation,
/// the atomic mixture, temperature, density, effective ionization potential.
///
/// 1998-2004 I. Smirnov
 
class MatterDef : public AtomMixDef {
  std::string nameh = "none";
  std::string notationh = "none";
  double temperatureh;
  double densityh;
  double I_effh;
  // I_eff is still not very reliable.
  // There are too many approximations for this.
  // Here is a simplest and probably not good.
  void calc_I_eff();
 
 public:
  MatterDef() = default;
  MatterDef(const std::string& fname, const std::string& fnotation, long fqatom,
            const std::vector<std::string>& fatom_not,
            const std::vector<double>& fweight_quan, double fdensity,
            double ftemperature);
  MatterDef(const std::string& fname, const std::string& fnotation,
            const std::string& fatom_not, double fdensity, double ftemperature);
  MatterDef(const std::string& fname, const std::string& fnotation,
            const std::string& fatom_not1, double fweight_quan1,
            const std::string& fatom_not2, double fweight_quan2,
            double fdensity, double ftemperature);
  MatterDef(const std::string& fname, const std::string& fnotation,
            const std::string& fatom_not1, double fweight_quan1,
            const std::string& fatom_not2, double fweight_quan2,
            const std::string& fatom_not3, double fweight_quan3,
            double fdensity, double ftemperature);
  virtual ~MatterDef() = default;
 
  const std::string& name() const { return nameh; }
  const std::string& notation() const { return notationh; }
  double density() const { return densityh; }
  double temperature() const { return temperatureh; }
  double I_eff() const { return I_effh; }
 
  void print(std::ostream& file, int l) const;
  MatterDef* copy() const { return new MatterDef(*this); }
};
std::ostream& operator<<(std::ostream& file, const MatterDef& f);
 
}
#endif