########################
# Simple shielding example
########################
d:Ge/World/HLX = 2 m
d:Ge/World/HLY = 2 m
d:Ge/World/HLZ = 2.1 m
s:Ge/wall/Type = "TsBox"
s:Ge/wall/Parent = "World"
s:Ge/wall/Material = "G4_CONCRETE"
d:Ge/wall/HLX = 1 m
d:Ge/wall/HLY = 1 m
d:Ge/wall/HLZ = 1 m
d:Ge/wall/TransX = 0 m
d:Ge/wall/TransY = 0 m
d:Ge/wall/TransZ = 0 m
d:Ge/wall/RotX = 0 deg
d:Ge/wall/RotY = 0 deg
d:Ge/wall/RotZ = 0 deg
s:Ge/wall/Color = "red"
i:Ge/wall/XBins = 1
i:Ge/wall/YBins = 1
i:Ge/wall/ZBins = 10
#########################
# Gemetry for vrt
#########################
d:Ge/subComponent/Thickness = 20 cm
d:Ge/subComponent/StartPosition = -1.0 * Ge/wall/HLZ cm
s:Ge/VrtParallelWorld/Type = "TsBox"
s:Ge/VrtParallelWorld/Parent = "World"
d:Ge/VrtParallelWorld/HLX = 1.1 m
d:Ge/VrtParallelWorld/HLY = 1.1 m
d:Ge/VrtParallelWorld/HLZ = 1.1 m
d:Ge/VrtParallelWorld/TransX = 0 cm
d:Ge/VrtParallelWorld/TransY = 0 cm
d:Ge/VrtParallelWorld/TransZ = 0 cm
d:Ge/VrtParallelWorld/RotX = 0 deg
d:Ge/VrtParallelWorld/RotY = 0 deg
d:Ge/VrtParallelWorld/RotZ = 0 deg
b:Ge/VrtParallelWorld/IsParallel = "true"
s:Ge/subComponent1/Type = "TsBox"
s:Ge/subComponent1/Parent = "VrtParallelWorld"
d:Ge/subComponent1/HLX = Ge/Wall/HLX m
d:Ge/subComponent1/HLY = Ge/Wall/HLY m
d:Ge/subComponent1/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent1/TransX = 0 cm
d:Ge/subComponent1/TransY = 0 cm
d:Ge/subComponent1/TransZ = Ge/subComponent/StartPosition + Ge/subComponent1/HLZ cm
d:Ge/subComponent1/RotX = 0 deg
d:Ge/subComponent1/RotY = 0 deg
d:Ge/subComponent1/RotZ = 0 deg
b:Ge/subComponent1/IsParallel = "true"
s:Ge/subComponent2/Type = "TsBox"
s:Ge/subComponent2/Parent = "VrtParallelWorld"
d:Ge/subComponent2/HLX = Ge/Wall/HLX m
d:Ge/subComponent2/HLY = Ge/Wall/HLY m
d:Ge/subComponent2/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent2/TransX = 0 cm
d:Ge/subComponent2/TransY = 0 cm
d:Ge/subComponent2/TransZTmp = Ge/subComponent1/TransZ + Ge/subComponent1/HLZ cm
d:Ge/subComponent2/TransZ = Ge/subComponent2/TransZTmp + Ge/subComponent2/HLZ cm
d:Ge/subComponent2/RotX = 0 deg
d:Ge/subComponent2/RotY = 0 deg
d:Ge/subComponent2/RotZ = 0 deg
b:Ge/subComponent2/IsParallel = "true"
s:Ge/subComponent3/Type = "TsBox"
s:Ge/subComponent3/Parent = "VrtParallelWorld"
d:Ge/subComponent3/HLX = Ge/Wall/HLX m
d:Ge/subComponent3/HLY = Ge/Wall/HLY m
d:Ge/subComponent3/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent3/TransX = 0 cm
d:Ge/subComponent3/TransY = 0 cm
d:Ge/subComponent3/TransZTmp = Ge/subComponent2/TransZ + Ge/subComponent2/HLZ cm
d:Ge/subComponent3/TransZ = Ge/subComponent3/TransZTmp + Ge/subComponent3/HLZ cm
d:Ge/subComponent3/RotX = 0 deg
d:Ge/subComponent3/RotY = 0 deg
d:Ge/subComponent3/RotZ = 0 deg
b:Ge/subComponent3/IsParallel = "true"
s:Ge/subComponent4/Type = "TsBox"
s:Ge/subComponent4/Parent = "VrtParallelWorld"
d:Ge/subComponent4/HLX = Ge/Wall/HLX m
d:Ge/subComponent4/HLY = Ge/Wall/HLY m
d:Ge/subComponent4/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent4/TransX = 0 cm
d:Ge/subComponent4/TransY = 0 cm
d:Ge/subComponent4/TransZTmp = Ge/subComponent3/TransZ + Ge/subComponent3/HLZ cm
d:Ge/subComponent4/TransZ = Ge/subComponent4/TransZTmp + Ge/subComponent4/HLZ cm
d:Ge/subComponent4/RotX = 0 deg
d:Ge/subComponent4/RotY = 0 deg
d:Ge/subComponent4/RotZ = 0 deg
b:Ge/subComponent4/IsParallel = "true"
s:Ge/subComponent5/Type = "TsBox"
s:Ge/subComponent5/Parent = "VrtParallelWorld"
d:Ge/subComponent5/HLX = Ge/Wall/HLX m
d:Ge/subComponent5/HLY = Ge/Wall/HLY m
d:Ge/subComponent5/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent5/TransX = 0 cm
d:Ge/subComponent5/TransY = 0 cm
d:Ge/subComponent5/TransZTmp = Ge/subComponent4/TransZ + Ge/subComponent4/HLZ cm
d:Ge/subComponent5/TransZ = Ge/subComponent5/TransZTmp + Ge/subComponent5/HLZ cm
d:Ge/subComponent5/RotX = 0 deg
d:Ge/subComponent5/RotY = 0 deg
d:Ge/subComponent5/RotZ = 0 deg
b:Ge/subComponent5/IsParallel = "true"
s:Ge/subComponent6/Type = "TsBox"
s:Ge/subComponent6/Parent = "VrtParallelWorld"
d:Ge/subComponent6/HLX = Ge/Wall/HLX m
d:Ge/subComponent6/HLY = Ge/Wall/HLY m
d:Ge/subComponent6/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent6/TransX = 0 cm
d:Ge/subComponent6/TransY = 0 cm
d:Ge/subComponent6/TransZTmp = Ge/subComponent5/TransZ + Ge/subComponent5/HLZ cm
d:Ge/subComponent6/TransZ = Ge/subComponent6/TransZTmp + Ge/subComponent6/HLZ cm
d:Ge/subComponent6/RotX = 0 deg
d:Ge/subComponent6/RotY = 0 deg
d:Ge/subComponent6/RotZ = 0 deg
b:Ge/subComponent6/IsParallel = "true"
s:Ge/subComponent7/Type = "TsBox"
s:Ge/subComponent7/Parent = "VrtParallelWorld"
d:Ge/subComponent7/HLX = Ge/Wall/HLX m
d:Ge/subComponent7/HLY = Ge/Wall/HLY m
d:Ge/subComponent7/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent7/TransX = 0 cm
d:Ge/subComponent7/TransY = 0 cm
d:Ge/subComponent7/TransZTmp = Ge/subComponent6/TransZ + Ge/subComponent6/HLZ cm
d:Ge/subComponent7/TransZ = Ge/subComponent7/TransZTmp + Ge/subComponent7/HLZ cm
d:Ge/subComponent7/RotX = 0 deg
d:Ge/subComponent7/RotY = 0 deg
d:Ge/subComponent7/RotZ = 0 deg
b:Ge/subComponent7/IsParallel = "true"
s:Ge/subComponent8/Type = "TsBox"
s:Ge/subComponent8/Parent = "VrtParallelWorld"
d:Ge/subComponent8/HLX = Ge/Wall/HLX m
d:Ge/subComponent8/HLY = Ge/Wall/HLY m
d:Ge/subComponent8/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent8/TransX = 0 cm
d:Ge/subComponent8/TransY = 0 cm
d:Ge/subComponent8/TransZTmp = Ge/subComponent7/TransZ + Ge/subComponent7/HLZ cm
d:Ge/subComponent8/TransZ = Ge/subComponent8/TransZTmp + Ge/subComponent8/HLZ cm
d:Ge/subComponent8/RotX = 0 deg
d:Ge/subComponent8/RotY = 0 deg
d:Ge/subComponent8/RotZ = 0 deg
b:Ge/subComponent8/IsParallel = "true"
s:Ge/subComponent9/Type = "TsBox"
s:Ge/subComponent9/Parent = "VrtParallelWorld"
d:Ge/subComponent9/HLX = Ge/Wall/HLX m
d:Ge/subComponent9/HLY = Ge/Wall/HLY m
d:Ge/subComponent9/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent9/TransX = 0 cm
d:Ge/subComponent9/TransY = 0 cm
d:Ge/subComponent9/TransZTmp = Ge/subComponent8/TransZ + Ge/subComponent8/HLZ cm
d:Ge/subComponent9/TransZ = Ge/subComponent9/TransZTmp + Ge/subComponent9/HLZ cm
d:Ge/subComponent9/RotX = 0 deg
d:Ge/subComponent9/RotY = 0 deg
d:Ge/subComponent9/RotZ = 0 deg
b:Ge/subComponent9/IsParallel = "true"
s:Ge/subComponent10/Type = "TsBox"
s:Ge/subComponent10/Parent = "VrtParallelWorld"
d:Ge/subComponent10/HLX = Ge/Wall/HLX m
d:Ge/subComponent10/HLY = Ge/Wall/HLY m
d:Ge/subComponent10/HLZ = 0.5 * Ge/subComponent/Thickness cm
d:Ge/subComponent10/TransX = 0 cm
d:Ge/subComponent10/TransY = 0 cm
d:Ge/subComponent10/TransZTmp = Ge/subComponent9/TransZ + Ge/subComponent9/HLZ cm
d:Ge/subComponent10/TransZ = Ge/subComponent10/TransZTmp + Ge/subComponent10/HLZ cm
d:Ge/subComponent10/RotX = 0 deg
d:Ge/subComponent10/RotY = 0 deg
d:Ge/subComponent10/RotZ = 0 deg
b:Ge/subComponent10/IsParallel = "true"
#########################
# VRT
#########################
b:Vr/UseVarianceReduction = "true"
b:Vr/ParticleSplit/Active = "true"
sv:Vr/ParticleSplit/ParticleName = 1 "neutron"
s:Vr/ParticleSplit/Component = "VrtParallelWorld"
sv:Vr/ParticleSplit/SubComponents = 10
"subComponent1" "subComponent2"
"subComponent3" "subComponent4"
"subComponent5" "subComponent6"
"subComponent7" "subComponent8"
"subComponent9" "subComponent10"
s:Vr/ParticleSplit/Type = "WeightWindow"
uv:Vr/ParticleSplit/WeightMap = 10 1.0 0.5 0.25 0.125 0.0625 0.03125 0.015625 0.0078125 0.00390625 0.001953125
dv:Vr/ParticleSplit/EnergyMap = 10 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 GeV
u:Vr/ParticleSplit/UpperLimitFactor = 1
u:Vr/ParticleSplit/SurvivalFactor = 1
i:Vr/ParticleSplit/MaximumSplitNumber = 100
s:Vr/ParticleSplit/PlaceOfAction = "OnBoundary"
#########################
# Scorer
#########################
s:Sc/scorer/Quantity = "DoseToMedium"
s:Sc/scorer/Component = "Wall"
s:Sc/scorer/OutputFile = "WeightWindow"
s:Sc/scorer/IfOutputFileAlreadyExists = "Overwrite"
sv:Sc/scorer/Report = 3 "Sum" "Mean" "Variance"
#########################
# Beam setting
#########################
s:So/Example/Type = "Beam"
s:So/Example/Component = "BeamPosition"
s:So/Example/BeamParticle = "neutron"
d:So/Example/BeamEnergy = 10 MeV
u:So/Example/BeamEnergySpread = 0.757504
s:So/Example/BeamPositionDistribution = "Gaussian"
s:So/Example/BeamPositionCutoffShape = "Ellipse"
d:So/Example/BeamPositionCutoffX = 10. cm
d:So/Example/BeamPositionCutoffY = 10. cm
d:So/Example/BeamPositionSpreadX = 0.65 cm
d:So/Example/BeamPositionSpreadY = 0.65 cm
s:So/Example/BeamAngularDistribution = "Gaussian"
d:So/Example/BeamAngularCutoffX = 90. deg
d:So/Example/BeamAngularCutoffY = 90. deg
d:So/Example/BeamAngularSpreadX = 0.0032 rad
d:So/Example/BeamAngularSpreadY = 0.0032 rad
i:So/Example/NumberOfHistoriesInRun = 10
d:Ge/BeamPosition/TransZ = -1.0 * Ge/World/HLZ m
d:Ge/BeamPosition/RotX = 0 deg
i:Ts/ShowHistoryCountAtInterval = 1
b:Ts/PauseBeforeQuit = "True"
s:Gr/view/Type = "OpenGl"
b:Gr/view/IncludeAxes = "True"