67 G4cout <<
" magneticMoment: " << magneticMoment/muN <<
G4endl;
71 G4double g_factor = 2 * magneticMoment/muN;
77 G4double anomaly = (g_factor - 2.)/2.;
80 anomaly = (g_factor * (mass_neutron/mass_proton) - 2.)/2.;
84 g_factor = 2 * magneticMoment/muN;
86 anomaly = (g_factor - 2.)/2.;
98 G4cout <<
" magneticMoment: " << magneticMoment/muB <<
G4endl;
102 g_factor = magneticMoment/muB;
108 anomaly = (g_factor - 2.)/2.;
140 if (phi >
pi/2 && phi <
pi) phi =
pi-
phi;
142 G4double z = std::sin(phi)*std::cos(theta);
143 G4double x = std::sin(phi)*std::sin(theta);