Project

General

Profile

Meeting 2015 09 25 » History » Version 1

Version 1/6 - Next » - Current version
Nicolas Ponieman, 09/21/2015 05:13 PM


Meeting 2015 09 25

Generic settings

/N02/det/trackMaxSteps 1000000000
/N02/det/trackMaxLen 10 km
/gun/energy 10.9181415106 MeV
/N02/prec/setDeltaOne 1.0e-2 mm
/N02/prec/setDeltaInt 1.0e-5 mm
/N02/prec/setEpsilonMin 1.0e-3
/N02/prec/setEpsilonMax 1.0e-3
/N02/det/stepMax 20. mm

  • All measured times correspond to executions which don't produce output at all (or produce a minimum logging).
  • MSE is the Mean Squared Error of the position
  • MSE = 1/N * sum_{i=1}^{N} (G4x_i - ANx_i)^2 + (G4y_i - ANy_i)^2
  • N = Number of output points
  • r_Error = sqrt(MSE)
  • max(X_Error) is the maximum of all the absolute errors in x, i.e. max(abs(G4x_i - ANx_i))

Experiment 4 - Vary DeltaChord in G4

1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm
Delta Chord [mm] Simulation Time RHS evaluation steps r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] r_Error [mm] Max(X_Error) [mm] Output Steps Output steps / RHS steps ratio
0.01 3.9s 6.60e6 50001 7.58e-3
0.05 2.3s 3.30e6 50001 1.51e-2
0.25 1.5s 1.65e6 50001 3.03e-2
0.5 1.3s 1.10e6 50001 4.55e-2
1.0 1.3s 1.10e6 50001 4.55e-2

Although all the errors look similar, there is a clear and counterintuitive tendency for the error to become smaller as stepMax is increased.
A few plots for different values of stepMax: