49{
50 G4int n_particle = 1;
52
54
56 G4double mass_proton = particle->GetPDGMass();
57
58 G4double muN = 0.5*CLHEP::eplus*CLHEP::hbar_Planck/
59 (mass_proton/CLHEP::c_squared);
60
61 G4cout << " *** Neutron *** " << G4endl;
62
63 particle = particleTable->FindParticle("neutron");
64 G4double mass_neutron = particle->GetPDGMass();
65
66 G4double magneticMoment = particle->GetPDGMagneticMoment();
67 G4cout << " magneticMoment: " << magneticMoment/muN << G4endl;
68
69
70
71 G4double g_factor = 2 * magneticMoment/muN;
72 G4cout << " g_factor: " << g_factor << G4endl;
73
74 G4double charge = particle->GetPDGCharge();
75 G4cout << " charge: " << charge << G4endl;
76
77 G4double anomaly = (g_factor - 2.)/2.;
78 G4cout << " anomaly: " << anomaly << G4endl;
79
80 anomaly = (g_factor * (mass_neutron/mass_proton) - 2.)/2.;
81 G4cout << " corrected anomaly: " << anomaly << G4endl;
82
83 muN = 0.5*CLHEP::eplus*CLHEP::hbar_Planck/(mass_neutron/CLHEP::c_squared);
84 g_factor = 2 * magneticMoment/muN;
85
86 anomaly = (g_factor - 2.)/2.;
87 G4cout << " *** anomaly: " << anomaly << G4endl;
88
89 G4cout << " *** MuonPlus *** " << G4endl;
90
91 particle = particleTable->FindParticle("mu+");
92 G4double mass_muon = particle->GetPDGMass();
93
94 G4double muB = 0.5*CLHEP::eplus*CLHEP::hbar_Planck/
95 (mass_muon/CLHEP::c_squared);
96
97 magneticMoment = particle->GetPDGMagneticMoment();
98 G4cout << " magneticMoment: " << magneticMoment/muB << G4endl;
99
100
101
102 g_factor = magneticMoment/muB;
103 G4cout << " g_factor: " << g_factor << G4endl;
104
105 charge = particle->GetPDGCharge();
106 G4cout << " charge: " << charge << G4endl;
107
108 anomaly = (g_factor - 2.)/2.;
109 G4cout << " anomaly: " << anomaly << G4endl;
110
111
112
113
114 particle = particleTable->FindParticle("neutron");
116}
G4ParticleGun * fParticleGun