#
i:Ts/NumberOfThreads = 2

# Secondary biasing occurs in the region
# of this component
s:Ge/target/Parent   = "World"
s:Ge/target/Material = "Lead"
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 cm
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             = "e-"
d:So/Example/BeamEnergy               = 6 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"

#Secondary biasing
b:Vr/UseVarianceReduction = "True"
b:Vr/BremsstrahlungSplit/Active = "True"
s:Vr/BremsstrahlungSplit/Type   = "secondarybiasing"
# The secondary particles produced in this process will
# be split  100 times if their energies are below of 6 MeV
sv:Vr/ParTicleSplit1/forregion/targetregion/processesNamed   = 1 "eBrem"
uv:Vr/BremsstrahlungSplit/ForRegion/targetregion/SplitNumber      = 1   100.
dv:Vr/BremsstrahlungSplit/ForRegion/targetregion/MaximumEnergies  = 1   6  MeV

