G4Material Class Reference

#include <G4Material.hh>

Collaboration diagram for G4Material:

Public Member Functions
	G4Material (const G4String &name, G4double z, G4double a, G4double density, G4State state=kStateUndefined, G4double temp=NTP_Temperature, G4double pressure=CLHEP::STP_Pressure)
	G4Material (const G4String &name, G4double density, G4int nComponents, G4State state=kStateUndefined, G4double temp=NTP_Temperature, G4double pressure=CLHEP::STP_Pressure)
	G4Material (const G4String &name, G4double density, const G4Material *baseMaterial, G4State state=kStateUndefined, G4double temp=NTP_Temperature, G4double pressure=CLHEP::STP_Pressure)
void	AddElement (G4Element *element, G4int nAtoms)
void	AddElement (G4Element *element, G4double fraction)
void	AddMaterial (G4Material *material, G4double fraction)
virtual	~G4Material ()
void	SetChemicalFormula (const G4String &chF)
const G4String &	GetName () const
const G4String &	GetChemicalFormula () const
G4double	GetDensity () const
G4State	GetState () const
G4double	GetTemperature () const
G4double	GetPressure () const
size_t	GetNumberOfElements () const
const G4ElementVector *	GetElementVector () const
const G4double *	GetFractionVector () const
const G4int *	GetAtomsVector () const
const G4Element *	GetElement (G4int iel) const
const G4double *	GetVecNbOfAtomsPerVolume () const
G4double	GetTotNbOfAtomsPerVolume () const
G4double	GetTotNbOfElectPerVolume () const
const G4double *	GetAtomicNumDensityVector () const
G4double	GetElectronDensity () const
G4double	GetRadlen () const
G4double	GetNuclearInterLength () const
G4IonisParamMat *	GetIonisation () const
G4SandiaTable *	GetSandiaTable () const
const G4Material *	GetBaseMaterial () const
const std::map< G4Material *, G4double > &	GetMatComponents () const
G4double	GetMassOfMolecule () const
G4double	GetZ () const
G4double	GetA () const
void	SetMaterialPropertiesTable (G4MaterialPropertiesTable *anMPT)
G4MaterialPropertiesTable *	GetMaterialPropertiesTable () const
size_t	GetIndex () const
G4int	operator== (const G4Material &) const
G4int	operator!= (const G4Material &) const
	G4Material (__void__ &)
void	SetName (const G4String &name)

Static Public Member Functions
static G4MaterialTable *	GetMaterialTable ()
static size_t	GetNumberOfMaterials ()
static G4Material *	GetMaterial (const G4String &name, G4bool warning=true)

Private Member Functions
	G4Material (const G4Material &)
const G4Material &	operator= (const G4Material &)
void	InitializePointers ()
void	ComputeDerivedQuantities ()
void	ComputeRadiationLength ()
void	ComputeNuclearInterLength ()
void	CopyPointersOfBaseMaterial ()

Private Attributes
G4String	fName
G4String	fChemicalFormula
G4double	fDensity
G4State	fState
G4double	fTemp
G4double	fPressure
G4int	maxNbComponents
G4int	fArrayLength
G4int	fNumberOfComponents
G4int	fNumberOfElements
G4ElementVector *	theElementVector
G4double *	fMassFractionVector
G4int *	fAtomsVector
G4MaterialPropertiesTable *	fMaterialPropertiesTable
size_t	fIndexInTable
G4double *	VecNbOfAtomsPerVolume
G4double	TotNbOfAtomsPerVolume
G4double	TotNbOfElectPerVolume
G4double	fRadlen
G4double	fNuclInterLen
G4IonisParamMat *	fIonisation
G4SandiaTable *	fSandiaTable
const G4Material *	fBaseMaterial
G4double	fMassOfMolecule
std::map< G4Material *, G4double >	fMatComponents

Static Private Attributes
static G4MaterialTable	theMaterialTable

Friends
std::ostream &	operator<< (std::ostream &, const G4Material *)
std::ostream &	operator<< (std::ostream &, const G4Material &)
std::ostream &	operator<< (std::ostream &, G4MaterialTable)

Detailed Description

Definition at line 120 of file G4Material.hh.

Constructor & Destructor Documentation

G4Material() [1/5]

G4Material::G4Material	(	const G4String &	name,
		G4double	z,
		G4double	a,
		G4double	density,
		G4State	state = kStateUndefined,
		G4double	temp = NTP_Temperature,
		G4double	pressure = CLHEP::STP_Pressure
	)

Definition at line 90 of file G4Material.cc.

  : fName(name)            
{
  InitializePointers();
    
  if (density < universe_mean_density)
    { 
      G4cout << " G4Material WARNING:"
         << " define a material with density=0 is not allowed. \n"
         << " The material " << name << " will be constructed with the"
         << " default minimal density: " << universe_mean_density/(g/cm3) 
         << "g/cm3" << G4endl;
      density = universe_mean_density;
    } 

  fDensity  = density;
  fState    = state;
  fTemp     = temp;
  fPressure = pressure;

  // Initialize theElementVector allocating one
  // element corresponding to this material
  maxNbComponents        = fNumberOfComponents = fNumberOfElements = 1;
  fArrayLength           = maxNbComponents;
  theElementVector       = new G4ElementVector();

  const std::vector<G4String> elmnames = 
    G4NistManager::Instance()->GetNistElementNames();
  G4String enam, snam;
  G4int iz = G4lrint(z);
  if(iz < (G4int)elmnames.size()) { 
    snam = elmnames[iz]; 
    enam = snam; 
  } else { 
    enam = "ELM_" + name; 
    snam = name;
  }
  theElementVector->push_back(new G4Element(enam, snam, z, a));  

  fMassFractionVector    = new G4double[1];
  fMassFractionVector[0] = 1. ;
  fMassOfMolecule        = a/Avogadro;
  
  if (fState == kStateUndefined)
    {
      if (fDensity > kGasThreshold) { fState = kStateSolid; }
      else                          { fState = kStateGas; }
    }

  ComputeDerivedQuantities();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4Material() [2/5]

G4Material::G4Material	(	const G4String &	name,
		G4double	density,
		G4int	nComponents,
		G4State	state = kStateUndefined,
		G4double	temp = NTP_Temperature,
		G4double	pressure = CLHEP::STP_Pressure
	)

Definition at line 149 of file G4Material.cc.

  : fName(name)            
{
  InitializePointers();
    
  if (density < universe_mean_density)
    {
      G4cout << "--- Warning from G4Material::G4Material()"
         << " define a material with density=0 is not allowed. \n"
         << " The material " << name << " will be constructed with the"
         << " default minimal density: " << universe_mean_density/(g/cm3) 
         << "g/cm3" << G4endl;
      density = universe_mean_density;
    }
        
  fDensity  = density;
  fState    = state;
  fTemp     = temp;
  fPressure = pressure;
    
  maxNbComponents     = nComponents;
  fArrayLength        = maxNbComponents;
  fNumberOfComponents = fNumberOfElements = 0;
  theElementVector    = new G4ElementVector();
  theElementVector->reserve(maxNbComponents);  
    
  if (fState == kStateUndefined) 
    {
      if (fDensity > kGasThreshold) { fState = kStateSolid; }
      else                          { fState = kStateGas; }
    }
}

Here is the call graph for this function:

G4Material() [3/5]

G4Material::G4Material	(	const G4String &	name,
		G4double	density,
		const G4Material *	baseMaterial,
		G4State	state = kStateUndefined,
		G4double	temp = NTP_Temperature,
		G4double	pressure = CLHEP::STP_Pressure
	)

Definition at line 188 of file G4Material.cc.

  : fName(name)            
{
  InitializePointers();
    
  if (density < universe_mean_density)
    {
      G4cout << "--- Warning from G4Material::G4Material()"
         << " define a material with density=0 is not allowed. \n"
         << " The material " << name << " will be constructed with the"
         << " default minimal density: " << universe_mean_density/(g/cm3) 
         << "g/cm3" << G4endl;
      density = universe_mean_density;
    }

  fDensity  = density;
  fState    = state;
  fTemp     = temp;
  fPressure = pressure;

  fBaseMaterial = bmat;
  fChemicalFormula = fBaseMaterial->GetChemicalFormula();
  fMassOfMolecule  = fBaseMaterial->GetMassOfMolecule();

  fNumberOfElements = fBaseMaterial->GetNumberOfElements();     
  maxNbComponents = fNumberOfElements;
  fNumberOfComponents = fNumberOfElements;

  fMaterialPropertiesTable = fBaseMaterial->GetMaterialPropertiesTable();

  CopyPointersOfBaseMaterial();
}

Here is the call graph for this function:

~G4Material()

G4Material::~G4Material ( )

virtual

Definition at line 241 of file G4Material.cc.

{
  //  G4cout << "### Destruction of material " << fName << " started" <<G4endl;
  if(!fBaseMaterial) {
    if (theElementVector)       { delete    theElementVector; }
    if (fMassFractionVector)    { delete [] fMassFractionVector; }
    if (fAtomsVector)           { delete [] fAtomsVector; }
    if (fSandiaTable)           { delete fSandiaTable; }
  }
  if (fIonisation)            { delete fIonisation; }
  if (VecNbOfAtomsPerVolume)  { delete [] VecNbOfAtomsPerVolume; }

  // Remove this material from theMaterialTable.
  //
  theMaterialTable[fIndexInTable] = 0;
}

G4Material() [4/5]

G4Material::G4Material

(

__void__ &

)

Definition at line 228 of file G4Material.cc.

  : fChemicalFormula(""), fDensity(0.0), fState(kStateUndefined), fTemp(0.0), 
    fPressure(0.0), maxNbComponents(0), fArrayLength(0), 
    fNumberOfComponents(0), fNumberOfElements(0), theElementVector(0), 
    fMassFractionVector(0), fAtomsVector(0), fMaterialPropertiesTable(0), 
    fIndexInTable(0), VecNbOfAtomsPerVolume(0), TotNbOfAtomsPerVolume(0), 
    TotNbOfElectPerVolume(0), fRadlen(0.0), fNuclInterLen(0.0), fIonisation(0), 
    fSandiaTable(0), fBaseMaterial(0), fMassOfMolecule(0.0)
{
}

G4Material() [5/5]

G4Material::G4Material ( const G4Material &  )

private

Member Function Documentation

AddElement() [1/2]

void G4Material::AddElement	(	G4Element *	element,
		G4int	nAtoms
	)

Definition at line 364 of file G4Material.cc.

{   
  // initialization
  if ( fNumberOfElements == 0 ) {
     fAtomsVector        = new G4int   [fArrayLength];
     fMassFractionVector = new G4double[fArrayLength];
  }

  // filling ...
  if ( fNumberOfElements < maxNbComponents ) {
     theElementVector->push_back(element);     
     fAtomsVector[fNumberOfElements] = nAtoms;
     fNumberOfComponents = ++fNumberOfElements;
  } else {
    G4cout << "G4Material::AddElement ERROR for " << fName << " nElement= " 
       <<  fNumberOfElements << G4endl;
    G4Exception ("G4Material::AddElement()", "mat031", FatalException, 
           "Attempt to add more than the declared number of elements.");
  } 
  // filled.
  if ( fNumberOfElements == maxNbComponents ) {     
    // compute proportion by mass
    G4int i=0;
    G4double Amol = 0.;
    for (i=0; i<fNumberOfElements; ++i) {
      G4double w = fAtomsVector[i]*(*theElementVector)[i]->GetA(); 
      Amol += w;
      fMassFractionVector[i] = w;
    }
    for (i=0; i<fNumberOfElements; ++i) {
      fMassFractionVector[i] /= Amol;
    }

    fMassOfMolecule = Amol/Avogadro;
    ComputeDerivedQuantities();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

AddElement() [2/2]

void G4Material::AddElement	(	G4Element *	element,
		G4double	fraction
	)

Definition at line 406 of file G4Material.cc.

{
  if(fraction < 0.0 || fraction > 1.0) {
    G4cout << "G4Material::AddElement ERROR for " << fName << " and " 
       << element->GetName() << "  mass fraction= " << fraction 
       << " is wrong " << G4endl;
    G4Exception ("G4Material::AddElement()", "mat032", FatalException,     
                 "Attempt to add element with wrong mass fraction");
  }
  // initialization
  if (fNumberOfComponents == 0) {
    fMassFractionVector = new G4double[fArrayLength];
    fAtomsVector        = new G4int   [fArrayLength];
  }
  // filling ...
  if (fNumberOfComponents < maxNbComponents) {
    G4int el = 0;
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
    while ((el<fNumberOfElements)&&(element!=(*theElementVector)[el])) { ++el; }
    if (el<fNumberOfElements) fMassFractionVector[el] += fraction;
    else {
      theElementVector->push_back(element); 
      fMassFractionVector[el] = fraction;
      ++fNumberOfElements;
    }
    ++fNumberOfComponents;  
  } else {
    G4cout << "G4Material::AddElement ERROR for " << fName << " nElement= " 
       <<  fNumberOfElements << G4endl;
    G4Exception ("G4Material::AddElement()", "mat033", FatalException, 
           "Attempt to add more than the declared number of elements.");
  }    

  // filled.
  if (fNumberOfComponents == maxNbComponents) {

    G4int i=0;
    G4double Zmol(0.), Amol(0.);
    // check sum of weights -- OK?
    G4double wtSum(0.0);
    for (i=0; i<fNumberOfElements; ++i) {
      wtSum += fMassFractionVector[i];
      Zmol +=  fMassFractionVector[i]*(*theElementVector)[i]->GetZ();
      Amol +=  fMassFractionVector[i]*(*theElementVector)[i]->GetA();
    }
    if (std::fabs(1.-wtSum) > perThousand) {
      G4cerr << "WARNING !! for " << fName << " sum of fractional masses "
         <<  wtSum << " is not 1 - results may be wrong" 
         << G4endl;
    }
    for (i=0; i<fNumberOfElements; ++i) {
      fAtomsVector[i] = 
    G4lrint(fMassFractionVector[i]*Amol/(*theElementVector)[i]->GetA());
    }
     
    ComputeDerivedQuantities();
  }
}

Here is the call graph for this function:

AddMaterial()

void G4Material::AddMaterial	(	G4Material *	material,
		G4double	fraction
	)

store massFraction of material component

Definition at line 469 of file G4Material.cc.

{
  if(fraction < 0.0 || fraction > 1.0) {
    G4cout << "G4Material::AddMaterial ERROR for " << fName << " and " 
       << material->GetName() << "  mass fraction= " << fraction 
       << " is wrong ";
    G4Exception ("G4Material::AddMaterial()", "mat034", FatalException,
                 "Attempt to add material with wrong mass fraction");      
  }
  // initialization
  if (fNumberOfComponents == 0) {
    fMassFractionVector = new G4double[fArrayLength];
    fAtomsVector        = new G4int   [fArrayLength];
  }

  G4int nelm = material->GetNumberOfElements();

  // arrays should be extended
  if(nelm > 1) {
    G4int nold    = fArrayLength;
    fArrayLength += nelm - 1;
    G4double* v1 = new G4double[fArrayLength];
    G4int* i1    = new G4int[fArrayLength];
    for(G4int i=0; i<nold; ++i) {
      v1[i] = fMassFractionVector[i];
      i1[i] = fAtomsVector[i];
    }
    delete [] fAtomsVector;
    delete [] fMassFractionVector;
    fMassFractionVector = v1;
    fAtomsVector = i1;
  }

  // filling ...
  if (fNumberOfComponents < maxNbComponents) {
    for (G4int elm=0; elm<nelm; ++elm)
      {
        G4Element* element = (*(material->GetElementVector()))[elm];
        G4int el = 0;
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
        while ((el<fNumberOfElements)&&(element!=(*theElementVector)[el])) el++;
        if (el < fNumberOfElements) fMassFractionVector[el] += fraction
                                          *(material->GetFractionVector())[elm];
        else {
      theElementVector->push_back(element); 
          fMassFractionVector[el] = fraction
                                      *(material->GetFractionVector())[elm];
          ++fNumberOfElements;
        }
      } 
    ++fNumberOfComponents;
    fMatComponents[material] = fraction;
      
  } else {
    G4cout << "G4Material::AddMaterial ERROR for " << fName << " nElement= " 
       <<  fNumberOfElements << G4endl;
    G4Exception ("G4Material::AddMaterial()", "mat035", FatalException, 
           "Attempt to add more than the declared number of components.");
  }    

  // filled.
  if (fNumberOfComponents == maxNbComponents) {
    G4int i=0;
    G4double Zmol(0.), Amol(0.);
    // check sum of weights -- OK?
    G4double wtSum(0.0);
    for (i=0; i<fNumberOfElements; ++i) {
      wtSum += fMassFractionVector[i];
      Zmol +=  fMassFractionVector[i]*(*theElementVector)[i]->GetZ();
      Amol +=  fMassFractionVector[i]*(*theElementVector)[i]->GetA();
    }
    if (std::fabs(1.-wtSum) > perThousand) {
      G4cout << "G4Material::AddMaterial WARNING !! for " << fName 
         << " sum of fractional masses "
         <<  wtSum << " is not 1 - results may be wrong" 
         << G4endl;
    }
    for (i=0; i<fNumberOfElements; ++i) {
      fAtomsVector[i] = 
    G4lrint(fMassFractionVector[i]*Amol/(*theElementVector)[i]->GetA());
    }
     
    ComputeDerivedQuantities();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeDerivedQuantities()

void G4Material::ComputeDerivedQuantities ( )

private

Definition at line 302 of file G4Material.cc.

{
  // Header routine to compute various properties of material.
  // 

  // Number of atoms per volume (per element), total nb of electrons per volume
  G4double Zi, Ai;
  TotNbOfAtomsPerVolume = 0.;
  if (VecNbOfAtomsPerVolume) { delete [] VecNbOfAtomsPerVolume; }
  VecNbOfAtomsPerVolume = new G4double[fNumberOfElements];
  TotNbOfElectPerVolume = 0.;
  for (G4int i=0; i<fNumberOfElements; ++i) {
     Zi = (*theElementVector)[i]->GetZ();
     Ai = (*theElementVector)[i]->GetA();
     VecNbOfAtomsPerVolume[i] = Avogadro*fDensity*fMassFractionVector[i]/Ai;
     TotNbOfAtomsPerVolume += VecNbOfAtomsPerVolume[i];
     TotNbOfElectPerVolume += VecNbOfAtomsPerVolume[i]*Zi;
  }
        
  ComputeRadiationLength();
  ComputeNuclearInterLength();

  if (fIonisation) { delete fIonisation; }
  fIonisation  = new G4IonisParamMat(this);
  if (fSandiaTable) { delete fSandiaTable; }
  fSandiaTable = new G4SandiaTable(this);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeNuclearInterLength()

void G4Material::ComputeNuclearInterLength ( )

private

Definition at line 569 of file G4Material.cc.

{
  static const G4double lambda0  = 35*CLHEP::g/CLHEP::cm2;
  static const G4double twothird = 2.0/3.0;
  G4double NILinv = 0.0;
  for (G4int i=0; i<fNumberOfElements; ++i) {
    G4int Z = G4lrint( (*theElementVector)[i]->GetZ());
    G4double A = (*theElementVector)[i]->GetN();
    if(1 == Z) {
      NILinv += VecNbOfAtomsPerVolume[i]*A;
    } else {
      NILinv += VecNbOfAtomsPerVolume[i]*G4Exp(twothird*G4Log(A));
    } 
  }
  NILinv *= amu/lambda0; 
  fNuclInterLen = (NILinv <= 0.0 ? DBL_MAX : 1./NILinv);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeRadiationLength()

void G4Material::ComputeRadiationLength ( )

private

Definition at line 558 of file G4Material.cc.

{
  G4double radinv = 0.0 ;
  for (G4int i=0;i<fNumberOfElements;++i) {
     radinv += VecNbOfAtomsPerVolume[i]*((*theElementVector)[i]->GetfRadTsai());
  }
  fRadlen = (radinv <= 0.0 ? DBL_MAX : 1./radinv);
}

Here is the caller graph for this function:

CopyPointersOfBaseMaterial()

void G4Material::CopyPointersOfBaseMaterial ( )

private

Definition at line 332 of file G4Material.cc.

{
  G4double factor = fDensity/fBaseMaterial->GetDensity();
  TotNbOfAtomsPerVolume = factor*fBaseMaterial->GetTotNbOfAtomsPerVolume();
  TotNbOfElectPerVolume = factor*fBaseMaterial->GetTotNbOfElectPerVolume();

  theElementVector = 
    const_cast<G4ElementVector*>(fBaseMaterial->GetElementVector());
  fMassFractionVector = 
    const_cast<G4double*>(fBaseMaterial->GetFractionVector());
  fAtomsVector = const_cast<G4int*>(fBaseMaterial->GetAtomsVector());

  const G4double* v = fBaseMaterial->GetVecNbOfAtomsPerVolume();
  if (VecNbOfAtomsPerVolume)  { delete [] VecNbOfAtomsPerVolume; }
  VecNbOfAtomsPerVolume = new G4double[fNumberOfElements];
  for (G4int i=0; i<fNumberOfElements; ++i) {
    VecNbOfAtomsPerVolume[i] = factor*v[i];
  }
  fRadlen = fBaseMaterial->GetRadlen()/factor;
  fNuclInterLen = fBaseMaterial->GetNuclearInterLength()/factor;

  if (fIonisation) { delete fIonisation; }
  fIonisation = new G4IonisParamMat(this);

  fSandiaTable = fBaseMaterial->GetSandiaTable();
  fIonisation->SetMeanExcitationEnergy(fBaseMaterial->GetIonisation()->GetMeanExcitationEnergy());
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetA()

G4double G4Material::GetA

(

)

const

Definition at line 638 of file G4Material.cc.

{ 
  if (fNumberOfElements > 1) { 
    G4cout << "G4Material ERROR in GetA. The material: " << fName 
       << " is a mixture.";
    G4Exception ("G4Material::GetA()", "mat037", FatalException,  
         "the Atomic mass is not well defined." );
  } 
  return  (*theElementVector)[0]->GetA();      
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetAtomicNumDensityVector()

const G4double* G4Material::GetAtomicNumDensityVector ( ) const

inline

Definition at line 216 of file G4Material.hh.

{return VecNbOfAtomsPerVolume;}

Here is the caller graph for this function:

GetAtomsVector()

const G4int* G4Material::GetAtomsVector ( ) const

inline

Definition at line 198 of file G4Material.hh.

{return fAtomsVector;}

Here is the caller graph for this function:

GetBaseMaterial()

const G4Material* G4Material::GetBaseMaterial ( ) const

inline

Definition at line 233 of file G4Material.hh.

{return fBaseMaterial;}

Here is the caller graph for this function:

GetChemicalFormula()

const G4String& G4Material::GetChemicalFormula ( ) const

inline

Definition at line 179 of file G4Material.hh.

{return fChemicalFormula;}

Here is the caller graph for this function:

GetDensity()

G4double G4Material::GetDensity ( ) const

inline

Definition at line 180 of file G4Material.hh.

{return fDensity;}

GetElectronDensity()

G4double G4Material::GetElectronDensity ( ) const

inline

Definition at line 217 of file G4Material.hh.

{return TotNbOfElectPerVolume;}

GetElement()

const G4Element* G4Material::GetElement ( G4int  iel ) const

inline

Definition at line 202 of file G4Material.hh.

{return (*theElementVector)[iel];}

Here is the caller graph for this function:

GetElementVector()

const G4ElementVector* G4Material::GetElementVector ( ) const

inline

Definition at line 190 of file G4Material.hh.

{return theElementVector;}

GetFractionVector()

const G4double* G4Material::GetFractionVector ( ) const

inline

Definition at line 194 of file G4Material.hh.

{return fMassFractionVector;}

Here is the caller graph for this function:

GetIndex()

size_t G4Material::GetIndex ( ) const

inline

Definition at line 262 of file G4Material.hh.

{return fIndexInTable;}

Here is the call graph for this function:

GetIonisation()

G4IonisParamMat* G4Material::GetIonisation ( ) const

inline

Definition at line 226 of file G4Material.hh.

{return fIonisation;}

Here is the caller graph for this function:

GetMassOfMolecule()

G4double G4Material::GetMassOfMolecule ( ) const

inline

Definition at line 242 of file G4Material.hh.

{return fMassOfMolecule;}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMatComponents()

const std::map<G4Material*,G4double>& G4Material::GetMatComponents ( ) const

inline

Definition at line 237 of file G4Material.hh.

                                                {return fMatComponents;}

Here is the caller graph for this function:

GetMaterial()

G4Material * G4Material::GetMaterial ( const G4String &  name, G4bool  warning = true  )

static

Definition at line 604 of file G4Material.cc.

{  
  // search the material by its name 
  for (size_t J=0 ; J<theMaterialTable.size() ; ++J)
    {
      if (theMaterialTable[J]->GetName() == materialName)
    { return theMaterialTable[J]; }
    }
   
  // the material does not exist in the table
  if (warning) {
    G4cout << "G4Material::GetMaterial() WARNING: The material: "
       << materialName 
       << " does not exist in the table. Return NULL pointer."
       << G4endl;
  }  
  return 0;          
}

Here is the call graph for this function:

GetMaterialPropertiesTable()

G4MaterialPropertiesTable* G4Material::GetMaterialPropertiesTable ( ) const

inline

Definition at line 252 of file G4Material.hh.

  {return fMaterialPropertiesTable;}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMaterialTable()

G4MaterialTable * G4Material::GetMaterialTable ( )

static

Definition at line 589 of file G4Material.cc.

{
  return &theMaterialTable;
}

GetName()

const G4String& G4Material::GetName ( void  ) const

inline

Definition at line 178 of file G4Material.hh.

{return fName;}

GetNuclearInterLength()

G4double G4Material::GetNuclearInterLength ( ) const

inline

Definition at line 223 of file G4Material.hh.

{return fNuclInterLen;}

Here is the caller graph for this function:

GetNumberOfElements()

size_t G4Material::GetNumberOfElements ( ) const

inline

Definition at line 186 of file G4Material.hh.

{return fNumberOfElements;}

GetNumberOfMaterials()

size_t G4Material::GetNumberOfMaterials ( )

static

Definition at line 596 of file G4Material.cc.

{
  return theMaterialTable.size();
}

Here is the caller graph for this function:

GetPressure()

G4double G4Material::GetPressure ( ) const

inline

Definition at line 183 of file G4Material.hh.

{return fPressure;}

Here is the caller graph for this function:

GetRadlen()

G4double G4Material::GetRadlen ( ) const

inline

Definition at line 220 of file G4Material.hh.

{return fRadlen;}

Here is the caller graph for this function:

GetSandiaTable()

G4SandiaTable* G4Material::GetSandiaTable ( ) const

inline

Definition at line 229 of file G4Material.hh.

{return fSandiaTable;}

Here is the caller graph for this function:

GetState()

G4State G4Material::GetState ( ) const

inline

Definition at line 181 of file G4Material.hh.

{return fState;}

Here is the caller graph for this function:

GetTemperature()

G4double G4Material::GetTemperature ( ) const

inline

Definition at line 182 of file G4Material.hh.

{return fTemp;}

Here is the caller graph for this function:

GetTotNbOfAtomsPerVolume()

G4double G4Material::GetTotNbOfAtomsPerVolume ( ) const

inline

Definition at line 209 of file G4Material.hh.

{return TotNbOfAtomsPerVolume;}

Here is the caller graph for this function:

GetTotNbOfElectPerVolume()

G4double G4Material::GetTotNbOfElectPerVolume ( ) const

inline

Definition at line 212 of file G4Material.hh.

{return TotNbOfElectPerVolume;}

Here is the caller graph for this function:

GetVecNbOfAtomsPerVolume()

const G4double* G4Material::GetVecNbOfAtomsPerVolume ( ) const

inline

Definition at line 206 of file G4Material.hh.

{return VecNbOfAtomsPerVolume;}

Here is the caller graph for this function:

GetZ()

G4double G4Material::GetZ

(

)

const

Definition at line 625 of file G4Material.cc.

{ 
  if (fNumberOfElements > 1) {
     G4cout << "G4Material ERROR in GetZ. The material: " << fName 
        << " is a mixture.";
     G4Exception ("G4Material::GetZ()", "mat036", FatalException, 
                  "the Atomic number is not well defined." );
  } 
  return (*theElementVector)[0]->GetZ();      
}

Here is the call graph for this function:

Here is the caller graph for this function:

InitializePointers()

void G4Material::InitializePointers ( )

private

Definition at line 260 of file G4Material.cc.

{
  theElementVector         = 0;
  fMassFractionVector      = 0;
  fAtomsVector             = 0;
  fMaterialPropertiesTable = 0;
    
  VecNbOfAtomsPerVolume    = 0;
  fBaseMaterial            = 0;

  fChemicalFormula         = "";

  // initilized data members
  fDensity  = 0.0;
  fState    = kStateUndefined;
  fTemp     = 0.0;
  fPressure = 0.0;
  maxNbComponents     = 0;
  fArrayLength        = 0;
  TotNbOfAtomsPerVolume = 0;
  TotNbOfElectPerVolume = 0; 
  fRadlen = 0.0;
  fNuclInterLen = 0.0;
  fMassOfMolecule = 0.0;

  fIonisation = 0;
  fSandiaTable = 0;

  // Store in the static Table of Materials
  fIndexInTable = theMaterialTable.size();
  for(size_t i=0; i<fIndexInTable; ++i) {
    if(theMaterialTable[i]->GetName() == fName) {
      G4cout << "G4Material WARNING: duplicate name of material "
         << fName << G4endl;
      break;
    }
  }
  theMaterialTable.push_back(this);
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

G4int G4Material::operator!=

(

const G4Material &

right

)

const

Definition at line 658 of file G4Material.cc.

{
  return (this != (G4Material *) &right);
}

Here is the caller graph for this function:

operator=()

const G4Material& G4Material::operator= ( const G4Material &  )

private

Here is the caller graph for this function:

operator==()

G4int G4Material::operator==

(

const G4Material &

right

)

const

Definition at line 651 of file G4Material.cc.

{
  return (this == (G4Material *) &right);
}

Here is the caller graph for this function:

SetChemicalFormula()

void G4Material::SetChemicalFormula ( const G4String &  chF )

inline

Definition at line 173 of file G4Material.hh.

{fChemicalFormula=chF;}

Here is the caller graph for this function:

SetMaterialPropertiesTable()

void G4Material::SetMaterialPropertiesTable ( G4MaterialPropertiesTable *  anMPT )

inline

Definition at line 249 of file G4Material.hh.

  {fMaterialPropertiesTable = anMPT;}

Here is the caller graph for this function:

SetName()

void G4Material::SetName ( const G4String &  name )

inline

Definition at line 283 of file G4Material.hh.

{fName=name;}

Here is the call graph for this function:

Here is the caller graph for this function:

Friends And Related Function Documentation

operator<< [1/3]

std::ostream& operator<< ( std::ostream &  flux, const G4Material *  material  )

friend

Definition at line 665 of file G4Material.cc.

{
  std::ios::fmtflags mode = flux.flags();
  flux.setf(std::ios::fixed,std::ios::floatfield);
  G4long prec = flux.precision(3);
  
  flux
    << " Material: "         << std::setw(8) <<  material->fName
    << " " << material->fChemicalFormula << " "
    << "  density: "         << std::setw(6) << std::setprecision(3)  
    << G4BestUnit(material->fDensity,"Volumic Mass") 
    << "  RadL: "            << std::setw(7)  << std::setprecision(3)  
    << G4BestUnit(material->fRadlen,"Length")
    << "  Nucl.Int.Length: " << std::setw(7)  << std::setprecision(3)  
    << G4BestUnit(material->fNuclInterLen,"Length") 
    << "\n" << std::setw(30)   
    << "  Imean: "           << std::setw(7)  << std::setprecision(3)  
    << G4BestUnit(material->GetIonisation()->GetMeanExcitationEnergy(),
          "Energy");
    
  if(material->fState == kStateGas) {
    flux
      << "  temperature: " << std::setw(6) << std::setprecision(2)  
      << (material->fTemp)/kelvin << " K"
      << "  pressure: "    << std::setw(6) << std::setprecision(2)   
      << (material->fPressure)/atmosphere << " atm";
  }
  flux << "\n";
  
  for (G4int i=0; i<material->fNumberOfElements; i++) {
    flux 
      << "\n   ---> " << (*(material->theElementVector))[i] 
      << "\n          ElmMassFraction: " 
      << std::setw(6)<< std::setprecision(2) 
      << (material->fMassFractionVector[i])/perCent << " %" 
      << "  ElmAbundance "     << std::setw(6)<< std::setprecision(2) 
      << 100*(material->VecNbOfAtomsPerVolume[i])
      /(material->TotNbOfAtomsPerVolume)
      << " % \n";
  }
  flux.precision(prec);    
  flux.setf(mode,std::ios::floatfield);
            
  return flux;
}

operator<< [2/3]

std::ostream& operator<< ( std::ostream &  flux, const G4Material &  material  )

friend

Definition at line 713 of file G4Material.cc.

{
  flux << &material;        
  return flux;
}

operator<< [3/3]

std::ostream& operator<< ( std::ostream &  flux, G4MaterialTable  MaterialTable  )

friend

Definition at line 721 of file G4Material.cc.

{
  //Dump info for all known materials
  flux << "\n***** Table : Nb of materials = " << MaterialTable.size() 
       << " *****\n" << G4endl;
        
  for (size_t i=0; i<MaterialTable.size(); ++i) { 
    flux << MaterialTable[i] << G4endl << G4endl; 
  }

  return flux;
}      

Member Data Documentation

fArrayLength

G4int G4Material::fArrayLength

private

Definition at line 319 of file G4Material.hh.

fAtomsVector

G4int* G4Material::fAtomsVector

private

Definition at line 325 of file G4Material.hh.

fBaseMaterial

const G4Material* G4Material::fBaseMaterial

private

Definition at line 349 of file G4Material.hh.

fChemicalFormula

G4String G4Material::fChemicalFormula

private

Definition at line 310 of file G4Material.hh.

fDensity

G4double G4Material::fDensity

private

Definition at line 311 of file G4Material.hh.

fIndexInTable

size_t G4Material::fIndexInTable

private

Definition at line 331 of file G4Material.hh.

fIonisation

G4IonisParamMat* G4Material::fIonisation

private

Definition at line 344 of file G4Material.hh.

fMassFractionVector

G4double* G4Material::fMassFractionVector

private

Definition at line 324 of file G4Material.hh.

fMassOfMolecule

G4double G4Material::fMassOfMolecule

private

Definition at line 350 of file G4Material.hh.

fMatComponents

std::map<G4Material*,G4double> G4Material::fMatComponents

private

Definition at line 351 of file G4Material.hh.

fMaterialPropertiesTable

G4MaterialPropertiesTable* G4Material::fMaterialPropertiesTable

private

Definition at line 327 of file G4Material.hh.

fName

G4String G4Material::fName

private

Definition at line 309 of file G4Material.hh.

fNuclInterLen

G4double G4Material::fNuclInterLen

private

Definition at line 342 of file G4Material.hh.

fNumberOfComponents

G4int G4Material::fNumberOfComponents

private

Definition at line 320 of file G4Material.hh.

fNumberOfElements

G4int G4Material::fNumberOfElements

private

Definition at line 322 of file G4Material.hh.

fPressure

G4double G4Material::fPressure

private

Definition at line 316 of file G4Material.hh.

fRadlen

G4double G4Material::fRadlen

private

Definition at line 341 of file G4Material.hh.

fSandiaTable

G4SandiaTable* G4Material::fSandiaTable

private

Definition at line 345 of file G4Material.hh.

fState

G4State G4Material::fState

private

Definition at line 313 of file G4Material.hh.

fTemp

G4double G4Material::fTemp

private

Definition at line 315 of file G4Material.hh.

maxNbComponents

G4int G4Material::maxNbComponents

private

Definition at line 318 of file G4Material.hh.

theElementVector

G4ElementVector* G4Material::theElementVector

private

Definition at line 323 of file G4Material.hh.

theMaterialTable

G4MaterialTable G4Material::theMaterialTable

staticprivate

Definition at line 330 of file G4Material.hh.

TotNbOfAtomsPerVolume

G4double G4Material::TotNbOfAtomsPerVolume

private

Definition at line 339 of file G4Material.hh.

TotNbOfElectPerVolume

G4double G4Material::TotNbOfElectPerVolume

private

Definition at line 340 of file G4Material.hh.

VecNbOfAtomsPerVolume

G4double* G4Material::VecNbOfAtomsPerVolume

private

Definition at line 338 of file G4Material.hh.

---

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

• G4Material.hh
• G4Material.cc