## Meeting 2015 10 30 » History » Version 5

Nicolas Ponieman, 10/26/2015 04:07 PM

 1 2 3 1 Nicolas Ponieman `{{toc}}` 1 Nicolas Ponieman 1 Nicolas Ponieman `h1. Meeting 2015 10 30` 1 Nicolas Ponieman 1 Nicolas Ponieman `h2. Objectives Proposed during last meeting` 1 Nicolas Ponieman 4 Nicolas Ponieman `# (P0) (DONE) Check the values and precision of the parameters involved in Geant4, QSS and analytical solution are exactly the same. ` 4 Nicolas Ponieman `# (P0) (DONE) Plot the difference between the simulated velocity and the expected value (which is a parameter in the model).` 2 Nicolas Ponieman `# (P1) For each simulation with a certain number of planes, calculate the error of the simulation.` 2 Nicolas Ponieman `# (P3) Improve last meeting plots which compare Geant4, QSS3 and DOPRI5 performance.` 1 Nicolas Ponieman 1 Nicolas Ponieman ` * Format issues (scales, line colors, etc.)` 1 Nicolas Ponieman ` * Make statistic on many (10?) simulations for the points and add {simulationTime, simulationTimeError} for each number of planes.` 1 Nicolas Ponieman ` * How does DOPRI5 find the intersection with the planes?` 2 Nicolas Ponieman `# (P1) New Experiment: varying stepMax (from 0.1mm to 100mm) calculate the error in the first steps (in the position and in the velocity, also the error which is compared to epsilon in Geant4)` 2 Nicolas Ponieman `# (P1) Repeat Experiment 1: debug what is the error in Geant4 stepper, which is to be compared with epsilon.` 2 Nicolas Ponieman `# (P2) New Experiment: study impact of deltaInt parameter in Geant4 when changing the number of planes.` 2 Nicolas Ponieman `# (P3) Make some experiments with a different magnetic field:` 1 Nicolas Ponieman `B(x,y) = (B/R) sqrt( x*x + (y+R)*(y+R)) zhat (where B= constant, zhat is the directional vector along the z axis, R is the radius of the circle)` 1 Nicolas Ponieman `This field will keep the problem same as it is, but will give a deviation if the particle is off the circle` 1 Nicolas Ponieman `(due to the error of the stepper).` 2 Nicolas Ponieman `# (P3) Introduce discrete stochastic events simulating an artificial decay.` 1 Nicolas Ponieman `# Introduce a different graphical representation (such as changing the color of the scatter points) for each different type of errors.` 1 Nicolas Ponieman ` ` 1 Nicolas Ponieman ` * What are the different possible kind of errors?` 1 Nicolas Ponieman `# (P1) Summarize all the information Nico learnt with Soon about the transportation chain in a slide, or a handwritten sketch.` 1 Nicolas Ponieman `# (P1) Group presentation on Thursday by Nico.` 4 Nicolas Ponieman 4 Nicolas Ponieman `h2. Constant Values` 4 Nicolas Ponieman 4 Nicolas Ponieman `After last meeting, we found a few inconsistencies in the constants we were using as parameters in our external scripts and the ones used internally by Geant4.` 4 Nicolas Ponieman 4 Nicolas Ponieman `* c is consistent` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** In Line 54: /opt/geant4/geant4.10.01.p01/source/externals/clhep/include/CLHEP/Units/PhysicalConstants.` 4 Nicolas Ponieman `static const double c_light = 2.99792458e+8 * m/s;` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** We used:` 4 Nicolas Ponieman `c = 299.792458` 4 Nicolas Ponieman 4 Nicolas Ponieman `* gamma is inconsistent` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** We used:` 4 Nicolas Ponieman `gamma = 22.366272` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** which is inconsistent with the value` 4 Nicolas Ponieman `gamma = np.sqrt(1/(1 - (v0/c)**2)) = 22.3662720421` 4 Nicolas Ponieman 4 Nicolas Ponieman `* Electron mass (M) is inconsistent` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** In Line 85 /opt/geant4/geant4.10.01.p01/source/externals/clhep/include/CLHEP/Units/PhysicalConstants.h` 4 Nicolas Ponieman `static const double electron_mass_c2 = 0.510998910 * MeV;` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** We used:` 4 Nicolas Ponieman `old_M = 0.5109989` 4 Nicolas Ponieman 4 Nicolas Ponieman `* Kinetic Energy is inconsistent:` 4 Nicolas Ponieman ` ** In Geant4, we were using a value for Kinetic Energt which was different to the one that can be calculated either with the old or new values of the constants:` 4 Nicolas Ponieman `/gun/energy 10.9181415106 MeV (this 'wrong' value was taken from the original wiki entry).` 4 Nicolas Ponieman 4 Nicolas Ponieman ` ** The new and correct value for Kinetic Energy should be:` 4 Nicolas Ponieman `/gun/energy 10.9181417243 MeV` 4 Nicolas Ponieman 4 Nicolas Ponieman `* A few experiments were reran and no difference in the plots nor in the errors was noticed.` 4 Nicolas Ponieman `This seems reasonable considering this corrections are of very low order, and our errors are larger.` 4 Nicolas Ponieman 4 Nicolas Ponieman `* As a conclusion, this is an error we have to correct but which doesn't modify our previous results and leaves open the same questions as before regarding the invariance of error when changing epsilon.` 4 Nicolas Ponieman 4 Nicolas Ponieman `h2. Velocity Error Plots` 5 Nicolas Ponieman 5 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/29964/VError_vs_dist.png!` 5 Nicolas Ponieman `Figure 1: VError (difference between simulated and expected value) vs distance` 5 Nicolas Ponieman 5 Nicolas Ponieman `!https://cdcvs.fnal.gov/redmine/attachments/download/29963/VError_hist.png!` 5 Nicolas Ponieman `Figure 2: VError (difference between simulated and expected value) histogram` 5 Nicolas Ponieman 5 Nicolas Ponieman `* We notice the simulated velocity is *always* smaller than the expected velocity. ` 5 Nicolas Ponieman `* Why is this happening?`