| 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 /// \file biasing/ReverseMC01/src/G4AdjointPhy 26 /// \file biasing/ReverseMC01/src/G4AdjointPhysicsMessenger.cc
27 /// \brief Implementation of the G4AdjointPhys 27 /// \brief Implementation of the G4AdjointPhysicsMessenger class
28 // 28 //
>> 29 // $Id: G4AdjointPhysicsMessenger.cc 71233 2013-06-12 13:16:05Z gcosmo $
29 // 30 //
30 ////////////////////////////////////////////// 31 //////////////////////////////////////////////////////////////
31 // Class Name: G4AdjointPhysicsMes 32 // Class Name: G4AdjointPhysicsMessenger
32 // Author: L. Desorgher 33 // Author: L. Desorgher
33 // Organisation: SpaceIT GmbH 34 // Organisation: SpaceIT GmbH
34 // Contract: ESA contract 21435/ 35 // Contract: ESA contract 21435/08/NL/AT
35 // Customer: ESA/ESTEC 36 // Customer: ESA/ESTEC
36 ////////////////////////////////////////////// 37 //////////////////////////////////////////////////////////////
37 38
38 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
39 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
40 41
41 #include "G4AdjointPhysicsMessenger.hh" 42 #include "G4AdjointPhysicsMessenger.hh"
42 << 43 #include "G4UIdirectory.hh"
43 #include "G4AdjointPhysicsList.hh" <<
44 #include "G4UIcmdWith3VectorAndUnit.hh" <<
45 #include "G4UIcmdWithABool.hh" 44 #include "G4UIcmdWithABool.hh"
46 #include "G4UIcmdWithADouble.hh" <<
47 #include "G4UIcmdWithADoubleAndUnit.hh" <<
48 #include "G4UIcmdWithAString.hh" <<
49 #include "G4UIcmdWithAnInteger.hh" 45 #include "G4UIcmdWithAnInteger.hh"
>> 46 #include "G4UIcmdWithADoubleAndUnit.hh"
>> 47 #include "G4UIcmdWithADouble.hh"
50 #include "G4UIcmdWithoutParameter.hh" 48 #include "G4UIcmdWithoutParameter.hh"
51 #include "G4UIdirectory.hh" << 49 #include "G4UIcmdWithAString.hh"
52 #include "G4UnitsTable.hh" 50 #include "G4UnitsTable.hh"
>> 51 #include "G4AdjointPhysicsList.hh"
>> 52 #include "G4UIcmdWith3VectorAndUnit.hh"
53 53
54 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55 55
56 G4AdjointPhysicsMessenger::G4AdjointPhysicsMes << 56 G4AdjointPhysicsMessenger::G4AdjointPhysicsMessenger(
57 : G4UImessenger(), << 57 G4AdjointPhysicsList* pPhysicsList)
58 fPhysicsList(pPhysicsList), << 58 : G4UImessenger(),
59 fPhysicsDir(0), << 59 fPhysicsList(pPhysicsList),
60 fUsepIonisationCmd(0), << 60 fPhysicsDir(0),
61 fUseBremCmd(0), << 61 fUsepIonisationCmd(0),
62 fUseComptonCmd(0), << 62 fUseBremCmd(0),
63 fUseMSCmd(0), << 63 fUseComptonCmd(0),
64 fUsePEEffectCmd(0), << 64 fUseMSCmd(0),
65 fUseGammaConversionCmd(0), << 65 fUsePEEffectCmd(0),
66 fUseEgainFluctuationCmd(0), << 66 fUseGammaConversionCmd(0),
67 fSetEminAdjModelsCmd(0), << 67 fUseEgainFluctuationCmd(0),
68 fSetEmaxAdjModelsCmd(0) << 68 fSetEminAdjModelsCmd(0),
69 { << 69 fSetEmaxAdjModelsCmd(0)
>> 70 {
70 fPhysicsDir = new G4UIdirectory("/adjoint_ph 71 fPhysicsDir = new G4UIdirectory("/adjoint_physics/");
71 72
72 fPhysicsDir->SetGuidance("Definition of adjo << 73 fPhysicsDir->SetGuidance(
>> 74 "Definition of adjoint and forward physics processes");
73 //------- 75 //-------
74 fUsepIonisationCmd = new G4UIcmdWithABool("/ << 76 fUsepIonisationCmd = new G4UIcmdWithABool(
75 fUsepIonisationCmd->SetGuidance("If true (fa << 77 "/adjoint_physics/UseProtonIonisation",this);
>> 78 fUsepIonisationCmd->SetGuidance(
>> 79 "If true (false) the proton ionisation is (not) considered");
76 fUsepIonisationCmd->AvailableForStates(G4Sta 80 fUsepIonisationCmd->AvailableForStates(G4State_PreInit);
77 << 81
78 fUseBremCmd = new G4UIcmdWithABool("/adjoint << 82 fUseBremCmd = new G4UIcmdWithABool("/adjoint_physics/UseBremsstrahlung",this);
79 fUseBremCmd->SetGuidance("If true (false) th << 83 fUseBremCmd->SetGuidance(
80 fUseBremCmd->AvailableForStates(G4State_PreI << 84 "If true (false) the bremsstrahlung process is (not) considered");
81 << 85 fUseBremCmd->AvailableForStates(G4State_PreInit);
82 fUseComptonCmd = new G4UIcmdWithABool("/adjo << 86
83 fUseComptonCmd->SetGuidance("If true (false) << 87 fUseComptonCmd = new G4UIcmdWithABool("/adjoint_physics/UseCompton",this);
84 fUseComptonCmd->AvailableForStates(G4State_P << 88 fUseComptonCmd->SetGuidance(
85 << 89 "If true (false) the Compton scattering is (not) considered");
86 fUseMSCmd = new G4UIcmdWithABool("/adjoint_p << 90 fUseComptonCmd->AvailableForStates(G4State_PreInit);
87 fUseMSCmd->SetGuidance("If true (false) the << 91
88 fUseMSCmd->AvailableForStates(G4State_PreIni << 92 fUseMSCmd = new G4UIcmdWithABool("/adjoint_physics/UseMS",this);
89 << 93 fUseMSCmd->SetGuidance(
90 fUseEgainFluctuationCmd = new G4UIcmdWithABo << 94 "If true (false) the continuous multiple scattering is (not) considered");
>> 95 fUseMSCmd->AvailableForStates(G4State_PreInit);
>> 96
>> 97 fUseEgainFluctuationCmd = new G4UIcmdWithABool(
>> 98 "/adjoint_physics/UseEgainElossFluctuation",this);
91 fUseEgainFluctuationCmd->SetGuidance( 99 fUseEgainFluctuationCmd->SetGuidance(
92 "Switch on/off the fluctation for continuo 100 "Switch on/off the fluctation for continuous energy gain/loss");
93 fUseEgainFluctuationCmd->AvailableForStates( << 101 fUseEgainFluctuationCmd->AvailableForStates(G4State_PreInit);
94 << 102
95 fUsePEEffectCmd = new G4UIcmdWithABool("/adj << 103 fUsePEEffectCmd = new G4UIcmdWithABool("/adjoint_physics/UsePEEffect",this);
96 fUsePEEffectCmd->AvailableForStates(G4State_ << 104 fUsePEEffectCmd->AvailableForStates(G4State_PreInit);
97 fUsePEEffectCmd->SetGuidance("If true (false << 105 fUsePEEffectCmd->SetGuidance(
98 << 106 "If true (false) the photo electric effect is (not) considered");
99 fUseGammaConversionCmd = new G4UIcmdWithABoo << 107
>> 108 fUseGammaConversionCmd = new G4UIcmdWithABool(
>> 109 "/adjoint_physics/UseGammaConversion",this);
100 fUseGammaConversionCmd->AvailableForStates(G 110 fUseGammaConversionCmd->AvailableForStates(G4State_PreInit);
101 fUseGammaConversionCmd->SetGuidance("If true << 111 fUseGammaConversionCmd->SetGuidance(
102 << 112 "If true the fwd gamma pair conversion is considered");
103 fSetEminAdjModelsCmd = << 113
104 new G4UIcmdWithADoubleAndUnit("/adjoint_ph << 114 fSetEminAdjModelsCmd = new G4UIcmdWithADoubleAndUnit(
105 fSetEminAdjModelsCmd->SetGuidance("Set the m << 115 "/adjoint_physics/SetEminForAdjointModels",this);
106 fSetEminAdjModelsCmd->SetParameterName("Emin << 116 fSetEminAdjModelsCmd->SetGuidance(
>> 117 "Set the minimum energy of the adjoint models");
>> 118 fSetEminAdjModelsCmd->SetParameterName("Emin",false);
107 fSetEminAdjModelsCmd->SetUnitCategory("Energ 119 fSetEminAdjModelsCmd->SetUnitCategory("Energy");
108 fSetEminAdjModelsCmd->AvailableForStates(G4S 120 fSetEminAdjModelsCmd->AvailableForStates(G4State_PreInit);
109 << 121
110 fSetEmaxAdjModelsCmd = << 122 fSetEmaxAdjModelsCmd = new G4UIcmdWithADoubleAndUnit(
111 new G4UIcmdWithADoubleAndUnit("/adjoint_ph << 123 "/adjoint_physics/SetEmaxForAdjointModels",this);
112 fSetEmaxAdjModelsCmd->SetGuidance("Set the m << 124 fSetEmaxAdjModelsCmd->SetGuidance(
113 fSetEmaxAdjModelsCmd->SetParameterName("Emax << 125 "Set the minimum energy of the adjoint models.");
>> 126 fSetEmaxAdjModelsCmd->SetParameterName("Emax",false);
114 fSetEmaxAdjModelsCmd->SetUnitCategory("Energ 127 fSetEmaxAdjModelsCmd->SetUnitCategory("Energy");
115 fSetEmaxAdjModelsCmd->AvailableForStates(G4S 128 fSetEmaxAdjModelsCmd->AvailableForStates(G4State_PreInit);
116 } 129 }
117 130
118 //....oooOO0OOooo........oooOO0OOooo........oo 131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
119 132
120 G4AdjointPhysicsMessenger::~G4AdjointPhysicsMe 133 G4AdjointPhysicsMessenger::~G4AdjointPhysicsMessenger()
121 { 134 {
122 delete fUsepIonisationCmd; 135 delete fUsepIonisationCmd;
123 delete fUseBremCmd; 136 delete fUseBremCmd;
124 delete fUseComptonCmd; 137 delete fUseComptonCmd;
125 delete fUseMSCmd; 138 delete fUseMSCmd;
126 delete fUsePEEffectCmd; 139 delete fUsePEEffectCmd;
127 delete fUseGammaConversionCmd; 140 delete fUseGammaConversionCmd;
128 delete fUseEgainFluctuationCmd; 141 delete fUseEgainFluctuationCmd;
129 delete fSetEminAdjModelsCmd; 142 delete fSetEminAdjModelsCmd;
130 delete fSetEmaxAdjModelsCmd; 143 delete fSetEmaxAdjModelsCmd;
131 } 144 }
132 145
133 //....oooOO0OOooo........oooOO0OOooo........oo 146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
134 147
135 void G4AdjointPhysicsMessenger::SetNewValue(G4 << 148 void G4AdjointPhysicsMessenger::SetNewValue(G4UIcommand* command,
>> 149 G4String newValue)
136 { 150 {
137 if (command == fUsepIonisationCmd) { << 151 if ( command==fUsepIonisationCmd){
138 fPhysicsList->SetUseProtonIonisation(fUsep << 152 fPhysicsList->SetUseProtonIonisation(
>> 153 fUsepIonisationCmd->GetNewBoolValue(newValue));
>> 154
139 } 155 }
140 else if (command == fUseBremCmd) { << 156 else if ( command==fUseBremCmd){
141 fPhysicsList->SetUseBrem(fUseBremCmd->GetN << 157 fPhysicsList->SetUseBrem(fUseBremCmd->GetNewBoolValue(newValue));
142 } 158 }
143 else if (command == fUseComptonCmd) { << 159 else if ( command==fUseComptonCmd){
144 fPhysicsList->SetUseCompton(fUseComptonCmd << 160 fPhysicsList->SetUseCompton(fUseComptonCmd->GetNewBoolValue(newValue));
145 } 161 }
146 else if (command == fUseMSCmd) { << 162 else if ( command==fUseMSCmd){
147 fPhysicsList->SetUseMS(fUseMSCmd->GetNewBo 163 fPhysicsList->SetUseMS(fUseMSCmd->GetNewBoolValue(newValue));
148 } 164 }
149 else if (command == fUsePEEffectCmd) { << 165 else if ( command==fUsePEEffectCmd){
150 fPhysicsList->SetUsePEEffect(fUsePEEffectC 166 fPhysicsList->SetUsePEEffect(fUsePEEffectCmd->GetNewBoolValue(newValue));
151 } 167 }
152 else if (command == fUseGammaConversionCmd) << 168 else if ( command==fUseGammaConversionCmd){
153 fPhysicsList->SetUseGammaConversion(fUseGa << 169 fPhysicsList->SetUseGammaConversion(
154 } << 170 fUseGammaConversionCmd->GetNewBoolValue(newValue));
155 else if (command == fUseEgainFluctuationCmd) << 171 }
156 fPhysicsList->SetUseEgainFluctuation(fUseE << 172 else if ( command==fUseEgainFluctuationCmd){
157 } << 173 fPhysicsList->SetUseEgainFluctuation(
158 << 174 fUseEgainFluctuationCmd->GetNewBoolValue(newValue));
159 else if (command == fSetEminAdjModelsCmd) { << 175 }
160 fPhysicsList->SetEminAdjModels(fSetEminAdj << 176
161 } << 177 else if ( command== fSetEminAdjModelsCmd){
162 else if (command == fSetEmaxAdjModelsCmd) { << 178 fPhysicsList->SetEminAdjModels(
163 fPhysicsList->SetEmaxAdjModels(fSetEmaxAdj << 179 fSetEminAdjModelsCmd->GetNewDoubleValue(newValue));
>> 180 }
>> 181 else if ( command== fSetEmaxAdjModelsCmd){
>> 182 fPhysicsList->SetEmaxAdjModels(
>> 183 fSetEmaxAdjModelsCmd->GetNewDoubleValue(newValue));
164 } 184 }
165 } 185 }
166 186
167 //....oooOO0OOooo........oooOO0OOooo........oo 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 188
168 189