| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // -------------------------------------------
27 // GEANT 4 class implementation file
28 //
29 // History: first implementation, based o
30 // 4th April 1996, G.Cosmo
31 // H.Kurashige 7 July
32 // add neutron life time Oct 17 2000
33 // *******************************************
34 // New impelemenataion as an utility class M
35 // -------------------------------------------
36
37 #include "G4Neutron.hh"
38
39 #include "G4DecayTable.hh"
40 #include "G4NeutronBetaDecayChannel.hh"
41 #include "G4ParticleTable.hh"
42 #include "G4PhysicalConstants.hh"
43 #include "G4String.hh"
44 #include "G4SystemOfUnits.hh"
45 #include "G4Types.hh"
46 #include "G4VDecayChannel.hh"
47
48 G4Neutron* G4Neutron::theInstance = nullptr;
49
50 G4Neutron* G4Neutron::Definition()
51 {
52 if (theInstance != nullptr) return theInstan
53 const G4String name = "neutron";
54 // search in particle table]
55 G4ParticleTable* pTable = G4ParticleTable::G
56 auto anInstance = static_cast<G4Ions*>(pTabl
57 if (anInstance == nullptr) {
58 // create particle
59 //
60 // Arguments for constructor are as fol
61 // name mass
62 // 2*spin parity C-
63 // 2*Isospin 2*Isospin3
64 // type lepton number ba
65 // stable lifetime
66 // shortlived subType
67 // use constants in CLHEP
68 // static const double neutron_mass_c2 =
69
70 // clang-format off
71 anInstance = new G4Ions(
72 name, neutron_mass_c2, 7.478e
73 1, +1, 0,
74 1, -1, 0,
75 "baryon", 0, +
76 false, 880.2*second, nullptr,
77 false, "nucleon", -2112,
78 0.0, 0
79 );
80 // clang-format on
81
82 // Magnetic Moment
83 G4double mN = eplus * hbar_Planck / 2. / (
84 anInstance->SetPDGMagneticMoment(-1.913042
85 // create Decay Table
86 auto table = new G4DecayTable();
87 // create a decay channel
88 G4VDecayChannel* mode = new G4NeutronBetaD
89 table->Insert(mode);
90 anInstance->SetDecayTable(table);
91 }
92 theInstance = static_cast<G4Neutron*>(anInst
93 return theInstance;
94 }
95
96 G4Neutron* G4Neutron::NeutronDefinition()
97 {
98 return Definition();
99 }
100
101 G4Neutron* G4Neutron::Neutron()
102 {
103 return Definition();
104 }
105