CCalMaterialFactory Class Reference

#include <CCalMaterialFactory.hh>

Collaboration diagram for CCalMaterialFactory:

Public Types
enum	MatDescription { byWeight, byVolume, byAtomic }

Public Member Functions
	~CCalMaterialFactory ()
G4Material *	findMaterial (const G4String &) const
G4Element *	findElement (const G4String &) const
G4Element *	addElement (const G4String &, const G4String &, G4double, G4double, G4double)
G4Material *	addMaterial (const G4String &nam, G4double density, G4int nconst, G4String mats[], G4double prop[], MatDescription md=byWeight)
void	readElements (const G4String &)
void	readMaterials (const G4String &)

Static Public Member Functions
static CCalMaterialFactory *	getInstance (const G4String &, const G4String &)
static CCalMaterialFactory *	getInstance (const G4String &)
static CCalMaterialFactory *	getInstance ()

Protected Member Functions
void	readElements (std::ifstream &)
void	readMaterials (std::ifstream &)

Private Member Functions
	CCalMaterialFactory ()
G4Material *	findG4Material (const G4String &) const
CCalMaterial *	findCCalMaterial (const G4String &) const
CCalAMaterial *	findCCalAMaterial (const G4String &) const
CCalMaterial *	addCCalMaterial (const G4String &nam, G4double density, G4int nconst, G4String mats[], G4double prop[], MatDescription md=byWeight)

Private Attributes
CCalMaterialTable	theCCalMaterials
CCalAMaterialTable	theCCalAMaterials

Static Private Attributes
static CCalMaterialFactory *	instance = 0
static G4String	elementfile = ""
static G4String	mixturefile = ""

Detailed Description

Definition at line 47 of file CCalMaterialFactory.hh.

Member Enumeration Documentation

MatDescription

enum CCalMaterialFactory::MatDescription

Enumerator
byWeight
byVolume
byAtomic

Definition at line 49 of file CCalMaterialFactory.hh.

{byWeight, byVolume, byAtomic};

Constructor & Destructor Documentation

~CCalMaterialFactory()

CCalMaterialFactory::~CCalMaterialFactory

(

)

Definition at line 96 of file CCalMaterialFactory.cc.

                                         {
  CCalMaterialTable::iterator ite;
  for(ite = theCCalMaterials.begin(); ite != theCCalMaterials.end(); ite++ ){
    delete *ite;
  }
  theCCalMaterials.clear();
  CCalAMaterialTable::iterator itea;
  for(itea = theCCalAMaterials.begin(); itea != theCCalAMaterials.end(); 
      itea++ ){
    delete *itea;
  }
  theCCalAMaterials.clear();
}

CCalMaterialFactory()

CCalMaterialFactory::CCalMaterialFactory ( )

private

Definition at line 437 of file CCalMaterialFactory.cc.

                                         {
  readElements (elementfile);
  readMaterials(mixturefile);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Function Documentation

addCCalMaterial()

CCalMaterial * CCalMaterialFactory::addCCalMaterial ( const G4String &  nam, G4double  density, G4int  nconst, G4String  mats[], G4double  prop[], MatDescription  md = byWeight  )

private

Definition at line 277 of file CCalMaterialFactory.cc.

                                             {
  ptrCCalMaterial* matcol=0;
  ptrCCalAMaterial* amatcol=0;

  if (md==byAtomic)
    amatcol = new ptrCCalAMaterial[nconst];
  else
    matcol = new ptrCCalMaterial[nconst];

  for (G4int i=0; i<nconst; i++){
    if (md==byAtomic) {
      CCalAMaterial* amat = findCCalAMaterial(mats[i]);
      if (amat)
    amatcol[i]=amat;
      else {
    G4cerr << "ERROR: Trying to build" << name << " out of unknown " 
         << mats[i] << "." << G4endl 
         << "Skiping this material!" << G4endl;
    delete[] amatcol;
    return 0;
      }
    } //by Atomic fractions
    else {
      CCalMaterial* mat = findCCalMaterial(mats[i]);
      if (mat)
    matcol[i]=mat;
      else {
    G4cerr << "ERROR: Trying to build" <<name << " out of unknown " 
         << mats[i] << "." << G4endl 
         << "Skiping this material!" << G4endl;
    delete[] matcol;
    return 0;
      }
    }
  } //for

  //Let's do the CCalMaterial!
  if (md==byAtomic) {
    CCalAMaterial* amaterial = new CCalAMaterial(name, density, nconst, 
                         amatcol, prop);
    delete[] amatcol;
    theCCalAMaterials.push_back(amaterial);
#ifdef ddebug
    G4cout << *amaterial << G4endl;
#endif
    return amaterial;
  } else {
    CCalMaterial::FractionType ft;
    if (md == byWeight)
      ft=CCalMaterial::FTWeight;
    else
      ft=CCalMaterial::FTVolume;
    CCalMaterial* material = new CCalMaterial(name, density, nconst, 
                        matcol, prop, ft);
    delete[] matcol;
    theCCalMaterials.push_back(material);
#ifdef ddebug
    G4cout << *material << G4endl;
#endif
    return material;
  }  
}

Here is the call graph for this function:

Here is the caller graph for this function:

addElement()

G4Element * CCalMaterialFactory::addElement	(	const G4String &	name,
		const G4String &	symbol,
		G4double	Z,
		G4double	A,
		G4double	density
	)

Definition at line 162 of file CCalMaterialFactory.cc.

                                         {

  G4Element* theEl = new G4Element(name, symbol, Z, A*g/mole);
  //Make it also as a material.
  CCalAMaterial* theMat = new CCalAMaterial(name,A,density);
  theCCalAMaterials.push_back(theMat);

#ifdef ddebug
  G4cout << "Element " << name << " created!" << G4endl;
#endif
  return theEl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

addMaterial()

G4Material * CCalMaterialFactory::addMaterial	(	const G4String &	nam,
		G4double	density,
		G4int	nconst,
		G4String	mats[],
		G4double	prop[],
		MatDescription	md = byWeight
	)

Definition at line 179 of file CCalMaterialFactory.cc.

                                           {
  addCCalMaterial(name, density, nconst, mats, prop, md);
  return findMaterial(name);
}

Here is the call graph for this function:

findCCalAMaterial()

CCalAMaterial * CCalMaterialFactory::findCCalAMaterial ( const G4String &  mat ) const

private

Definition at line 264 of file CCalMaterialFactory.cc.

        {
  for (unsigned int i=0; i<theCCalAMaterials.size(); i++)
    if (theCCalAMaterials[i]->Name()==mat){
#ifdef ddebug
      G4cout << "CCalMaterial " << mat << " found!" << G4endl;
#endif
      return theCCalAMaterials[i];
    }
  return 0;
}

Here is the caller graph for this function:

findCCalMaterial()

CCalMaterial * CCalMaterialFactory::findCCalMaterial ( const G4String &  mat ) const

private

Definition at line 251 of file CCalMaterialFactory.cc.

        {
  for (unsigned int i=0; i<theCCalMaterials.size(); i++)
    if (theCCalMaterials[i]->Name()==mat){
#ifdef ddebug
      G4cout << "CCalMaterial " << mat << " found!" << G4endl;
#endif
      return theCCalMaterials[i];
    }
  return (CCalMaterial*) findCCalAMaterial(mat);
}

Here is the call graph for this function:

Here is the caller graph for this function:

findElement()

G4Element * CCalMaterialFactory::findElement

(

const G4String &

mat

)

const

Definition at line 149 of file CCalMaterialFactory.cc.

                                                                      {
  const G4ElementTable  theElements = *(G4Element::GetElementTable());
  for (unsigned int i=0; i<theElements.size(); i++)
    if (theElements[i]->GetName()==mat){
#ifdef ddebug
      G4cout << "Element " << mat << " found!" << G4endl;
#endif
      return theElements[i];
    }
  return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

findG4Material()

G4Material * CCalMaterialFactory::findG4Material ( const G4String &  mat ) const

private

Definition at line 240 of file CCalMaterialFactory.cc.

                                                                          {
  const G4MaterialTable theG4Materials = *(G4Material::GetMaterialTable());
  for (unsigned int i=0; i<theG4Materials.size(); i++) {
    if (theG4Materials[i]->GetName()==mat){
      return theG4Materials[i];
    }
  }
  return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

findMaterial()

G4Material * CCalMaterialFactory::findMaterial

(

const G4String &

mat

)

const

Definition at line 111 of file CCalMaterialFactory.cc.

                                                                        {
  G4Material* theMat=findG4Material(mat);

  if (theMat) {
#ifdef ddebug
    G4cout << "Material " << mat << " already defined. Returning previous "
     << "instance." << G4endl;
#endif
    return theMat;
  } else { 
    CCalMaterial* CCalmat=findCCalMaterial(mat);
    if (CCalmat){
      G4Material* G4Mat = new G4Material(CCalmat->Name(),
                     CCalmat->Density()*g/cm3, 
                     CCalmat->NElements());
      for(G4int i=0; i<CCalmat->NElements(); i++) {
    G4Element* elem = findElement(CCalmat->Element(i));
    if (!elem) {
      G4ExceptionDescription ed;
      ed << "   Could not build material " << mat << "." << G4endl;
      G4Exception("CCalMaterialFactory::findMaterial()","ccal001",
              FatalException,ed);
    }
    G4Mat->AddElement(elem, CCalmat->Weight(i));
      }
#ifdef ddebug
    G4cout << "Material " << mat << " has been built successfully." << G4endl;
#endif
      return G4Mat;
    } else {
      G4cerr << "ERROR: Material " << mat << " not found in CCal database!!!" 
       << G4endl;
      return 0;
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

getInstance() [1/3]

CCalMaterialFactory * CCalMaterialFactory::getInstance ( const G4String &  matfile, const G4String &  mixfile  )

static

Definition at line 53 of file CCalMaterialFactory.cc.

                                                      {
  if ((matfile=="" || matfile==elementfile) &&
      (mixfile=="" || mixfile==mixturefile))
    return getInstance();
  else if ((matfile != "" && elementfile != "" && matfile != elementfile) ||
           (mixfile != "" && mixturefile != "" && mixfile != mixturefile)) {
    G4cerr << "ERROR: Trying to get materials from " << matfile << " and " 
         << mixfile << " while previously were retrieved from " 
         << elementfile << " and " << mixturefile << "." << G4endl;
    return 0;
  } else {
    if (elementfile == "") 
      elementfile=matfile;
    if (mixturefile == "") 
      mixturefile=mixfile;
    return getInstance();
  }
}

Here is the call graph for this function:

getInstance() [2/3]

CCalMaterialFactory * CCalMaterialFactory::getInstance ( const G4String &  matfile )

static

Definition at line 74 of file CCalMaterialFactory.cc.

                                                                            {
  return getInstance(matfile,matfile);
}

Here is the call graph for this function:

getInstance() [3/3]

CCalMaterialFactory * CCalMaterialFactory::getInstance ( )

static

Definition at line 79 of file CCalMaterialFactory.cc.

                                                     {
  if (elementfile=="" || mixturefile=="") {
    G4cerr << "ERROR: You haven't defined files to be used for materials in "
     << "CCalMaterialFactory::getInstance(const G4String&,const G4String&)"
     << G4endl;
    return 0;
  }

  if (instance==0) {
    instance = new CCalMaterialFactory;
    return instance;
  }
  else
    return instance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

readElements() [1/2]

void CCalMaterialFactory::readElements

(

const G4String &

matfile

)

Definition at line 190 of file CCalMaterialFactory.cc.

                                                              {

  G4String path = "NULL";
  if (getenv("CCAL_GLOBALPATH"))
    path = getenv("CCAL_GLOBALPATH");

  G4cout << " ==> Opening file " << matfile << " to read elements..." << G4endl;
  std::ifstream is;
  bool ok = openGeomFile(is, path, matfile);
  if (!ok) {
    G4cerr << "ERROR: Could not open file " << matfile << G4endl;
    return;
  }

  // Find *DO GMAT
  findDO(is, G4String("GMAT"));
  
  readElements(is);
  
  is.close();
}

Here is the call graph for this function:

Here is the caller graph for this function:

readElements() [2/2]

void CCalMaterialFactory::readElements ( std::ifstream &  is )

protected

Definition at line 346 of file CCalMaterialFactory.cc.

                                                     {
  G4String name, symbol;
  
  G4cout << "     ==> Reading elements... " << G4endl;
#ifdef debug
  G4cout << "       Element    \tsymbol\tA\tZ\tdensity\tX_0          abs_l"<< G4endl;
#endif
  //There should come the list of materials. #. Defines a comment
  //*DO defines the beguining of the Mixes block.

  readName(is,name);
  while (name != "*ENDDO") {
    //It should be an element definition
    G4double A, Z, density;
    is >> symbol >> A >> Z >> density >> jump;    
#ifdef debug
    G4cout << "       " << name << "    \t" << symbol << "\t" 
     << A << "\t" << Z << "\t" << density << G4endl;
#endif    
    addElement(name, symbol, Z, A, density);
    readName(is,name);
  };
  G4cout << "     " << G4Element::GetElementTable()->size() 
       << " elements read from file" << G4endl << G4endl;
}

Here is the call graph for this function:

readMaterials() [1/2]

void CCalMaterialFactory::readMaterials

(

const G4String &

matfile

)

Definition at line 213 of file CCalMaterialFactory.cc.

                                                               {

  G4String path = "NULL";
  if (getenv("CCAL_GLOBALPATH"))
    path = getenv("CCAL_GLOBALPATH");

  G4cout << " ==> Opening file " << matfile << " to read materials..." << G4endl;
  std::ifstream is;
  bool ok = openGeomFile(is, path, matfile);
  if (!ok) {
    G4cerr << "ERROR: Could not open file " << matfile << G4endl;
    return;
  }

  // Find *DO GMIX
  findDO(is, G4String("GMIX"));

  readMaterials(is);

  is.close();
}

Here is the call graph for this function:

Here is the caller graph for this function:

readMaterials() [2/2]

void CCalMaterialFactory::readMaterials ( std::ifstream &  is )

protected

Definition at line 373 of file CCalMaterialFactory.cc.

                                                      {
  G4String name, matname;

  G4cout << "     ==> Reading materials... " << G4endl;

  //Take into account the special case of vacuum...
#ifdef debug
  G4cout <<"       \"Vacuum\"" << G4endl;
#endif
  G4double density  = universe_mean_density;   //from PhysicalConstants.h
  G4double pressure = 1.E-19*pascal;
  G4double temperature = 0.1*kelvin;
  new G4Material("Vacuum", /*Z=*/ 1., /*A=*/ 1.01*g/mole,
         density, kStateGas, temperature, pressure);

  //There should come the list of materials. #. Defines a comment
  //*ENDDO defines the block.
  readName(is,name);
  while (name != "*ENDDO") {
    //It should be a material definition
    matname=name;
    G4int nElem;
    G4double dens;
    is >> nElem >> dens >> jump;

#ifdef debug
    G4cout <<"       " << matname
       << " made of " << nElem 
       << " elements. Density=" << dens
       << G4endl;
#endif

    G4int absnelem = std::abs(nElem);

    G4String* mats    = new G4String[absnelem];
    G4double* weights = new G4double[absnelem];
    
    G4double prop;
    for(int i=0; i<absnelem; i++) {
      readName(is, name);
      is >> prop >> jump;
      mats[i]=name;
      weights[i]=std::abs(prop);
    } //for...
    MatDescription md;
    if (nElem>0 && prop<0)
      md = byAtomic;
    else if (nElem>0)
      md = byWeight;
    else
      md = byVolume;

    addCCalMaterial(matname, dens, absnelem, mats, weights, md);
    delete[] mats;
    delete[] weights;
    
    readName(is,name);
  };  //while

  G4cout << "     " << theCCalMaterials.size() << " materials read from " 
       << mixturefile << G4endl << G4endl;
}

Here is the call graph for this function:

Member Data Documentation

elementfile

G4String CCalMaterialFactory::elementfile = ""

staticprivate

Definition at line 92 of file CCalMaterialFactory.hh.

instance

CCalMaterialFactory * CCalMaterialFactory::instance = 0

staticprivate

Definition at line 90 of file CCalMaterialFactory.hh.

mixturefile

G4String CCalMaterialFactory::mixturefile = ""

staticprivate

Definition at line 93 of file CCalMaterialFactory.hh.

theCCalAMaterials

CCalAMaterialTable CCalMaterialFactory::theCCalAMaterials

private

Definition at line 96 of file CCalMaterialFactory.hh.

theCCalMaterials

CCalMaterialTable CCalMaterialFactory::theCCalMaterials

private

Definition at line 95 of file CCalMaterialFactory.hh.

---

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

• CCalMaterialFactory.hh
• CCalMaterialFactory.cc