BDSMaterials.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2023.
 
This file is part of BDSIM.
 
BDSIM 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 version 3 of the License.
 
BDSIM 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 BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSMATERIALS_H
#define BDSMATERIALS_H
 
#include <list>
#include <map>
#include <set>
#include <vector>
 
#include "globals.hh"
#include "G4Material.hh" // for G4State
 
class G4Element;
class G4MaterialPropertiesTable;
 
class BDSMaterials
{
public:
 
  static BDSMaterials* Instance();
  ~BDSMaterials(); 
 
  void PrepareRequiredMaterials(G4bool verbose = false);
 
  void ListMaterials()const;
 
  G4Material* GetMaterial(G4String material) const;
  G4Element*  GetElement(const G4String& symbol) const;
 
  void CacheMaterialsFromGDML(const std::map<G4String, G4Material*>& materialsGDML);
 
  void CheckForConflictingMaterialsAfterLoad(const G4String& geometryFileName,
                                             const G4String& componentName) const;
 
protected:
  BDSMaterials();
  void AddMaterial(G4Material* material, G4String name);
  void AddExistingMaterialAlias(const G4String &existingMaterialName,
                                G4String alias);
 
  void AddMaterial(G4String name,
                   G4double Z,
                   G4double A,
                   G4double density,
                   G4State  state, 
                   G4double temperature, 
                   G4double pressure);
 
  template <typename Type>
  void AddMaterial(G4String name, 
                   G4double density, 
                   G4State  state, 
                   G4double temperature, 
                   G4double pressure,
                   const std::list<G4String>& components,
                   const std::list<Type>&     componentFractions);
 
  G4Element* CheckElement(const G4String& symbol) const;
 
private:
  static BDSMaterials* instance;
 
  void DefineMetals();
  void DefineSuperconductors();
  void DefineNonMetalSolids();
  void DefineScintillators();
  void DefineLHCComponents();
  void DefineLiquids();
  void DefineGases();
  void DefinePlasmas();
  void DefineVacuums();
 
  void AddElement(G4Element* element, const G4String& symbol);
  void AddElement(const G4String& name, const G4String& symbol, G4double itsZ, G4double itsA);
 
  void DensityCheck(G4double  density,
                    const G4String& materialName) const;
 
  void PrintBasicMaterialMassFraction(G4Material* material) const;
 
  std::set<G4String> materialNames;
  std::map<G4String, G4Material*> materials;
  std::set<G4String> aliasNames;
  std::map<G4String, G4Material*> aliases;
  std::map<G4Material*, std::vector<G4String>> materialToAliases;
  std::set<G4String> externalMaterialNames;
  std::map<G4String, G4Material*> externalMaterials;
  std::map<G4String,G4Element*>  elements;
  std::vector<G4MaterialPropertiesTable*> propertiesTables;
  G4MaterialPropertiesTable* CreatePropertiesTable();
};
 
#endif