#
i:Ts/NumberOfThreads = 1
# Secondary biasing occurs in the region
# of this component
s:Ge/target/Parent = "World"
s:Ge/target/Material = "G4_WATER"
s:Ge/target/Type = "TsBox"
d:Ge/target/HLX = 1.5 cm
d:Ge/target/HLY = 1.5 cm
d:Ge/target/HLZ = 0.5 mm
s:Ge/target/Color = "yellow"
# Set the region
s:Ge/target/AssignToRegionNaMed = "targetregion"
# Component for scoring
s:Ge/scorer/Parent = "World"
s:Ge/scorer/Material = "G4_WATER"
s:Ge/scorer/Type = "TsCylinder"
d:Ge/scorer/RMax = 5 cm
d:Ge/scorer/HL = 0.05 mm
d:Ge/scorer/SPhi = 0 deg
d:Ge/scorer/DPhi = 360 deg
d:Ge/scorer/TransZ = 1 m + Ge/scorer/HL
s:Ge/scorer/Color = "yellow"
i:Ge/scorer/RBins = 20
i:Ge/scorer/PhiBins = 1
i:Ge/scorer/ZBins = 1
s:Sc/scorer/Quantity = "EnergyFluence"
s:Sc/scorer/Component = "scorer"
s:Sc/scorer/OutputFile = "secondaryBiasing"
s:Sc/scorer/IfOutputFileAlreadyExists = "Overwrite"
# VISUALIZATION
s:Gr/view/Type = "OpenGl"
i:Gr/view/WindowSizeX = 600
i:Gr/view/WindowSizeY = 600
d:Gr/view/Phi = 0 deg
d:Gr/view/Theta = 270 deg
u:Gr/view/Zoom = 5
b:Ts/PauseBeforeQuit = "True"
# SOURCE
s:So/Example/Type = "Beam"
s:So/Example/Component = "BeamPosition"
s:So/Example/BeamParticle = "gamma"
d:So/Example/BeamEnergy = 10 MeV
u:So/Example/BeamEnergySpread = 0.1
s:So/Example/BeamPositionDistribution = "None"
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
Ge/BeamPosition/Parent = "World"
Ge/BeamPosition/TransX = 0.0 cm
Ge/BeamPosition/TransZ = -1.0 * Ge/target/HLZ cm
Ge/BeamPosition/TransY = 0.0 mm
Ge/BeamPosition/RotX = 0 deg
# GEOMETRY
d:Ge/World/HLX = 16.1 cm
d:Ge/World/HLY = 16.1 cm
d:Ge/World/HLZ = 1.5 m
s:Ge/World/Material = "Air"
sv:Ph/Default/Modules = 1 "g4em-standard_opt3"
i:So/Example/NumberOfHistoriesInRun = 1000
i:Ts/ShowHistoryCountAtInterval = 100
b:Ts/ShowCPUTime = "True"
#Forced interaction in a region
#A distance value shorter than the thickness of the volume
#of interest must be set.
#By setting that distance of interaction large enough,
#then the true physic length is taken into account. This
#will mimic activation deactivation of this technique.
s:Tf/ActiveVRT/Function = "step"
dv:Tf/ActiveVRT/Times = 5 1 2 3 4 5 ms
dv:Tf/ActiveVRT/Values = 5 0.1 Ge/World/HLZ 0.1 Ge/World/HLZ 0.1 mm
Tf/TimeLineEnd = 5 ms
Tf/NumberOfSequentialTimes = 5
b:Vr/UseVarianceReduction = "true"
s:Vr/ForcedInteraction/Type = "ForcedInteraction"
sv:Vr/ForcedInteraction/forregion/targetregion/processesNamed = 2 "compt" "conv"
dv:Vr/ForcedInteraction/ForRegion/targetregion/ForcedDistances = 2 Tf/ActiveVRT/Value Tf/ActiveVRT/Value mm
b:Vr/ForcedInteraction/ForRegion/targetregion/CorrectByWeight = "True"