Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 /// \file biasing/ReverseMC01/src/G4AdjointPhysicsMessenger.cc 27 /// \brief Implementation of the G4AdjointPhysicsMessenger class 28 // 29 // 30 ////////////////////////////////////////////////////////////// 31 // Class Name: G4AdjointPhysicsMessenger 32 // Author: L. Desorgher 33 // Organisation: SpaceIT GmbH 34 // Contract: ESA contract 21435/08/NL/AT 35 // Customer: ESA/ESTEC 36 ////////////////////////////////////////////////////////////// 37 38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 40 41 #include "G4AdjointPhysicsMessenger.hh" 42 43 #include "G4AdjointPhysicsList.hh" 44 #include "G4UIcmdWith3VectorAndUnit.hh" 45 #include "G4UIcmdWithABool.hh" 46 #include "G4UIcmdWithADouble.hh" 47 #include "G4UIcmdWithADoubleAndUnit.hh" 48 #include "G4UIcmdWithAString.hh" 49 #include "G4UIcmdWithAnInteger.hh" 50 #include "G4UIcmdWithoutParameter.hh" 51 #include "G4UIdirectory.hh" 52 #include "G4UnitsTable.hh" 53 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 55 56 G4AdjointPhysicsMessenger::G4AdjointPhysicsMessenger(G4AdjointPhysicsList* pPhysicsList) 57 : G4UImessenger(), 58 fPhysicsList(pPhysicsList), 59 fPhysicsDir(0), 60 fUsepIonisationCmd(0), 61 fUseBremCmd(0), 62 fUseComptonCmd(0), 63 fUseMSCmd(0), 64 fUsePEEffectCmd(0), 65 fUseGammaConversionCmd(0), 66 fUseEgainFluctuationCmd(0), 67 fSetEminAdjModelsCmd(0), 68 fSetEmaxAdjModelsCmd(0) 69 { 70 fPhysicsDir = new G4UIdirectory("/adjoint_physics/"); 71 72 fPhysicsDir->SetGuidance("Definition of adjoint and forward physics processes"); 73 //------- 74 fUsepIonisationCmd = new G4UIcmdWithABool("/adjoint_physics/UseProtonIonisation", this); 75 fUsepIonisationCmd->SetGuidance("If true (false) the proton ionisation is (not) considered"); 76 fUsepIonisationCmd->AvailableForStates(G4State_PreInit); 77 78 fUseBremCmd = new G4UIcmdWithABool("/adjoint_physics/UseBremsstrahlung", this); 79 fUseBremCmd->SetGuidance("If true (false) the bremsstrahlung process is (not) considered"); 80 fUseBremCmd->AvailableForStates(G4State_PreInit); 81 82 fUseComptonCmd = new G4UIcmdWithABool("/adjoint_physics/UseCompton", this); 83 fUseComptonCmd->SetGuidance("If true (false) the Compton scattering is (not) considered"); 84 fUseComptonCmd->AvailableForStates(G4State_PreInit); 85 86 fUseMSCmd = new G4UIcmdWithABool("/adjoint_physics/UseMS", this); 87 fUseMSCmd->SetGuidance("If true (false) the continuous multiple scattering is (not) considered"); 88 fUseMSCmd->AvailableForStates(G4State_PreInit); 89 90 fUseEgainFluctuationCmd = new G4UIcmdWithABool("/adjoint_physics/UseEgainElossFluctuation", this); 91 fUseEgainFluctuationCmd->SetGuidance( 92 "Switch on/off the fluctation for continuous energy gain/loss"); 93 fUseEgainFluctuationCmd->AvailableForStates(G4State_PreInit); 94 95 fUsePEEffectCmd = new G4UIcmdWithABool("/adjoint_physics/UsePEEffect", this); 96 fUsePEEffectCmd->AvailableForStates(G4State_PreInit); 97 fUsePEEffectCmd->SetGuidance("If true (false) the photo electric effect is (not) considered"); 98 99 fUseGammaConversionCmd = new G4UIcmdWithABool("/adjoint_physics/UseGammaConversion", this); 100 fUseGammaConversionCmd->AvailableForStates(G4State_PreInit); 101 fUseGammaConversionCmd->SetGuidance("If true the fwd gamma pair conversion is considered"); 102 103 fSetEminAdjModelsCmd = 104 new G4UIcmdWithADoubleAndUnit("/adjoint_physics/SetEminForAdjointModels", this); 105 fSetEminAdjModelsCmd->SetGuidance("Set the minimum energy of the adjoint models"); 106 fSetEminAdjModelsCmd->SetParameterName("Emin", false); 107 fSetEminAdjModelsCmd->SetUnitCategory("Energy"); 108 fSetEminAdjModelsCmd->AvailableForStates(G4State_PreInit); 109 110 fSetEmaxAdjModelsCmd = 111 new G4UIcmdWithADoubleAndUnit("/adjoint_physics/SetEmaxForAdjointModels", this); 112 fSetEmaxAdjModelsCmd->SetGuidance("Set the minimum energy of the adjoint models."); 113 fSetEmaxAdjModelsCmd->SetParameterName("Emax", false); 114 fSetEmaxAdjModelsCmd->SetUnitCategory("Energy"); 115 fSetEmaxAdjModelsCmd->AvailableForStates(G4State_PreInit); 116 } 117 118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 119 120 G4AdjointPhysicsMessenger::~G4AdjointPhysicsMessenger() 121 { 122 delete fUsepIonisationCmd; 123 delete fUseBremCmd; 124 delete fUseComptonCmd; 125 delete fUseMSCmd; 126 delete fUsePEEffectCmd; 127 delete fUseGammaConversionCmd; 128 delete fUseEgainFluctuationCmd; 129 delete fSetEminAdjModelsCmd; 130 delete fSetEmaxAdjModelsCmd; 131 } 132 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 135 void G4AdjointPhysicsMessenger::SetNewValue(G4UIcommand* command, G4String newValue) 136 { 137 if (command == fUsepIonisationCmd) { 138 fPhysicsList->SetUseProtonIonisation(fUsepIonisationCmd->GetNewBoolValue(newValue)); 139 } 140 else if (command == fUseBremCmd) { 141 fPhysicsList->SetUseBrem(fUseBremCmd->GetNewBoolValue(newValue)); 142 } 143 else if (command == fUseComptonCmd) { 144 fPhysicsList->SetUseCompton(fUseComptonCmd->GetNewBoolValue(newValue)); 145 } 146 else if (command == fUseMSCmd) { 147 fPhysicsList->SetUseMS(fUseMSCmd->GetNewBoolValue(newValue)); 148 } 149 else if (command == fUsePEEffectCmd) { 150 fPhysicsList->SetUsePEEffect(fUsePEEffectCmd->GetNewBoolValue(newValue)); 151 } 152 else if (command == fUseGammaConversionCmd) { 153 fPhysicsList->SetUseGammaConversion(fUseGammaConversionCmd->GetNewBoolValue(newValue)); 154 } 155 else if (command == fUseEgainFluctuationCmd) { 156 fPhysicsList->SetUseEgainFluctuation(fUseEgainFluctuationCmd->GetNewBoolValue(newValue)); 157 } 158 159 else if (command == fSetEminAdjModelsCmd) { 160 fPhysicsList->SetEminAdjModels(fSetEminAdjModelsCmd->GetNewDoubleValue(newValue)); 161 } 162 else if (command == fSetEmaxAdjModelsCmd) { 163 fPhysicsList->SetEmaxAdjModels(fSetEmaxAdjModelsCmd->GetNewDoubleValue(newValue)); 164 } 165 } 166 167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 168