| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // G4DecayTableMessenger class implementation 26 // G4DecayTableMessenger class implementation
27 // 27 //
28 // Author: H.Kurashige, 13 June 1997 28 // Author: H.Kurashige, 13 June 1997
29 //-------------------------------------------- 29 //---------------------------------------------------------------------
30 30
31 #include "G4DecayTableMessenger.hh" 31 #include "G4DecayTableMessenger.hh"
32 << 32 #include "G4UImanager.hh"
>> 33 #include "G4UIdirectory.hh"
>> 34 #include "G4UIcmdWithoutParameter.hh"
>> 35 #include "G4UIcmdWithAnInteger.hh"
>> 36 #include "G4UIcmdWithADouble.hh"
>> 37 #include "G4VDecayChannel.hh"
33 #include "G4DecayTable.hh" 38 #include "G4DecayTable.hh"
34 #include "G4ParticleDefinition.hh" 39 #include "G4ParticleDefinition.hh"
35 #include "G4ParticleTable.hh" 40 #include "G4ParticleTable.hh"
36 #include "G4UIcmdWithADouble.hh" << 41 #include "G4ios.hh" // Include from 'system'
37 #include "G4UIcmdWithAnInteger.hh" << 42 #include <iomanip> // Include from 'system'
38 #include "G4UIcmdWithoutParameter.hh" <<
39 #include "G4UIdirectory.hh" <<
40 #include "G4UImanager.hh" <<
41 #include "G4VDecayChannel.hh" <<
42 #include "G4ios.hh" // Include from 'system' <<
43 43
44 #include <iomanip> // Include from 'system' << 44 G4DecayTableMessenger::G4DecayTableMessenger(G4ParticleTable* pTable)
45 << 45 : theParticleTable(pTable)
46 G4DecayTableMessenger::G4DecayTableMessenger(G <<
47 { 46 {
48 if (theParticleTable == nullptr) { << 47 if ( theParticleTable == nullptr )
>> 48 {
49 theParticleTable = G4ParticleTable::GetPar 49 theParticleTable = G4ParticleTable::GetParticleTable();
50 } 50 }
51 currentParticle = nullptr; 51 currentParticle = nullptr;
52 52
53 // Command /particle/property/decay/ 53 // Command /particle/property/decay/
54 thisDirectory = new G4UIdirectory("/particle 54 thisDirectory = new G4UIdirectory("/particle/property/decay/");
55 thisDirectory->SetGuidance("Decay Table cont 55 thisDirectory->SetGuidance("Decay Table control commands.");
56 56
57 // Command /particle/property/decay/select 57 // Command /particle/property/decay/select
58 selectCmd = new G4UIcmdWithAnInteger("/parti << 58 selectCmd = new G4UIcmdWithAnInteger("/particle/property/decay/select",this);
59 selectCmd->SetGuidance("Enter index of decay 59 selectCmd->SetGuidance("Enter index of decay mode.");
60 selectCmd->SetParameterName("mode", true); 60 selectCmd->SetParameterName("mode", true);
61 selectCmd->SetDefaultValue(0); 61 selectCmd->SetDefaultValue(0);
62 selectCmd->SetRange("mode >=0"); 62 selectCmd->SetRange("mode >=0");
63 currentChannel = nullptr; 63 currentChannel = nullptr;
64 64
65 // Command /particle/property/decay/dump 65 // Command /particle/property/decay/dump
66 dumpCmd = new G4UIcmdWithoutParameter("/part << 66 dumpCmd = new G4UIcmdWithoutParameter("/particle/property/decay/dump",this);
67 dumpCmd->SetGuidance("Dump decay mode inform 67 dumpCmd->SetGuidance("Dump decay mode information.");
68 68
69 // Command /particle/property/decay/br 69 // Command /particle/property/decay/br
70 brCmd = new G4UIcmdWithADouble("/particle/pr << 70 brCmd = new G4UIcmdWithADouble("/particle/property/decay/br",this);
71 brCmd->SetGuidance("Set branching ratio. [0< 71 brCmd->SetGuidance("Set branching ratio. [0< BR <1.0]");
72 brCmd->SetParameterName("br", false); << 72 brCmd->SetParameterName("br",false);
73 brCmd->SetRange("(br >=0.0) && (br <=1.0)"); 73 brCmd->SetRange("(br >=0.0) && (br <=1.0)");
74 } 74 }
75 75
76 G4DecayTableMessenger::~G4DecayTableMessenger( 76 G4DecayTableMessenger::~G4DecayTableMessenger()
77 { 77 {
78 delete dumpCmd; << 78 if (dumpCmd != nullptr) delete dumpCmd;
79 delete selectCmd; << 79 if (selectCmd != nullptr) delete selectCmd;
80 delete brCmd; << 80 if (brCmd != nullptr) delete brCmd;
81 delete thisDirectory; << 81 if (thisDirectory != nullptr) delete thisDirectory;
82 } 82 }
83 83
84 void G4DecayTableMessenger::SetNewValue(G4UIco 84 void G4DecayTableMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
85 { 85 {
86 if (SetCurrentParticle() == nullptr) { << 86 if (SetCurrentParticle()== nullptr)
>> 87 {
87 G4cout << "Particle is not selected yet !! 88 G4cout << "Particle is not selected yet !! Command ignored." << G4endl;
88 return; 89 return;
89 } 90 }
90 if (currentDecayTable == nullptr) { << 91 if (currentDecayTable== nullptr)
>> 92 {
91 G4cout << "The particle has no decay table 93 G4cout << "The particle has no decay table !! Command ignored." << G4endl;
92 return; 94 return;
93 } 95 }
94 96
95 if (command == dumpCmd) { << 97 if( command == dumpCmd )
>> 98 {
96 // Command /particle/property/decay/dump 99 // Command /particle/property/decay/dump
97 currentDecayTable->DumpInfo(); 100 currentDecayTable->DumpInfo();
>> 101
98 } 102 }
99 else if (command == selectCmd) { << 103 else if ( command == selectCmd )
>> 104 {
100 // Command /particle/property/decay/select 105 // Command /particle/property/decay/select
101 G4int index = selectCmd->GetNewIntValue(ne << 106 G4int index = selectCmd->GetNewIntValue(newValue) ;
102 currentChannel = currentDecayTable->GetDec 107 currentChannel = currentDecayTable->GetDecayChannel(index);
103 if (currentChannel == nullptr) { << 108 if ( currentChannel == nullptr )
>> 109 {
104 G4cout << "Invalid index. Command ignore 110 G4cout << "Invalid index. Command ignored." << G4endl;
105 } 111 }
106 else { << 112 else
>> 113 {
107 idxCurrentChannel = index; 114 idxCurrentChannel = index;
108 } 115 }
>> 116
109 } 117 }
110 else { << 118 else
111 if (currentChannel == nullptr) { << 119 {
>> 120 if ( currentChannel == nullptr )
>> 121 {
112 G4cout << "Select a decay channel. Comma 122 G4cout << "Select a decay channel. Command ignored." << G4endl;
113 return; 123 return;
114 } 124 }
115 if (command == brCmd) { << 125 if (command == brCmd)
>> 126 {
116 // Command /particle/property/decay/br 127 // Command /particle/property/decay/br
117 G4double br = brCmd->GetNewDoubleValue(n << 128 G4double br = brCmd->GetNewDoubleValue(newValue);
118 if ((br < 0.0) || (br > 1.0)) { << 129 if( (br<0.0) || (br>1.0) )
>> 130 {
119 G4cout << "Invalid brancing ratio. Com 131 G4cout << "Invalid brancing ratio. Command ignored." << G4endl;
120 } 132 }
121 else { << 133 else
>> 134 {
122 currentChannel->SetBR(br); 135 currentChannel->SetBR(br);
123 } 136 }
124 } 137 }
125 } 138 }
126 } 139 }
127 140
>> 141
128 G4ParticleDefinition* G4DecayTableMessenger::S 142 G4ParticleDefinition* G4DecayTableMessenger::SetCurrentParticle()
129 { 143 {
130 // set currentParticle pointer << 144 // set currentParticle pointer
131 // get particle name by asking G4ParticleMes 145 // get particle name by asking G4ParticleMessenger via UImanager
132 146
133 G4String particleName = G4UImanager::GetUIpo << 147 G4String particleName
>> 148 = G4UImanager::GetUIpointer()->GetCurrentStringValue("/particle/select");
134 149
135 if (currentParticle != nullptr) { << 150 if (currentParticle != nullptr )
136 // check whether selection is changed << 151 {
137 if (currentParticle->GetParticleName() != << 152 // check whether selection is changed
>> 153 if (currentParticle->GetParticleName() != particleName)
>> 154 {
138 currentParticle = theParticleTable->Find 155 currentParticle = theParticleTable->FindParticle(particleName);
139 idxCurrentChannel = -1; 156 idxCurrentChannel = -1;
140 currentDecayTable = nullptr; 157 currentDecayTable = nullptr;
141 } 158 }
142 else { << 159 else
143 // no change << 160 {
>> 161 // no change
144 return currentParticle; 162 return currentParticle;
145 } 163 }
>> 164
146 } 165 }
147 else { << 166 else
>> 167 {
148 currentParticle = theParticleTable->FindPa 168 currentParticle = theParticleTable->FindParticle(particleName);
149 idxCurrentChannel = -1; 169 idxCurrentChannel = -1;
150 currentDecayTable = nullptr; 170 currentDecayTable = nullptr;
151 } 171 }
152 172
153 if (currentParticle != nullptr) { << 173 if (currentParticle != nullptr )
>> 174 {
154 currentDecayTable = currentParticle->GetDe 175 currentDecayTable = currentParticle->GetDecayTable();
155 if ((currentDecayTable != nullptr) && (idx << 176 if ( (currentDecayTable != nullptr ) && (idxCurrentChannel >0) )
>> 177 {
156 currentChannel = currentDecayTable->GetD 178 currentChannel = currentDecayTable->GetDecayChannel(idxCurrentChannel);
157 } 179 }
158 else { << 180 else
>> 181 {
159 idxCurrentChannel = -1; 182 idxCurrentChannel = -1;
160 currentChannel = nullptr; 183 currentChannel = nullptr;
161 } 184 }
162 } 185 }
163 186
164 return currentParticle; 187 return currentParticle;
165 } 188 }
166 189
167 G4String G4DecayTableMessenger::GetCurrentValu 190 G4String G4DecayTableMessenger::GetCurrentValue(G4UIcommand* command)
168 { 191 {
169 G4String returnValue(1, '\0'); << 192 G4String returnValue('\0');
170 193
171 if (SetCurrentParticle() == nullptr) { << 194 if (SetCurrentParticle() == nullptr)
172 // no particle is selected. return null << 195 {
>> 196 // no particle is selected. return null
173 return returnValue; 197 return returnValue;
174 } 198 }
175 199
176 if (command == selectCmd) { << 200 if( command == selectCmd )
>> 201 {
177 // Command /particle/property/decay/select 202 // Command /particle/property/decay/select
178 returnValue = selectCmd->ConvertToString(i 203 returnValue = selectCmd->ConvertToString(idxCurrentChannel);
>> 204
179 } 205 }
180 else if (command == brCmd) { << 206 else if( command == brCmd )
181 if (currentChannel != nullptr) { << 207 {
>> 208 if ( currentChannel != nullptr)
>> 209 {
182 returnValue = brCmd->ConvertToString(cur 210 returnValue = brCmd->ConvertToString(currentChannel->GetBR());
183 } << 211 }
184 } 212 }
185 return returnValue; 213 return returnValue;
186 } 214 }
187 215