Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/materials/src/G4NistMessenger.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 /materials/src/G4NistMessenger.cc (Version 11.3.0) and /materials/src/G4NistMessenger.cc (Version 10.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 // $Id: G4NistMessenger.cc 72057 2013-07-04 13:07:29Z gcosmo $
                                                   >>  27 //
                                                   >>  28 //
 26 // File name:     G4NistMessenger                  29 // File name:     G4NistMessenger
 27 //                                                 30 //
 28 // Author:        Vladimir Ivanchenko              31 // Author:        Vladimir Ivanchenko
 29 //                                                 32 //
 30 // Creation date: 23.12.2004                       33 // Creation date: 23.12.2004
 31 //                                                 34 //
 32 // Modifications:                                  35 // Modifications:
 33 //                                                 36 //
 34 //                                                 37 //
 35 // -------------------------------------------     38 // -------------------------------------------------------------------
                                                   >>  39 //
                                                   >>  40 
                                                   >>  41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >>  42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 36                                                    43 
 37 #include "G4NistMessenger.hh"                      44 #include "G4NistMessenger.hh"
 38                                                    45 
 39 #include "G4ApplicationState.hh"               << 
 40 #include "G4DensityEffectData.hh"              << 
 41 #include "G4IonisParamMat.hh"                  << 
 42 #include "G4NistManager.hh"                        46 #include "G4NistManager.hh"
 43 #include "G4UIcmdWithAString.hh"               << 
 44 #include "G4UIcmdWithAnInteger.hh"             << 
 45 #include "G4UIdirectory.hh"                    << 
 46 #include "G4UIcommand.hh"                      << 
 47 #include "G4UIparameter.hh"                    << 
 48                                                    47 
 49 #include <sstream>                             <<  48 #include "G4UIdirectory.hh"
                                                   >>  49 #include "G4UIcmdWithAnInteger.hh"
                                                   >>  50 #include "G4UIcmdWithAString.hh"
                                                   >>  51 #include "G4IonisParamMat.hh"
                                                   >>  52 #include "G4DensityEffectData.hh"
 50                                                    53 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 52                                                    55 
 53 G4NistMessenger::G4NistMessenger(G4NistManager <<  56 G4NistMessenger::G4NistMessenger(G4NistManager* man)
                                                   >>  57   : manager(man)
 54 {                                                  58 {
 55   matDir = new G4UIdirectory("/material/");        59   matDir = new G4UIdirectory("/material/");
 56   matDir->SetGuidance("Commands for materials"     60   matDir->SetGuidance("Commands for materials");
 57                                                    61 
 58   verCmd = new G4UIcmdWithAnInteger("/material <<  62   verCmd = new G4UIcmdWithAnInteger("/material/verbose",this);
 59   verCmd->SetGuidance("Set verbose level.");       63   verCmd->SetGuidance("Set verbose level.");
 60   verCmd->AvailableForStates(G4State_PreInit,  <<  64   
 61   verCmd->SetToBeBroadcasted(false);           << 
 62                                                << 
 63   nistDir = new G4UIdirectory("/material/nist/     65   nistDir = new G4UIdirectory("/material/nist/");
 64   nistDir->SetGuidance("Commands for the nist      66   nistDir->SetGuidance("Commands for the nist dataBase");
 65                                                <<  67     
 66   prtElmCmd = new G4UIcmdWithAString("/materia <<  68   prtElmCmd = new G4UIcmdWithAString("/material/nist/printElement",this);
 67   prtElmCmd->SetGuidance("print element(s) in      69   prtElmCmd->SetGuidance("print element(s) in dataBase.");
 68   prtElmCmd->SetGuidance("symbol = element.");     70   prtElmCmd->SetGuidance("symbol = element.");
 69   prtElmCmd->SetGuidance("all    = all element     71   prtElmCmd->SetGuidance("all    = all elements.");
 70   prtElmCmd->SetParameterName("symbol", true);     72   prtElmCmd->SetParameterName("symbol", true);
 71   prtElmCmd->SetDefaultValue("all");               73   prtElmCmd->SetDefaultValue("all");
 72   prtElmCmd->AvailableForStates(G4State_PreIni <<  74   
 73   prtElmCmd->SetToBeBroadcasted(false);        <<  75   przElmCmd = new G4UIcmdWithAnInteger("/material/nist/printElementZ",this);
 74                                                << 
 75   przElmCmd = new G4UIcmdWithAnInteger("/mater << 
 76   przElmCmd->SetGuidance("print element Z in d     76   przElmCmd->SetGuidance("print element Z in dataBase.");
 77   przElmCmd->SetGuidance("0 = all elements.");     77   przElmCmd->SetGuidance("0 = all elements.");
 78   przElmCmd->SetParameterName("Z", true);          78   przElmCmd->SetParameterName("Z", true);
 79   przElmCmd->SetDefaultValue(0);                   79   przElmCmd->SetDefaultValue(0);
 80   przElmCmd->SetRange("0<=Z && Z<108");            80   przElmCmd->SetRange("0<=Z && Z<108");
 81   przElmCmd->AvailableForStates(G4State_PreIni <<  81    
 82   przElmCmd->SetToBeBroadcasted(false);        <<  82   lisMatCmd = new G4UIcmdWithAString("/material/nist/listMaterials",this);
 83                                                << 
 84   lisMatCmd = new G4UIcmdWithAString("/materia << 
 85   lisMatCmd->SetGuidance("Materials in Geant4      83   lisMatCmd->SetGuidance("Materials in Geant4 dataBase.");
 86   lisMatCmd->SetGuidance("simple - simple NIST     84   lisMatCmd->SetGuidance("simple - simple NIST materials.");
 87   lisMatCmd->SetGuidance("compound - compound      85   lisMatCmd->SetGuidance("compound - compound NIST materials.");
 88   lisMatCmd->SetGuidance("hep - HEP materials.     86   lisMatCmd->SetGuidance("hep - HEP materials.");
 89   lisMatCmd->SetGuidance("bio - biomedical mat     87   lisMatCmd->SetGuidance("bio - biomedical materials.");
 90   lisMatCmd->SetGuidance("all - list of all Ge     88   lisMatCmd->SetGuidance("all - list of all Geant4 materials.");
 91   lisMatCmd->SetParameterName("matlist", true)     89   lisMatCmd->SetParameterName("matlist", true);
 92   // lisMatCmd->SetCandidates("simple compound     90   // lisMatCmd->SetCandidates("simple compound hep bio all");
 93   lisMatCmd->SetDefaultValue("all");               91   lisMatCmd->SetDefaultValue("all");
 94   lisMatCmd->AvailableForStates(G4State_PreIni <<  92   
 95   lisMatCmd->SetToBeBroadcasted(false);        << 
 96                                                << 
 97   g4Dir = new G4UIdirectory("/material/g4/");      93   g4Dir = new G4UIdirectory("/material/g4/");
 98   g4Dir->SetGuidance("Commands for G4MaterialT     94   g4Dir->SetGuidance("Commands for G4MaterialTable");
 99                                                <<  95   
100   g4ElmCmd = new G4UIcmdWithAString("/material <<  96   g4ElmCmd = new G4UIcmdWithAString("/material/g4/printElement",this);
101   g4ElmCmd->SetGuidance("print Element from G4     97   g4ElmCmd->SetGuidance("print Element from G4ElementTable.");
102   g4ElmCmd->SetGuidance("all - all elements.")     98   g4ElmCmd->SetGuidance("all - all elements.");
103   g4ElmCmd->SetParameterName("elm", true);         99   g4ElmCmd->SetParameterName("elm", true);
104   g4ElmCmd->SetDefaultValue("all");               100   g4ElmCmd->SetDefaultValue("all");
105                                                << 101     
106   g4MatCmd = new G4UIcmdWithAString("/material << 102   g4MatCmd = new G4UIcmdWithAString("/material/g4/printMaterial",this);
107   g4MatCmd->SetGuidance("print Material from G    103   g4MatCmd->SetGuidance("print Material from G4MaterialTable.");
108   g4MatCmd->SetGuidance("all - all materials")    104   g4MatCmd->SetGuidance("all - all materials");
109   g4MatCmd->SetParameterName("pmat", true);       105   g4MatCmd->SetParameterName("pmat", true);
110   g4MatCmd->SetDefaultValue("all");               106   g4MatCmd->SetDefaultValue("all");
111   g4MatCmd->AvailableForStates(G4State_PreInit << 
112   g4MatCmd->SetToBeBroadcasted(false);         << 
113                                                   107 
114   g4DensCmd = new G4UIcmdWithAString("/materia << 108   g4DensCmd = new G4UIcmdWithAString("/material/g4/printDensityEffParam",this);
115   g4DensCmd->SetGuidance("print Material from     109   g4DensCmd->SetGuidance("print Material from G4DensityEffectData.");
116   g4DensCmd->SetGuidance("all - all materials"    110   g4DensCmd->SetGuidance("all - all materials");
117   g4DensCmd->SetParameterName("dmat", true);      111   g4DensCmd->SetParameterName("dmat", true);
118   g4DensCmd->SetDefaultValue("all");              112   g4DensCmd->SetDefaultValue("all");
119   g4DensCmd->AvailableForStates(G4State_PreIni << 
120   g4DensCmd->SetToBeBroadcasted(false);        << 
121                                                << 
122   densCmd = new G4UIcmdWithAString("/material/ << 
123   densCmd->SetGuidance("enable accurate comput << 
124   densCmd->SetGuidance("all - all materials.") << 
125   densCmd->SetParameterName("dens", true);     << 
126   densCmd->SetDefaultValue("all");             << 
127   densCmd->AvailableForStates(G4State_PreInit) << 
128   lisMatCmd->SetToBeBroadcasted(false);        << 
129                                                << 
130   adensCmd = new G4UIcmdWithAString("/material << 
131   adensCmd->SetGuidance("disable accurate comp << 
132   adensCmd->SetGuidance("all - all materials." << 
133   adensCmd->SetParameterName("dens", true);    << 
134   adensCmd->SetDefaultValue("all");            << 
135   adensCmd->AvailableForStates(G4State_PreInit << 
136   adensCmd->SetToBeBroadcasted(false);         << 
137                                                << 
138   fPosiCmd = new G4UIcommand("/material/g4/ort << 
139   fPosiCmd->SetGuidance("defined orto-positron << 
140   fPosiCmd->SetGuidance("via material name, al << 
141   fPosiCmd->AvailableForStates(G4State_PreInit << 
142   fPosiCmd->SetToBeBroadcasted(false);         << 
143                                                << 
144   auto p1 = new G4UIparameter("matname",'s', f << 
145   fPosiCmd->SetParameter(p1);                  << 
146                                                << 
147   auto p2 = new G4UIparameter("fraction", 'd', << 
148   p2->SetParameterRange("fraction>=0. && fract << 
149   fPosiCmd->SetParameter(p2);                  << 
150 }                                                 113 }
151                                                   114 
152 //....oooOO0OOooo........oooOO0OOooo........oo    115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
153                                                   116 
154 G4NistMessenger::~G4NistMessenger()               117 G4NistMessenger::~G4NistMessenger()
155 {                                                 118 {
156   delete verCmd;                               << 119   delete verCmd;  
157   delete prtElmCmd;                               120   delete prtElmCmd;
158   delete przElmCmd;                               121   delete przElmCmd;
159   delete lisMatCmd;                               122   delete lisMatCmd;
160   delete nistDir;                                 123   delete nistDir;
161                                                << 124   
162   delete g4ElmCmd;                             << 125   delete g4ElmCmd;   
163   delete g4MatCmd;                                126   delete g4MatCmd;
164   delete g4DensCmd;                               127   delete g4DensCmd;
165   delete densCmd;                              << 
166   delete adensCmd;                             << 
167   delete fPosiCmd;                             << 
168                                                << 
169   delete g4Dir;                                   128   delete g4Dir;
170   delete matDir;                               << 129   delete matDir;  
171 }                                                 130 }
172                                                   131 
173 //....oooOO0OOooo........oooOO0OOooo........oo    132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
174                                                   133 
175 void G4NistMessenger::SetNewValue(G4UIcommand*    134 void G4NistMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
176 {                                                 135 {
177   // G4cout << "G4NistMessenger::SetNewValue < << 136   //G4cout << "G4NistMessenger::SetNewValue <" << newValue << ">" << G4endl;
178   if (command == verCmd) {                     << 137   if (command == verCmd) 
179     manager->SetVerbose(verCmd->GetNewIntValue << 138    { manager->SetVerbose(verCmd->GetNewIntValue(newValue));}
180   }                                            << 139     
181   else if (command == prtElmCmd) {             << 140   else if (command == prtElmCmd)
182     manager->PrintElement(newValue);           << 141    { manager->PrintElement(newValue); }
183   }                                            << 142     
184   else if (command == przElmCmd) {                143   else if (command == przElmCmd) {
185     G4int Z = przElmCmd->GetNewIntValue(newVal    144     G4int Z = przElmCmd->GetNewIntValue(newValue);
186     if (Z >= 0 && Z < 108) {                   << 145     if(Z >= 0 && Z < 108) { manager->PrintElement(Z); }
187       manager->PrintElement(Z);                << 146   }   
188     }                                          << 147   else if (command == lisMatCmd) 
189   }                                            << 148    { manager->ListMaterials(newValue); }
190   else if (command == lisMatCmd) {             << 149 
191     manager->ListMaterials(newValue);          << 150   else if (command == g4ElmCmd)
192   }                                            << 151    { manager->PrintG4Element(newValue); }
193   else if (command == g4ElmCmd) {              << 152    
194     manager->PrintG4Element(newValue);         << 153   else if (command == g4MatCmd)
195   }                                            << 154    { manager->PrintG4Material(newValue); }
196   else if (command == g4MatCmd) {              << 155 
197     manager->PrintG4Material(newValue);        << 156   else if (command == g4DensCmd)
198   }                                            << 157     { G4IonisParamMat::GetDensityEffectData()->PrintData(newValue); }
199   else if (command == g4DensCmd) {             << 
200     G4IonisParamMat::GetDensityEffectData()->P << 
201   }                                            << 
202   else if (command == densCmd) {               << 
203     manager->SetDensityEffectCalculatorFlag(ne << 
204   }                                            << 
205   else if (command == adensCmd) {              << 
206     manager->SetDensityEffectCalculatorFlag(ne << 
207   }                                            << 
208   else if (command == fPosiCmd) {              << 
209     G4String mnam{""};                         << 
210     G4double f{0.0};                           << 
211     std::istringstream ss(newValue);           << 
212     ss >> mnam >> f;                           << 
213     // set fraction for all materials          << 
214     if (mnam == "all" || mnam == "none") {     << 
215       if (mnam == "none" || f < 0.0) { f = 0.0 << 
216       auto mtable = G4Material::GetMaterialTab << 
217       for ( auto const & mat : *mtable ) {     << 
218   mat->GetIonisation()->SetOrtoPositroniumFrac << 
219       }                                        << 
220     } else {                                   << 
221       // set fraction for one material         << 
222       auto mat = manager->FindOrBuildMaterial( << 
223       if (nullptr != mat) {                    << 
224   mat->GetIonisation()->SetOrtoPositroniumFrac << 
225       }                                        << 
226     }                                          << 
227   }                                            << 
228 }                                                 158 }
                                                   >> 159 
                                                   >> 160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
229                                                   161