Project

General

Profile

Meeting 2015 09 25 » History » Version 1

Nicolas Ponieman, 09/21/2015 05:13 PM

1 1 Nicolas Ponieman
{{toc}}
2 1 Nicolas Ponieman
3 1 Nicolas Ponieman
h1. Meeting 2015 09 25
4 1 Nicolas Ponieman
5 1 Nicolas Ponieman
h2. Generic settings
6 1 Nicolas Ponieman
7 1 Nicolas Ponieman
/N02/det/trackMaxSteps 1000000000
8 1 Nicolas Ponieman
/N02/det/trackMaxLen   10 km
9 1 Nicolas Ponieman
/gun/energy 10.9181415106 MeV
10 1 Nicolas Ponieman
/N02/prec/setDeltaOne   1.0e-2 mm  
11 1 Nicolas Ponieman
/N02/prec/setDeltaInt   1.0e-5 mm
12 1 Nicolas Ponieman
/N02/prec/setEpsilonMin 1.0e-3
13 1 Nicolas Ponieman
/N02/prec/setEpsilonMax 1.0e-3
14 1 Nicolas Ponieman
/N02/det/stepMax 20. mm
15 1 Nicolas Ponieman
16 1 Nicolas Ponieman
* All measured times correspond to executions which don't produce output at all (or produce a minimum logging). 
17 1 Nicolas Ponieman
* MSE is the Mean Squared Error of the position
18 1 Nicolas Ponieman
* MSE = 1/N * sum_{i=1}^{N} (G4x_i - ANx_i)^2 + (G4y_i - ANy_i)^2
19 1 Nicolas Ponieman
* N = Number of output points
20 1 Nicolas Ponieman
* r_Error = sqrt(MSE)
21 1 Nicolas Ponieman
* max(X_Error) is the maximum of all the absolute errors in x, i.e. max(abs(G4x_i - ANx_i))
22 1 Nicolas Ponieman
23 1 Nicolas Ponieman
h2. Experiment 4 - Vary DeltaChord in G4
24 1 Nicolas Ponieman
25 1 Nicolas Ponieman
26 1 Nicolas Ponieman
|\3=. |\2=. 1 km |\2=. 100 m |\2=. 10 m |\2=. 1 m |\2=. 100 mm |\2=. 10 mm |\2=. 1 mm |\2=. |
27 1 Nicolas Ponieman
|_. 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 |
28 1 Nicolas Ponieman
| 0.01   | 3.9s | 6.60e6 |  |  |  |  |  |  |  |  |  |  |  |  |  |  | 50001 | 7.58e-3 |
29 1 Nicolas Ponieman
| 0.05   | 2.3s | 3.30e6 |  |  |  |  |  |  |  |  |  |  |  |  |  |  | 50001 | 1.51e-2 |
30 1 Nicolas Ponieman
| 0.25   | 1.5s | 1.65e6 |  |  |  |  |  |  |  |  |  |  |  |  |  |  | 50001 | 3.03e-2 |
31 1 Nicolas Ponieman
| 0.5    | 1.3s | 1.10e6 |  |  |  |  |  |  |  |  |  |  |  |  |  |  | 50001 | 4.55e-2 |
32 1 Nicolas Ponieman
| 1.0    | 1.3s | 1.10e6 |  |  |  |  |  |  |  |  |  |  |  |  |  |  | 50001 | 4.55e-2 |
33 1 Nicolas Ponieman
34 1 Nicolas Ponieman
35 1 Nicolas Ponieman
Although all the errors look similar, there is a clear and counterintuitive tendency for the error to become smaller as stepMax is increased. 
36 1 Nicolas Ponieman
A few plots for different values of stepMax: