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 << 27 // 26 // 28 // Author: H.Kurashige, 13 June 1997 << 27 // $Id$ 29 //-------------------------------------------- << 28 // >> 29 // >> 30 //--------------------------------------------------------------- >> 31 // >> 32 // G4DecayTableMessenger.cc >> 33 // >> 34 // Description: >> 35 // This is a messenger class to interface to exchange information >> 36 // between ParticleDefinition and UI. >> 37 // >> 38 // History: >> 39 // 13 June 1997, H. Kurashige : The 1st version created. >> 40 // 10 Nov. 1997 H. Kurashige : fixed bugs >> 41 // 08 jan. 1998 H. Kurashige : new UIcommnds >> 42 // >> 43 //--------------------------------------------------------------- 30 44 31 #include "G4DecayTableMessenger.hh" 45 #include "G4DecayTableMessenger.hh" 32 << 46 #include "G4UImanager.hh" >> 47 #include "G4UIdirectory.hh" >> 48 #include "G4UIcmdWithoutParameter.hh" >> 49 #include "G4UIcmdWithAnInteger.hh" >> 50 #include "G4UIcmdWithADouble.hh" >> 51 #include "G4VDecayChannel.hh" 33 #include "G4DecayTable.hh" 52 #include "G4DecayTable.hh" 34 #include "G4ParticleDefinition.hh" 53 #include "G4ParticleDefinition.hh" 35 #include "G4ParticleTable.hh" 54 #include "G4ParticleTable.hh" 36 #include "G4UIcmdWithADouble.hh" << 55 #include "G4ios.hh" // Include from 'system' 37 #include "G4UIcmdWithAnInteger.hh" << 56 #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 57 44 #include <iomanip> // Include from 'system' << 58 G4DecayTableMessenger::G4DecayTableMessenger(G4ParticleTable* pTable) 45 << 59 :theParticleTable(pTable), 46 G4DecayTableMessenger::G4DecayTableMessenger(G << 60 currentParticle(nullptr), >> 61 currentDecayTable(nullptr), >> 62 idxCurrentChannel(-1), >> 63 currentChannel(nullptr), >> 64 thisDirectory(nullptr), >> 65 dumpCmd(nullptr), >> 66 selectCmd(nullptr), >> 67 brCmd(nullptr) 47 { 68 { 48 if (theParticleTable == nullptr) { << 69 if ( theParticleTable == nullptr ) { 49 theParticleTable = G4ParticleTable::GetPar 70 theParticleTable = G4ParticleTable::GetParticleTable(); 50 } 71 } 51 currentParticle = nullptr; 72 currentParticle = nullptr; 52 73 53 // Command /particle/property/decay/ << 74 //Commnad /particle/property/decay/ 54 thisDirectory = new G4UIdirectory("/particle 75 thisDirectory = new G4UIdirectory("/particle/property/decay/"); 55 thisDirectory->SetGuidance("Decay Table cont 76 thisDirectory->SetGuidance("Decay Table control commands."); 56 77 57 // Command /particle/property/decay/select << 78 //Commnad /particle/property/decay/select 58 selectCmd = new G4UIcmdWithAnInteger("/parti << 79 selectCmd = new G4UIcmdWithAnInteger("/particle/property/decay/select",this); 59 selectCmd->SetGuidance("Enter index of decay 80 selectCmd->SetGuidance("Enter index of decay mode."); 60 selectCmd->SetParameterName("mode", true); 81 selectCmd->SetParameterName("mode", true); 61 selectCmd->SetDefaultValue(0); 82 selectCmd->SetDefaultValue(0); 62 selectCmd->SetRange("mode >=0"); 83 selectCmd->SetRange("mode >=0"); 63 currentChannel = nullptr; 84 currentChannel = nullptr; 64 85 65 // Command /particle/property/decay/dump << 86 //Commnad /particle/property/decay/dump 66 dumpCmd = new G4UIcmdWithoutParameter("/part << 87 dumpCmd = new G4UIcmdWithoutParameter("/particle/property/decay/dump",this); 67 dumpCmd->SetGuidance("Dump decay mode inform 88 dumpCmd->SetGuidance("Dump decay mode information."); 68 89 69 // Command /particle/property/decay/br << 90 //Command /particle/property/decay/br 70 brCmd = new G4UIcmdWithADouble("/particle/pr << 91 brCmd = new G4UIcmdWithADouble("/particle/property/decay/br",this); 71 brCmd->SetGuidance("Set branching ratio. [0< 92 brCmd->SetGuidance("Set branching ratio. [0< BR <1.0]"); 72 brCmd->SetParameterName("br", false); << 93 brCmd->SetParameterName("br",false); 73 brCmd->SetRange("(br >=0.0) && (br <=1.0)"); 94 brCmd->SetRange("(br >=0.0) && (br <=1.0)"); >> 95 74 } 96 } 75 97 76 G4DecayTableMessenger::~G4DecayTableMessenger( 98 G4DecayTableMessenger::~G4DecayTableMessenger() 77 { 99 { 78 delete dumpCmd; << 100 if (dumpCmd != nullptr) delete dumpCmd; 79 delete selectCmd; << 101 if (selectCmd != nullptr) delete selectCmd; 80 delete brCmd; << 102 if (brCmd != nullptr) delete brCmd; 81 delete thisDirectory; << 103 if (thisDirectory != nullptr) delete thisDirectory; 82 } 104 } 83 105 84 void G4DecayTableMessenger::SetNewValue(G4UIco << 106 void G4DecayTableMessenger::SetNewValue(G4UIcommand * command,G4String newValue) 85 { 107 { 86 if (SetCurrentParticle() == nullptr) { << 108 if (SetCurrentParticle()== nullptr) { 87 G4cout << "Particle is not selected yet !! 109 G4cout << "Particle is not selected yet !! Command ignored." << G4endl; 88 return; 110 return; 89 } 111 } 90 if (currentDecayTable == nullptr) { << 112 if (currentDecayTable== nullptr) { 91 G4cout << "The particle has no decay table 113 G4cout << "The particle has no decay table !! Command ignored." << G4endl; 92 return; 114 return; 93 } 115 } 94 116 95 if (command == dumpCmd) { << 117 if( command == dumpCmd ){ 96 // Command /particle/property/decay/dump << 118 //Commnad /particle/property/decay/dump 97 currentDecayTable->DumpInfo(); 119 currentDecayTable->DumpInfo(); 98 } << 120 99 else if (command == selectCmd) { << 121 } else if ( command == selectCmd ){ 100 // Command /particle/property/decay/select << 122 //Commnad /particle/property/decay/select 101 G4int index = selectCmd->GetNewIntValue(ne << 123 G4int index = selectCmd->GetNewIntValue(newValue) ; 102 currentChannel = currentDecayTable->GetDec 124 currentChannel = currentDecayTable->GetDecayChannel(index); 103 if (currentChannel == nullptr) { << 125 if ( currentChannel == nullptr ) { 104 G4cout << "Invalid index. Command ignore 126 G4cout << "Invalid index. Command ignored." << G4endl; 105 } << 127 } else { 106 else { << 107 idxCurrentChannel = index; 128 idxCurrentChannel = index; 108 } 129 } 109 } << 130 110 else { << 131 } else { 111 if (currentChannel == nullptr) { << 132 if ( currentChannel == nullptr ) { 112 G4cout << "Select a decay channel. Comma 133 G4cout << "Select a decay channel. Command ignored." << G4endl; 113 return; 134 return; 114 } 135 } 115 if (command == brCmd) { 136 if (command == brCmd) { 116 // Command /particle/property/decay/br << 137 //Commnad /particle/property/decay/br 117 G4double br = brCmd->GetNewDoubleValue(n << 138 G4double br = brCmd->GetNewDoubleValue(newValue); 118 if ((br < 0.0) || (br > 1.0)) { << 139 if( (br<0.0) || (br>1.0) ) { 119 G4cout << "Invalid brancing ratio. Com << 140 G4cout << "Invalid brancing ratio. Command ignored." << G4endl; 120 } << 141 } else { 121 else { << 142 currentChannel->SetBR(br); 122 currentChannel->SetBR(br); << 123 } 143 } 124 } 144 } 125 } 145 } 126 } 146 } 127 147 >> 148 128 G4ParticleDefinition* G4DecayTableMessenger::S 149 G4ParticleDefinition* G4DecayTableMessenger::SetCurrentParticle() 129 { 150 { 130 // set currentParticle pointer << 151 // set currentParticle pointer 131 // get particle name by asking G4ParticleMes 152 // get particle name by asking G4ParticleMessenger via UImanager 132 153 133 G4String particleName = G4UImanager::GetUIpo 154 G4String particleName = G4UImanager::GetUIpointer()->GetCurrentStringValue("/particle/select"); 134 155 135 if (currentParticle != nullptr) { << 156 if (currentParticle != nullptr ){ 136 // check whether selection is changed << 157 // check whether selection is changed 137 if (currentParticle->GetParticleName() != 158 if (currentParticle->GetParticleName() != particleName) { 138 currentParticle = theParticleTable->Find 159 currentParticle = theParticleTable->FindParticle(particleName); 139 idxCurrentChannel = -1; 160 idxCurrentChannel = -1; 140 currentDecayTable = nullptr; 161 currentDecayTable = nullptr; 141 } << 162 } else { 142 else { << 163 // no change 143 // no change << 144 return currentParticle; 164 return currentParticle; 145 } 165 } 146 } << 166 147 else { << 167 } else { 148 currentParticle = theParticleTable->FindPa 168 currentParticle = theParticleTable->FindParticle(particleName); 149 idxCurrentChannel = -1; 169 idxCurrentChannel = -1; 150 currentDecayTable = nullptr; 170 currentDecayTable = nullptr; 151 } << 171 if (currentParticle != nullptr ){ 152 << 172 currentDecayTable = currentParticle->GetDecayTable(); 153 if (currentParticle != nullptr) { << 173 if ((currentDecayTable != nullptr ) && (idxCurrentChannel >0) ) { 154 currentDecayTable = currentParticle->GetDe << 174 currentChannel = currentDecayTable->GetDecayChannel(idxCurrentChannel); 155 if ((currentDecayTable != nullptr) && (idx << 175 } else { 156 currentChannel = currentDecayTable->GetD << 176 idxCurrentChannel = -1; 157 } << 177 currentChannel = nullptr; 158 else { << 178 } 159 idxCurrentChannel = -1; << 160 currentChannel = nullptr; << 161 } 179 } 162 } << 163 180 >> 181 } 164 return currentParticle; 182 return currentParticle; 165 } 183 } 166 184 167 G4String G4DecayTableMessenger::GetCurrentValu << 185 G4String G4DecayTableMessenger::GetCurrentValue(G4UIcommand * command) 168 { 186 { 169 G4String returnValue(1, '\0'); << 187 G4String returnValue('\0'); 170 188 171 if (SetCurrentParticle() == nullptr) { << 189 if (SetCurrentParticle()==nullptr) { 172 // no particle is selected. return null << 190 // no particle is selected. return null 173 return returnValue; 191 return returnValue; 174 } 192 } 175 193 176 if (command == selectCmd) { << 194 if( command == selectCmd ){ 177 // Command /particle/property/decay/select << 195 //Commnad /particle/property/decay/select 178 returnValue = selectCmd->ConvertToString(i 196 returnValue = selectCmd->ConvertToString(idxCurrentChannel); 179 } << 197 180 else if (command == brCmd) { << 198 } else if( command == brCmd ){ 181 if (currentChannel != nullptr) { << 199 if ( currentChannel != nullptr) { 182 returnValue = brCmd->ConvertToString(cur 200 returnValue = brCmd->ConvertToString(currentChannel->GetBR()); 183 } << 201 } 184 } 202 } 185 return returnValue; 203 return returnValue; 186 } 204 } >> 205 >> 206 >> 207 >> 208 >> 209 >> 210 >> 211 >> 212 187 213