source: trunk/source/physics_lists/builders/src/G4OpticalPhysicsMessenger.cc@ 1344

Last change on this file since 1344 was 1340, checked in by garnier, 15 years ago

update ti head
File size: 12.1 KB
Line 
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#include "G4OpticalPhysicsMessenger.hh"
27#include "G4OpticalPhysics.hh"
28
29#include "G4UIcommand.hh"
30#include "G4UIdirectory.hh"
31
32#include "G4UIcommand.hh"
33#include "G4UIdirectory.hh"
34#include "G4UIcmdWithABool.hh"
35#include "G4UIcmdWithAString.hh"
36#include "G4UIcmdWithADouble.hh"
37#include "G4UIcmdWithAnInteger.hh"
38#include "G4UIcmdWithADoubleAndUnit.hh"
39
40#include "G4ParticleTable.hh"
41#include "G4ProcessVector.hh"
42#include "G4ProcessManager.hh"
43#include "G4ParticleDefinition.hh"
44
45//_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
46G4OpticalPhysicsMessenger::G4OpticalPhysicsMessenger(
47 G4OpticalPhysics* opticalPhysics)
48 : G4UImessenger(),
49 fOpticalPhysics(opticalPhysics),
50 fSelectedProcess(0),
51 fSelectOpProcessCmd(0),
52 fSetOpProcessActivationCmd(0),
53 fSetOpProcessVerboseCmd(0),
54 fSetCerenkovMaxPhotonsCmd(0),
55 fSetCerenkovMaxBetaChangeCmd(0),
56 fSetScintillationYieldFactorCmd(0),
57 fSetScintillationByParticleTypeCmd(0),
58 fSetOpticalSurfaceModelCmd(0),
59 fSetWLSTimeProfileCmd(0),
60 fSetTrackSecondariesFirstCmd(0)
61{
62/// Standard constructor
63
64 fDir = new G4UIdirectory("/optics_engine/");
65 fDir->SetGuidance("Commands related to the optical physics simulation engine.");
66
67 fSelectOpProcessCmd
68 = new G4UIcmdWithAString("/optics_engine/selectOpProcess", this);
69 fSelectOpProcessCmd
70 ->SetGuidance("Select optical process for applying verbose/activation/trackfirst commands");
71 fSelectOpProcessCmd->SetParameterName("OpProcess", false);
72 fSelectOpProcessCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
73
74 fSetOpProcessActivationCmd
75 = new G4UIcmdWithABool("/optics_engine/setOpProcessActivation", this);
76 fSetOpProcessActivationCmd->SetGuidance("Activate/Inactivate selected optical process");
77 fSetOpProcessActivationCmd->SetParameterName("OpProcessActivation", false);
78 fSetOpProcessActivationCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
79
80 fSetOpProcessVerboseCmd
81 = new G4UIcmdWithAnInteger("/optics_engine/setOpProcessVerbose", this);
82 fSetOpProcessVerboseCmd->SetGuidance("Set verbosity level for selected optical process");
83 fSetOpProcessVerboseCmd->SetParameterName("OpProcessVerbose", false);
84 fSetOpProcessVerboseCmd->SetRange("OpProcessVerbose>=0");
85 fSetOpProcessVerboseCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
86
87 fSetCerenkovMaxPhotonsCmd
88 = new G4UIcmdWithAnInteger("/optics_engine/setCerenkovMaxPhotons", this);
89 fSetCerenkovMaxPhotonsCmd->SetGuidance("Set maximum number of photons per step");
90 fSetCerenkovMaxPhotonsCmd->SetParameterName("CerenkovMaxPhotons", false);
91 fSetCerenkovMaxPhotonsCmd->SetRange("CerenkovMaxPhotons>=0");
92 fSetCerenkovMaxPhotonsCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
93
94 fSetCerenkovMaxBetaChangeCmd
95 = new G4UIcmdWithADouble("/optics_engine/setCerenkovMaxBetaChange", this);
96 fSetCerenkovMaxBetaChangeCmd
97 ->SetGuidance("Set maximum change of beta of parent particle per step");
98 fSetCerenkovMaxBetaChangeCmd->SetParameterName("CerenkovMaxBetaChange", false);
99 fSetCerenkovMaxBetaChangeCmd->SetRange("CerenkovMaxBetaChange>=0");
100 fSetCerenkovMaxBetaChangeCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
101
102 fSetScintillationYieldFactorCmd
103 = new G4UIcmdWithADouble("/optics_engine/setScintillationYieldFactor", this);
104 fSetScintillationYieldFactorCmd->SetGuidance("Set scintillation yield factor");
105 fSetScintillationYieldFactorCmd->SetParameterName("ScintillationYieldFactor", false);
106 fSetScintillationYieldFactorCmd->SetRange("ScintillationYieldFactor>=0");
107 fSetScintillationYieldFactorCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
108
109 fSetScintillationByParticleTypeCmd
110 = new G4UIcmdWithABool("/optics_engine/setScintillationByParticleType", this);
111 fSetScintillationByParticleTypeCmd->SetGuidance("Activate/Inactivate scintillation process by particle type");
112 fSetScintillationByParticleTypeCmd->SetParameterName("ScintillationByParticleTypeActivation", false);
113 fSetScintillationByParticleTypeCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
114
115 fSetOpticalSurfaceModelCmd
116 = new G4UIcmdWithAString("/optics_engine/setOpticalSurfaceModel", this);
117 fSetOpticalSurfaceModelCmd
118 ->SetGuidance("Set optical surface model (glisur or unified)");
119 fSetOpticalSurfaceModelCmd->SetParameterName("OpticalSurfaceModel", false);
120 fSetOpticalSurfaceModelCmd->SetCandidates("glisur unified");
121 fSetOpticalSurfaceModelCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
122
123 fSetWLSTimeProfileCmd
124 = new G4UIcmdWithAString("/optics_engine/setWLSTimeProfile", this);
125 fSetWLSTimeProfileCmd
126 ->SetGuidance("Set the WLS time profile (delta or exponential)");
127 fSetWLSTimeProfileCmd->SetParameterName("WLSTimeProfile", false);
128 fSetWLSTimeProfileCmd->SetCandidates("delta exponential");
129 fSetWLSTimeProfileCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
130
131 fSetTrackSecondariesFirstCmd
132 = new G4UIcmdWithABool("/optics_engine/setTrackSecondariesFirst", this);
133 fSetTrackSecondariesFirstCmd
134 ->SetGuidance("Set option to track secondaries before finishing their parent track");
135 fSetTrackSecondariesFirstCmd->SetParameterName("TrackSecondariesFirst", false);
136 fSetTrackSecondariesFirstCmd->AvailableForStates(G4State_Idle,G4State_GeomClosed,G4State_EventProc);
137}
138
139//_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
140G4OpticalPhysicsMessenger::~G4OpticalPhysicsMessenger()
141{
142/// Destructor
143 delete fDir;
144 delete fSelectOpProcessCmd;
145 delete fSetOpProcessActivationCmd;
146 delete fSetOpProcessVerboseCmd;
147 delete fSetCerenkovMaxPhotonsCmd;
148 delete fSetCerenkovMaxBetaChangeCmd;
149 delete fSetScintillationYieldFactorCmd;
150 delete fSetScintillationByParticleTypeCmd;
151 delete fSetOpticalSurfaceModelCmd;
152 delete fSetWLSTimeProfileCmd;
153 delete fSetTrackSecondariesFirstCmd;
154}
155
156//_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
157void G4OpticalPhysicsMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
158{
159/// Apply command to the associated object.
160
161 if (command == fSelectOpProcessCmd) {
162 if ( newValue == "Cerenkov" ) fSelectedProcess = fOpticalPhysics->
163 GetCerenkovProcess();
164 else if ( newValue == "Scintillation" ) fSelectedProcess = fOpticalPhysics->
165 GetScintillationProcess();
166 else if ( newValue == "OpAbsorption" ) fSelectedProcess = fOpticalPhysics->
167 GetOpAbsorptionProcess();
168 else if ( newValue == "OpRayleigh" ) fSelectedProcess = fOpticalPhysics->
169 GetOpRayleighProcess();
170 else if ( newValue == "OpMieHG" ) fSelectedProcess = fOpticalPhysics->
171 GetOpMieHGProcess();
172 else if ( newValue == "OpBoundary" ) fSelectedProcess = fOpticalPhysics->
173 GetOpBoundaryProcess();
174 else if ( newValue == "OpWLS" ) fSelectedProcess = fOpticalPhysics->
175 GetOpWLSProcess();
176 }
177 else if (command == fSetOpProcessActivationCmd) {
178 G4ParticleTable* theParticleTable = G4ParticleTable::GetParticleTable();
179 G4ParticleTable::G4PTblDicIterator* theParticleIterator =
180 theParticleTable->GetIterator();
181 if ( fSelectedProcess ) {
182 theParticleIterator->reset();
183 while( (*theParticleIterator)() ) {
184 G4ParticleDefinition* particle = theParticleIterator->value();
185 G4ProcessManager* processManager = particle->GetProcessManager();
186 G4ProcessVector* processVector = processManager->GetProcessList();
187
188 // activate or inactivate processes according to
189 for (G4int i=0; i<processVector->length(); i++) {
190 if ( (*processVector)[i] != fSelectedProcess ) continue;
191 G4cout << "Setting process activation: "
192 << fSelectedProcess->GetProcessName() << " "
193 << fSetOpProcessActivationCmd->GetNewBoolValue(newValue)
194 << G4endl;
195 processManager->SetProcessActivation(fSelectedProcess,
196 fSetOpProcessActivationCmd->GetNewBoolValue(newValue));
197 }
198 }
199 }
200 }
201 else if (command == fSetOpProcessVerboseCmd) {
202 if ( fSelectedProcess )
203 fSelectedProcess->SetVerboseLevel
204 (fSetOpProcessVerboseCmd->GetNewIntValue(newValue));
205 }
206 else if (command == fSetCerenkovMaxPhotonsCmd) {
207 fOpticalPhysics
208 ->SetMaxNumPhotonsPerStep(
209 fSetCerenkovMaxPhotonsCmd->GetNewIntValue(newValue));
210 }
211 else if (command == fSetCerenkovMaxBetaChangeCmd) {
212 fOpticalPhysics
213 ->SetMaxBetaChangePerStep(
214 fSetCerenkovMaxBetaChangeCmd->GetNewDoubleValue(newValue));
215 }
216 else if (command == fSetScintillationYieldFactorCmd) {
217 fOpticalPhysics
218 ->SetScintillationYieldFactor(
219 fSetScintillationYieldFactorCmd->GetNewDoubleValue(newValue));
220 }
221 else if (command == fSetScintillationByParticleTypeCmd) {
222 fOpticalPhysics
223 ->SetScintillationByParticleType(
224 fSetScintillationByParticleTypeCmd->GetNewBoolValue(newValue));
225 }
226 else if (command == fSetOpticalSurfaceModelCmd) {
227 if ( newValue == "glisur" ) {
228 fOpticalPhysics
229 ->SetOpticalSurfaceModel(glisur);
230 }
231 if ( newValue == "unified" ) {
232 fOpticalPhysics
233 ->SetOpticalSurfaceModel(unified);
234 }
235 }
236 else if (command == fSetWLSTimeProfileCmd) {
237 if ( newValue == "delta" ) {
238 fOpticalPhysics
239 ->SetWLSTimeProfile("delta"); }
240 if ( newValue == "exponential" ) {
241 fOpticalPhysics
242 ->SetWLSTimeProfile("exponential");
243 }
244 }
245 else if (command == fSetTrackSecondariesFirstCmd) {
246 if ( fSelectedProcess ) {
247 G4Scintillation* scintillation =
248 static_cast<G4Scintillation*>(fSelectedProcess);
249 if (scintillation) scintillation ->
250 SetTrackSecondariesFirst(
251 fSetTrackSecondariesFirstCmd->GetNewBoolValue(newValue));
252 G4Cerenkov* cerenkov =
253 static_cast<G4Cerenkov*>(fSelectedProcess);
254 if (cerenkov) cerenkov ->
255 SetTrackSecondariesFirst(
256 fSetTrackSecondariesFirstCmd->GetNewBoolValue(newValue));
257 }
258 }
259}
Note: See TracBrowser for help on using the repository browser.