Project

General

Profile

Meeting 2015 09 25 » History » Version 3

Nicolas Ponieman, 09/21/2015 06:18 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 2 Nicolas Ponieman
|\3=. |\2=. 1 km |\2=. 100 m |\2=. 10 m |\2=. 1 m |\2=. 100 mm |\2=. |
27 2 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] |_. Output Steps |_. Output steps / RHS steps ratio |
28 2 Nicolas Ponieman
| 0.01   | 3.9s | 6.60e6 | 1.74E+00 | 3.13E+00 | 1.74E-01 | 3.09E-01 | 1.74E-02 | 2.98E-02 | 1.74E-03 | 2.80E-03 | 1.30E-04 | 2.40E-04 | 50001 | 7.58e-3 |
29 2 Nicolas Ponieman
| 0.05   | 2.3s | 3.30e6 |  1.74E+00 | 3.13E+00 | 1.74E-01 | 3.09E-01 | 1.74E-02 | 2.97E-02 | 1.74E-03 | 2.79E-03 | 1.30E-04 | 2.40E-04 | 50001 | 1.51e-2 |
30 2 Nicolas Ponieman
| 0.25   | 1.5s | 1.65e6 | 1.73E+00 | 3.12E+00 | 1.73E-01 | 3.07E-01 | 1.73E-02 | 2.96E-02 | 1.72E-03 | 2.78E-03 | 1.26E-04 | 2.30E-04 | 50001 | 3.03e-2 |
31 2 Nicolas Ponieman
| 0.5    | 1.3s | 1.10e6 | 1.70E+00 | 3.06E+00 | 1.70E-01 | 3.02E-01 | 1.70E-02 | 2.91E-02 | 1.68E-03 | 2.73E-03 | 1.14E-04 | 2.10E-04 | 50001 | 4.55e-2 |
32 2 Nicolas Ponieman
| 1.0    | 1.3s | 1.10e6 | 1.60E+00 | 2.88E+00 | 1.60E-01 | 2.84E-01 | 1.59E-02 | 2.73E-02 | 1.52E-03 | 2.55E-03 | 5.27E-05 | 7.00E-05  |50001 | 4.55e-2 |
33 1 Nicolas Ponieman
34 1 Nicolas Ponieman
35 2 Nicolas Ponieman
There is a *clear tendence* for the error to become smaller as deltaChord is increased.
36 2 Nicolas Ponieman
A few plots of the error for different values of deltaChord:
37 3 Nicolas Ponieman
!https://cdcvs.fnal.gov/redmine/attachments/download/29100/G4_vs_AN_error_no_cross_deltaChord0.01mm.png!
38 3 Nicolas Ponieman
!https://cdcvs.fnal.gov/redmine/attachments/download/29098/G4_vs_AN_error_no_cross_deltaChord0.25mm.png!
39 3 Nicolas Ponieman
!https://cdcvs.fnal.gov/redmine/attachments/download/29099/G4_vs_AN_error_no_cross_deltaChord1.0mm.png!