G4gsmixt.cc File Reference

#include <iomanip>
#include "globals.hh"
#include "G4SystemOfUnits.hh"
#include "G3toG4.hh"
#include "G3EleTable.hh"
#include "G3MatTable.hh"
#include "G4Material.hh"
#include "G4Isotope.hh"

Go to the source code of this file.

Functions
void	PG4gsmixt (G4String *tokens)
void	G4gsmixt (G4int imate, G4String name, G4double *a, G4double *z, G4double dens, G4int nlmat, G4double *wmat)

Function Documentation

G4gsmixt()

void G4gsmixt	(	G4int	imate,
		G4String	name,
		G4double *	a,
		G4double *	z,
		G4double	dens,
		G4int	nlmat,
		G4double *	wmat )

Definition at line 72 of file G4gsmixt.cc.

{
  // in Geant3:
  // After a call with ratios by number (negative number of elements), 
  // the ratio array is changed to the ratio by weight, so all successive 
  // calls with the same array must specify the number of elements as 
  // positive 
  G4int i=0;
  if (nlmat<0) {
    // in case of proportions given in atom counts (nlmat<0),
    // the wmat[i] are converted to weight fractions
    G4double aMol = 0.;
    for (i=0; i<std::abs(nlmat); i++) { 
      // total molecular weight 
      aMol += wmat[i]*a[i]; 
    }  
    if (aMol == 0.) {
      G4String text = "G4mixt: Total molecular weight in " + name + " = 0.";       
      G4Exception("G4gsmixt()", "G3toG40016", FatalException, text);
      return;
    }
    for (i=0; i<std::abs(nlmat); i++) {
      // weight fractions
      wmat[i] = wmat[i]*a[i]/aMol;
    }
  }
 
  // create material with given number of components
  // (elements)
 
  G4Material* material 
    = new G4Material(name, dens*g/cm3, std::abs(nlmat));
  for (i=0; i<std::abs(nlmat); i++) {
    // add units
    // G4Element* element = G4Element(z[i], a[i]*g/mole, name);
    G4Element* element = G3Ele.GetEle(z[i]);
    material->AddElement(element, wmat[i]);    
  }
 
  // add the material to the List
  G3Mat.put(imate, material);
}

Referenced by PG4gsmixt().

PG4gsmixt()

void PG4gsmixt

(

G4String *

tokens

)

Definition at line 42 of file G4gsmixt.cc.

{
    // fill the parameter containers
    G3fillParams(tokens,PTgsmixt);
 
    // interpret the parameters
    G4String name = Spar[0].data();
    G4int imate = Ipar[0];
    G4int nlmat = Ipar[1];
    //G4double dens = Rpar[0]*g/cm3;
    G4double dens = Rpar[0];
    G4double *a = Rpar + 1;
    G4double *z = Rpar + 1+std::abs(nlmat);
    G4double *wmat = Rpar + 1 + 2*std::abs(nlmat);
/*
    for (int i=0; i<std::abs(nlmat); i++){
      //Rpar[i]=Rpar[i]*g/mole;
      Rpar[i]=Rpar[i];
    };
*/
    G4gsmixt(imate,name,a,z,dens,nlmat,wmat);
}

Referenced by G3CLEval().