Stephan/geometryfactor
1
0
Fork 2
Code Issues Pull requests Projects Releases Packages Wiki Activity

Compare commits

base: Stephan:1d26575ae7dc71a63a709a654da348145f4e88b2
Branches Tags
Stephan:master
..
compare: Stephan:9c4a219f9bdb3ac53a47e759dd702856b32b4abe
Branches Tags
Stephan:master

2 commits
1d26575ae7 ... 9c4a219f9b

Author SHA1 Message Date
Stephan I. Böttcher
9c4a219f9b chaos-jr: fixes 2023-11-20 13:25:02 +01:00
Stephan I. Böttcher
80e23279de geometryfactor: provide atan(), atan2() from math for arguments 2023-11-20 13:24:47 +01:00
2 changed files with 7 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

10
chaos-jr.gf
View file

@ -2,23 +2,25 @@
# CHAOS Jr. HET without A-detectors # CHAOS Jr. HET without A-detectors
φ1=0 φ1=0
φ2=π/12 φ2=π/6
Δr=0.1 Δr=0.1
Δθ=0.5° Δθ=0.5°
radius=26
#BGO #BGO
define=W=20 define=W=20
define=D=40 define=D=40
prism= 6, ${D}/2, ${W}/2 prism= 6, ${D}/2*sqrt(3/4), ${W}/2
# SSDs # SSDs
define=a=24.3 define=a=24.3
define=w=0.3 define=w=0.3
define=d=36.5 define=d=36.5
define=n=32 define=n=36
prism= ${n}, ${d}/2, ${w}/2 prism= ${n}, ${d}/2, ${w}/2
move= 0,0,-${a}/2 move= 0,0,-${a}/2
copy=-1 copy=-1
move= 0,0,${a} move= 0,0,${a}
θ₂=atan2(${d},${a}+${w})
radius=sqrt(${d}**2 + (${a}+${w})**2)/2

1
geometryfactor.py
View file

@ -2,6 +2,7 @@
from sys import argv, stdout, stderr from sys import argv, stdout, stderr
from math import pi as π from math import pi as π
from math import atan, atan2
import numpy import numpy
from numpy import sqrt, sin, cos, tan, array, arange, empty, cross, newaxis from numpy import sqrt, sin, cos, tan, array, arange, empty, cross, newaxis
from numpy.linalg import solve from numpy.linalg import solve