Message: Creating a mixture of materials using G4Material::AddMaterial() Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Creating a mixture of materials using G4Material::AddMaterial() 

Forum: Geometry
Date: 19 Jun, 2017
From: Sam de Jong <Sam de Jong>

Hello GEANT4 forums,

I'm having some problems with a material definition. My goal is to create graphite with some boron impurities. As a check, I created Graphite using the nist database, then created a new material that contained only the graphite:

	//basic graphite with custom density
	G4Material *baseMat=  nist->BuildMaterialWithNewDensity("Graphite_1.63","G4_GRAPHITE",1.630*g/cm3);

	//new material with one component:
	pile_mat=new G4Material("impure", 1.630*g/cm3, 1);

	//add only graphite
	pile_mat->AddMaterial(baseMat, 100.*perCent);

	//print out the material tables
	G4cout << *(G4Material::GetMaterialTable()) << endl;

When I look at the material tables, I get this:

 Material: G4_GRAPHITE Graphite   density:  2.210 g/cm3   RadL:  19.320 cm   Nucl.Int.Length:  36.269 cm 
                       Imean:  78.000 eV 

   --->  Element: C (C)   Z =  6.0   N =    12   A = 12.011 g/mole
         --->  Isotope:   C12   Z =  6   N =  12   A =  12.00 g/mole   abundance: 98.930 %
         --->  Isotope:   C13   Z =  6   N =  13   A =  13.00 g/mole   abundance:  1.070 %
          ElmMassFraction: 100.00 %  ElmAbundance 100.00 % 

 Material: Graphite_1.63 Graphite   density:  1.630 g/cm3   RadL:  26.195 cm   Nucl.Int.Length:  49.174 cm 
                       Imean:  78.000 eV 

   --->  Element: C (C)   Z =  6.0   N =    12   A = 12.011 g/mole
         --->  Isotope:   C12   Z =  6   N =  12   A =  12.00 g/mole   abundance: 98.930 %
         --->  Isotope:   C13   Z =  6   N =  13   A =  13.00 g/mole   abundance:  1.070 %
          ElmMassFraction: 100.00 %  ElmAbundance 100.00 % 

 Material:   impure    density:  1.630 g/cm3   RadL:  26.195 cm   Nucl.Int.Length:  49.174 cm 
                       Imean:  81.000 eV 

   --->  Element: C (C)   Z =  6.0   N =    12   A = 12.011 g/mole
         --->  Isotope:   C12   Z =  6   N =  12   A =  12.00 g/mole   abundance: 98.930 %
         --->  Isotope:   C13   Z =  6   N =  13   A =  13.00 g/mole   abundance:  1.070 %
          ElmMassFraction: 100.00 %  ElmAbundance 100.00 % 

The impure material has a different Imean than the Graphite_1.63 material, even though 'impure' is composed only of 'Graphite_1.63'. It seems the mixture is using the imean for carbon, not graphite. Is there something I'm doing wrong here? How can I solve this problem?

Thank you for your assistance,

-Sam

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Creating a mixture of materials using G4Material::AddMaterial()   (michel maire - 21 Jun, 2017)
 Add Message Add Message
to: "Creating a mixture of materials using G4Material::AddMaterial()"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]