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 // ------------------------------------------- 26 // ------------------------------------------------------------------- 27 // 27 // 28 // GEANT4 Class file 28 // GEANT4 Class file 29 // 29 // 30 // File name: G4EmParametersMessenger 30 // File name: G4EmParametersMessenger 31 // 31 // 32 // Author: Vladimir Ivanchenko created 32 // Author: Vladimir Ivanchenko created from G4EnergyLossMessenger 33 // 33 // 34 // Creation date: 22-05-2013 34 // Creation date: 22-05-2013 35 // 35 // 36 // ------------------------------------------- 36 // ------------------------------------------------------------------- 37 // 37 // 38 38 39 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 40 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 41 41 42 #include "G4EmParametersMessenger.hh" 42 #include "G4EmParametersMessenger.hh" 43 #include "G4UIdirectory.hh" 43 #include "G4UIdirectory.hh" 44 #include "G4UIcommand.hh" 44 #include "G4UIcommand.hh" 45 #include "G4UIparameter.hh" 45 #include "G4UIparameter.hh" 46 #include "G4UIcmdWithABool.hh" 46 #include "G4UIcmdWithABool.hh" 47 #include "G4UIcmdWithAnInteger.hh" 47 #include "G4UIcmdWithAnInteger.hh" 48 #include "G4UIcmdWithADouble.hh" 48 #include "G4UIcmdWithADouble.hh" 49 #include "G4UIcmdWithADoubleAndUnit.hh" 49 #include "G4UIcmdWithADoubleAndUnit.hh" 50 #include "G4UIcmdWithAString.hh" 50 #include "G4UIcmdWithAString.hh" 51 #include "G4UIcmdWith3VectorAndUnit.hh" 51 #include "G4UIcmdWith3VectorAndUnit.hh" 52 #include "G4UImanager.hh" 52 #include "G4UImanager.hh" 53 #include "G4MscStepLimitType.hh" 53 #include "G4MscStepLimitType.hh" 54 #include "G4NuclearFormfactorType.hh" 54 #include "G4NuclearFormfactorType.hh" 55 #include "G4EmParameters.hh" 55 #include "G4EmParameters.hh" 56 56 57 #include <sstream> 57 #include <sstream> 58 58 59 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 60 60 61 G4EmParametersMessenger::G4EmParametersMesseng 61 G4EmParametersMessenger::G4EmParametersMessenger(G4EmParameters* ptr) 62 : theParameters(ptr) 62 : theParameters(ptr) 63 { 63 { 64 emDirectory = new G4UIdirectory("/process/em << 64 gconvDirectory = new G4UIdirectory("/process/gconv/"); 65 emDirectory->SetGuidance("General commands f << 66 eLossDirectory = new G4UIdirectory("/process << 67 eLossDirectory->SetGuidance("Commands for en << 68 mscDirectory = new G4UIdirectory("/process/m << 69 mscDirectory->SetGuidance("Commands for EM s << 70 gconvDirectory = new G4UIdirectory("/process << 71 gconvDirectory->SetGuidance("Commands for EM 65 gconvDirectory->SetGuidance("Commands for EM gamma conversion BH5D model."); 72 dnaDirectory = new G4UIdirectory("/process/d << 66 eLossDirectory = new G4UIdirectory("/process/eLoss/"); >> 67 eLossDirectory->SetGuidance("Commands for EM processes."); >> 68 mscDirectory = new G4UIdirectory("/process/msc/"); >> 69 mscDirectory->SetGuidance("Commands for EM scattering processes."); >> 70 emDirectory = new G4UIdirectory("/process/em/"); >> 71 emDirectory->SetGuidance("General commands for EM processes."); >> 72 dnaDirectory = new G4UIdirectory("/process/dna/"); 73 dnaDirectory->SetGuidance("Commands for DNA 73 dnaDirectory->SetGuidance("Commands for DNA processes."); 74 74 75 flucCmd = new G4UIcmdWithABool("/process/eLo 75 flucCmd = new G4UIcmdWithABool("/process/eLoss/fluct",this); 76 flucCmd->SetGuidance("Enable/disable energy 76 flucCmd->SetGuidance("Enable/disable energy loss fluctuations."); 77 flucCmd->SetParameterName("choice",true); 77 flucCmd->SetParameterName("choice",true); 78 flucCmd->SetDefaultValue(true); 78 flucCmd->SetDefaultValue(true); 79 flucCmd->AvailableForStates(G4State_PreInit, 79 flucCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 80 flucCmd->SetToBeBroadcasted(false); 80 flucCmd->SetToBeBroadcasted(false); 81 81 82 rangeCmd = new G4UIcmdWithABool("/process/eL 82 rangeCmd = new G4UIcmdWithABool("/process/eLoss/CSDARange",this); 83 rangeCmd->SetGuidance("Enable/disable CSDA r 83 rangeCmd->SetGuidance("Enable/disable CSDA range calculation"); 84 rangeCmd->SetParameterName("range",true); 84 rangeCmd->SetParameterName("range",true); 85 rangeCmd->SetDefaultValue(false); 85 rangeCmd->SetDefaultValue(false); 86 rangeCmd->AvailableForStates(G4State_PreInit 86 rangeCmd->AvailableForStates(G4State_PreInit); 87 rangeCmd->SetToBeBroadcasted(false); 87 rangeCmd->SetToBeBroadcasted(false); 88 88 89 lpmCmd = new G4UIcmdWithABool("/process/eLos 89 lpmCmd = new G4UIcmdWithABool("/process/eLoss/LPM",this); 90 lpmCmd->SetGuidance("Enable/disable LPM effe 90 lpmCmd->SetGuidance("Enable/disable LPM effect calculation"); 91 lpmCmd->SetParameterName("lpm",true); 91 lpmCmd->SetParameterName("lpm",true); 92 lpmCmd->SetDefaultValue(true); 92 lpmCmd->SetDefaultValue(true); 93 lpmCmd->AvailableForStates(G4State_PreInit,G 93 lpmCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 94 lpmCmd->SetToBeBroadcasted(false); 94 lpmCmd->SetToBeBroadcasted(false); 95 95 96 rsCmd = new G4UIcmdWithABool("/process/eLoss 96 rsCmd = new G4UIcmdWithABool("/process/eLoss/useCutAsFinalRange",this); 97 rsCmd->SetGuidance("Enable/disable use of cu 97 rsCmd->SetGuidance("Enable/disable use of cut in range as a final range"); 98 rsCmd->SetParameterName("choice",true); 98 rsCmd->SetParameterName("choice",true); 99 rsCmd->SetDefaultValue(false); 99 rsCmd->SetDefaultValue(false); 100 rsCmd->AvailableForStates(G4State_PreInit,G4 100 rsCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 101 rsCmd->SetToBeBroadcasted(false); 101 rsCmd->SetToBeBroadcasted(false); 102 102 103 aplCmd = new G4UIcmdWithABool("/process/em/a 103 aplCmd = new G4UIcmdWithABool("/process/em/applyCuts",this); 104 aplCmd->SetGuidance("Enable/disable applying 104 aplCmd->SetGuidance("Enable/disable applying cuts for gamma processes"); 105 aplCmd->SetParameterName("apl",true); 105 aplCmd->SetParameterName("apl",true); 106 aplCmd->SetDefaultValue(false); 106 aplCmd->SetDefaultValue(false); 107 aplCmd->AvailableForStates(G4State_PreInit,G 107 aplCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 108 aplCmd->SetToBeBroadcasted(false); 108 aplCmd->SetToBeBroadcasted(false); 109 109 110 intCmd = new G4UIcmdWithABool("/process/em/i << 111 intCmd->SetGuidance("Enable/disable integral << 112 intCmd->SetParameterName("choice",true); << 113 intCmd->SetDefaultValue(true); << 114 intCmd->AvailableForStates(G4State_PreInit); << 115 intCmd->SetToBeBroadcasted(false); << 116 << 117 latCmd = new G4UIcmdWithABool("/process/msc/ 110 latCmd = new G4UIcmdWithABool("/process/msc/LateralDisplacement",this); 118 latCmd->SetGuidance("Enable/disable sampling 111 latCmd->SetGuidance("Enable/disable sampling of lateral displacement"); 119 latCmd->SetParameterName("lat",true); 112 latCmd->SetParameterName("lat",true); 120 latCmd->SetDefaultValue(true); 113 latCmd->SetDefaultValue(true); 121 latCmd->AvailableForStates(G4State_PreInit,G 114 latCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 122 latCmd->SetToBeBroadcasted(false); 115 latCmd->SetToBeBroadcasted(false); 123 116 124 lat96Cmd = new G4UIcmdWithABool("/process/ms 117 lat96Cmd = new G4UIcmdWithABool("/process/msc/LateralDisplacementAlg96",this); 125 lat96Cmd->SetGuidance("Enable/disable sampli 118 lat96Cmd->SetGuidance("Enable/disable sampling of lateral displacement"); 126 lat96Cmd->SetParameterName("lat96",true); 119 lat96Cmd->SetParameterName("lat96",true); 127 lat96Cmd->SetDefaultValue(false); 120 lat96Cmd->SetDefaultValue(false); 128 lat96Cmd->AvailableForStates(G4State_PreInit 121 lat96Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 129 lat96Cmd->SetToBeBroadcasted(false); 122 lat96Cmd->SetToBeBroadcasted(false); 130 123 131 mulatCmd = new G4UIcmdWithABool("/process/ms 124 mulatCmd = new G4UIcmdWithABool("/process/msc/MuHadLateralDisplacement",this); 132 mulatCmd->SetGuidance("Enable/disable sampli 125 mulatCmd->SetGuidance("Enable/disable sampling of lateral displacement for muons and hadrons"); 133 mulatCmd->SetParameterName("mulat",true); 126 mulatCmd->SetParameterName("mulat",true); 134 mulatCmd->SetDefaultValue(true); 127 mulatCmd->SetDefaultValue(true); 135 mulatCmd->AvailableForStates(G4State_PreInit 128 mulatCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 136 mulatCmd->SetToBeBroadcasted(false); 129 mulatCmd->SetToBeBroadcasted(false); 137 130 138 delCmd = new G4UIcmdWithABool("/process/eLos 131 delCmd = new G4UIcmdWithABool("/process/eLoss/UseAngularGenerator",this); 139 delCmd->SetGuidance("Enable usage of angular 132 delCmd->SetGuidance("Enable usage of angular generator for ionisation"); 140 delCmd->SetParameterName("del",true); 133 delCmd->SetParameterName("del",true); 141 delCmd->SetDefaultValue(false); 134 delCmd->SetDefaultValue(false); 142 delCmd->AvailableForStates(G4State_PreInit); 135 delCmd->AvailableForStates(G4State_PreInit); 143 delCmd->SetToBeBroadcasted(false); 136 delCmd->SetToBeBroadcasted(false); 144 137 145 mottCmd = new G4UIcmdWithABool("/process/msc 138 mottCmd = new G4UIcmdWithABool("/process/msc/UseMottCorrection",this); 146 mottCmd->SetGuidance("Enable usage of Mott c 139 mottCmd->SetGuidance("Enable usage of Mott corrections for e- elastic scattering"); 147 mottCmd->SetParameterName("mott",true); 140 mottCmd->SetParameterName("mott",true); 148 mottCmd->SetDefaultValue(false); 141 mottCmd->SetDefaultValue(false); 149 mottCmd->AvailableForStates(G4State_PreInit) 142 mottCmd->AvailableForStates(G4State_PreInit); 150 mottCmd->SetToBeBroadcasted(false); 143 mottCmd->SetToBeBroadcasted(false); 151 144 152 birksCmd = new G4UIcmdWithABool("/process/em << 145 birksCmd = new G4UIcmdWithABool("/process/msc/UseG4EmSaturation",this); 153 birksCmd->SetGuidance("Enable usage of built 146 birksCmd->SetGuidance("Enable usage of built-in Birks saturation"); 154 birksCmd->SetParameterName("birks",true); 147 birksCmd->SetParameterName("birks",true); 155 birksCmd->SetDefaultValue(false); 148 birksCmd->SetDefaultValue(false); 156 birksCmd->AvailableForStates(G4State_PreInit 149 birksCmd->AvailableForStates(G4State_PreInit,G4State_Init); 157 birksCmd->SetToBeBroadcasted(false); 150 birksCmd->SetToBeBroadcasted(false); 158 151 159 sharkCmd = new G4UIcmdWithABool("/process/em 152 sharkCmd = new G4UIcmdWithABool("/process/em/UseGeneralProcess",this); 160 sharkCmd->SetGuidance("Enable gamma, e+- gen 153 sharkCmd->SetGuidance("Enable gamma, e+- general process"); 161 sharkCmd->SetParameterName("gen",true); 154 sharkCmd->SetParameterName("gen",true); 162 sharkCmd->SetDefaultValue(false); 155 sharkCmd->SetDefaultValue(false); 163 sharkCmd->AvailableForStates(G4State_PreInit 156 sharkCmd->AvailableForStates(G4State_PreInit); 164 sharkCmd->SetToBeBroadcasted(false); 157 sharkCmd->SetToBeBroadcasted(false); 165 158 166 poCmd = new G4UIcmdWithABool("/process/em/Po << 167 poCmd->SetGuidance("Enable polarisation"); << 168 poCmd->AvailableForStates(G4State_PreInit); << 169 poCmd->SetToBeBroadcasted(false); << 170 << 171 sampleTCmd = new G4UIcmdWithABool("/process/ 159 sampleTCmd = new G4UIcmdWithABool("/process/em/enableSamplingTable",this); 172 sampleTCmd->SetGuidance("Enable usage of sam 160 sampleTCmd->SetGuidance("Enable usage of sampling table for secondary generation"); 173 sampleTCmd->SetParameterName("sampleT",true) 161 sampleTCmd->SetParameterName("sampleT",true); 174 sampleTCmd->SetDefaultValue(false); 162 sampleTCmd->SetDefaultValue(false); 175 sampleTCmd->AvailableForStates(G4State_PreIn 163 sampleTCmd->AvailableForStates(G4State_PreInit); 176 sampleTCmd->SetToBeBroadcasted(false); 164 sampleTCmd->SetToBeBroadcasted(false); 177 165 178 icru90Cmd = new G4UIcmdWithABool("/process/e 166 icru90Cmd = new G4UIcmdWithABool("/process/eLoss/UseICRU90",this); 179 icru90Cmd->SetGuidance("Enable usage of ICRU 167 icru90Cmd->SetGuidance("Enable usage of ICRU90 stopping powers"); 180 icru90Cmd->SetParameterName("icru90",true); 168 icru90Cmd->SetParameterName("icru90",true); 181 icru90Cmd->SetDefaultValue(false); 169 icru90Cmd->SetDefaultValue(false); 182 icru90Cmd->AvailableForStates(G4State_PreIni 170 icru90Cmd->AvailableForStates(G4State_PreInit); 183 icru90Cmd->SetToBeBroadcasted(false); 171 icru90Cmd->SetToBeBroadcasted(false); 184 172 185 mudatCmd = new G4UIcmdWithABool("/process/em 173 mudatCmd = new G4UIcmdWithABool("/process/em/MuDataFromFile",this); 186 mudatCmd->SetGuidance("Enable usage of muon 174 mudatCmd->SetGuidance("Enable usage of muon data from file"); 187 mudatCmd->SetParameterName("mudat",true); 175 mudatCmd->SetParameterName("mudat",true); 188 mudatCmd->SetDefaultValue(false); 176 mudatCmd->SetDefaultValue(false); 189 mudatCmd->AvailableForStates(G4State_PreInit 177 mudatCmd->AvailableForStates(G4State_PreInit); 190 mudatCmd->SetToBeBroadcasted(false); 178 mudatCmd->SetToBeBroadcasted(false); 191 179 192 peKCmd = new G4UIcmdWithABool("/process/em/P << 193 peKCmd->SetGuidance("Enable sampling of phot << 194 peKCmd->SetParameterName("peK",true); << 195 peKCmd->SetDefaultValue(true); << 196 peKCmd->AvailableForStates(G4State_PreInit); << 197 peKCmd->SetToBeBroadcasted(false); << 198 << 199 mscPCmd = new G4UIcmdWithABool("/process/msc << 200 mscPCmd->SetGuidance("Enable msc positron co << 201 mscPCmd->SetParameterName("mscPC",true); << 202 mscPCmd->SetDefaultValue(true); << 203 mscPCmd->AvailableForStates(G4State_PreInit, << 204 mscPCmd->SetToBeBroadcasted(false); << 205 << 206 pepicsCmd = new G4UIcmdWithABool("/process/e << 207 pepicsCmd->SetGuidance("Use EPICS2017 data f << 208 pepicsCmd->SetParameterName("pepics",true); << 209 pepicsCmd->SetDefaultValue(false); << 210 pepicsCmd->AvailableForStates(G4State_PreIni << 211 pepicsCmd->SetToBeBroadcasted(false); << 212 << 213 f3gCmd = new G4UIcmdWithABool("/process/em/3 << 214 f3gCmd->SetGuidance("Enable/disable 3 gamma << 215 f3gCmd->SetParameterName("f3gamma",true); << 216 f3gCmd->SetDefaultValue(false); << 217 f3gCmd->AvailableForStates(G4State_PreInit); << 218 f3gCmd->SetToBeBroadcasted(false); << 219 << 220 minEnCmd = new G4UIcmdWithADoubleAndUnit("/p 180 minEnCmd = new G4UIcmdWithADoubleAndUnit("/process/eLoss/minKinEnergy",this); 221 minEnCmd->SetGuidance("Set the min kinetic e 181 minEnCmd->SetGuidance("Set the min kinetic energy for EM tables"); 222 minEnCmd->SetParameterName("emin",true); 182 minEnCmd->SetParameterName("emin",true); 223 minEnCmd->SetUnitCategory("Energy"); 183 minEnCmd->SetUnitCategory("Energy"); 224 minEnCmd->AvailableForStates(G4State_PreInit 184 minEnCmd->AvailableForStates(G4State_PreInit); 225 minEnCmd->SetToBeBroadcasted(false); 185 minEnCmd->SetToBeBroadcasted(false); 226 186 227 maxEnCmd = new G4UIcmdWithADoubleAndUnit("/p 187 maxEnCmd = new G4UIcmdWithADoubleAndUnit("/process/eLoss/maxKinEnergy",this); 228 maxEnCmd->SetGuidance("Set the max kinetic e 188 maxEnCmd->SetGuidance("Set the max kinetic energy for EM tables"); 229 maxEnCmd->SetParameterName("emax",true); 189 maxEnCmd->SetParameterName("emax",true); 230 maxEnCmd->SetUnitCategory("Energy"); 190 maxEnCmd->SetUnitCategory("Energy"); 231 maxEnCmd->AvailableForStates(G4State_PreInit 191 maxEnCmd->AvailableForStates(G4State_PreInit); 232 maxEnCmd->SetToBeBroadcasted(false); 192 maxEnCmd->SetToBeBroadcasted(false); 233 193 234 cenCmd = new G4UIcmdWithADoubleAndUnit("/pro 194 cenCmd = new G4UIcmdWithADoubleAndUnit("/process/eLoss/maxKinEnergyCSDA",this); 235 cenCmd->SetGuidance("Set the max kinetic ene 195 cenCmd->SetGuidance("Set the max kinetic energy for CSDA table"); 236 cenCmd->SetParameterName("emaxCSDA",true); 196 cenCmd->SetParameterName("emaxCSDA",true); 237 cenCmd->SetUnitCategory("Energy"); 197 cenCmd->SetUnitCategory("Energy"); 238 cenCmd->AvailableForStates(G4State_PreInit); 198 cenCmd->AvailableForStates(G4State_PreInit); 239 cenCmd->SetToBeBroadcasted(false); 199 cenCmd->SetToBeBroadcasted(false); 240 200 241 max5DCmd = new G4UIcmdWithADoubleAndUnit("/p 201 max5DCmd = new G4UIcmdWithADoubleAndUnit("/process/em/max5DMuPairEnergy",this); 242 max5DCmd->SetGuidance("Set the max kinetic e 202 max5DCmd->SetGuidance("Set the max kinetic energy for 5D muon pair production"); 243 max5DCmd->SetParameterName("emax5D",true); 203 max5DCmd->SetParameterName("emax5D",true); 244 max5DCmd->SetUnitCategory("Energy"); 204 max5DCmd->SetUnitCategory("Energy"); 245 max5DCmd->AvailableForStates(G4State_PreInit 205 max5DCmd->AvailableForStates(G4State_PreInit); 246 max5DCmd->SetToBeBroadcasted(false); 206 max5DCmd->SetToBeBroadcasted(false); 247 207 248 lowEnCmd = new G4UIcmdWithADoubleAndUnit("/p 208 lowEnCmd = new G4UIcmdWithADoubleAndUnit("/process/em/lowestElectronEnergy",this); 249 lowEnCmd->SetGuidance("Set the lowest kineti 209 lowEnCmd->SetGuidance("Set the lowest kinetic energy for e+-"); 250 lowEnCmd->SetParameterName("elow",true); 210 lowEnCmd->SetParameterName("elow",true); 251 lowEnCmd->SetUnitCategory("Energy"); 211 lowEnCmd->SetUnitCategory("Energy"); 252 lowEnCmd->AvailableForStates(G4State_PreInit 212 lowEnCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 253 lowEnCmd->SetToBeBroadcasted(false); 213 lowEnCmd->SetToBeBroadcasted(false); 254 214 255 lowhEnCmd = new G4UIcmdWithADoubleAndUnit("/ 215 lowhEnCmd = new G4UIcmdWithADoubleAndUnit("/process/em/lowestMuHadEnergy",this); 256 lowhEnCmd->SetGuidance("Set the lowest kinet 216 lowhEnCmd->SetGuidance("Set the lowest kinetic energy for muons and hadrons"); 257 lowhEnCmd->SetParameterName("elowh",true); 217 lowhEnCmd->SetParameterName("elowh",true); 258 lowhEnCmd->SetUnitCategory("Energy"); 218 lowhEnCmd->SetUnitCategory("Energy"); 259 lowhEnCmd->AvailableForStates(G4State_PreIni 219 lowhEnCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 260 lowhEnCmd->SetToBeBroadcasted(false); 220 lowhEnCmd->SetToBeBroadcasted(false); 261 221 262 lowEn3Cmd = new G4UIcmdWithADoubleAndUnit("/ 222 lowEn3Cmd = new G4UIcmdWithADoubleAndUnit("/process/em/lowestTripletEnergy",this); 263 lowEn3Cmd->SetGuidance("Set the lowest kinet 223 lowEn3Cmd->SetGuidance("Set the lowest kinetic energy for triplet production"); 264 lowEn3Cmd->SetParameterName("elow3",true); 224 lowEn3Cmd->SetParameterName("elow3",true); 265 lowEn3Cmd->SetUnitCategory("Energy"); 225 lowEn3Cmd->SetUnitCategory("Energy"); 266 lowEn3Cmd->AvailableForStates(G4State_PreIni 226 lowEn3Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 267 lowEn3Cmd->SetToBeBroadcasted(false); 227 lowEn3Cmd->SetToBeBroadcasted(false); 268 228 269 lllCmd = new G4UIcmdWithADouble("/process/eL 229 lllCmd = new G4UIcmdWithADouble("/process/eLoss/linLossLimit",this); 270 lllCmd->SetGuidance("Set linearLossLimit par 230 lllCmd->SetGuidance("Set linearLossLimit parameter"); 271 lllCmd->SetParameterName("linlim",true); 231 lllCmd->SetParameterName("linlim",true); 272 lllCmd->AvailableForStates(G4State_PreInit,G 232 lllCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 273 lllCmd->SetToBeBroadcasted(false); 233 lllCmd->SetToBeBroadcasted(false); 274 234 275 brCmd = new G4UIcmdWithADoubleAndUnit("/proc 235 brCmd = new G4UIcmdWithADoubleAndUnit("/process/eLoss/bremThreshold",this); 276 brCmd->SetGuidance("Set e+- bremsstrahlung e 236 brCmd->SetGuidance("Set e+- bremsstrahlung energy threshold"); 277 brCmd->SetParameterName("emaxBrem",true); 237 brCmd->SetParameterName("emaxBrem",true); 278 brCmd->SetUnitCategory("Energy"); 238 brCmd->SetUnitCategory("Energy"); 279 brCmd->AvailableForStates(G4State_PreInit); 239 brCmd->AvailableForStates(G4State_PreInit); 280 brCmd->SetToBeBroadcasted(false); 240 brCmd->SetToBeBroadcasted(false); 281 241 282 br1Cmd = new G4UIcmdWithADoubleAndUnit("/pro 242 br1Cmd = new G4UIcmdWithADoubleAndUnit("/process/eLoss/bremMuHadThreshold",this); 283 br1Cmd->SetGuidance("Set muon/hadron bremsst 243 br1Cmd->SetGuidance("Set muon/hadron bremsstrahlung energy threshold"); 284 br1Cmd->SetParameterName("emaxMuHadBrem",tru 244 br1Cmd->SetParameterName("emaxMuHadBrem",true); 285 br1Cmd->SetUnitCategory("Energy"); 245 br1Cmd->SetUnitCategory("Energy"); 286 br1Cmd->AvailableForStates(G4State_PreInit); 246 br1Cmd->AvailableForStates(G4State_PreInit); 287 br1Cmd->SetToBeBroadcasted(false); 247 br1Cmd->SetToBeBroadcasted(false); 288 248 289 labCmd = new G4UIcmdWithADouble("/process/eL 249 labCmd = new G4UIcmdWithADouble("/process/eLoss/LambdaFactor",this); 290 labCmd->SetGuidance("Set lambdaFactor parame 250 labCmd->SetGuidance("Set lambdaFactor parameter for integral option"); 291 labCmd->SetParameterName("Fl",true); 251 labCmd->SetParameterName("Fl",true); 292 labCmd->AvailableForStates(G4State_PreInit,G 252 labCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 293 labCmd->SetToBeBroadcasted(false); 253 labCmd->SetToBeBroadcasted(false); 294 254 295 mscfCmd = new G4UIcmdWithADouble("/process/m 255 mscfCmd = new G4UIcmdWithADouble("/process/msc/FactorForAngleLimit",this); 296 mscfCmd->SetGuidance("Set factor for computa 256 mscfCmd->SetGuidance("Set factor for computation of a limit for -t (invariant transfer)"); 297 mscfCmd->SetParameterName("Fact",true); 257 mscfCmd->SetParameterName("Fact",true); 298 mscfCmd->SetRange("Fact>0"); 258 mscfCmd->SetRange("Fact>0"); 299 mscfCmd->SetDefaultValue(1.); 259 mscfCmd->SetDefaultValue(1.); 300 mscfCmd->AvailableForStates(G4State_PreInit) 260 mscfCmd->AvailableForStates(G4State_PreInit); 301 mscfCmd->SetToBeBroadcasted(false); 261 mscfCmd->SetToBeBroadcasted(false); 302 262 303 angCmd = new G4UIcmdWithADoubleAndUnit("/pro 263 angCmd = new G4UIcmdWithADoubleAndUnit("/process/msc/ThetaLimit",this); 304 angCmd->SetGuidance("Set the limit on the po 264 angCmd->SetGuidance("Set the limit on the polar angle for msc and single scattering"); 305 angCmd->SetParameterName("theta",true); 265 angCmd->SetParameterName("theta",true); 306 angCmd->SetUnitCategory("Angle"); 266 angCmd->SetUnitCategory("Angle"); 307 angCmd->AvailableForStates(G4State_PreInit); 267 angCmd->AvailableForStates(G4State_PreInit); 308 angCmd->SetToBeBroadcasted(false); 268 angCmd->SetToBeBroadcasted(false); 309 269 310 msceCmd = new G4UIcmdWithADoubleAndUnit("/pr 270 msceCmd = new G4UIcmdWithADoubleAndUnit("/process/msc/EnergyLimit",this); 311 msceCmd->SetGuidance("Set the upper energy l 271 msceCmd->SetGuidance("Set the upper energy limit for msc"); 312 msceCmd->SetParameterName("mscE",true); 272 msceCmd->SetParameterName("mscE",true); 313 msceCmd->SetUnitCategory("Energy"); 273 msceCmd->SetUnitCategory("Energy"); 314 msceCmd->AvailableForStates(G4State_PreInit) 274 msceCmd->AvailableForStates(G4State_PreInit); 315 msceCmd->SetToBeBroadcasted(false); 275 msceCmd->SetToBeBroadcasted(false); 316 276 317 nielCmd = new G4UIcmdWithADoubleAndUnit("/pr 277 nielCmd = new G4UIcmdWithADoubleAndUnit("/process/em/MaxEnergyNIEL",this); 318 nielCmd->SetGuidance("Set the upper energy l 278 nielCmd->SetGuidance("Set the upper energy limit for NIEL"); 319 nielCmd->SetParameterName("niel",true); 279 nielCmd->SetParameterName("niel",true); 320 nielCmd->SetUnitCategory("Energy"); 280 nielCmd->SetUnitCategory("Energy"); 321 nielCmd->AvailableForStates(G4State_PreInit) 281 nielCmd->AvailableForStates(G4State_PreInit); 322 nielCmd->SetToBeBroadcasted(false); 282 nielCmd->SetToBeBroadcasted(false); 323 283 324 frCmd = new G4UIcmdWithADouble("/process/msc 284 frCmd = new G4UIcmdWithADouble("/process/msc/RangeFactor",this); 325 frCmd->SetGuidance("Set RangeFactor for msc 285 frCmd->SetGuidance("Set RangeFactor for msc processes of e+-"); 326 frCmd->SetParameterName("Fr",true); 286 frCmd->SetParameterName("Fr",true); 327 frCmd->SetRange("Fr>0"); 287 frCmd->SetRange("Fr>0"); 328 frCmd->SetDefaultValue(0.04); 288 frCmd->SetDefaultValue(0.04); 329 frCmd->AvailableForStates(G4State_PreInit,G4 289 frCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 330 frCmd->SetToBeBroadcasted(false); 290 frCmd->SetToBeBroadcasted(false); 331 291 332 fr1Cmd = new G4UIcmdWithADouble("/process/ms 292 fr1Cmd = new G4UIcmdWithADouble("/process/msc/RangeFactorMuHad",this); 333 fr1Cmd->SetGuidance("Set RangeFactor for msc 293 fr1Cmd->SetGuidance("Set RangeFactor for msc processes of muons/hadrons"); 334 fr1Cmd->SetParameterName("Fr1",true); 294 fr1Cmd->SetParameterName("Fr1",true); 335 fr1Cmd->SetRange("Fr1>0"); 295 fr1Cmd->SetRange("Fr1>0"); 336 fr1Cmd->SetDefaultValue(0.2); 296 fr1Cmd->SetDefaultValue(0.2); 337 fr1Cmd->AvailableForStates(G4State_PreInit,G 297 fr1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 338 fr1Cmd->SetToBeBroadcasted(false); 298 fr1Cmd->SetToBeBroadcasted(false); 339 299 340 fgCmd = new G4UIcmdWithADouble("/process/msc 300 fgCmd = new G4UIcmdWithADouble("/process/msc/GeomFactor",this); 341 fgCmd->SetGuidance("Set GeomFactor parameter 301 fgCmd->SetGuidance("Set GeomFactor parameter for msc processes"); 342 fgCmd->SetParameterName("Fg",true); 302 fgCmd->SetParameterName("Fg",true); 343 fgCmd->SetRange("Fg>0"); 303 fgCmd->SetRange("Fg>0"); 344 fgCmd->SetDefaultValue(2.5); 304 fgCmd->SetDefaultValue(2.5); 345 fgCmd->AvailableForStates(G4State_PreInit,G4 305 fgCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 346 fgCmd->SetToBeBroadcasted(false); 306 fgCmd->SetToBeBroadcasted(false); 347 307 348 skinCmd = new G4UIcmdWithADouble("/process/m 308 skinCmd = new G4UIcmdWithADouble("/process/msc/Skin",this); 349 skinCmd->SetGuidance("Set skin parameter for 309 skinCmd->SetGuidance("Set skin parameter for msc processes"); 350 skinCmd->SetParameterName("skin",true); 310 skinCmd->SetParameterName("skin",true); 351 skinCmd->AvailableForStates(G4State_PreInit, 311 skinCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 352 skinCmd->SetToBeBroadcasted(false); 312 skinCmd->SetToBeBroadcasted(false); 353 313 354 screCmd = new G4UIcmdWithADouble("/process/m 314 screCmd = new G4UIcmdWithADouble("/process/msc/ScreeningFactor",this); 355 screCmd->SetGuidance("Set screening factor") 315 screCmd->SetGuidance("Set screening factor"); 356 screCmd->SetParameterName("screen",true); 316 screCmd->SetParameterName("screen",true); 357 screCmd->AvailableForStates(G4State_PreInit) 317 screCmd->AvailableForStates(G4State_PreInit); 358 screCmd->SetToBeBroadcasted(false); 318 screCmd->SetToBeBroadcasted(false); 359 319 360 safCmd = new G4UIcmdWithADouble("/process/ms 320 safCmd = new G4UIcmdWithADouble("/process/msc/SafetyFactor",this); 361 safCmd->SetGuidance("Set safety factor"); 321 safCmd->SetGuidance("Set safety factor"); 362 safCmd->SetParameterName("fsafe",true); 322 safCmd->SetParameterName("fsafe",true); 363 safCmd->AvailableForStates(G4State_PreInit); 323 safCmd->AvailableForStates(G4State_PreInit); 364 safCmd->SetToBeBroadcasted(false); 324 safCmd->SetToBeBroadcasted(false); 365 325 366 llimCmd = new G4UIcmdWithADoubleAndUnit("/pr 326 llimCmd = new G4UIcmdWithADoubleAndUnit("/process/msc/LambdaLimit",this); 367 llimCmd->SetGuidance("Set the upper energy l 327 llimCmd->SetGuidance("Set the upper energy limit for NIEL"); 368 llimCmd->SetParameterName("ll",true); 328 llimCmd->SetParameterName("ll",true); 369 llimCmd->SetUnitCategory("Length"); 329 llimCmd->SetUnitCategory("Length"); 370 llimCmd->AvailableForStates(G4State_PreInit) 330 llimCmd->AvailableForStates(G4State_PreInit); 371 llimCmd->SetToBeBroadcasted(false); 331 llimCmd->SetToBeBroadcasted(false); 372 332 373 amCmd = new G4UIcmdWithAnInteger("/process/e << 333 amCmd = new G4UIcmdWithAnInteger("/process/eLoss/binsPerDecade",this); 374 amCmd->SetGuidance("Set number of bins per d 334 amCmd->SetGuidance("Set number of bins per decade for EM tables"); 375 amCmd->SetParameterName("bins",true); 335 amCmd->SetParameterName("bins",true); 376 amCmd->SetDefaultValue(7); 336 amCmd->SetDefaultValue(7); 377 amCmd->AvailableForStates(G4State_PreInit); 337 amCmd->AvailableForStates(G4State_PreInit); 378 amCmd->SetToBeBroadcasted(false); 338 amCmd->SetToBeBroadcasted(false); 379 339 380 verCmd = new G4UIcmdWithAnInteger("/process/ 340 verCmd = new G4UIcmdWithAnInteger("/process/eLoss/verbose",this); 381 verCmd->SetGuidance("Set verbose level for E 341 verCmd->SetGuidance("Set verbose level for EM physics"); 382 verCmd->SetParameterName("verb",true); 342 verCmd->SetParameterName("verb",true); 383 verCmd->SetDefaultValue(1); 343 verCmd->SetDefaultValue(1); 384 verCmd->AvailableForStates(G4State_PreInit,G 344 verCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 385 verCmd->SetToBeBroadcasted(false); 345 verCmd->SetToBeBroadcasted(false); 386 346 387 ver1Cmd = new G4UIcmdWithAnInteger("/process 347 ver1Cmd = new G4UIcmdWithAnInteger("/process/em/verbose",this); 388 ver1Cmd->SetGuidance("Set verbose level for 348 ver1Cmd->SetGuidance("Set verbose level for EM physics"); 389 ver1Cmd->SetParameterName("verb1",true); 349 ver1Cmd->SetParameterName("verb1",true); 390 ver1Cmd->SetDefaultValue(1); 350 ver1Cmd->SetDefaultValue(1); 391 ver1Cmd->AvailableForStates(G4State_PreInit, 351 ver1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 392 ver1Cmd->SetToBeBroadcasted(false); 352 ver1Cmd->SetToBeBroadcasted(false); 393 353 394 ver2Cmd = new G4UIcmdWithAnInteger("/process 354 ver2Cmd = new G4UIcmdWithAnInteger("/process/em/workerVerbose",this); 395 ver2Cmd->SetGuidance("Set worker verbose lev 355 ver2Cmd->SetGuidance("Set worker verbose level for EM physics"); 396 ver2Cmd->SetParameterName("verb2",true); 356 ver2Cmd->SetParameterName("verb2",true); 397 ver2Cmd->SetDefaultValue(0); 357 ver2Cmd->SetDefaultValue(0); 398 ver2Cmd->AvailableForStates(G4State_PreInit, 358 ver2Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 399 ver2Cmd->SetToBeBroadcasted(false); 359 ver2Cmd->SetToBeBroadcasted(false); 400 360 401 nFreeCmd = new G4UIcmdWithAnInteger("/proces << 402 nFreeCmd->SetGuidance("Set number for logari << 403 nFreeCmd->SetParameterName("nFree",true); << 404 nFreeCmd->SetDefaultValue(2); << 405 nFreeCmd->AvailableForStates(G4State_PreInit << 406 nFreeCmd->SetToBeBroadcasted(false); << 407 << 408 transWithMscCmd = new G4UIcmdWithAString("/p << 409 transWithMscCmd->SetGuidance("Enable/disable << 410 transWithMscCmd->SetParameterName("trans",tr << 411 transWithMscCmd->SetCandidates("Disabled Ena << 412 transWithMscCmd->AvailableForStates(G4State_ << 413 transWithMscCmd->SetToBeBroadcasted(false); << 414 << 415 mscCmd = new G4UIcmdWithAString("/process/ms 361 mscCmd = new G4UIcmdWithAString("/process/msc/StepLimit",this); 416 mscCmd->SetGuidance("Set msc step limitation 362 mscCmd->SetGuidance("Set msc step limitation type"); 417 mscCmd->SetParameterName("StepLim",true); 363 mscCmd->SetParameterName("StepLim",true); 418 mscCmd->SetCandidates("Minimal UseSafety Use 364 mscCmd->SetCandidates("Minimal UseSafety UseSafetyPlus UseDistanceToBoundary"); 419 mscCmd->AvailableForStates(G4State_PreInit,G 365 mscCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 420 mscCmd->SetToBeBroadcasted(false); 366 mscCmd->SetToBeBroadcasted(false); 421 367 422 msc1Cmd = new G4UIcmdWithAString("/process/m 368 msc1Cmd = new G4UIcmdWithAString("/process/msc/StepLimitMuHad",this); 423 msc1Cmd->SetGuidance("Set msc step limitatio 369 msc1Cmd->SetGuidance("Set msc step limitation type for muons/hadrons"); 424 msc1Cmd->SetParameterName("StepLim1",true); 370 msc1Cmd->SetParameterName("StepLim1",true); 425 msc1Cmd->SetCandidates("Minimal UseSafety Us 371 msc1Cmd->SetCandidates("Minimal UseSafety UseSafetyPlus UseDistanceToBoundary"); 426 msc1Cmd->AvailableForStates(G4State_PreInit, 372 msc1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 427 msc1Cmd->SetToBeBroadcasted(false); 373 msc1Cmd->SetToBeBroadcasted(false); 428 374 429 dumpCmd = new G4UIcommand("/process/em/print 375 dumpCmd = new G4UIcommand("/process/em/printParameters",this); 430 dumpCmd->SetGuidance("Print all EM parameter 376 dumpCmd->SetGuidance("Print all EM parameters."); 431 dumpCmd->AvailableForStates(G4State_PreInit, 377 dumpCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 432 dumpCmd->SetToBeBroadcasted(false); 378 dumpCmd->SetToBeBroadcasted(false); 433 379 434 nffCmd = new G4UIcmdWithAString("/process/em 380 nffCmd = new G4UIcmdWithAString("/process/em/setNuclearFormFactor",this); 435 nffCmd->SetGuidance("Define type of nuclear 381 nffCmd->SetGuidance("Define type of nuclear form-factor"); 436 nffCmd->SetParameterName("NucFF",true); 382 nffCmd->SetParameterName("NucFF",true); 437 nffCmd->SetCandidates("None Exponential Gaus 383 nffCmd->SetCandidates("None Exponential Gaussian Flat"); 438 nffCmd->AvailableForStates(G4State_PreInit); 384 nffCmd->AvailableForStates(G4State_PreInit); 439 nffCmd->SetToBeBroadcasted(false); 385 nffCmd->SetToBeBroadcasted(false); 440 386 441 ssCmd = new G4UIcmdWithAString("/process/em/ 387 ssCmd = new G4UIcmdWithAString("/process/em/setSingleScattering",this); 442 ssCmd->SetGuidance("Define type of e+- singl 388 ssCmd->SetGuidance("Define type of e+- single scattering model"); 443 ssCmd->SetParameterName("SS",true); 389 ssCmd->SetParameterName("SS",true); 444 ssCmd->SetCandidates("WVI Mott DPWA"); 390 ssCmd->SetCandidates("WVI Mott DPWA"); 445 ssCmd->AvailableForStates(G4State_PreInit); 391 ssCmd->AvailableForStates(G4State_PreInit); 446 ssCmd->SetToBeBroadcasted(false); 392 ssCmd->SetToBeBroadcasted(false); 447 393 448 fluc1Cmd = new G4UIcmdWithAString("/process/ << 449 fluc1Cmd->SetGuidance("Define type of energy << 450 fluc1Cmd->SetParameterName("Fluc1",true); << 451 fluc1Cmd->SetCandidates("Dummy Universal Urb << 452 fluc1Cmd->AvailableForStates(G4State_PreInit << 453 fluc1Cmd->SetToBeBroadcasted(false); << 454 << 455 posiCmd = new G4UIcmdWithAString("/process/e << 456 posiCmd->SetGuidance("Define model of positr << 457 posiCmd->SetParameterName("Posi",true); << 458 posiCmd->SetCandidates("Simple Allison OrePa << 459 posiCmd->AvailableForStates(G4State_PreInit) << 460 posiCmd->SetToBeBroadcasted(false); << 461 << 462 tripletCmd = new G4UIcmdWithAnInteger("/proc 394 tripletCmd = new G4UIcmdWithAnInteger("/process/gconv/conversionType",this); 463 tripletCmd->SetGuidance("gamma conversion tr 395 tripletCmd->SetGuidance("gamma conversion triplet/nuclear generation type:"); 464 tripletCmd->SetGuidance("0 - (default) both 396 tripletCmd->SetGuidance("0 - (default) both triplet and nuclear"); 465 tripletCmd->SetGuidance("1 - force nuclear") 397 tripletCmd->SetGuidance("1 - force nuclear"); 466 tripletCmd->SetGuidance("2 - force triplet") 398 tripletCmd->SetGuidance("2 - force triplet"); 467 tripletCmd->SetParameterName("type",false); 399 tripletCmd->SetParameterName("type",false); 468 tripletCmd->SetRange("type >= 0 && type <= 2 400 tripletCmd->SetRange("type >= 0 && type <= 2"); 469 tripletCmd->SetDefaultValue(0); 401 tripletCmd->SetDefaultValue(0); 470 tripletCmd->AvailableForStates(G4State_PreIn 402 tripletCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 471 tripletCmd->SetToBeBroadcasted(false); 403 tripletCmd->SetToBeBroadcasted(false); 472 404 473 onIsolatedCmd = new G4UIcmdWithABool("/proce 405 onIsolatedCmd = new G4UIcmdWithABool("/process/gconv/onIsolated",this); 474 onIsolatedCmd->SetGuidance("Conversion on is 406 onIsolatedCmd->SetGuidance("Conversion on isolated charged particles"); 475 onIsolatedCmd->SetGuidance("false (default) 407 onIsolatedCmd->SetGuidance("false (default) : atomic electron screening"); 476 onIsolatedCmd->SetGuidance("true : conversio 408 onIsolatedCmd->SetGuidance("true : conversion on isolated particles."); 477 onIsolatedCmd->SetParameterName("flag",false 409 onIsolatedCmd->SetParameterName("flag",false); 478 onIsolatedCmd->SetDefaultValue(false); 410 onIsolatedCmd->SetDefaultValue(false); 479 onIsolatedCmd->AvailableForStates(G4State_Pr 411 onIsolatedCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 480 onIsolatedCmd->SetToBeBroadcasted(false); 412 onIsolatedCmd->SetToBeBroadcasted(false); 481 } 413 } 482 414 483 //....oooOO0OOooo........oooOO0OOooo........oo 415 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 484 416 485 G4EmParametersMessenger::~G4EmParametersMessen 417 G4EmParametersMessenger::~G4EmParametersMessenger() 486 { 418 { 487 delete gconvDirectory; 419 delete gconvDirectory; 488 delete eLossDirectory; 420 delete eLossDirectory; 489 delete mscDirectory; 421 delete mscDirectory; 490 delete emDirectory; 422 delete emDirectory; 491 delete dnaDirectory; 423 delete dnaDirectory; 492 424 493 delete flucCmd; 425 delete flucCmd; 494 delete rangeCmd; 426 delete rangeCmd; 495 delete lpmCmd; 427 delete lpmCmd; 496 delete rsCmd; 428 delete rsCmd; 497 delete aplCmd; 429 delete aplCmd; 498 delete intCmd; << 499 delete latCmd; 430 delete latCmd; 500 delete lat96Cmd; 431 delete lat96Cmd; 501 delete mulatCmd; 432 delete mulatCmd; 502 delete delCmd; 433 delete delCmd; 503 delete mottCmd; 434 delete mottCmd; 504 delete birksCmd; 435 delete birksCmd; 505 delete sharkCmd; 436 delete sharkCmd; 506 delete onIsolatedCmd; 437 delete onIsolatedCmd; 507 delete sampleTCmd; 438 delete sampleTCmd; 508 delete poCmd; << 509 delete icru90Cmd; 439 delete icru90Cmd; 510 delete mudatCmd; 440 delete mudatCmd; 511 delete peKCmd; << 512 delete f3gCmd; << 513 delete mscPCmd; << 514 delete pepicsCmd; << 515 441 516 delete minEnCmd; 442 delete minEnCmd; 517 delete maxEnCmd; 443 delete maxEnCmd; 518 delete max5DCmd; 444 delete max5DCmd; 519 delete cenCmd; 445 delete cenCmd; 520 delete lowEnCmd; 446 delete lowEnCmd; 521 delete lowhEnCmd; 447 delete lowhEnCmd; 522 delete lowEn3Cmd; 448 delete lowEn3Cmd; 523 delete lllCmd; 449 delete lllCmd; 524 delete brCmd; 450 delete brCmd; 525 delete br1Cmd; 451 delete br1Cmd; 526 delete labCmd; 452 delete labCmd; 527 delete mscfCmd; 453 delete mscfCmd; 528 delete angCmd; 454 delete angCmd; 529 delete msceCmd; 455 delete msceCmd; 530 delete nielCmd; 456 delete nielCmd; 531 delete frCmd; 457 delete frCmd; 532 delete fr1Cmd; 458 delete fr1Cmd; 533 delete fgCmd; 459 delete fgCmd; 534 delete skinCmd; 460 delete skinCmd; 535 delete safCmd; 461 delete safCmd; 536 delete llimCmd; 462 delete llimCmd; 537 delete screCmd; 463 delete screCmd; 538 464 539 delete amCmd; 465 delete amCmd; 540 delete verCmd; 466 delete verCmd; 541 delete ver1Cmd; 467 delete ver1Cmd; 542 delete ver2Cmd; 468 delete ver2Cmd; 543 delete transWithMscCmd; << 544 delete nFreeCmd; << 545 delete tripletCmd; 469 delete tripletCmd; 546 470 547 delete mscCmd; 471 delete mscCmd; 548 delete msc1Cmd; 472 delete msc1Cmd; 549 delete nffCmd; 473 delete nffCmd; 550 delete ssCmd; 474 delete ssCmd; 551 delete fluc1Cmd; << 552 delete posiCmd; << 553 475 554 delete dumpCmd; 476 delete dumpCmd; 555 } 477 } 556 478 557 //....oooOO0OOooo........oooOO0OOooo........oo 479 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 558 480 559 void G4EmParametersMessenger::SetNewValue(G4UI 481 void G4EmParametersMessenger::SetNewValue(G4UIcommand* command, 560 G4St 482 G4String newValue) 561 { 483 { 562 G4bool physicsModified = false; 484 G4bool physicsModified = false; 563 if (command == flucCmd) { 485 if (command == flucCmd) { 564 theParameters->SetLossFluctuations(flucCmd 486 theParameters->SetLossFluctuations(flucCmd->GetNewBoolValue(newValue)); 565 physicsModified = true; 487 physicsModified = true; 566 } else if (command == rangeCmd) { 488 } else if (command == rangeCmd) { 567 theParameters->SetBuildCSDARange(rangeCmd- 489 theParameters->SetBuildCSDARange(rangeCmd->GetNewBoolValue(newValue)); 568 } else if (command == lpmCmd) { 490 } else if (command == lpmCmd) { 569 theParameters->SetLPM(lpmCmd->GetNewBoolVa 491 theParameters->SetLPM(lpmCmd->GetNewBoolValue(newValue)); 570 physicsModified = true; 492 physicsModified = true; 571 } else if (command == rsCmd) { 493 } else if (command == rsCmd) { 572 theParameters->SetUseCutAsFinalRange(rsCmd 494 theParameters->SetUseCutAsFinalRange(rsCmd->GetNewBoolValue(newValue)); 573 physicsModified = true; 495 physicsModified = true; 574 } else if (command == aplCmd) { 496 } else if (command == aplCmd) { 575 theParameters->SetApplyCuts(aplCmd->GetNew 497 theParameters->SetApplyCuts(aplCmd->GetNewBoolValue(newValue)); 576 physicsModified = true; 498 physicsModified = true; 577 } else if (command == intCmd) { << 578 theParameters->SetIntegral(intCmd->GetNewB << 579 } else if (command == latCmd) { 499 } else if (command == latCmd) { 580 theParameters->SetLateralDisplacement(latC 500 theParameters->SetLateralDisplacement(latCmd->GetNewBoolValue(newValue)); 581 physicsModified = true; 501 physicsModified = true; 582 } else if (command == lat96Cmd) { 502 } else if (command == lat96Cmd) { 583 theParameters->SetLateralDisplacementAlg96 503 theParameters->SetLateralDisplacementAlg96(lat96Cmd->GetNewBoolValue(newValue)); 584 physicsModified = true; 504 physicsModified = true; 585 } else if (command == mulatCmd) { 505 } else if (command == mulatCmd) { 586 theParameters->SetMuHadLateralDisplacement 506 theParameters->SetMuHadLateralDisplacement(mulatCmd->GetNewBoolValue(newValue)); 587 physicsModified = true; 507 physicsModified = true; 588 } else if (command == delCmd) { 508 } else if (command == delCmd) { 589 theParameters->ActivateAngularGeneratorFor 509 theParameters->ActivateAngularGeneratorForIonisation(delCmd->GetNewBoolValue(newValue)); 590 } else if (command == mottCmd) { 510 } else if (command == mottCmd) { 591 theParameters->SetUseMottCorrection(mottCm 511 theParameters->SetUseMottCorrection(mottCmd->GetNewBoolValue(newValue)); 592 } else if (command == birksCmd) { 512 } else if (command == birksCmd) { 593 theParameters->SetBirksActive(birksCmd->Ge 513 theParameters->SetBirksActive(birksCmd->GetNewBoolValue(newValue)); 594 } else if (command == icru90Cmd) { 514 } else if (command == icru90Cmd) { 595 theParameters->SetUseICRU90Data(icru90Cmd- 515 theParameters->SetUseICRU90Data(icru90Cmd->GetNewBoolValue(newValue)); 596 } else if (command == sharkCmd) { 516 } else if (command == sharkCmd) { 597 theParameters->SetGeneralProcessActive(sha 517 theParameters->SetGeneralProcessActive(sharkCmd->GetNewBoolValue(newValue)); 598 } else if (command == poCmd) { << 599 theParameters->SetEnablePolarisation(poCmd << 600 } else if (command == sampleTCmd) { 518 } else if (command == sampleTCmd) { 601 theParameters->SetEnableSamplingTable(samp 519 theParameters->SetEnableSamplingTable(sampleTCmd->GetNewBoolValue(newValue)); 602 } else if (command == mudatCmd) { 520 } else if (command == mudatCmd) { 603 theParameters->SetRetrieveMuDataFromFile(m 521 theParameters->SetRetrieveMuDataFromFile(mudatCmd->GetNewBoolValue(newValue)); 604 } else if (command == peKCmd) { << 605 theParameters->SetPhotoeffectBelowKShell(p << 606 } else if (command == f3gCmd) { << 607 theParameters->Set3GammaAnnihilationOnFly( << 608 } else if (command == mscPCmd) { << 609 theParameters->SetMscPositronCorrection(ms << 610 } else if (command == pepicsCmd) { << 611 theParameters->SetUseEPICS2017XS(pepicsCmd << 612 522 613 } else if (command == minEnCmd) { 523 } else if (command == minEnCmd) { 614 theParameters->SetMinEnergy(minEnCmd->GetN 524 theParameters->SetMinEnergy(minEnCmd->GetNewDoubleValue(newValue)); 615 } else if (command == maxEnCmd) { 525 } else if (command == maxEnCmd) { 616 theParameters->SetMaxEnergy(maxEnCmd->GetN 526 theParameters->SetMaxEnergy(maxEnCmd->GetNewDoubleValue(newValue)); 617 } else if (command == max5DCmd) { 527 } else if (command == max5DCmd) { 618 theParameters->SetMaxEnergyFor5DMuPair(max 528 theParameters->SetMaxEnergyFor5DMuPair(max5DCmd->GetNewDoubleValue(newValue)); 619 } else if (command == cenCmd) { 529 } else if (command == cenCmd) { 620 theParameters->SetMaxEnergyForCSDARange(ce 530 theParameters->SetMaxEnergyForCSDARange(cenCmd->GetNewDoubleValue(newValue)); 621 physicsModified = true; 531 physicsModified = true; 622 } else if (command == lowEnCmd) { 532 } else if (command == lowEnCmd) { 623 theParameters->SetLowestElectronEnergy(low 533 theParameters->SetLowestElectronEnergy(lowEnCmd->GetNewDoubleValue(newValue)); 624 physicsModified = true; 534 physicsModified = true; 625 } else if (command == lowEn3Cmd) { 535 } else if (command == lowEn3Cmd) { 626 theParameters->SetLowestTripletEnergy(lowE 536 theParameters->SetLowestTripletEnergy(lowEn3Cmd->GetNewDoubleValue(newValue)); 627 physicsModified = true; 537 physicsModified = true; 628 } else if (command == lowhEnCmd) { 538 } else if (command == lowhEnCmd) { 629 theParameters->SetLowestMuHadEnergy(lowhEn 539 theParameters->SetLowestMuHadEnergy(lowhEnCmd->GetNewDoubleValue(newValue)); 630 physicsModified = true; 540 physicsModified = true; 631 } else if (command == lllCmd) { 541 } else if (command == lllCmd) { 632 theParameters->SetLinearLossLimit(lllCmd-> 542 theParameters->SetLinearLossLimit(lllCmd->GetNewDoubleValue(newValue)); 633 physicsModified = true; 543 physicsModified = true; 634 } else if (command == brCmd) { 544 } else if (command == brCmd) { 635 theParameters->SetBremsstrahlungTh(brCmd-> 545 theParameters->SetBremsstrahlungTh(brCmd->GetNewDoubleValue(newValue)); 636 physicsModified = true; 546 physicsModified = true; 637 } else if (command == br1Cmd) { 547 } else if (command == br1Cmd) { 638 theParameters->SetMuHadBremsstrahlungTh(br 548 theParameters->SetMuHadBremsstrahlungTh(br1Cmd->GetNewDoubleValue(newValue)); 639 physicsModified = true; 549 physicsModified = true; 640 } else if (command == labCmd) { 550 } else if (command == labCmd) { 641 theParameters->SetLambdaFactor(labCmd->Get 551 theParameters->SetLambdaFactor(labCmd->GetNewDoubleValue(newValue)); 642 physicsModified = true; 552 physicsModified = true; 643 } else if (command == mscfCmd) { 553 } else if (command == mscfCmd) { 644 theParameters->SetFactorForAngleLimit(mscf 554 theParameters->SetFactorForAngleLimit(mscfCmd->GetNewDoubleValue(newValue)); 645 } else if (command == angCmd) { 555 } else if (command == angCmd) { 646 theParameters->SetMscThetaLimit(angCmd->Ge 556 theParameters->SetMscThetaLimit(angCmd->GetNewDoubleValue(newValue)); 647 } else if (command == msceCmd) { 557 } else if (command == msceCmd) { 648 theParameters->SetMscEnergyLimit(msceCmd-> 558 theParameters->SetMscEnergyLimit(msceCmd->GetNewDoubleValue(newValue)); 649 } else if (command == nielCmd) { 559 } else if (command == nielCmd) { 650 theParameters->SetMaxNIELEnergy(nielCmd->G 560 theParameters->SetMaxNIELEnergy(nielCmd->GetNewDoubleValue(newValue)); 651 } else if (command == frCmd) { 561 } else if (command == frCmd) { 652 theParameters->SetMscRangeFactor(frCmd->Ge 562 theParameters->SetMscRangeFactor(frCmd->GetNewDoubleValue(newValue)); 653 physicsModified = true; 563 physicsModified = true; 654 } else if (command == fr1Cmd) { 564 } else if (command == fr1Cmd) { 655 theParameters->SetMscMuHadRangeFactor(fr1C 565 theParameters->SetMscMuHadRangeFactor(fr1Cmd->GetNewDoubleValue(newValue)); 656 physicsModified = true; 566 physicsModified = true; 657 } else if (command == fgCmd) { 567 } else if (command == fgCmd) { 658 theParameters->SetMscGeomFactor(fgCmd->Get 568 theParameters->SetMscGeomFactor(fgCmd->GetNewDoubleValue(newValue)); 659 physicsModified = true; 569 physicsModified = true; 660 } else if (command == skinCmd) { 570 } else if (command == skinCmd) { 661 theParameters->SetMscSkin(skinCmd->GetNewD 571 theParameters->SetMscSkin(skinCmd->GetNewDoubleValue(newValue)); 662 physicsModified = true; 572 physicsModified = true; 663 } else if (command == safCmd) { 573 } else if (command == safCmd) { 664 theParameters->SetMscSafetyFactor(safCmd-> 574 theParameters->SetMscSafetyFactor(safCmd->GetNewDoubleValue(newValue)); 665 } else if (command == llimCmd) { 575 } else if (command == llimCmd) { 666 theParameters->SetMscLambdaLimit(llimCmd-> 576 theParameters->SetMscLambdaLimit(llimCmd->GetNewDoubleValue(newValue)); 667 } else if (command == screCmd) { 577 } else if (command == screCmd) { 668 theParameters->SetScreeningFactor(screCmd- 578 theParameters->SetScreeningFactor(screCmd->GetNewDoubleValue(newValue)); 669 } else if (command == amCmd) { << 670 theParameters->SetNumberOfBinsPerDecade(am << 671 } else if (command == verCmd) { 579 } else if (command == verCmd) { 672 theParameters->SetVerbose(verCmd->GetNewIn 580 theParameters->SetVerbose(verCmd->GetNewIntValue(newValue)); 673 } else if (command == ver1Cmd) { 581 } else if (command == ver1Cmd) { 674 theParameters->SetVerbose(ver1Cmd->GetNewI 582 theParameters->SetVerbose(ver1Cmd->GetNewIntValue(newValue)); 675 } else if (command == ver2Cmd) { 583 } else if (command == ver2Cmd) { 676 theParameters->SetWorkerVerbose(ver2Cmd->G 584 theParameters->SetWorkerVerbose(ver2Cmd->GetNewIntValue(newValue)); 677 } else if (command == nFreeCmd) { << 678 theParameters->SetNumberForFreeVector(nFre << 679 } else if (command == dumpCmd) { 585 } else if (command == dumpCmd) { 680 theParameters->SetIsPrintedFlag(false); 586 theParameters->SetIsPrintedFlag(false); 681 theParameters->Dump(); 587 theParameters->Dump(); 682 } else if (command == transWithMscCmd) { << 683 G4TransportationWithMscType type = G4Trans << 684 if(newValue == "Disabled") { << 685 type = G4TransportationWithMscType::fDis << 686 } else if(newValue == "Enabled") { << 687 type = G4TransportationWithMscType::fEna << 688 } else if(newValue == "MultipleSteps") { << 689 type = G4TransportationWithMscType::fMul << 690 } else { << 691 G4ExceptionDescription ed; << 692 ed << " TransportationWithMsc type <" << << 693 G4Exception("G4EmParametersMessenger", " << 694 } << 695 theParameters->SetTransportationWithMsc(ty << 696 } else if (command == mscCmd || command == m 588 } else if (command == mscCmd || command == msc1Cmd) { 697 G4MscStepLimitType msctype = fUseSafety; 589 G4MscStepLimitType msctype = fUseSafety; 698 if(newValue == "Minimal") { 590 if(newValue == "Minimal") { 699 msctype = fMinimal; 591 msctype = fMinimal; 700 } else if(newValue == "UseDistanceToBounda 592 } else if(newValue == "UseDistanceToBoundary") { 701 msctype = fUseDistanceToBoundary; 593 msctype = fUseDistanceToBoundary; 702 } else if(newValue == "UseSafety") { 594 } else if(newValue == "UseSafety") { 703 msctype = fUseSafety; 595 msctype = fUseSafety; 704 } else if(newValue == "UseSafetyPlus") { 596 } else if(newValue == "UseSafetyPlus") { 705 msctype = fUseSafetyPlus; 597 msctype = fUseSafetyPlus; 706 } else { 598 } else { 707 G4ExceptionDescription ed; 599 G4ExceptionDescription ed; 708 ed << " StepLimit type <" << newValue << 600 ed << " StepLimit type <" << newValue << "> unknown!"; 709 G4Exception("G4EmParametersMessenger", " 601 G4Exception("G4EmParametersMessenger", "em0044", JustWarning, ed); 710 return; 602 return; 711 } 603 } 712 if (command == mscCmd) { 604 if (command == mscCmd) { 713 theParameters->SetMscStepLimitType(mscty 605 theParameters->SetMscStepLimitType(msctype); 714 } else { 606 } else { 715 theParameters->SetMscMuHadStepLimitType( 607 theParameters->SetMscMuHadStepLimitType(msctype); 716 } 608 } 717 physicsModified = true; 609 physicsModified = true; 718 } else if (command == nffCmd) { 610 } else if (command == nffCmd) { 719 G4NuclearFormfactorType x = fNoneNF; 611 G4NuclearFormfactorType x = fNoneNF; 720 if(newValue == "Exponential") { x = fExpon 612 if(newValue == "Exponential") { x = fExponentialNF; } 721 else if(newValue == "Gaussian") { x = fGau 613 else if(newValue == "Gaussian") { x = fGaussianNF; } 722 else if(newValue == "Flat") { x = fFlatNF; 614 else if(newValue == "Flat") { x = fFlatNF; } 723 else if(newValue != "None") { 615 else if(newValue != "None") { 724 G4ExceptionDescription ed; 616 G4ExceptionDescription ed; 725 ed << " NuclearFormFactor type <" << new 617 ed << " NuclearFormFactor type <" << newValue << "> unknown!"; 726 G4Exception("G4EmParametersMessenger", " 618 G4Exception("G4EmParametersMessenger", "em0044", JustWarning, ed); 727 return; 619 return; 728 } 620 } 729 theParameters->SetNuclearFormfactorType(x) 621 theParameters->SetNuclearFormfactorType(x); 730 } else if (command == ssCmd) { 622 } else if (command == ssCmd) { 731 G4eSingleScatteringType x = fWVI; 623 G4eSingleScatteringType x = fWVI; 732 if(newValue == "DPWA") { x = fDPWA; } 624 if(newValue == "DPWA") { x = fDPWA; } 733 else if(newValue == "Mott") { x = fMott; } 625 else if(newValue == "Mott") { x = fMott; } 734 else if(newValue != "WVI") { 626 else if(newValue != "WVI") { 735 G4ExceptionDescription ed; 627 G4ExceptionDescription ed; 736 ed << " G4eSingleScatteringType type <" 628 ed << " G4eSingleScatteringType type <" << newValue << "> unknown!"; 737 G4Exception("G4EmParametersMessenger", " 629 G4Exception("G4EmParametersMessenger", "em0044", JustWarning, ed); 738 return; 630 return; 739 } 631 } 740 theParameters->SetSingleScatteringType(x); 632 theParameters->SetSingleScatteringType(x); 741 } else if (command == fluc1Cmd) { << 742 G4EmFluctuationType x = fUniversalFluctuat << 743 if(newValue == "Dummy") { x = fDummyFluctu << 744 else if(newValue == "Urban") { x = fUrbanF << 745 theParameters->SetFluctuationType(x); << 746 } else if (command == posiCmd) { << 747 G4PositronAtRestModelType x = fSimplePosit << 748 if (newValue == "Allison") { x = fAllisonP << 749 else if (newValue == "OrePowell") { x = fO << 750 else if (newValue == "OrePowellPolar") { x << 751 theParameters->SetPositronAtRestModelType( << 752 } else if ( command==tripletCmd ) { 633 } else if ( command==tripletCmd ) { 753 theParameters->SetConversionType(tripletCm 634 theParameters->SetConversionType(tripletCmd->GetNewIntValue(newValue)); 754 } else if ( command==onIsolatedCmd ) { 635 } else if ( command==onIsolatedCmd ) { 755 theParameters->SetOnIsolated(onIsolatedCmd 636 theParameters->SetOnIsolated(onIsolatedCmd->GetNewBoolValue(newValue)); 756 physicsModified = true; 637 physicsModified = true; 757 } 638 } 758 639 759 if(physicsModified) { 640 if(physicsModified) { 760 G4UImanager::GetUIpointer()->ApplyCommand( 641 G4UImanager::GetUIpointer()->ApplyCommand("/run/physicsModified"); 761 } 642 } 762 } 643 } 763 644 764 //....oooOO0OOooo........oooOO0OOooo........oo 645 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 765 646