Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/particles/management/src/G4DecayTableMessenger.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /particles/management/src/G4DecayTableMessenger.cc (Version 11.3.0) and /particles/management/src/G4DecayTableMessenger.cc (Version 8.2)


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