## Meeting 2015 09 04 » History » Version 2

Nicolas Ponieman, 09/02/2015 09:16 AM

 1 2 3 1 Nicolas Ponieman `{{toc}}` 1 Nicolas Ponieman 1 Nicolas Ponieman `h1. Meeting 2015 09 04` 1 Nicolas Ponieman 1 Nicolas Ponieman `h2. Generic settings` 1 Nicolas Ponieman 1 Nicolas Ponieman `/N02/det/trackMaxSteps 1000000000` 1 Nicolas Ponieman `/N02/det/trackMaxLen 10 km` 1 Nicolas Ponieman `/gun/energy 10.9181415106 MeV` 1 Nicolas Ponieman `/N02/prec/setDeltaOne 1.0e-2 mm ` 1 Nicolas Ponieman `/N02/prec/setDeltaInt 1.0e-5 mm` 1 Nicolas Ponieman `/N02/prec/setEpsilonMin 1.0e-3` 1 Nicolas Ponieman `/N02/prec/setEpsilonMax 1.0e-3` 1 Nicolas Ponieman `/N02/prec/setDeltaChord 0.25 mm` 1 Nicolas Ponieman 1 Nicolas Ponieman `MSE is the Mean Squared Error` 1 Nicolas Ponieman `MSE = 1/N * sum_{i=1}^{n} (G4_i - AN_i)^2` 1 Nicolas Ponieman `N = Number of output points` 2 Nicolas Ponieman `r_err` 1 Nicolas Ponieman 1 Nicolas Ponieman `h2. Experiment 1 - Vary epsilon in G4` 1 Nicolas Ponieman 1 Nicolas Ponieman `(see [[Meeting_2015_08_21|August 21 meeting]])` 1 Nicolas Ponieman 1 Nicolas Ponieman `/N02/det/stepMax 20. mm (i.e., stepMax = 20mm)` 1 Nicolas Ponieman `epsilonMax = epsilonMin = epsilon` 1 Nicolas Ponieman `The mean square error doesn't seem to vary.` 1 Nicolas Ponieman `The error, however, keeps growing with distance.` 1 Nicolas Ponieman 1 Nicolas Ponieman `|_. epsilon |_. Simulation Time |_. RHS evaluation steps |_. Mean Square Error [mm^2] |_. sqrt(MSE) [mm] |_. Max(r_err) [mm] |` 1 Nicolas Ponieman `| 1.0e-3 | 0m13.725s | 16500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-4 | 0m13.725s | 16500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-5 | 0m13.725s | 16500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-6 | 0m21.815s | 38500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-7 | 0m27.980s | 60500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-8 | 0m41.054s | 104499999 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-9 | 0m55.830s | 159500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-10 | 1m25.948s | 269499999 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-11 | 2m17.996s | 456500010 | 294 | 17.14 | |` 1 Nicolas Ponieman `| 1.0e-12 | 3m50.224s | 786518324 | 294 | 17.14 | |` 1 Nicolas Ponieman 1 Nicolas Ponieman `A few plots for different values of epsilon:` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28621/G4_vs_AN_error_eps1e-3.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28620/G4_vs_AN_error_eps1e-6.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28619/G4_vs_AN_error_eps1e-9.png!` 1 Nicolas Ponieman 1 Nicolas Ponieman `The 3 plots look exactly the same!` 1 Nicolas Ponieman `The output for those epsilon values is different, though.` 1 Nicolas Ponieman 1 Nicolas Ponieman `h2. Experiment 2` 1 Nicolas Ponieman 1 Nicolas Ponieman `/N02/prec/setEpsilonMin 1.0e-3 (i.e, epsilon = 1.0e-3)` 1 Nicolas Ponieman `/N02/prec/setEpsilonMax 1.0e-3` 1 Nicolas Ponieman `Mean Square Error measured for first 15000 points (final time differs a lot, maybe we should consider some other measurment)` 1 Nicolas Ponieman 1 Nicolas Ponieman `|_. stepMax [mm] |_. Simulation Time |_. RHS evaluation steps |_. Mean Square Error [mm^2] |_. sqrt(MSE)[mm] |_. Max(r_err) [mm] |` 1 Nicolas Ponieman `| 0.2 | 13m34.816s | 550000011 | 0.00002 | | |` 1 Nicolas Ponieman `| 1 | 2m51.026s | 110000000 | 0.0006 | | |` 1 Nicolas Ponieman `| 2 | 1m28.192s | 55000000 | 0.002 | | |` 1 Nicolas Ponieman `| 5 | 0m33.121s | 22000000 | 0.017 | | |` 1 Nicolas Ponieman `| 10 | 0m21.907s | 22000010 | 0.068 | | |` 1 Nicolas Ponieman `| 20 | 0m13.797s | 16500010 | 0.27 | | |` 1 Nicolas Ponieman `| 50 | 0m8.319s | 13200010 | 1.65 | | |` 1 Nicolas Ponieman `| 100 | 0m7.179s | 13200020 | 6.34 | | |` 1 Nicolas Ponieman `| 150 | 0m6.848s | 13200020 | 13.68 | | |` 1 Nicolas Ponieman `| 200 | 0m6.382s | 12650020 | 23.31 | | |` 1 Nicolas Ponieman 1 Nicolas Ponieman `h2. Experiment 3 - Vary DeltaQ in PD` 1 Nicolas Ponieman 1 Nicolas Ponieman `Simulated with QSS3 method` 1 Nicolas Ponieman `DeltaQmin = DeltaQ * 1e-3` 1 Nicolas Ponieman `*Several values of the table were left blank deliberately as they weren't calculated yet*` 1 Nicolas Ponieman 1 Nicolas Ponieman `|_. DeltaQ |_. Simulation Time |_. Equivalente in PD to RHS eval steps |_. Mean Square Error [mm^2] |_. sqrt(MSE)[mm] |_. Max(r_err) [mm] |` 1 Nicolas Ponieman `| 1.0e-1 | 0.41s | 5.60e5 | | | |` 1 Nicolas Ponieman `| 1.0e-2 | 0.42s | 6.40e5 | | | |` 1 Nicolas Ponieman `| 1.0e-3 | 1.33s | 2.40e6 | 704.43 | 26.54 | |` 1 Nicolas Ponieman `| 1.0e-4 | 8.89s | 1.42e7 | 568.32 | 23.83 | |` 1 Nicolas Ponieman `| 1.0e-5 | 16.3s | 3.08e7 | 86.89 | 9.32 | |` 1 Nicolas Ponieman `| 1.0e-6 | 35.6s | 6.64e7 | 56.79 | 7.54 | |` 1 Nicolas Ponieman `| 1.0e-7 | 77.6s | 1.43e8 | 57.45 | 7.58 | |` 1 Nicolas Ponieman `| 1.0e-8 | 165s | 3.09e8 | | | |` 1 Nicolas Ponieman `| 1.0e-9 | 354s | 6.64e8 | | | |` 1 Nicolas Ponieman `| 1.0e-10 | 740s | 1.43e9 | | | |` 1 Nicolas Ponieman `| 1.0e-11 | 1630s | 3.09e9 | | | |` 1 Nicolas Ponieman `| 1.0e-12 | | | | | |` 1 Nicolas Ponieman 1 Nicolas Ponieman `A few plots for different values of DeltaQ:` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28625/PD_vs_AN_error_DeltaQ1e-3.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28624/PD_vs_AN_error_DeltaQ1e-4.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28623/PD_vs_AN_error_DeltaQ1e-5.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28622/PD_vs_AN_error_DeltaQ1e-6.png!` 1 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/28626/PD_vs_AN_error_DeltaQ1e-7.png!` 1 Nicolas Ponieman `*The last two plots look exactly the same*`