CCalMaterial Class Reference

#include <CCalMaterial.hh>

Inheritance diagram for CCalMaterial:

Collaboration diagram for CCalMaterial:

Public Types
enum	FractionType { FTWeight, FTVolume }

Public Member Functions
	CCalMaterial (G4String mat, double dens, int nelem, CCalMaterial **constituents, double *weights, FractionType=FTWeight)
	CCalMaterial (const CCalMaterial &)
virtual	~CCalMaterial ()
G4String	Name () const
double	Density () const
int	NElements () const
G4String	Element (int i) const
double	Weight (int i) const
G4bool	operator== (const CCalMaterial &) const
G4bool	operator!= (const CCalMaterial &) const
CCalMaterial &	operator= (const CCalMaterial &)

Protected Member Functions
	CCalMaterial ()
void	computeDensity (int nconst, CCalMaterial **constituents, double *weights, FractionType ft)
void	closeMaterial ()

Protected Attributes
G4String	name
double	density
int	nElem
G4String *	theElements
double *	theWeights

Friends
std::ostream &	operator<< (std::ostream &, const CCalMaterial &)

Detailed Description

Definition at line 36 of file CCalMaterial.hh.

Member Enumeration Documentation

enum CCalMaterial::FractionType

Enumerator
FTWeight
FTVolume

Definition at line 41 of file CCalMaterial.hh.

{FTWeight, FTVolume};

Constructor & Destructor Documentation

CCalMaterial::CCalMaterial	(	G4String	mat,
		double	dens,
		int	nelem,
		CCalMaterial **	constituents,
		double *	weights,
		FractionType	ft = FTWeight
	)

Definition at line 35 of file CCalMaterial.cc.

                               : name(mat), density(dens) {
  nElem = 0;
  
  int i=0;
  for (i=0; i<nconst; i++)
    nElem += constituents[i]->NElements();

  theElements = new G4String[nElem];
  theWeights  = new double[nElem];

  double factor;
  int nelem=0;
  for (i=0; i<nconst; i++) {
    if (ft==FTWeight)
      factor=1.0;
    else
      factor=constituents[i]->Density();
    for (int j=0; j<constituents[i]->NElements(); j++) {
      theElements[nelem] = constituents[i]->Element(j);
      theWeights[nelem]  = constituents[i]->Weight(j)* weights[i] * factor;
      nelem++;
    }
  }

  if (density<0) { //Let's compute density
    computeDensity(nconst, constituents, weights, ft);
  }
  closeMaterial();
}

Here is the call graph for this function:

CCalMaterial::CCalMaterial

(

const CCalMaterial &

mat

)

Definition at line 67 of file CCalMaterial.cc.

                                                 :
  name(mat.name), density(mat.density), nElem(mat.nElem) {
  theElements = new G4String[nElem];
  theWeights  = new double[nElem];
  for (int i=0; i<nElem; i++){
    theElements[i]=mat.theElements[i];
    theWeights[i]=mat.theWeights[i];
  }
}

CCalMaterial::~CCalMaterial ( )

virtual

Definition at line 77 of file CCalMaterial.cc.

                            {
  if (theElements)
    delete[] theElements;
  if (theWeights)
    delete[] theWeights;
}

CCalMaterial::CCalMaterial ( )

inlineprotected

Definition at line 63 of file CCalMaterial.hh.

{} //Default constructor

Member Function Documentation

void CCalMaterial::closeMaterial ( )

protected

Definition at line 129 of file CCalMaterial.cc.

                                 {
  int trueConst=0;

  double norm=0;

  for (int i=0; i<nElem; i++) {
    norm+=theWeights[i];
    if (theElements[i]!="") {
      trueConst++;
      for (int j=i+1; j<nElem; j++) {
    if(theElements[i]==theElements[j]){
      theWeights[i]+=theWeights[j];
      theElements[j]="";
    }
      }//for j
    } //if
  }//for i

  if (trueConst != nElem) {
    G4String* newConst = new G4String[trueConst];
    double* newWeight = new double[trueConst];
    
    int newi=0;
    for(int i=0; i<nElem; i++){
      if (theElements[i]!="") {
    newConst[newi]  = theElements[i];
    newWeight[newi] = theWeights[i]/norm;
    newi++;
      }
    }

#ifdef debug    
    G4cout << "\tGoing from " << nElem <<" constituents to " << trueConst << G4endl;
#endif
    nElem=trueConst;
    
    delete[] theElements;
    delete[] theWeights;

    theElements=newConst;
    theWeights=newWeight;
  }
  else { //Let's normalize the weights
    for (int i=0; i<nElem; i++)
      theWeights[i] = theWeights[i]/norm;
  }
}

Here is the caller graph for this function:

void CCalMaterial::computeDensity ( int  nconst, CCalMaterial **  constituents, double *  weights, FractionType  ft  )

protected

Definition at line 84 of file CCalMaterial.cc.

                                                    {
  double mass=0;
  double volume=0;
  for (int i=0; i<nconst; i++) {
    if (ft==FTWeight) {
      mass+=weights[i];
      volume+=(weights[i]/constituents[i]->Density());
    }
    else { //by volume
      mass+=(weights[i]*constituents[i]->Density());
      volume+=weights[i];
    }
  }
  density=mass/volume;
}

Here is the call graph for this function:

Here is the caller graph for this function:

double CCalMaterial::Density ( ) const

inline

Definition at line 52 of file CCalMaterial.hh.

{return density;}        //Density in g/cm3.

Here is the caller graph for this function:

G4String CCalMaterial::Element ( int  i ) const

inline

Definition at line 54 of file CCalMaterial.hh.

{return theElements[i];} //Should be protected.

Here is the caller graph for this function:

G4String CCalMaterial::Name ( ) const

inline

Definition at line 51 of file CCalMaterial.hh.

{return name;}           //Material name.

Here is the caller graph for this function:

int CCalMaterial::NElements ( ) const

inline

Definition at line 53 of file CCalMaterial.hh.

{return nElem;}          //Number of Elements.

Here is the caller graph for this function:

G4bool CCalMaterial::operator!=

(

const CCalMaterial &

mat

)

const

Definition at line 125 of file CCalMaterial.cc.

                                                            {
  return (name!=mat.name);
}

CCalMaterial & CCalMaterial::operator=

(

const CCalMaterial &

mat

)

Definition at line 102 of file CCalMaterial.cc.

                                                            {
  if(theElements)
    delete[] theElements;
  if(theWeights)
    delete[] theWeights;

  name=mat.name;
  density=mat.density;
  nElem=mat.nElem;
  
  theElements = new G4String[nElem];
  theWeights  = new double[nElem];
  for (int i=0; i<nElem; i++){
    theElements[i]=mat.theElements[i];
    theWeights[i]=mat.theWeights[i];
  }
  return *this;
}

G4bool CCalMaterial::operator==

(

const CCalMaterial &

mat

)

const

Definition at line 121 of file CCalMaterial.cc.

                                                            {
  return (name==mat.name);
}

double CCalMaterial::Weight ( int  i ) const

inline

Definition at line 55 of file CCalMaterial.hh.

{return theWeights[i];}  //Should be protected.

Here is the caller graph for this function:

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  os, const CCalMaterial &  mat  )

friend

Definition at line 177 of file CCalMaterial.cc.

                                                                {
  os << mat.name << G4endl;
  os << "Density= " << mat.density << " g/cm3. Number of Elements: "
     << mat.nElem << G4endl;
  for (int i=0; i<mat.nElem; i++)
    os << '\t' << mat.theElements[i] << '\t' << mat.theWeights[i] << G4endl;
  return os;
}

Member Data Documentation

double CCalMaterial::density

protected

Definition at line 71 of file CCalMaterial.hh.

G4String CCalMaterial::name

protected

Definition at line 70 of file CCalMaterial.hh.

int CCalMaterial::nElem

protected

Definition at line 72 of file CCalMaterial.hh.

G4String* CCalMaterial::theElements

protected

Definition at line 73 of file CCalMaterial.hh.

double* CCalMaterial::theWeights

protected

Definition at line 74 of file CCalMaterial.hh.

---

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

• source/geant4.10.03.p03/examples/advanced/composite_calorimeter/include/CCalMaterial.hh
• source/geant4.10.03.p03/examples/advanced/composite_calorimeter/src/CCalMaterial.cc