47 {
48
49#define gr UNITS::gram
50#define mole UNITS::mole
51#define cm3 UNITS::cm3
52
53
54 GeoElement* Nitrogen =
55 new GeoElement(
"Nitrogen",
"N", 7.0, 14.0067 *
gr /
mole);
56 GeoElement* Oxygen =
new GeoElement(
"Oxygen",
"O", 8.0, 15.9995 *
gr /
mole);
57 GeoElement* Argon =
new GeoElement(
"Argon",
"Ar", 18.0, 39.948 *
gr /
mole);
58 GeoElement* Hydrogen =
59 new GeoElement(
"Hydrogen",
"H", 1.0, 1.00797 *
gr /
mole);
60
61
62 double densityOfAir = 0.001214 *
gr /
cm3;
63 GeoMaterial* air = new GeoMaterial("Air", densityOfAir);
64 air->add(Nitrogen, 0.7494);
65 air->add(Oxygen, 0.2369);
66 air->add(Argon, 0.0129);
67 air->add(Hydrogen, 0.0008);
68 air->lock();
69
70
71
72
73 const GeoBox* worldBox =
74 new GeoBox(1000 * UNITS::cm, 1000 * UNITS::cm, 1000 * UNITS::cm);
75 const GeoLogVol* worldLog = new GeoLogVol("WorldLog", worldBox, air);
76 GeoPhysVol* world = new GeoPhysVol(worldLog);
77 return world;
78}