#=============================================================================
#       Allpix: simulation of the LUCID experiment (with visualisation)
#=============================================================================

# The LUCID TPX0 detector 
/allpix/det/setId        10000
/allpix/det/setPosition   -3.01   -7.600   0.0 mm
/allpix/det/setRotation   90.0     0.0    90.0   deg
/allpix/det/setFoilThickness 0.0174 mm
/allpix/det/setLowTHL 13. keV

# The LUCID TPX1 detector 
/allpix/det/setId        10001
/allpix/det/setPosition   0.0   18.3  18.3 mm
/allpix/det/setRotation   0.0    0.0   0.0 deg
/allpix/det/setFoilThickness 0.0174 mm
/allpix/det/setLowTHL 13. keV

# The LUCID TPX2 detector 
/allpix/det/setId        10002
/allpix/det/setPosition  18.3   18.3   0.0 mm
/allpix/det/setRotation   0.0  -90.0   0.0 deg
/allpix/det/setFoilThickness 0.0174 mm
/allpix/det/setLowTHL 13. keV

# The LUCID TPX3 detector 
/allpix/det/setId        10003
/allpix/det/setPosition   0.0   18.3 -18.3 mm
/allpix/det/setRotation   0.0  180.0   0.0 deg
/allpix/det/setFoilThickness 0.0174 mm
/allpix/det/setLowTHL 13. keV

# The LUCID TPX4 detector.
/allpix/det/setId        10004
/allpix/det/setPosition -18.3   18.3   0.0 mm
/allpix/det/setRotation   0.0   90.0   0.0 deg
/allpix/det/setFoilThickness 0.0174 mm
/allpix/det/setLowTHL 13. keV

# The test structure - the LUCID experiment.
/allpix/extras/setTestStructure LUCID
/allpix/extras/setTestStructurePosition 0.0 0.0 0.0 mm
/allpix/extras/setTestStructureRotation 0.0 0.0 0.0 deg
/allpix/extras/setTestStructureShieldingThickness 0.68 mm

# We're in space, of course.
/allpix/extras/setWorldMaterial Vacuum

# The maximum step length (applies only to the sensor).
/allpix/det/setMaxStepLengthSensor 10 um

# Initialize the simulation.
/run/initialize

# Build the detectors.
/allpix/det/update

##############################################################################
#                         Visualisation configuration
##############################################################################
/vis/scene/create
#/vis/scene/add/axes 0. 0. 0. 10. cm
/vis/scene/add/volume World -1 2
/vis/viewer/set/style s

/run/verbose 0
/control/verbose 0
/control/saveHistory
/tracking/verbose 0
/allpix/phys/verbose 0

/vis/open OIX 1024x768-100+100
#/vis/open OGLIXm 1024x768-100+100
#/vis/open RayTracer
#/vis/open OGLIQt
#/vis/open OGLSQt
#/vis/open OIX
#/vis/open OGLIX 1024x768-100+100
#/vis/open OGLSXm
#/vis/open DAWNFILE
#/vis/open OGLSX
#/vis/open OGL 600x600-0+0

#/vis/viewer/set/background 0.4 0.5 0.6
/vis/viewer/set/background 1.0 1.0 1.0
/vis/viewer/set/viewpointThetaPhi 20 50
#/vis/viewer/set/background 0 0 0
/vis/viewer/zoom 2.0
/vis/viewer/flush

# Draw trajectories at end of event, showing trajectory points as
# markers of size 2 pixels.
/vis/scene/add/trajectories
/vis/modeling/trajectories/create/drawByCharge
/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
/vis/scene/endOfEventAction accumulate

##############################################################################

/tracking/verbose 0

################################
# General Particle Source (GPS)
################################

# For testing the individual detectors.
#/gps/particle proton
#/gps/ene/type Mono
#/gps/energy 40 MeV
#/gps/pos/type Point
#
## Test TPX0
#/gps/pos/centre -6.0 25.0 2.0 mm
#/gps/direction 0 -1 0
#
## Test TPX1
#/gps/pos/centre -6.0 18.3 25.0 mm
#/gps/direction 0 0 -1
#
## Test TPX2
#/gps/pos/centre 25.0 18.3 6.0 mm
#/gps/direction  -1 0 0
#
## Test TPX3
#/gps/pos/centre 6.0 24.0 -25.0 mm
#/gps/direction 0 0 1
#
## Test TPX4
#/gps/pos/centre -25.0 18.3 -6.0 mm
#/gps/direction 1 0 0
#
## Test TPX1 and TPX3
#/gps/pos/centre -6.0 18.3 25.0 mm
#/gps/direction 0.2 0 -1

#
# Spherical source for outer space
#----------------------------------
/gps/ene/type Mono
/gps/energy 40 MeV
/gps/particle proton
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. mm
/gps/pos/radius 49.9 mm
/gps/pos/confine PseudoDetector_phys
/gps/ang/type cos

/gps/source/list

##############################################################################

# Run the simulation with the number of events specified.
/run/beamOn 50

##############################################################################