ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4tgbMaterialMixtureByWeight.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4tgbMaterialMixtureByWeight.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 //
29 // class G4tgbMaterialMixtureByWeight
30 
31 // History:
32 // - Created. P.Arce, CIEMAT (November 2007)
33 // -------------------------------------------------------------------------
34 
36 
37 #include "G4SystemOfUnits.hh"
38 #include "G4tgbMaterial.hh"
39 #include "G4tgbMaterialMgr.hh"
40 #include "G4tgrMessenger.hh"
41 
42 // -------------------------------------------------------------------------
44 {
45 }
46 
47 
48 // -------------------------------------------------------------------------
50 {
51 }
52 
53 
54 // -------------------------------------------------------------------------
56 {
57  theTgrMate = hg;
58 }
59 
60 
61 // -------------------------------------------------------------------------
63 {
64  //----- construct new G4Material with component materials (a mixture)
65  G4Material* mate = new G4Material( theTgrMate->GetName(),
71 #ifdef G4VERBOSE
73  {
74  G4cout << " G4tgbMaterialMixtureByWeight::BuildG4Material() -"
75  << " Constructing new G4Material:"
76  << " " << theTgrMate->GetName()
77  << " " << theTgrMate->GetDensity()/g*cm3
79  << " " << theTgrMate->GetState()
80  << " " << theTgrMate->GetTemperature()
81  << " " << theTgrMate->GetPressure() << G4endl;
82  }
83 #endif
84 
85  //--- Add components
86  G4Element* compElem;
87  G4Material* compMate;
89 
90  for( G4int ii = 0; ii < theTgrMate->GetNumberOfComponents(); ii++)
91  {
92  // Look if this component is an element
93  compElem = mf->FindOrBuildG4Element( GetComponent(ii), false );
94  if( compElem != 0 )
95  {
96  mate->AddElement( compElem, GetFraction( ii ) );
97  }
98  else // If it is not an element, it must be a material
99  {
100  compMate = mf->FindOrBuildG4Material( GetComponent(ii) );
101  if( compMate != 0 )
102  {
103  // If it is a material add it by weight fraction
104  mate->AddMaterial( compMate, GetFraction( ii ) );
105  }
106  else
107  {
108  G4String ErrMessage = "Component " + GetComponent(ii)
109  + " of material " + theTgrMate->GetName()
110  + "\n" + "is not an element nor a material !";
111  G4Exception("G4tgbMaterialMixtureByWeight::BuildG4Material()",
112  "InvalidSetup", FatalException, ErrMessage);
113  }
114  }
115  }
116 
117 #ifdef G4VERBOSE
119  {
120  G4cout << " Constructing new G4Material by weight: " << *mate << G4endl;
121  G4cout << " TEMPERATURE " << mate->GetTemperature()/kelvin
122  << " kelvin " << G4endl;
123  }
124 #endif
125 
126  return mate;
127 }