ADD chaos.{gf,gpt}

chaos.gf

@@ -0,0 +1,56 @@
+#! /usr/bin/env ./geometryfactor.py
+# CHAOS Jr.  HET without A-detectors
+
+φ1=0
+φ2=π/6
+Δr=0.1
+Δθ=0.5°
+
+#BGO
+define=W=20
+define=D=2*51.962
+prism= 6, ${D}/2*sqrt(3/4), ${W}/2
+
+#Aerogel
+#box = 62,62,40
+#move = 0,0,${W}/2+10
+prism= 4,31,20
+move= 0,0,40
+
+# SSDs
+define=w=0.3
+define=d=36.5
+define=n=36
+define=d1=17.4
+
+prism= ${n}, ${d}/2, ${w}/2
+#A
+move= 0,0,${W}/2+55
+#C
+copy=-1
+move= 0,0,-50
+#E
+copy=-1
+move= 0,0,-30
+
+#E123
+copy=-1
+move=-18.5,-10.53,-4
+copy=-1
+move=2*18.5,0,-4
+copy=-1
+move=-18.25,31.6,-4
+
+#A1
+prism= ${n}, ${d1}/2, ${w}/2
+move= 0,0,${W}/2+55
+#C1
+copy=-1
+move= 0,0,-50
+
+#?
+define=a=24.7
+θ₂=atan2(${d},${a}+${w})
+radius=sqrt(${d}**2 + (${a}+${w})**2)/2
+
+