add my simulation setup for my experiment

toms_erste_welt.gf

@@ -0,0 +1,147 @@
+#! /usr/bin/env ./geometryfactor.py
+# Toms position dependency. My quick first model
+
+# all units in mm
+#BGO but Big 
+define=W=20
+# side length
+define=D=51.962
+#apothem
+define=apo=${D}*sqrt(3/4)
+
+# number of sides, apothem of hexagon, thickness but here divided by 2
+prism= 6, ${apo},${W}/2
+
+# define SSDs variables. photosensitive area is 1cmx1cm and thickness is 0.007cm
+define=wx=5
+define=wy=5
+define=wz=0.07
+
+# define 
+define=a=18.85
+define=b=36.65
+define=c=19.5
+define=d=37.5
+define=e=18.85
+define=f=27.5
+
+# sum thickness of teflon and millipor
+define=dc=0.5
+
+# box in the middle
+box= ${wx},${wx},${wz}
+move=0,0,(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# box on the sides
+box= ${wx},${wx},${wz}
+move=${b},0,(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# box on the sides
+box= ${wx},${wx},${wz}
+move=-${b},0,(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# box on the corners
+box= ${wx},${wx},${wz}
+move=0,${d},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# box on the corners
+box= ${wx},${wx},${wz}
+move=0,-${d},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# boxes between the edges
+box= ${wx},${wx},${wz}
+move= ${a},0, (${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# boxes between the edges
+define=dm=17
+box= ${wx},${wx},${wz}
+move= -${a},0,(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# boxes between the edges
+define=dm=17
+box= ${wx},${wx},${wz}
+move= 0,${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# boxes between the edges
+define=dm=17
+box= ${wx},${wx},${wz}
+move= 0,-${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# x,y,z: x:photodiodes, y: corner, z:height
+
+#boxes in on the sides of the edges
+box= ${wx},${wx},${wz}
+move=${b},${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+#boxes in on the sides of the edges
+box= ${wx},${wx},${wz}
+move=${b},-${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+#boxes in on the sides of the edges
+box= ${wx},${wx},${wz}
+move=-${b},${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+#boxes in on the sides of the edges
+box= ${wx},${wx},${wz}
+move=-${b},-${c},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# sidies
+define=dm=17
+box= ${wx},${wx},${wz}
+move= ${e},${f},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# sidie
+define=dm=17
+box= ${wx},${wx},${wz}
+move= -${e},${f},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# sidies
+define=dm=17
+box= ${wx},${wx},${wz}
+move= ${e},-${f},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+# sidie
+define=dm=17
+box= ${wx},${wx},${wz}
+move= -${e},-${f},(${W}+${dc})/2
+copy=-1
+move=0,0, -${W}-${dc}
+
+
+
+
+
+
+