gar::garg4::MaterialPropertyLoader Class Reference

#include <MaterialPropertyLoader.h>

Public Member Functions
	MaterialPropertyLoader ()
	~MaterialPropertyLoader ()
std::map< double, double >	GetMaterialProperty (std::string const &Material, std::string const &Property)
double	GetMaterialConstProperty (std::string const &Material, std::string const &Property)
std::map< std::string, double >	GetMaterialConstProperties (std::string const &Material)
std::map< std::string, std::map< double, double > >	GetMaterialProperties (std::string const &Material)
void	SetMaterialProperty (std::string const &Material, std::string const &Property, std::map< double, double > const &Values, double Unit)
void	SetMaterialConstProperty (std::string const &Material, std::string const &Property, double Value, double Unit)
void	SetBirksConstant (std::string const &, double, double)
void	SetReflectances (std::string const &, std::map< std::string, std::map< double, double > > const &, std::map< std::string, std::map< double, double > > const &)
void	GetPropertiesFromServices ()
void	UpdateGeometry (G4LogicalVolumeStore *)

Private Attributes
std::map< std::string, std::map< std::string, double > >	fConstPropertyList
std::map< std::string, std::map< std::string, std::map< double, double > > >	fPropertyList
std::map< std::string, double >	fBirksConstants

Detailed Description

Definition at line 20 of file MaterialPropertyLoader.h.

Constructor & Destructor Documentation

gar::garg4::MaterialPropertyLoader::MaterialPropertyLoader ( )

inline

Definition at line 24 of file MaterialPropertyLoader.h.

{}

gar::garg4::MaterialPropertyLoader::~MaterialPropertyLoader ( )

inline

Definition at line 25 of file MaterialPropertyLoader.h.

{}

Member Function Documentation

std::map<std::string,double> gar::garg4::MaterialPropertyLoader::GetMaterialConstProperties ( std::string const &  Material )

inline

Definition at line 40 of file MaterialPropertyLoader.h.

      {return fConstPropertyList[Material];}

double gar::garg4::MaterialPropertyLoader::GetMaterialConstProperty ( std::string const &  Material, std::string const &  Property  )

inline

Definition at line 36 of file MaterialPropertyLoader.h.

      {return fConstPropertyList[Material][Property];}

std::map<std::string,std::map<double,double> > gar::garg4::MaterialPropertyLoader::GetMaterialProperties ( std::string const &  Material )

inline

Definition at line 43 of file MaterialPropertyLoader.h.

      {return fPropertyList[Material];}

std::map<double,double> gar::garg4::MaterialPropertyLoader::GetMaterialProperty ( std::string const &  Material, std::string const &  Property  )

inline

Definition at line 32 of file MaterialPropertyLoader.h.

      {return fPropertyList[Material][Property];}

void gar::garg4::MaterialPropertyLoader::GetPropertiesFromServices

(

)

Definition at line 193 of file MaterialPropertyLoader.cxx.

    {
//      const detinfo::GArProperties*      GarProp = gar::providerFrom<detinfo::GArPropertiesService>();
//      const detinfo::DetectorProperties* DetProp = gar::providerFrom<detinfo::DetectorPropertiesService>();

        const detinfo::ECALProperties*     EcalProp = gar::providerFrom<detinfo::ECALPropertiesService>();

        SetMaterialConstProperty("Scintillator", "DENSITY", 1.06, CLHEP::g/CLHEP::cm3);
        SetBirksConstant("Scintillator", EcalProp->ScintBirksConstant(), CLHEP::mm/CLHEP::MeV);
    }

void gar::garg4::MaterialPropertyLoader::SetBirksConstant	(	std::string const &	Material,
		double	PropertyValue,
		double	Unit
	)

Definition at line 64 of file MaterialPropertyLoader.cxx.

    {
      fBirksConstants[Material] = PropertyValue * Unit;
      MF_LOG_INFO("MaterialPropertyLoader")
      << "Set Birks constant "
      << Material;
    }

void gar::garg4::MaterialPropertyLoader::SetMaterialConstProperty	(	std::string const &	Material,
		std::string const &	Property,
		double	Value,
		double	Unit
	)

Definition at line 47 of file MaterialPropertyLoader.cxx.

    {
      fConstPropertyList[Material][Property]=PropertyValue*Unit;

      MF_LOG_INFO("MaterialPropertyLoader")
      << "Added const property "
      << Material
      << "  "
      << Property
      << " = "
      << PropertyValue;
    }

void gar::garg4::MaterialPropertyLoader::SetMaterialProperty	(	std::string const &	Material,
		std::string const &	Property,
		std::map< double, double > const &	Values,
		double	Unit
	)

Definition at line 28 of file MaterialPropertyLoader.cxx.

    {
      std::map<double,double> PropVectorWithUnit;
      for(auto const& itr : PropertyVector){
        PropVectorWithUnit[itr.first * CLHEP::eV] = itr.second * Unit;
      }

      fPropertyList[Material][Property] = PropVectorWithUnit;

      MF_LOG_INFO("MaterialPropertyLoader")
      <<"Added property "
      << Material<< "  "
                  << Property;
    }

void gar::garg4::MaterialPropertyLoader::SetReflectances	(	std::string const &	,
		std::map< std::string, std::map< double, double > > const &	Reflectances,
		std::map< std::string, std::map< double, double > > const &	DiffuseFractions
	)

Definition at line 162 of file MaterialPropertyLoader.cxx.

    {
      std::map<double, double> ReflectanceToStore;
      std::map<double, double> DiffuseToStore;

      for(auto const& itMat : Reflectances){
        std::string ReflectancePropName = std::string("REFLECTANCE_") + itMat.first;
        ReflectanceToStore.clear();
        for(auto const& itEn : itMat.second){
          ReflectanceToStore[itEn.first] = itEn.second;
        }
        SetMaterialProperty("GAr", ReflectancePropName, ReflectanceToStore,1);
      }

      for(auto const& itMat : DiffuseFractions){
        std::string DiffusePropName = std::string("DIFFUSE_REFLECTANCE_FRACTION_") + itMat.first;
        DiffuseToStore.clear();
        for(auto const& itEn : itMat.second){
          DiffuseToStore[itEn.first] = itEn.second;
        }
        SetMaterialProperty("GAr",
                            DiffusePropName,
                            DiffuseToStore,
                            1);
      }

    }

void gar::garg4::MaterialPropertyLoader::UpdateGeometry

(

G4LogicalVolumeStore *

lvs

)

Definition at line 75 of file MaterialPropertyLoader.cxx.

    {
      std::map<std::string,G4MaterialPropertiesTable*> MaterialTables;
      std::map<std::string,bool> MaterialsSet;

      MF_LOG_INFO("MaterialPropertyLoader")
      << "UPDATING GEOMETRY";

      // Loop over each material with a property vector and create a new material table for it
      for(auto const& itr : fPropertyList) {
        std::string Material     = itr.first;
        MaterialsSet[Material]   = true;
        MaterialTables[Material] = new G4MaterialPropertiesTable;
      }

      // Loop over each material with a const property,
      // if material table does not exist, create one
      for(auto const& itr : fConstPropertyList){
        std::string Material = itr.first;
        if(!MaterialsSet[Material]){
          MaterialsSet[Material]   = true;
          MaterialTables[Material] = new G4MaterialPropertiesTable;
        }
      }

      // For each property vector, convert to an array of g4doubles and
      // feed to materials table Lots of firsts and seconds!  See annotation
      // in MaterialPropertyLoader.h to follow what each element is

      for(auto const& itr : fPropertyList){
        std::string Material = itr.first;
        for(auto const& jitr : itr.second){
          std::string Property = jitr.first;
          std::vector<G4double> g4MomentumVector;
          std::vector<G4double> g4PropertyVector;

          for(auto const& kitr : jitr.second){
            g4MomentumVector.push_back(kitr.first);
            g4PropertyVector.push_back(kitr.second);
          }
          int NoOfElements = g4MomentumVector.size();
          MaterialTables[Material]->AddProperty(Property.c_str(),
                                                &g4MomentumVector[0],
                                                &g4PropertyVector[0],
                                                NoOfElements);
          MF_LOG_INFO("MaterialPropertyLoader")
          << "Added property "
          << Property
          << " to material table "
          << Material;
        }
      }

      //Add each const property element
      for(auto const& itr : fConstPropertyList){
        std::string Material = itr.first;
        for(auto const& jitr : itr.second){
          std::string Property   = jitr.first;
          G4double PropertyValue = jitr.second;
          MaterialTables[Material]->AddConstProperty(Property.c_str(), PropertyValue);

          MF_LOG_INFO("MaterialPropertyLoader")
          << "Added const property "
          << Property
          << " to material table "
          << Material;
        }
      }

      //Loop through geometry elements and apply relevant material table where materials match
      for ( G4LogicalVolumeStore::iterator i = lvs->begin(); i != lvs->end(); ++i ){
        G4LogicalVolume* volume = (*i);
        G4Material* TheMaterial = volume->GetMaterial();
        std::string Material = TheMaterial->GetName();
        for(auto const& jitr : MaterialTables){
          if(Material == jitr.first){
            TheMaterial->SetMaterialPropertiesTable(jitr.second);
            //Birks Constant, for some reason, must be set separately
            if(fBirksConstants[Material]!=0)
              TheMaterial->GetIonisation()->SetBirksConstant(fBirksConstants[Material]);
            volume->SetMaterial(TheMaterial);
          }
        }
      }
    }

Member Data Documentation

std::map<std::string, double> gar::garg4::MaterialPropertyLoader::fBirksConstants

private

Definition at line 82 of file MaterialPropertyLoader.h.

std::map< std::string , std::map < std::string,double> > gar::garg4::MaterialPropertyLoader::fConstPropertyList

private

Definition at line 77 of file MaterialPropertyLoader.h.

std::map< std::string , std::map < std::string , std::map < double , double > > > gar::garg4::MaterialPropertyLoader::fPropertyList

private

Definition at line 80 of file MaterialPropertyLoader.h.

---

The documentation for this class was generated from the following files:

• garsoft/GArG4/MaterialPropertyLoader.h
• garsoft/GArG4/MaterialPropertyLoader.cxx