G4MaterialPropertiesTable Class Reference

#include <G4MaterialPropertiesTable.hh>

Inheritance diagram for G4MaterialPropertiesTable:

Collaboration diagram for G4MaterialPropertiesTable:

Public Member Functions
	G4MaterialPropertiesTable ()
	~G4MaterialPropertiesTable ()
void	AddConstProperty (const char *key, G4double PropertyValue)
G4MaterialPropertyVector *	AddProperty (const char *key, G4double *PhotonEnergies, G4double *PropertyValues, G4int NumEntries)
void	AddProperty (const char *key, G4MaterialPropertyVector *opv)
void	RemoveConstProperty (const char *key)
void	RemoveProperty (const char *key)
G4double	GetConstProperty (const char *key)
G4bool	ConstPropertyExists (const char *key)
G4MaterialPropertyVector *	GetProperty (const char *key)
void	AddEntry (const char *key, G4double aPhotonEnergy, G4double aPropertyValue)
void	DumpTable ()
const std::map< G4String, G4MaterialPropertyVector *, std::less< G4String > > *	GetPropertiesMap () const
const std::map< G4String, G4double, std::less< G4String > > *	GetPropertiesCMap () const

Private Types
typedef std::map< G4String, G4MaterialPropertyVector *, std::less< G4String > >::iterator	MPTiterator
typedef std::map< G4String, G4double, std::less< G4String > >::iterator	MPTCiterator

Private Member Functions
G4MaterialPropertyVector *	SetGROUPVEL ()

Private Attributes
std::map< G4String, G4MaterialPropertyVector *, std::less< G4String > >	MPT
std::map< G4String, G4double, std::less< G4String > >	MPTC

Detailed Description

Definition at line 70 of file G4MaterialPropertiesTable.hh.

Member Typedef Documentation

MPTCiterator

typedef std::map< G4String, G4double, std::less<G4String> >::iterator G4MaterialPropertiesTable::MPTCiterator

private

Definition at line 135 of file G4MaterialPropertiesTable.hh.

MPTiterator

typedef std::map< G4String, G4MaterialPropertyVector*, std::less<G4String> >::iterator G4MaterialPropertiesTable::MPTiterator

private

Definition at line 131 of file G4MaterialPropertiesTable.hh.

Constructor & Destructor Documentation

G4MaterialPropertiesTable()

G4MaterialPropertiesTable::G4MaterialPropertiesTable

(

)

Definition at line 57 of file G4MaterialPropertiesTable.cc.

{
}

~G4MaterialPropertiesTable()

G4MaterialPropertiesTable::~G4MaterialPropertiesTable

(

)

Definition at line 65 of file G4MaterialPropertiesTable.cc.

{
  MPTiterator i;
  for (i = MPT.begin(); i != MPT.end(); ++i)
  {
    delete (*i).second;
  }
  MPT.clear();
  MPTC.clear();
}

Member Function Documentation

AddConstProperty()

void G4MaterialPropertiesTable::AddConstProperty ( const char *  key, G4double  PropertyValue  )

inline

Here is the caller graph for this function:

AddEntry()

void G4MaterialPropertiesTable::AddEntry ( const char *  key, G4double  aPhotonEnergy, G4double  aPropertyValue  )

inline

Here is the caller graph for this function:

AddProperty() [1/2]

G4MaterialPropertyVector* G4MaterialPropertiesTable::AddProperty ( const char *  key, G4double *  PhotonEnergies, G4double *  PropertyValues, G4int  NumEntries  )

inline

Here is the caller graph for this function:

AddProperty() [2/2]

void G4MaterialPropertiesTable::AddProperty ( const char *  key, G4MaterialPropertyVector *  opv  )

inline

ConstPropertyExists()

G4bool G4MaterialPropertiesTable::ConstPropertyExists ( const char *  key )

inline

Here is the caller graph for this function:

DumpTable()

void G4MaterialPropertiesTable::DumpTable

(

)

Definition at line 80 of file G4MaterialPropertiesTable.cc.

{
  MPTiterator i;
  for (i = MPT.begin(); i != MPT.end(); ++i)
  {
    G4cout << (*i).first << G4endl;
    if ( (*i).second != 0 )
    {
      (*i).second->DumpValues();
    }
    else
    {
      G4Exception("G4MaterialPropertiesTable::DumpTable()", "mat204",
                  JustWarning, "NULL Material Property Vector Pointer.");
    }
  }
  MPTCiterator j;
  for (j = MPTC.begin(); j != MPTC.end(); ++j)
  {
    G4cout << j->first << G4endl;
    if ( j->second != 0 )
    {
      G4cout << j->second << G4endl;
    }
    else
    {
      G4Exception("G4MaterialPropertiesTable::DumpTable()", "mat202",
                  JustWarning, "No Material Constant Property.");
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetConstProperty()

G4double G4MaterialPropertiesTable::GetConstProperty ( const char *  key )

inline

Here is the caller graph for this function:

GetPropertiesCMap()

const std::map< G4String, G4double, std::less<G4String> >* G4MaterialPropertiesTable::GetPropertiesCMap ( ) const

inline

Definition at line 120 of file G4MaterialPropertiesTable.hh.

{ return &MPTC; }

Here is the call graph for this function:

Here is the caller graph for this function:

GetPropertiesMap()

const std::map< G4String, G4MaterialPropertyVector*, std::less<G4String> >* G4MaterialPropertiesTable::GetPropertiesMap ( ) const

inline

Definition at line 118 of file G4MaterialPropertiesTable.hh.

{ return &MPT; }

Here is the caller graph for this function:

GetProperty()

G4MaterialPropertyVector* G4MaterialPropertiesTable::GetProperty ( const char *  key )

inline

Here is the caller graph for this function:

RemoveConstProperty()

void G4MaterialPropertiesTable::RemoveConstProperty ( const char *  key )

inline

Here is the caller graph for this function:

RemoveProperty()

void G4MaterialPropertiesTable::RemoveProperty ( const char *  key )

inline

SetGROUPVEL()

G4MaterialPropertyVector * G4MaterialPropertiesTable::SetGROUPVEL ( )

private

Definition at line 119 of file G4MaterialPropertiesTable.cc.

{
#ifdef G4MULTITHREADED
  G4AutoLock mptm(&materialPropertyTableMutex);
#endif

  // check if "GROUPVEL" already exists
  MPTiterator itr;
  itr = MPT.find("GROUPVEL");
  if(itr != MPT.end()) return itr->second;

  // fetch RINDEX data, give up if unavailable
  //
  G4MaterialPropertyVector *rindex = this->GetProperty("RINDEX");
  if (rindex==0)  { return 0; }

  // RINDEX exists but has no entries, give up
  //
  if ( rindex->GetVectorLength() == 0 ) { return 0; }

  // add GROUPVEL vector
  //
  G4MaterialPropertyVector* groupvel = new G4MaterialPropertyVector();

  // fill GROUPVEL vector using RINDEX values
  // rindex built-in "iterator" was advanced to first entry above
  //
  G4double E0 = rindex->Energy(0);
  G4double n0 = (*rindex)[0];

  if (E0 <= 0.)
  {
    G4Exception("G4MaterialPropertiesTable::SetGROUPVEL()", "mat205",
                FatalException, "Optical Photon Energy <= 0");
  }
                                                                                
  if ( rindex->GetVectorLength() >= 2 )
  {
    // good, we have at least two entries in RINDEX
    // get next energy/value pair

    G4double E1 = rindex->Energy(1);
    G4double n1 = (*rindex)[1];

    if (E1 <= 0.)
    {
      G4Exception("G4MaterialPropertiesTable::SetGROUPVEL()", "mat205",
                  FatalException, "Optical Photon Energy <= 0");
    }

    G4double vg;

    // add entry at first photon energy
    //
    vg = c_light/(n0+(n1-n0)/std::log(E1/E0));

    // allow only for 'normal dispersion' -> dn/d(logE) > 0
    //
    if((vg<0) || (vg>c_light/n0))  { vg = c_light/n0; }

    groupvel->InsertValues( E0, vg );

    // add entries at midpoints between remaining photon energies
    //

    for (size_t i = 2; i < rindex->GetVectorLength(); i++)
    {
      vg = c_light/( 0.5*(n0+n1)+(n1-n0)/std::log(E1/E0));

      // allow only for 'normal dispersion' -> dn/d(logE) > 0
      //
      if((vg<0) || (vg>c_light/(0.5*(n0+n1))))  { vg = c_light/(0.5*(n0+n1)); }
      groupvel->InsertValues( 0.5*(E0+E1), vg );

      // get next energy/value pair, or exit loop
      //
      E0 = E1;
      n0 = n1;
      E1 = rindex->Energy(i);
      n1 = (*rindex)[i];

      if (E1 <= 0.)
      {
        G4Exception("G4MaterialPropertiesTable::SetGROUPVEL()", "mat205",
                    FatalException, "Optical Photon Energy <= 0");
      }
    }

    // add entry at last photon energy
    //
    vg = c_light/(n1+(n1-n0)/std::log(E1/E0));

    // allow only for 'normal dispersion' -> dn/d(logE) > 0
    //
    if((vg<0) || (vg>c_light/n1))  { vg = c_light/n1; }
    groupvel->InsertValues( E1, vg );
  }
  else // only one entry in RINDEX -- weird!
  {
    groupvel->InsertValues( E0, c_light/n0 );
  }
                                                                                
  this->AddProperty( "GROUPVEL", groupvel );
                                                                                
  return groupvel;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

MPT

std::map<G4String, G4MaterialPropertyVector*, std::less<G4String> > G4MaterialPropertiesTable::MPT

private

Definition at line 129 of file G4MaterialPropertiesTable.hh.

MPTC

std::map< G4String, G4double, std::less<G4String> > G4MaterialPropertiesTable::MPTC

private

Definition at line 133 of file G4MaterialPropertiesTable.hh.

---

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

• G4MaterialPropertiesTable.hh
• G4MaterialPropertiesTable.cc