Gabriele Cosmo
135 lines
4.3 KiB
Plaintext
135 lines
4.3 KiB
Plaintext
#
|
|
# Macro file for the test of field01 example
|
|
#
|
|
# It creates the default geometry (simple absorber cylinder )
|
|
#
|
|
/tracking/verbose 1
|
|
/run/verbose 1
|
|
/control/verbose 2
|
|
#/control/cout/ignoreThreadsExcept 0
|
|
#
|
|
# ****** Start of initialisation of classes for field propagation
|
|
# Change the type of stepper used to integrate the ODE of motion
|
|
#
|
|
# Recommended and default since Geant4 10.4 (Dec 2017) - Dormand Prince (7) 4/5
|
|
/field/stepperType DormandPrince745
|
|
|
|
# High efficiency 5th order embedded stepper - Tsitouras
|
|
# /field/stepperType TsitourasRK45
|
|
|
|
# Higher order Dormand Prince stepper - 6th order (for smoothest fields & long steps)
|
|
# /field/stepperType DormandPrinceRK56
|
|
#
|
|
#
|
|
# Different types of steppers can be chosen for integration in B-field:
|
|
#
|
|
# 745 - Dormand Prince 745 : established, efficient embedded method
|
|
# recommended in literature, e.g. Hairer, Numerical Recipes
|
|
# Used in several RK programs (e.g. DOPRI5 code)
|
|
#
|
|
# 45 - Bogacki Shampine 4/5 : more efficient 4/5 embedded stepper
|
|
# 13 - Nystrom RK4 : stepper with few field calls & analytic estimate of error
|
|
#
|
|
# Good choices for reasonably smooth fields - available since Geant4 1.0
|
|
# 8 - Cash Karp RKF 45 : 'embedded' RK method - 4th/5th faster, robust
|
|
# ( uses difference of 4th & 5th order for error estimate )
|
|
#
|
|
# Default - good choice for unknown fields
|
|
# 4 - ClassicalRK4 : original Runge-Kutta method, very robust but slower )
|
|
# ( obtains error estimate by doing 2 half steps )
|
|
#
|
|
# Good choices for non-smooth fields:
|
|
# 3 - SimpleHeum : low order, with error obtained from half-steps
|
|
# 23 - BogackiShampine23 : lower order embedded method (new in 10.3-beta)
|
|
#
|
|
#
|
|
# Control the mininum step size for integration in field
|
|
# with a value of 10 um ==> accept any amount of integration error in smaller steps!
|
|
#
|
|
/field/setMinimumStep 0.01 mm
|
|
#
|
|
# Control the (relative) accuracy of the integration.
|
|
# ---------------------------------------------------
|
|
# 1. DeltaOneStep:
|
|
## Errors of this size in an integration sub-step are acceptable except when
|
|
# it is limited by the relative integration error limits (epsilon_min/max)
|
|
#
|
|
/field/setDeltaOneStep 2.0e-05 mm
|
|
#
|
|
# DeltaOneStep is used to obtain a relative accuracy for the step:
|
|
# relative_eps = DeltaOneStep / step-length
|
|
# To ensure robustness the value of this is kept within the following limits:
|
|
# 2. Epsilon Min: the minimum value of the relative integration error (EpsilonMin) -
|
|
#
|
|
/field/setMinimumEpsilonStep 5.0e-06
|
|
#
|
|
# This Limit to ensure that large steps do NOT result in a very low value of relative_eps
|
|
# (ie very high accuracy) that would result in integration requiring many steps and CPU cycles.
|
|
# Note: it is a dimensionless number
|
|
#
|
|
# 3. Epsilon Max:
|
|
#
|
|
/field/setMaximumEpsilonStep 5.0e-04
|
|
#
|
|
# A minimum value of the relative integration error (EpsilonMax). Also a dimensionless number");
|
|
# This Limit serves to ensure that a very small step does NOT result in a large relative error
|
|
# (ie low accuracy) and unreliable results.");
|
|
#
|
|
|
|
#
|
|
# ****** Initialisation of classes for field propagation is done
|
|
#
|
|
/run/initialize
|
|
/field/printParameters
|
|
#
|
|
# /run/particle/dumpCutValues
|
|
#
|
|
/gun/particle e-
|
|
#
|
|
# /gun/particle proton
|
|
# /gun/particle chargedgeantino
|
|
#
|
|
/gun/energy 500.0 MeV
|
|
/run/beamOn 1
|
|
#
|
|
/gun/energy 250.0 MeV
|
|
/run/beamOn 1
|
|
#
|
|
/gun/energy 200.0 MeV
|
|
/run/beamOn 1
|
|
#
|
|
/gun/energy 100.0 MeV
|
|
/run/beamOn 1
|
|
#
|
|
/gun/energy 50.0 MeV
|
|
/run/beamOn 1
|
|
#
|
|
# Change the value of the B-field
|
|
#
|
|
/field01/setField 0.1 0 0 tesla
|
|
/run/beamOn 1
|
|
#
|
|
/gun/energy 500.0 MeV
|
|
/tracking/verbose 0
|
|
/run/printProgress 10
|
|
/run/beamOn 100
|
|
#
|
|
# Test commands defined in this example
|
|
/control/verbose 2
|
|
#
|
|
/calor/setAbsMat Xe20CO2 # default Air
|
|
/calor/setWorldMat Kr20CO2 # default Air
|
|
/calor/setAbsThick 0.8 mm # default 1 mm
|
|
/calor/setAbsRad 1900 mm # default 20000.*mm
|
|
/calor/setAbsZpos 20990 mm # default 21990.0*mm
|
|
/calor/setWorldZ 45000 mm # default 44000.*mm
|
|
/calor/setWorldR 23000 mm # default 22000.*mm
|
|
/field01/setField 0 0 1.0 kG # default field 3*tesla in Z-direction
|
|
#
|
|
/gun/random on # default "off"
|
|
/gun/xvertex 100 mm # default 0
|
|
/gun/yvertex 100 mm # default 0
|
|
/gun/zvertex 100 mm # default 0
|
|
#
|
|
/run/beamOn 1
|