legend-pygeom-hades 0.1.0a0__py3-none-any.whl

pygeomhades/dimensions.py

@@ -0,0 +1,283 @@
+from __future__ import annotations
+
+from dbetto import AttrsDict
+
+
+def get_bottom_plate_metadata() -> AttrsDict:
+    """Extract the metadata describing the bottom plate."""
+
+    return AttrsDict(
+        {
+            "width": 750,
+            "depth": 750,
+            "height": 15,
+            "cavity": {
+                "width": 120,
+                "depth": 940,  # <!--475*2-->
+                "height": 20,
+            },
+        }
+    )
+
+
+def get_cryostat_metadata(det_type: str, order: int, xtal_slice: str) -> AttrsDict:
+    """Extract the metadata corresponding to the cryostat
+
+    In future this will be moved into external metadata.
+
+    Parameters
+    ----------
+    det_type
+        The detector type (should be icpc or bege).
+    order
+        The order number.
+    xtal_slice
+        The slice of the crystal (typically A or B).
+    """
+    cryostat = {
+        "width": 101.6,
+        "height": 122.2,
+        "thickness": 1.5,
+        "position_cavity_from_top": 1.5,
+        "position_cavity_from_bottom": 0.8,
+        "position_from_bottom": 250.0,
+    }
+    xl_orders = [3, 8, 9, 10]
+
+    if det_type == "bege":
+        cryostat["height"] = 122.2
+        cryostat["width"] = 101.6
+    elif (det_type == "icpc") and (order in xl_orders):
+        cryostat["width"] = 114.3
+        cryostat["height"] = 171.0
+    elif det_type == "icpc":
+        cryostat["width"] = 101.6
+        cryostat["height"] = 171.0
+    else:
+        msg = "Only detector type icpc or bege are supported."
+        raise ValueError(msg)
+
+    # override batch 9
+    if order == 9 and xtal_slice == "B":
+        cryostat["width"] = 107.95
+
+    return AttrsDict(cryostat)
+
+
+def get_castle_dimensions(table_num: int) -> AttrsDict:
+    """Extract the lead castle dimensions for a given table.
+
+    Parameters
+    ----------
+    table_num
+        The number of the table to use, can be 1 or 2.
+    """
+
+    if table_num == 1:
+        lead_castle = {
+            "base": {
+                "width": 480,
+                "depth": 450,
+                "height": 500,
+            },
+            "inner_cavity": {
+                "width": 300,
+                "depth": 250,
+                "height": 500,
+            },
+            "cavity": {
+                "width": 120,
+                "depth": 100,
+                "height": 400,
+            },
+            "top": {
+                "width": 300,
+                "depth": 300,
+                "height": 90,
+            },
+            "front": {
+                "width": 160,
+                "depth": 100,
+                "height": 400,
+            },
+        }
+    elif table_num == 2:
+        lead_castle = {
+            "base": {
+                "width": 350,
+                "depth": 350,
+                "height": 400,
+            },
+            "inner_cavity": {
+                "width": 250,
+                "depth": 250,
+                "height": 400,
+            },
+            "top": {
+                "width": 200,
+                "depth": 200,
+                "height": 50,
+            },
+            "copper_plate": {
+                "width": 350,
+                "depth": 350,
+                "height": 10,
+            },
+        }
+
+    else:
+        msg = "Table number must be 1 or 2"
+        raise ValueError(msg)
+
+    return AttrsDict(lead_castle)
+
+
+def get_source_metadata(source_type: str, meas_type: str = "") -> AttrsDict:
+    """Get the dimensions of the source and collimator.
+
+    Parameters
+    ----------
+    source_type
+        The type of source (am_collimated, am, ba, co or th)
+    meas_type
+        The measurement (for th only) either lat or top.
+    """
+    if source_type == "am_collimated":
+        source = {
+            "height": 2.0,
+            "width": 1.0,
+            "capsule": {
+                "width": 20,
+                "depth": None,
+                "height": 10,
+            },
+            "collimator": {
+                "width": 30,
+                "depth": 30,
+                "height": 65,
+                "beam_height": 25.6,
+                "beam_width": 1.0,
+                "window": 0.2,
+            },
+        }
+    elif source_type == "am":
+        source = {
+            "height": 0.1,
+            "width": 1.0,
+            "capsule": {
+                "width": 11.08,
+                "depth": 23.08,
+                "height": 2.02,
+            },
+        }
+    elif source_type == "co":
+        source = {
+            "height": 0.1,
+            "width": 5.0,
+            "foil": {"width": 20, "height": 0.5},
+            "al_ring": {"height": 3.0, "width_max": 30, "width_min": 20},
+        }
+    elif source_type == "ba":
+        source = {
+            "height": 0.1,
+            "width": 5.0,
+            "foil": {
+                "width": 26.0,
+                "height": 0.5,
+            },
+            "al_ring": {"height": 3.0, "width_max": 30, "width_min": 26},
+        }
+
+    elif source_type == "th":
+        source = {
+            "height": 1.0,
+            "width": 1.0,
+            "capsule": {
+                "height": 7.0,
+                "width": 2.0,
+            },
+            "epoxy": {"height": 2.2, "width": 1.6},
+            "plates": {"height": 2.0, "width": 8.0, "cavity_width": 2.0},
+            "collimator": {
+                "height": 30.0,
+                "depth": 30.0,
+                "width": 30.0,
+                "beam_height": 15.0,
+                "beam_width": 1.0,
+            },
+        }
+
+        if meas_type == "top":
+            source["offset_height"] = 0.0
+        elif meas_type == "lat":
+            source["offset_height"] = 18.0
+        else:
+            msg = "can only have top or lat measurements"
+            raise RuntimeError(msg)
+    else:
+        msg = f"source type can only be am_collimated, ba, co, am or th not {source_type}"
+        raise RuntimeError(msg)
+
+    return AttrsDict(source)
+
+
+def get_source_holder_metadata(source_type: str, meas_type: str = "lat") -> AttrsDict:
+    """Get the dimensions of the source holder.
+
+    Parameters
+    ----------
+    source_type
+        The type of source (am_collimated, am, ba, co or th)
+    meas_type
+        The measurement (for th only) either lat or top.
+    """
+
+    if source_type in ["co", "ba", "am_collimated"]:
+        source_holder = {
+            "source": {
+                "top_plate_height": 3.0,
+                "top_plate_width": 30.0,
+                "top_height": 10.0,
+                "top_inner_width": 20.0,
+                "top_bottom_height": 6.1,
+                "bottom_inner_width": 102.0,
+            },
+            "outer_width": 108.0,
+            "inner_width": 87.0,
+        }
+
+    elif source_type == "am":
+        source_holder = {
+            "source": {
+                "top_height": 10.0,
+                "top_inner_width": 7.39,
+                "top_inner_depth": 15.39,
+                "bottom_inner_width": 102.0,
+                "top_bottom_height": 5.6,
+                "top_plate_width": 11.08,
+                "top_plate_depth": 23.08,
+                "top_plate_height": 2.0,
+            },
+            "outer_width": 108.0,
+            "inner_width": 87.0,
+        }
+
+    elif source_type == "th":
+        source_holder = {
+            "source": {
+                "height": 30.0,
+                "cavity_width": 3.0,
+                "bottom_height": 3.0,
+                "bottom_width": 50.0,
+            },
+        }
+
+        if meas_type == "lat":
+            source_holder["outer_width"] = 181.6
+            source_holder["inner_width"] = 101.6
+            source_holder["lat"] = {"height": 65.0, "cavity_height": 60.0, "cavity_width": 50.0}
+    else:
+        msg = f"Source must be co, ba, am_collimated, am or th not {source_type}"
+        raise RuntimeError(msg)
+
+    return AttrsDict(source_holder)

pygeomhades/metadata.py

@@ -0,0 +1,28 @@
+from __future__ import annotations
+
+import copy
+from importlib import resources
+
+from dbetto import AttrsDict, TextDB
+
+
+class PublicMetadataProxy:
+    """Provides proxies to transparently replace legend hardware metadata with sample data."""
+
+    def __init__(self):
+        dummy = TextDB(resources.files("pygeomhades") / "configs" / "dummy_geom")
+
+        self.hardware = AttrsDict({"detectors": {"germanium": {"diodes": _DiodeProxy(dummy)}}})
+
+
+class _DiodeProxy:
+    def __init__(self, dummy_detectors: TextDB):
+        self.dummy_detectors = dummy_detectors
+
+    def __getitem__(self, det_name: str) -> AttrsDict:
+        det = self.dummy_detectors[det_name[0] + "99000A"]
+        m = copy.copy(det)
+        m.name = det_name
+        m.production.order = 0
+        m.production.slice = "A"
+        return m

pygeomhades/models/dummy/bottom_plate_dummy.gdml

@@ -0,0 +1,44 @@
+<?xml version="1.0" ?>
+<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">
+
+ 	<define>
+		<quantity name="bottom_plate_x" type="length" value="bottom_plate_width" unit="mm"/>
+		<quantity name="bottom_plate_y" type="length" value="bottom_plate_depth" unit="mm"/>
+		<quantity name="bottom_plate_z" type="length" value="bottom_plate_height" unit="mm"/>
+		<quantity name="cavity_bottom_plate_x" type="length" value="bottom_cavity_plate_width" unit="mm"/>
+		<quantity name="cavity_bottom_plate_y" type="length" value="bottom_cavity_plate_depth" unit="mm"/>
+		<quantity name="cavity_bottom_plate_z" type="length" value="bottom_cavity_plate_height" unit="mm"/>
+		<position name="pos_cavity_bottom_plate" y="(bottom_plate_y)/2"/>
+	</define>
+
+	<materials>
+		<element Z="13" formula="Al" name="Aluminium" >
+			<atom value="26.98" />
+		</element>
+		<material name="Al" formula="Al" >
+			<D value="2.7" unit="g/cm3" />
+			<fraction n="1" ref="Aluminium"/>
+		</material>
+	</materials>
+
+	<solids>
+		<box name="bottom_plate" x="bottom_plate_x" y="bottom_plate_y" z="bottom_plate_z"/>
+		<box name="cavity_bottom_plate" x="cavity_bottom_plate_x" y="cavity_bottom_plate_y" z="cavity_bottom_plate_z"/>
+		<subtraction name="final_bottom_plate" >
+			<first ref="bottom_plate"/> <second ref="cavity_bottom_plate"/>
+			<positionref ref="pos_cavity_bottom_plate" />
+		</subtraction>
+	</solids>
+
+	<structure>
+		<volume name="Bottom_plate">
+			<materialref ref="Al"/>
+			<solidref ref="final_bottom_plate"/>
+		</volume>
+	</structure>
+
+	<setup name="Default" version="1.0">
+		<world ref="Bottom_plate"/>
+	</setup>
+
+</gdml>

pygeomhades/models/dummy/cryostat_dummy.gdml

@@ -0,0 +1,57 @@
+<?xml version="1.0" ?>
+<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">
+
+    <define>
+        <constant name="TWOPI" value="2.*pi"/>
+        <quantity name="cryostat_z" type="length" value="cryostat_height" unit="mm"/>
+        <quantity name="cryostat_radius" type="length" value="cryostat_width/2" unit="mm"/>
+        <quantity name="cryostat_cavity_radius" type="length" value="(cryostat_width-2*cryostat_thickness)/2" unit="mm"/>
+        <quantity name="start_cavity_z" type="length" value="position_cryostat_cavity_fromTop" unit="mm"/>
+        <quantity name="end_cavity_z" type="length" value="(cryostat_height-position_cryostat_cavity_fromBottom)" unit="mm"/>
+    </define>
+
+	<materials>
+		<element Z="13" formula="Al" name="Aluminium" >
+			<atom value="26.98" />
+		</element>
+		<element Z="29" formula="Cu" name="Copper" >
+			<atom value="63.546" />
+		</element>
+		<element Z="82" formula="Pb" name="Lead" >
+			<atom value="207.2" />
+		</element>
+		<element Z="83" formula="Bi" name="Bismuth" >
+			<atom value="208.98" />
+		</element>
+		<material name="EN_AW-2011T8" formula="Al_Cu6BiPb" >
+			<D value="2.84" unit="g/cm3" />
+			<fraction n="0.06" ref="Copper" />
+			<fraction n="0.004" ref="Lead" />
+			<fraction n="0.004" ref="Bismuth" />
+			<fraction n="0.932" ref="Aluminium" />
+		</material>
+	</materials>
+
+	<solids>
+		<polycone name="cryostat" startphi="0" deltaphi="TWOPI" aunit="rad">
+			<zplane rmin="0" rmax="cryostat_radius" z="0" />
+			<zplane rmin="0" rmax="cryostat_radius" z="start_cavity_z" />
+			<zplane rmin="cryostat_cavity_radius" rmax="cryostat_radius" z="start_cavity_z" />
+			<zplane rmin="cryostat_cavity_radius" rmax="cryostat_radius" z="end_cavity_z" />
+			<zplane rmin="0" rmax="cryostat_radius" z="end_cavity_z" />
+			<zplane rmin="0" rmax="cryostat_radius" z="cryostat_z" />
+		</polycone>
+	</solids>
+
+	<structure>
+		<volume name="Cryostat">
+			<materialref ref="EN_AW-2011T8"/>
+			<solidref ref="cryostat"/>
+		</volume>
+	</structure>
+
+	<setup name="Default" version="1.0">
+		<world ref="Cryostat"/>
+	</setup>
+
+</gdml>

pygeomhades/models/dummy/holder_bege_dummy.gdml

@@ -0,0 +1,61 @@
+<?xml version="1.0" ?>
+<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">
+
+	<define>
+	  	<constant name="TWOPI" value="2.*pi"/>
+	  	<quantity name="max_radius" type="length"  value="max_radius_in_mm" unit="mm"/>
+	  	<quantity name="outer_height" type="length" value="outer_height_in_mm" unit="mm"/>
+	  	<quantity name="inner_height" type="length" value="inner_height_in_mm" unit="mm"/>
+	  	<quantity name="outer_radius" type="length" value="outer_radius_in_mm" unit="mm"/>
+	  	<quantity name="inner_radius" type="length" value="inner_radius_in_mm" unit="mm"/>
+	  	<quantity name="pos_top_ring_z" type="length"  value="position_top_ring_in_mm" unit="mm"/>
+	  	<quantity name="end_top_ring_z" type="length" value="end_top_ring_in_mm"  unit="mm"/>
+	</define>
+
+	<materials>
+	    <element Z="29" formula="Cu" name="Copper" >
+	        <atom value="63.546" />
+	    </element>
+	    <element Z="82" formula="Pb" name="Lead" >
+	        <atom value="207.2" />
+	    </element>
+	    <element Z="83" formula="Bi" name="Bismuth" >
+	        <atom value="208.98" />
+	    </element>
+	    <element Z="13" formula="Al" name="Aluminium" >
+	        <atom value="26.98" />
+	    </element>
+	    <material name="EN_AW-2011T8" formula="Al_Cu6BiPb" >
+	        <D value="2.84" unit="g/cm3" />
+	        <fraction n="0.06" ref="Copper" />
+	        <fraction n="0.004" ref="Lead" />
+	        <fraction n="0.004" ref="Bismuth" />
+	        <fraction n="0.932" ref="Aluminium" />
+	    </material>
+	</materials>
+
+	<solids>
+		<polycone name="holder" startphi="0" deltaphi="TWOPI" aunit="rad">
+	        <zplane rmin="inner_radius" rmax="outer_radius" z="0" />
+	        <zplane rmin="inner_radius" rmax="outer_radius" z="pos_top_ring_z" />
+	        <zplane rmin="inner_radius" rmax="max_radius" z="pos_top_ring_z" />
+	        <zplane rmin="inner_radius" rmax="max_radius" z="end_top_ring_z" />
+	        <zplane rmin="inner_radius" rmax="outer_radius" z=" end_top_ring_z" />
+	        <zplane rmin="inner_radius" rmax="outer_radius" z="inner_height" />
+	        <zplane rmin="0" rmax="outer_radius" z="inner_height" />
+	        <zplane rmin="0" rmax="outer_radius" z="outer_height" />
+	    </polycone>
+	</solids>
+
+	<structure>
+		<volume name="Holder">
+		 	<materialref ref="EN_AW-2011T8"/>
+		 	<solidref ref="holder"/>
+		</volume>
+	</structure>
+
+	<setup name="Default" version="1.0">
+		<world ref="Holder"/>
+	</setup>
+
+</gdml>

pygeomhades/models/dummy/holder_icpc_dummy.gdml

@@ -0,0 +1,78 @@
+<?xml version="1.0" ?>
+<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">
+
+	<define>
+  		<constant name="TWOPI" value="2.*pi"/>
+  		<quantity name="max_radius" type="length"  value="max_radius_in_mm" unit="mm"/>
+  		<quantity name="outer_height" type="length" value="outer_height_in_mm" unit="mm"/>
+  		<quantity name="inner_height" type="length" value="inner_height_in_mm" unit="mm"/>
+  		<quantity name="outer_radius" type="length" value="outer_radius_in_mm" unit="mm"/>
+  		<quantity name="inner_radius" type="length" value="inner_radius_in_mm" unit="mm"/>
+  		<quantity name="bottom_cyl_inner_radius" type="length" value="inner_bottom_cyl_radius_in_mm" unit="mm"/>
+  		<quantity name="bottom_cyl_outer_radius" type="length" value="outer_bottom_cyl_radius_in_mm" unit="mm"/>
+  		<quantity name="edge_height" type="length" value="edge_height_in_mm" unit="mm"/>
+  		<quantity name="pos_bottom_ring_z" type="length"  value="pos_bottom_ring_in_mm" unit="mm"/>
+  		<quantity name="pos_top_ring_z" type="length"  value="pos_top_ring_in_mm" unit="mm"/>
+  		<quantity name="end_top_ring_z" type="length" value="end_top_ring_in_mm"  unit="mm"/>
+  		<quantity name="end_bottom_ring_z" type="length" value="end_bottom_ring_in_mm"  unit="mm"/>
+  		<quantity name="end_bottom_cyl_outer_z" type="length" value="end_bottom_cyl_outer_in_mm" unit="mm" />
+  		<quantity name="end_bottom_cyl_inner_z" type="length" value="end_bottom_cyl_inner_in_mm" unit="mm" />
+ 	</define>
+
+ 	<materials>
+	  	<element Z="13" formula="Al" name="Aluminium" >
+		 	<atom value="26.98" />
+	 	</element>
+	 	<element Z="29" formula="Cu" name="Copper" >
+		 	<atom value="63.546" />
+	 	</element>
+	 	<element Z="82" formula="Pb" name="Lead" >
+		 	<atom value="207.2" />
+	 	</element>
+	 	<element Z="83" formula="Bi" name="Bismuth" >
+		 	<atom value="208.98" />
+	 	</element>
+	 	<material name="EN_AW-2011T8" formula="Al_Cu6BiPb" >
+		 	<D value="2.84" unit="g/cm3" />
+		 	<fraction n="0.06" ref="Copper" />
+		 	<fraction n="0.004" ref="Lead" />
+		 	<fraction n="0.004" ref="Bismuth" />
+		 	<fraction n="0.932" ref="Aluminium" />
+	 	</material>
+ 	</materials>
+
+ 	<solids>
+	 	<polycone name="holder" startphi="0" deltaphi="TWOPI" aunit="rad">
+			<zplane rmin="inner_radius" rmax="max_radius" z="0" />
+			<zplane rmin="inner_radius" rmax="max_radius" z="edge_height" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z="edge_height" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z="pos_top_ring_z" />
+			<zplane rmin="inner_radius" rmax="max_radius" z="pos_top_ring_z" />
+			<zplane rmin="inner_radius" rmax="max_radius" z="end_top_ring_z" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z=" end_top_ring_z" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z="pos_bottom_ring_z" />
+			<zplane rmin="inner_radius" rmax="max_radius" z="pos_bottom_ring_z" />
+			<zplane rmin="inner_radius" rmax="max_radius" z="end_bottom_ring_z" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z="end_bottom_ring_z" />
+			<zplane rmin="inner_radius" rmax="outer_radius" z="inner_height" />
+			<zplane rmin="bottom_cyl_inner_radius" rmax="outer_radius" z="inner_height" />
+			<zplane rmin="bottom_cyl_inner_radius" rmax="outer_radius" z="outer_height" />
+			<zplane rmin="bottom_cyl_inner_radius" rmax="bottom_cyl_outer_radius" z="outer_height" />
+			<zplane rmin="bottom_cyl_inner_radius" rmax="bottom_cyl_outer_radius" z="end_bottom_cyl_inner_z" />
+			<zplane rmin="0" rmax="bottom_cyl_outer_radius" z="end_bottom_cyl_inner_z" />
+			<zplane rmin="0" rmax="bottom_cyl_outer_radius" z="end_bottom_cyl_outer_z" />
+	 	</polycone>
+ 	</solids>
+
+ 	<structure>
+		<volume name="Holder">
+			<materialref ref="EN_AW-2011T8"/>
+			<solidref ref="holder"/>
+		</volume>
+ 	</structure>
+
+ 	<setup name="Default" version="1.0">
+		<world ref="Holder"/>
+ 	</setup>
+
+</gdml>

pygeomhades/models/dummy/lead_castle_table1_dummy.gdml

@@ -0,0 +1,71 @@
+<?xml version="1.0" ?>
+<gdml xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://service-spi.web.cern.ch/service-spi/app/releases/GDML/schema/gdml.xsd">
+
+	<define>
+	  	<quantity name="base_x_1" type="length" value="base_width_1" unit="mm" />
+	  	<quantity name="base_y_1" type="length" value="base_depth_1" unit="mm"/>
+	  	<quantity name="base_z_1" type="length" value="base_height_1" unit="mm"/>
+	  	<quantity name="inner_cavity_x_1" type="length" value="inner_cavity_width_1" unit="mm"/>
+	  	<quantity name="inner_cavity_y_1" type="length" value="inner_cavity_depth_1" unit="mm"/>
+	  	<quantity name="inner_cavity_z_1" type="length" value="inner_cavity_height_1" unit="mm"/>
+	  	<quantity name="cavity_x_1" type="length" value="cavity_width_1" unit="mm"/>
+	  	<quantity name="cavity_y_1" type="length" value="cavity_depth_1" unit="mm"/>
+	  	<quantity name="cavity_z_1" type="length" value="cavity_height_1" unit="mm"/>
+	  	<quantity name="top_x_1" type="length" value="top_width_1" unit="mm"/>
+	  	<quantity name="top_y_1" type="length" value="top_depth_1" unit="mm"/>
+	  	<quantity name="top_z_1" type="length" value="top_height_1" unit="mm"/>
+	  	<quantity name="front_x_1" type="length" value="front_width_1" unit="mm"/>
+	  	<quantity name="front_y_1" type="length" value="front_depth_1" unit="mm"/>
+	  	<quantity name="front_z_1" type="length" value="front_height_1" unit="mm"/>
+	  	<position name="pos_cavity_base_1" y="(inner_cavity_y_1/2+(base_y_1-inner_cavity_y_1)/4)" z="(inner_cavity_z_1-cavity_z_1)/2" />
+	  	<position name="pos_top_1" z="-(base_z_1+top_z_1)/2 - 0.01" />
+	  	<position name="pos_front_1" y="(base_y_1+front_y_1)/2 - 0.01" z="(base_z_1-front_z_1)/2" />
+	  	<position name="center" x="0" y="0" z="0" />
+ 	</define>
+
+	<materials>
+		<element Z="82" formula="Pb" name="Lead" >
+			<atom value="207.2" />
+		</element>
+		<material name="Pb" formula="Pb" >
+			<D value="2.84" unit="g/cm3" />
+			<fraction n="1" ref="Lead"/>
+		</material>
+	</materials>
+
+	<solids>
+		<box name="base_lead_castle_1" x="base_x_1" y="base_y_1" z="base_z_1"/>
+		<box name="inner_cavity_base_1" x="inner_cavity_x_1" y="inner_cavity_y_1" z="inner_cavity_z_1"/>
+		<box name="cavity_base_1" x="cavity_x_1" y="cavity_y_1" z="cavity_z_1"/>
+		<box name="top_lead_castle_1" x="top_x_1" y="top_y_1" z="top_z_1"/>
+		<box name="front_1" x="front_x_1" y="front_y_1" z="front_z_1"/>
+		<union name="total_cavity_1" >
+			<first ref="inner_cavity_base_1"/> <second ref="cavity_base_1"/>
+			<positionref ref="pos_cavity_base_1" />
+		</union>
+		<subtraction name="base_cavity_1" >
+			<first ref="base_lead_castle_1"/> <second ref="total_cavity_1"/>
+			<positionref ref="center" />
+		</subtraction>
+		<union name="top_base_1" >
+			<first ref="base_cavity_1"/> <second ref="top_lead_castle_1"/>
+			<positionref ref="pos_top_1" />
+		</union>
+		<union name="final_lead_castle_1" >
+			<first ref="top_base_1"/> <second ref="front_1"/>
+			<positionref ref="pos_front_1" />
+		</union>
+	</solids>
+
+	<structure>
+		<volume name="Lead_castle">
+			<materialref ref="G4_Pb"/>
+			<solidref ref="final_lead_castle_1"/>
+		</volume>
+	</structure>
+
+	<setup name="Default" version="1.0">
+		<world ref="Lead_castle"/>
+	</setup>
+
+</gdml>