Project

General

Profile

Meeting 2015 07 17

Theoretical radius: ~ 38 [mm]
Radius = (M/c^2) * gamma * v/ B = 0.5109989 * 22.366272 *0.999 * 1000 / 299.7922458 ~ 38.0854 [mm]

The relative error in each dimension should be: RelativeError = RelativeTolerance * Radius
For G4 we assume temporarily that RelativeTolerance = Epsilon (maybe we should convert it from 3D to 1D for each axis)
For PD we know that RelativeError = DeltaQrel

All PD simulations with QSS4 method.

Plots that can be useful for discussion during the meeting:

DeltaQrel = 1e-3

Simulation Time (with output): 497 ms
Number of points: 1550





DeltaQrel = 1e-4

Simulation Time (with output): 547 ms
Number of points: 3450





DeltaQrel = 1e-5

Simulation Time (with output): 610 ms
Number of points: 7600





Geant 4

Simulation Time (verbosity = 0, has some logging anyway): 1693 ms
Number of points: 100000 (trackMaxLen/stepMax = 10m / 1e-1 mm = 10^5, is this casual?)