source: trunk/source/processes/electromagnetic/lowenergy/src/G4DNAGenericIonsManager.cc @ 1315

Last change on this file since 1315 was 1315, checked in by garnier, 14 years ago

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File size: 7.4 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26// $Id: G4DNAGenericIonsManager.cc,v 1.6 2009/06/10 13:32:36 mantero Exp $
27// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
28
29#include "G4DNAGenericIonsManager.hh"
30#include "G4Alpha.hh"
31#include "G4Ions.hh"
32
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
34
35G4DNAGenericIonsManager * G4DNAGenericIonsManager :: Instance(void)
36{
37 if (!theInstance)
38  theInstance=new G4DNAGenericIonsManager;
39 
40 return theInstance;
41}
42
43//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
44
45G4ParticleDefinition * G4DNAGenericIonsManager :: GetIon(const G4String & name)
46{
47 IonsMap::const_iterator i(map.find(name));
48 
49 if (i==map.end())
50  return 0;
51 
52 return i->second;
53}
54
55//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
56
57G4DNAGenericIonsManager :: G4DNAGenericIonsManager()
58{
59  //               name             mass          width         charge
60  //             2*spin           parity  C-conjugation
61  //          2*Isospin       2*Isospin3       G-parity
62  //               type    lepton number  baryon number   PDG encoding
63  //             stable         lifetime    decay table
64  //         shortlived          subType  anti_encoding
65  //         excitation   
66 
67 G4Ions *helium;
68 G4Ions *hydrogen;
69 G4Ions *alphaPlus;
70 G4Ions *positronium1s;
71 G4Ions *positronium2s;
72 
73 
74 helium=     new G4Ions(
75                        "helium",    3.727417*GeV,       0.0*MeV,  +0.0*eplus,
76                        0,              +1,             0,
77                        0,               0,             0,
78                        "nucleus",              +2,            +4,           0,
79                        true,                -1.0,             0,       
80                        false,                "",               0,             
81                        0.0);
82
83 alphaPlus=  new G4Ions("alpha+",    3.727417*GeV,       0.0*MeV,  +1.0*eplus,
84                               1,              +1,             0,
85                               0,               0,             0,
86                       "nucleus",              +1,            +4,           0,
87                        true,            -1.0,             0, false,
88                              "",               0,             0.0);
89
90 hydrogen= new G4Ions("hydrogen",   0.9382723*GeV,       0.0*MeV,  +0.0*eplus,
91                               0,              +1,             0,
92                               0,               0,             0,
93                       "nucleus",              +1,            +1,           0,
94                        true,            -1.0,             0, false,
95                              "",               0,             0.0);
96
97 positronium1s= new G4Ions("Ps-1s",   2*electron_mass_c2,      0.0*MeV,  +0.0*eplus,
98                               0,               0,             0,
99                               0,               0,             0,
100                       "nucleus",               0,             0,           0,
101                            true,            -1.0,             0, false,
102                              "",               0,             0.0);
103
104 positronium2s= new G4Ions("Ps-2s",   2*electron_mass_c2,      0.0*MeV,  +0.0*eplus,
105                               0,               0,             0,
106                               0,               0,             0,
107                       "nucleus",               0,             0,           0,
108                            true,            -1.0,             0, false,
109                              "",               0,             0.0);
110
111
112 /*
113 // molechules construction
114
115 G4Ions* oxonium; // H3O -- it will become H3O+
116 G4Ions* hydroxyl; // OH -- it will produce OH- too
117 G4Ions* molHydrogen; // H2
118 //G4Ions* hydroxide; // OH-
119 G4Ions* hydroPeroxide; // H2O2
120 G4Ions* water; // H2O -- it will become also H2O+
121
122
123 G4double mass = 19.02*g/Avogadro - 11*electron_mass_c2;
124
125 oxonium = new G4Ions("H3O",        mass,             0,  +11.0*eplus,
126                      0,               0,             0,
127                      0,               0,             0,
128                      "molecule",      0,             0,           0,
129                      true,         -1.0,             0,
130                      false,          "",             0,             
131                      0.0);
132 
133 mass = 17.00734*g/Avogadro - 9*electron_mass_c2;
134 
135 hydroxyl = new G4Ions("OH",        mass,             0,  +9.0*eplus,
136                      0,               0,             0,
137                      0,               0,             0,
138                      "molecule",      0,             0,           0,
139                      true,         -1.0,             0,
140                      false,          "",             0,             
141                      0.0);
142
143 mass = 2.01588*g/Avogadro - 2*electron_mass_c2;
144
145 molHydrogen = new G4Ions("H2",         mass,             0,  +2.0*eplus,
146                          0,               0,             0,
147                          0,               0,             0,
148                          "molecule",      0,             0,           0,
149                          true,         -1.0,             0,
150                          false,          "",             0,             
151                          0.0);
152
153 mass = 34.01468*g/Avogadro - 18*electron_mass_c2;
154
155 hydroPeroxide = new G4Ions("H2O2",      mass,             0, +18.0*eplus,
156                           0,               0,             0,
157                           0,               0,             0,
158                           "molecule",      0,             0,           0,
159                           true,         -1.0,             0,
160                           false,          "",             0,             
161                           0.0);
162 
163 mass = 18.015*g/Avogadro - 10*electron_mass_c2;
164
165 water = new G4Ions("H2O",        mass,             0,  +10.0*eplus,
166                    0,               0,             0,
167                    0,               0,             0,
168                    "molecule",      0,             0,           0,
169                    true,         -1.0,             0,
170                    false,          "",             0,             
171                    0.0);
172
173 map["H3O" ]  =oxonium;
174 map["OH"  ]  =hydroxyl;
175 map["H2"  ]  =molHydrogen;
176 map["H2O2"]  =hydroPeroxide;
177 map["H2O" ]  =water;
178 */
179
180
181 map["helium"  ]=helium;
182 map["hydrogen"]=hydrogen;
183 map["alpha+"  ]=alphaPlus;
184 map["alpha++" ]=G4Alpha::Alpha();
185 map["Ps-1s"   ]=positronium1s;
186 map["Ps-2s"   ]=positronium2s;
187
188}
189
190//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
191
192G4DNAGenericIonsManager * G4DNAGenericIonsManager::theInstance(0);
193   
Note: See TracBrowser for help on using the repository browser.