Meeting 2015 10 30 » History » Version 4

« Previous - Version 4/10 (diff) - Next » - Current version
Nicolas Ponieman, 10/26/2015 04:05 PM

Meeting 2015 10 30

Objectives Proposed during last meeting

  1. (P0) (DONE) Check the values and precision of the parameters involved in Geant4, QSS and analytical solution are exactly the same.
  2. (P0) (DONE) Plot the difference between the simulated velocity and the expected value (which is a parameter in the model).
  3. (P1) For each simulation with a certain number of planes, calculate the error of the simulation.
  4. (P3) Improve last meeting plots which compare Geant4, QSS3 and DOPRI5 performance.
    • Format issues (scales, line colors, etc.) * Make statistic on many (10?) simulations for the points and add {simulationTime, simulationTimeError} for each number of planes. * How does DOPRI5 find the intersection with the planes?
  5. (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)
  6. (P1) Repeat Experiment 1: debug what is the error in Geant4 stepper, which is to be compared with epsilon.
  7. (P2) New Experiment: study impact of deltaInt parameter in Geant4 when changing the number of planes.
  8. (P3) Make some experiments with a different magnetic field:
    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)
    This field will keep the problem same as it is, but will give a deviation if the particle is off the circle
    (due to the error of the stepper).
  9. (P3) Introduce discrete stochastic events simulating an artificial decay.
  10. Introduce a different graphical representation (such as changing the color of the scatter points) for each different type of errors.
    • What are the different possible kind of errors?
  11. (P1) Summarize all the information Nico learnt with Soon about the transportation chain in a slide, or a handwritten sketch.
  12. (P1) Group presentation on Thursday by Nico.

Constant Values

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.

  • c is consistent
    • In Line 54: /opt/geant4/geant4.10.01.p01/source/externals/clhep/include/CLHEP/Units/PhysicalConstants.
      static const double c_light = 2.99792458e+8 * m/s;
    • We used:
      c = 299.792458
  • gamma is inconsistent
    • We used:
      gamma = 22.366272
    • which is inconsistent with the value
      gamma = np.sqrt(1/(1 - (v0/c)**2)) = 22.3662720421
  • Electron mass (M) is inconsistent
    • In Line 85 /opt/geant4/geant4.10.01.p01/source/externals/clhep/include/CLHEP/Units/PhysicalConstants.h
      static const double electron_mass_c2 = 0.510998910 * MeV;
    • We used:
      old_M = 0.5109989
  • Kinetic Energy is inconsistent: ** 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:
    /gun/energy 10.9181415106 MeV (this 'wrong' value was taken from the original wiki entry).
    • The new and correct value for Kinetic Energy should be:
      /gun/energy 10.9181417243 MeV
  • A few experiments were reran and no difference in the plots nor in the errors was noticed.
    This seems reasonable considering this corrections are of very low order, and our errors are larger.
  • 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.

Velocity Error Plots