| 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 hadronic/Hadr01/src/DetectorMessenge 26 /// \file hadronic/Hadr01/src/DetectorMessenger.cc
27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class
28 // 28 //
>> 29 // $Id$
29 // 30 //
30 ////////////////////////////////////////////// 31 /////////////////////////////////////////////////////////////////////////
31 // 32 //
32 // DetectorMessenger 33 // DetectorMessenger
33 // 34 //
34 // Created: 31.01.03 V.Ivanchenko 35 // Created: 31.01.03 V.Ivanchenko
35 // 36 //
36 // Modified: 37 // Modified:
37 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko)
38 // 16.11.2006 Add beamCmd (V.Ivanchenko) 39 // 16.11.2006 Add beamCmd (V.Ivanchenko)
39 // 40 //
40 ////////////////////////////////////////////// 41 ////////////////////////////////////////////////////////////////////////
41 // 42 //
42 43
43 #include "DetectorMessenger.hh" 44 #include "DetectorMessenger.hh"
44 45
45 #include "DetectorConstruction.hh" 46 #include "DetectorConstruction.hh"
46 #include "HistoManager.hh" << 47 #include "G4UIdirectory.hh"
47 <<
48 #include "G4HadronicParameters.hh" <<
49 #include "G4UIcmdWith3Vector.hh" <<
50 #include "G4UIcmdWithABool.hh" 48 #include "G4UIcmdWithABool.hh"
51 #include "G4UIcmdWithADoubleAndUnit.hh" <<
52 #include "G4UIcmdWithAString.hh" 49 #include "G4UIcmdWithAString.hh"
53 #include "G4UIcmdWithAnInteger.hh" 50 #include "G4UIcmdWithAnInteger.hh"
>> 51 #include "G4UIcmdWith3Vector.hh"
>> 52 #include "G4UIcmdWithADoubleAndUnit.hh"
54 #include "G4UIcmdWithoutParameter.hh" 53 #include "G4UIcmdWithoutParameter.hh"
55 #include "G4UIdirectory.hh" << 54 #include "HistoManager.hh"
56 55
57 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58 57
59 DetectorMessenger::DetectorMessenger(DetectorC << 58 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det)
60 : G4UImessenger(), << 59 :fDetector(Det)
61 fDetector(Det), <<
62 fTestDir(0), <<
63 fMatCmd(0), <<
64 fMat1Cmd(0), <<
65 fRCmd(0), <<
66 fLCmd(0), <<
67 fEdepCmd(0), <<
68 fBinCmd(0), <<
69 fNOfAbsCmd(0), <<
70 fVerbCmd(0), <<
71 fBeamCmd(0) <<
72 { 60 {
73 fTestDir = new G4UIdirectory("/testhadr/"); 61 fTestDir = new G4UIdirectory("/testhadr/");
74 fTestDir->SetGuidance(" Hadronic Extended Ex 62 fTestDir->SetGuidance(" Hadronic Extended Example.");
75 63
76 fMatCmd = new G4UIcmdWithAString("/testhadr/ << 64 fMatCmd = new G4UIcmdWithAString("/testhadr/TargetMat",this);
77 fMatCmd->SetGuidance("Select Material for th 65 fMatCmd->SetGuidance("Select Material for the target");
78 fMatCmd->SetParameterName("tMaterial", false << 66 fMatCmd->SetParameterName("tMaterial",false);
79 fMatCmd->AvailableForStates(G4State_PreInit, << 67 fMatCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
80 68
81 fMat1Cmd = new G4UIcmdWithAString("/testhadr << 69 fMat1Cmd = new G4UIcmdWithAString("/testhadr/WorldMat",this);
82 fMat1Cmd->SetGuidance("Select Material for w 70 fMat1Cmd->SetGuidance("Select Material for world");
83 fMat1Cmd->SetParameterName("wMaterial", fals << 71 fMat1Cmd->SetParameterName("wMaterial",false);
84 fMat1Cmd->AvailableForStates(G4State_PreInit << 72 fMat1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
85 73
86 fRCmd = new G4UIcmdWithADoubleAndUnit("/test << 74 fRCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetRadius",this);
87 fRCmd->SetGuidance("Set radius of the target 75 fRCmd->SetGuidance("Set radius of the target");
88 fRCmd->SetParameterName("radius", false); << 76 fRCmd->SetParameterName("radius",false);
89 fRCmd->SetUnitCategory("Length"); 77 fRCmd->SetUnitCategory("Length");
90 fRCmd->SetRange("radius>0"); 78 fRCmd->SetRange("radius>0");
91 fRCmd->AvailableForStates(G4State_PreInit); << 79 fRCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
92 80
93 fLCmd = new G4UIcmdWithADoubleAndUnit("/test << 81 fLCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetLength",this);
94 fLCmd->SetGuidance("Set length of the target 82 fLCmd->SetGuidance("Set length of the target");
95 fLCmd->SetParameterName("length", false); << 83 fLCmd->SetParameterName("length",false);
96 fLCmd->SetUnitCategory("Length"); 84 fLCmd->SetUnitCategory("Length");
97 fLCmd->SetRange("length>0"); 85 fLCmd->SetRange("length>0");
98 fLCmd->AvailableForStates(G4State_PreInit); << 86 fLCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
99 87
100 fBinCmd = new G4UIcmdWithAnInteger("/testhad << 88 fBinCmd = new G4UIcmdWithAnInteger("/testhadr/NumberOfBinsE",this);
101 fBinCmd->SetGuidance("Set number of bins for 89 fBinCmd->SetGuidance("Set number of bins for Energy");
102 fBinCmd->SetParameterName("NEbins", false); << 90 fBinCmd->SetParameterName("NEbins",false);
103 fBinCmd->AvailableForStates(G4State_PreInit) << 91 fBinCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
104 92
105 fNOfAbsCmd = new G4UIcmdWithAnInteger("/test << 93 fNOfAbsCmd = new G4UIcmdWithAnInteger("/testhadr/NumberDivZ",this);
106 fNOfAbsCmd->SetGuidance("Set number of slice 94 fNOfAbsCmd->SetGuidance("Set number of slices");
107 fNOfAbsCmd->SetParameterName("NZ", false); << 95 fNOfAbsCmd->SetParameterName("NZ",false);
108 fNOfAbsCmd->AvailableForStates(G4State_PreIn << 96 fNOfAbsCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
109 97
110 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/t << 98 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/MaxEdep",this);
111 fEdepCmd->SetGuidance("Set max energy in his 99 fEdepCmd->SetGuidance("Set max energy in histogram");
112 fEdepCmd->SetParameterName("edep", false); << 100 fEdepCmd->SetParameterName("edep",false);
113 fEdepCmd->SetUnitCategory("Energy"); 101 fEdepCmd->SetUnitCategory("Energy");
114 fEdepCmd->SetRange("edep>0"); 102 fEdepCmd->SetRange("edep>0");
115 fEdepCmd->AvailableForStates(G4State_PreInit << 103 fEdepCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
116 104
117 fBeamCmd = new G4UIcmdWithABool("/testhadr/D << 105 fUpdateCmd = new G4UIcmdWithoutParameter("/testhadr/Update",this);
118 fBeamCmd->SetGuidance("Use default beam posi << 106 fUpdateCmd->SetGuidance("Update geometry.");
119 fBeamCmd->AvailableForStates(G4State_PreInit << 107 fUpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
>> 108 fUpdateCmd->SetGuidance("if you changed geometrical value(s)");
>> 109 fUpdateCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
120 110
121 fVerbCmd = new G4UIcmdWithAnInteger("/testha << 111 fBeamCmd = new G4UIcmdWithABool("/testhadr/DefaultBeamPosition",this);
122 fVerbCmd->SetGuidance("Set verbose for "); << 112 fBeamCmd->SetGuidance("show inelastic and elastic cross sections");
123 fVerbCmd->SetParameterName("verb", false); <<
124 fVerbCmd->AvailableForStates(G4State_PreInit <<
125 113
126 fBertCmd = new G4UIcmdWithAnInteger("/testha << 114 fVerbCmd = new G4UIcmdWithAnInteger("/testhadr/Verbose",this);
127 fBertCmd->SetGuidance("Set printout of Berti << 115 fVerbCmd->SetGuidance("Set verbose for ");
128 fBertCmd->SetParameterName("bert", false); << 116 fVerbCmd->SetParameterName("verb",false);
129 fBertCmd->AvailableForStates(G4State_PreInit << 117 fVerbCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
130 <<
131 fBCCmd = new G4UIcmdWithABool("/testhadr/BCP <<
132 fBCCmd->SetGuidance("Enable b-, c- hadronic <<
133 fBCCmd->SetParameterName("BCparticles", fals <<
134 fBCCmd->AvailableForStates(G4State_PreInit); <<
135 <<
136 fHNCmd = new G4UIcmdWithABool("/testhadr/Hyp <<
137 fHNCmd->SetGuidance("Enable hyper-nuclei had <<
138 fHNCmd->SetParameterName("HyperNuclei", fals <<
139 fHNCmd->AvailableForStates(G4State_PreInit); <<
140 <<
141 fNGPCmd = new G4UIcmdWithABool("/testhadr/Ne <<
142 fNGPCmd->SetGuidance("Enable neutron general <<
143 fNGPCmd->SetParameterName("NeutronGeneralPro <<
144 fNGPCmd->AvailableForStates(G4State_PreInit) <<
145 } 118 }
146 119
147 //....oooOO0OOooo........oooOO0OOooo........oo 120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
148 121
149 DetectorMessenger::~DetectorMessenger() 122 DetectorMessenger::~DetectorMessenger()
150 { 123 {
151 delete fMatCmd; 124 delete fMatCmd;
152 delete fMat1Cmd; 125 delete fMat1Cmd;
153 delete fRCmd; 126 delete fRCmd;
154 delete fLCmd; 127 delete fLCmd;
155 delete fBinCmd; 128 delete fBinCmd;
156 delete fNOfAbsCmd; 129 delete fNOfAbsCmd;
>> 130 delete fUpdateCmd;
157 delete fTestDir; 131 delete fTestDir;
158 delete fBeamCmd; 132 delete fBeamCmd;
159 delete fVerbCmd; 133 delete fVerbCmd;
160 delete fBertCmd; <<
161 delete fEdepCmd; 134 delete fEdepCmd;
162 delete fBCCmd; <<
163 delete fHNCmd; <<
164 delete fNGPCmd; <<
165 } 135 }
166 136
167 //....oooOO0OOooo........oooOO0OOooo........oo 137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
168 138
169 void DetectorMessenger::SetNewValue(G4UIcomman 139 void DetectorMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
170 { 140 {
171 HistoManager* h = HistoManager::GetPointer() 141 HistoManager* h = HistoManager::GetPointer();
172 if (command == fMatCmd) { << 142 if( command == fMatCmd ) {
173 fDetector->SetTargetMaterial(newValue); << 143 fDetector->SetTargetMaterial(newValue);
174 } << 144 } else if( command == fMat1Cmd ) {
175 else if (command == fMat1Cmd) { << 145 fDetector->SetWorldMaterial(newValue);
176 fDetector->SetWorldMaterial(newValue); << 146 } else if( command == fRCmd ) {
177 } << 147 fDetector->SetTargetRadius(fRCmd->GetNewDoubleValue(newValue));
178 else if (command == fRCmd) { << 148 } else if( command == fLCmd ) {
179 h->SetTargetRadius(fRCmd->GetNewDoubleValu <<
180 } <<
181 else if (command == fLCmd) { <<
182 h->SetTargetLength(fLCmd->GetNewDoubleValu 149 h->SetTargetLength(fLCmd->GetNewDoubleValue(newValue));
183 } << 150 } else if( command == fNOfAbsCmd ) {
184 else if (command == fNOfAbsCmd) { <<
185 h->SetNumberOfSlices(fNOfAbsCmd->GetNewInt 151 h->SetNumberOfSlices(fNOfAbsCmd->GetNewIntValue(newValue));
186 } << 152 } else if( command == fBinCmd ) {
187 else if (command == fBinCmd) { <<
188 h->SetNumberOfBinsE(fBinCmd->GetNewIntValu 153 h->SetNumberOfBinsE(fBinCmd->GetNewIntValue(newValue));
189 } << 154 } else if( command == fVerbCmd ) {
190 else if (command == fVerbCmd) { <<
191 h->SetVerbose(fVerbCmd->GetNewIntValue(new 155 h->SetVerbose(fVerbCmd->GetNewIntValue(newValue));
192 } << 156 } else if (command == fBeamCmd) {
193 else if (command == fBertCmd) { <<
194 h->SetPrintBertiniXS(fBertCmd->GetNewIntVa <<
195 } <<
196 else if (command == fBeamCmd) { <<
197 h->SetDefaultBeamPositionFlag(fBeamCmd->Ge 157 h->SetDefaultBeamPositionFlag(fBeamCmd->GetNewBoolValue(newValue));
198 } << 158 } else if (command == fEdepCmd) {
199 else if (command == fEdepCmd) { <<
200 h->SetMaxEnergyDeposit(fEdepCmd->GetNewDou 159 h->SetMaxEnergyDeposit(fEdepCmd->GetNewDoubleValue(newValue));
201 } << 160 } else if( command == fUpdateCmd ) {
202 else if (command == fBCCmd) { << 161 fDetector->UpdateGeometry();
203 G4HadronicParameters::Instance()->SetEnabl <<
204 } <<
205 else if (command == fHNCmd) { <<
206 G4HadronicParameters::Instance()->SetEnabl <<
207 } <<
208 else if (command == fNGPCmd) { <<
209 G4HadronicParameters::Instance()->SetEnabl <<
210 fNGPCmd->GetNewBoolValue(newValue)); <<
211 } 162 }
212 } 163 }
213 164
214 //....oooOO0OOooo........oooOO0OOooo........oo 165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 166
215 167