| 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 electromagnetic/TestEm9/src/Detector 26 /// \file electromagnetic/TestEm9/src/DetectorMessenger.cc
27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class
28 // 28 //
29 // 29 //
30 // 30 //
31 ////////////////////////////////////////////// 31 /////////////////////////////////////////////////////////////////////////
32 // 32 //
33 // TestEm9: Crystal calorimeter 33 // TestEm9: Crystal calorimeter
34 // 34 //
35 // Created: 31.01.03 V.Ivanchenko 35 // Created: 31.01.03 V.Ivanchenko
36 // 36 //
37 // Modified: 37 // Modified:
38 // 38 //
39 ////////////////////////////////////////////// 39 ////////////////////////////////////////////////////////////////////////
40 // 40 //
41 41
42 //....oooOO0OOooo........oooOO0OOooo........oo 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 44
45 #include "DetectorMessenger.hh" 45 #include "DetectorMessenger.hh"
46 46
47 #include "DetectorConstruction.hh" 47 #include "DetectorConstruction.hh"
48 #include "HistoManager.hh" << 48 #include "G4UIdirectory.hh"
49 << 49 #include "G4UIcmdWithAString.hh"
50 #include "G4UIcmdWith3Vector.hh" 50 #include "G4UIcmdWith3Vector.hh"
51 #include "G4UIcmdWithADoubleAndUnit.hh" 51 #include "G4UIcmdWithADoubleAndUnit.hh"
52 #include "G4UIcmdWithAString.hh" <<
53 #include "G4UIcmdWithoutParameter.hh" 52 #include "G4UIcmdWithoutParameter.hh"
54 #include "G4UIdirectory.hh" << 53 #include "HistoManager.hh"
55 54
56 //....oooOO0OOooo........oooOO0OOooo........oo 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
57 56
58 DetectorMessenger::DetectorMessenger(DetectorC << 57 DetectorMessenger::DetectorMessenger(DetectorConstruction * det)
59 : G4UImessenger(), << 58 :G4UImessenger(),fDetector(det),
60 fDetector(det), << 59 fAtestemDir(0),
61 fAtestemDir(0), << 60 fAMaterCmd(0),
62 fAMaterCmd(0), << 61 fALBinCmd(0),
63 fALBinCmd(0), << 62 fAl1Cmd(0),
64 fAl1Cmd(0), << 63 fAl2Cmd(0),
65 fAl2Cmd(0), << 64 fAl3Cmd(0),
66 fAl3Cmd(0), << 65 fAl4Cmd(0),
67 fAl4Cmd(0), << 66 fAl5Cmd(0),
68 fAl5Cmd(0), << 67 fAl6Cmd(0),
69 fAl6Cmd(0), << 68 fAUpdateCmd(0),
70 fAUpdateCmd(0), << 69 fAaccCmd1(0),
71 fAaccCmd1(0), << 70 fAaccCmd2(0),
72 fAaccCmd2(0), << 71 fAaccCmd3(0)
73 fAaccCmd3(0) <<
74 { 72 {
75 fAtestemDir = new G4UIdirectory("/testem/"); 73 fAtestemDir = new G4UIdirectory("/testem/");
76 fAtestemDir->SetGuidance(" detector control. 74 fAtestemDir->SetGuidance(" detector control.");
77 75
78 fAMaterCmd = new G4UIcmdWithAString("/testem << 76 fAMaterCmd = new G4UIcmdWithAString("/testem/det/CalMat",this);
79 fAMaterCmd->SetGuidance("Select Material for 77 fAMaterCmd->SetGuidance("Select Material for calorimeter");
80 fAMaterCmd->SetParameterName("calMaterial", << 78 fAMaterCmd->SetParameterName("calMaterial",false);
81 fAMaterCmd->AvailableForStates(G4State_PreIn << 79 fAMaterCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
82 80
83 fALBinCmd = new G4UIcmdWithAString("/testem/ << 81 fALBinCmd = new G4UIcmdWithAString("/testem/det/AbsMat",this);
84 fALBinCmd->SetGuidance("Select Material for 82 fALBinCmd->SetGuidance("Select Material for absorber");
85 fALBinCmd->SetParameterName("absMarerial", f << 83 fALBinCmd->SetParameterName("absMarerial",false);
86 fALBinCmd->AvailableForStates(G4State_PreIni << 84 fALBinCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
87 85
88 fAl1Cmd = new G4UIcmdWithADoubleAndUnit("/te << 86 fAl1Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/EcalLength",this);
89 fAl1Cmd->SetGuidance("Set length of Ecal"); 87 fAl1Cmd->SetGuidance("Set length of Ecal");
90 fAl1Cmd->SetParameterName("lEcal", false); << 88 fAl1Cmd->SetParameterName("lEcal",false);
91 fAl1Cmd->SetUnitCategory("Length"); 89 fAl1Cmd->SetUnitCategory("Length");
92 fAl1Cmd->SetRange("lEcal>0"); 90 fAl1Cmd->SetRange("lEcal>0");
93 fAl1Cmd->AvailableForStates(G4State_PreInit, << 91 fAl1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
94 92
95 fAl2Cmd = new G4UIcmdWithADoubleAndUnit("/te << 93 fAl2Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/EcalWidth",this);
96 fAl2Cmd->SetGuidance("Set width of Ecal crys 94 fAl2Cmd->SetGuidance("Set width of Ecal crystal");
97 fAl2Cmd->SetParameterName("wEcal", false); << 95 fAl2Cmd->SetParameterName("wEcal",false);
98 fAl2Cmd->SetUnitCategory("Length"); 96 fAl2Cmd->SetUnitCategory("Length");
99 fAl2Cmd->AvailableForStates(G4State_PreInit, << 97 fAl2Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
100 98
101 fAl3Cmd = new G4UIcmdWithADoubleAndUnit("/te << 99 fAl3Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/AbsLength",this);
102 fAl3Cmd->SetGuidance("Set length of the abso 100 fAl3Cmd->SetGuidance("Set length of the absorber");
103 fAl3Cmd->SetParameterName("lAbs", false); << 101 fAl3Cmd->SetParameterName("lAbs",false);
104 fAl3Cmd->SetUnitCategory("Length"); 102 fAl3Cmd->SetUnitCategory("Length");
105 fAl3Cmd->AvailableForStates(G4State_PreInit, << 103 fAl3Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
106 104
107 fAl4Cmd = new G4UIcmdWithADoubleAndUnit("/te << 105 fAl4Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/VertexLength",this);
108 fAl4Cmd->SetGuidance("Set length of the vert 106 fAl4Cmd->SetGuidance("Set length of the vertex region");
109 fAl4Cmd->SetParameterName("lVert", false); << 107 fAl4Cmd->SetParameterName("lVert",false);
110 fAl4Cmd->SetUnitCategory("Length"); 108 fAl4Cmd->SetUnitCategory("Length");
111 fAl4Cmd->AvailableForStates(G4State_PreInit, << 109 fAl4Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
112 110
113 fAl5Cmd = new G4UIcmdWithADoubleAndUnit("/te << 111 fAl5Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/PadLength",this);
114 fAl5Cmd->SetGuidance("Set length of vertex d 112 fAl5Cmd->SetGuidance("Set length of vertex detector");
115 fAl5Cmd->SetParameterName("lPad", false); << 113 fAl5Cmd->SetParameterName("lPad",false);
116 fAl5Cmd->SetUnitCategory("Length"); 114 fAl5Cmd->SetUnitCategory("Length");
117 fAl5Cmd->AvailableForStates(G4State_PreInit, << 115 fAl5Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
118 116
119 fAl6Cmd = new G4UIcmdWithADoubleAndUnit("/te << 117 fAl6Cmd = new G4UIcmdWithADoubleAndUnit("/testem/det/PadWidth",this);
120 fAl6Cmd->SetGuidance("Set width of a vertex 118 fAl6Cmd->SetGuidance("Set width of a vertex pad");
121 fAl6Cmd->SetParameterName("wPad", false); << 119 fAl6Cmd->SetParameterName("wPad",false);
122 fAl6Cmd->SetUnitCategory("Length"); 120 fAl6Cmd->SetUnitCategory("Length");
123 fAl6Cmd->AvailableForStates(G4State_PreInit, << 121 fAl6Cmd->AvailableForStates(G4State_PreInit,G4State_Idle);
124 122
125 fAUpdateCmd = new G4UIcmdWithoutParameter("/ << 123 fAUpdateCmd = new G4UIcmdWithoutParameter("/testem/det/update",this);
126 fAUpdateCmd->SetGuidance("Update geometry.") 124 fAUpdateCmd->SetGuidance("Update geometry.");
127 fAUpdateCmd->SetGuidance("This command MUST 125 fAUpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
128 fAUpdateCmd->SetGuidance("if you changed geo 126 fAUpdateCmd->SetGuidance("if you changed geometrical value(s)");
129 fAUpdateCmd->AvailableForStates(G4State_PreI << 127 fAUpdateCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
130 128
131 fAaccCmd1 = new G4UIcmdWith3Vector("/testem/ << 129 fAaccCmd1 = new G4UIcmdWith3Vector("/testem/det/acceptance1",this);
132 fAaccCmd1->SetGuidance("set Edep and RMS"); 130 fAaccCmd1->SetGuidance("set Edep and RMS");
133 fAaccCmd1->SetGuidance("acceptance values fo 131 fAaccCmd1->SetGuidance("acceptance values for central cell");
134 fAaccCmd1->SetParameterName("edep", "rms", " << 132 fAaccCmd1->SetParameterName("edep","rms","limit",true);
135 fAaccCmd1->SetRange("edep>0 && edep<1 && rms 133 fAaccCmd1->SetRange("edep>0 && edep<1 && rms>0");
136 fAaccCmd1->AvailableForStates(G4State_PreIni << 134 fAaccCmd1->AvailableForStates(G4State_PreInit,G4State_Idle);
137 135
138 fAaccCmd2 = new G4UIcmdWith3Vector("/testem/ << 136 fAaccCmd2 = new G4UIcmdWith3Vector("/testem/det/acceptance9",this);
139 fAaccCmd2->SetGuidance("set Edep and RMS"); 137 fAaccCmd2->SetGuidance("set Edep and RMS");
140 fAaccCmd2->SetGuidance("acceptance values fo 138 fAaccCmd2->SetGuidance("acceptance values for 3x3 matrix");
141 fAaccCmd2->SetParameterName("edep", "rms", " << 139 fAaccCmd2->SetParameterName("edep","rms","limit",true);
142 fAaccCmd2->SetRange("edep>0 && edep<1 && rms 140 fAaccCmd2->SetRange("edep>0 && edep<1 && rms>0");
143 fAaccCmd2->AvailableForStates(G4State_PreIni << 141 fAaccCmd2->AvailableForStates(G4State_PreInit,G4State_Idle);
144 142
145 fAaccCmd3 = new G4UIcmdWith3Vector("/testem/ << 143 fAaccCmd3 = new G4UIcmdWith3Vector("/testem/det/acceptance25",this);
146 fAaccCmd3->SetGuidance("set Edep and RMS"); 144 fAaccCmd3->SetGuidance("set Edep and RMS");
147 fAaccCmd3->SetGuidance("acceptance values fo 145 fAaccCmd3->SetGuidance("acceptance values for 5x5 matrix");
148 fAaccCmd3->SetParameterName("edep", "rms", " << 146 fAaccCmd3->SetParameterName("edep","rms","limit",true);
149 fAaccCmd3->SetRange("edep>0 && edep<1 && rms 147 fAaccCmd3->SetRange("edep>0 && edep<1 && rms>0");
150 fAaccCmd3->AvailableForStates(G4State_PreIni << 148 fAaccCmd3->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 149
151 } 150 }
152 151
153 //....oooOO0OOooo........oooOO0OOooo........oo 152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
154 153
155 DetectorMessenger::~DetectorMessenger() 154 DetectorMessenger::~DetectorMessenger()
156 { 155 {
157 delete fAMaterCmd; 156 delete fAMaterCmd;
158 delete fALBinCmd; 157 delete fALBinCmd;
159 delete fAl1Cmd; 158 delete fAl1Cmd;
160 delete fAl2Cmd; 159 delete fAl2Cmd;
161 delete fAl3Cmd; 160 delete fAl3Cmd;
162 delete fAl4Cmd; 161 delete fAl4Cmd;
163 delete fAl5Cmd; 162 delete fAl5Cmd;
164 delete fAl6Cmd; 163 delete fAl6Cmd;
165 delete fAUpdateCmd; 164 delete fAUpdateCmd;
166 delete fAtestemDir; 165 delete fAtestemDir;
167 delete fAaccCmd1; 166 delete fAaccCmd1;
168 delete fAaccCmd2; 167 delete fAaccCmd2;
169 delete fAaccCmd3; 168 delete fAaccCmd3;
170 } 169 }
171 170
172 //....oooOO0OOooo........oooOO0OOooo........oo 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
173 172
174 void DetectorMessenger::SetNewValue(G4UIcomman << 173 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
175 { 174 {
176 if (command == fAMaterCmd) { << 175 if( command == fAMaterCmd )
177 fDetector->SetEcalMaterial(newValue); << 176 { fDetector->SetEcalMaterial(newValue);}
178 } << 177
179 << 178 if( command == fALBinCmd )
180 if (command == fALBinCmd) { << 179 { fDetector->SetAbsMaterial(newValue);}
181 fDetector->SetAbsMaterial(newValue); << 180
182 } << 181 if( command == fAl1Cmd )
183 << 182 { fDetector->SetEcalLength(fAl1Cmd->GetNewDoubleValue(newValue));}
184 if (command == fAl1Cmd) { << 183
185 fDetector->SetEcalLength(fAl1Cmd->GetNewDo << 184 if( command == fAl2Cmd )
186 } << 185 { fDetector->SetEcalWidth(fAl2Cmd->GetNewDoubleValue(newValue));}
187 << 186
188 if (command == fAl2Cmd) { << 187 if( command == fAl3Cmd )
189 fDetector->SetEcalWidth(fAl2Cmd->GetNewDou << 188 { fDetector->SetAbsLength(fAl3Cmd->GetNewDoubleValue(newValue));}
190 } << 189
191 << 190 if( command == fAl4Cmd )
192 if (command == fAl3Cmd) { << 191 { fDetector->SetVertexLength(fAl4Cmd->GetNewDoubleValue(newValue));}
193 fDetector->SetAbsLength(fAl3Cmd->GetNewDou << 192
194 } << 193 if( command == fAl5Cmd )
195 << 194 { fDetector->SetPadLength(fAl5Cmd->GetNewDoubleValue(newValue));}
196 if (command == fAl4Cmd) { << 195
197 fDetector->SetVertexLength(fAl4Cmd->GetNew << 196 if( command == fAl6Cmd )
198 } << 197 { fDetector->SetPadWidth(fAl6Cmd->GetNewDoubleValue(newValue));}
199 << 198
200 if (command == fAl5Cmd) { << 199 if( command == fAUpdateCmd )
201 fDetector->SetPadLength(fAl5Cmd->GetNewDou << 200 { fDetector->UpdateGeometry();}
202 } <<
203 <<
204 if (command == fAl6Cmd) { <<
205 fDetector->SetPadWidth(fAl6Cmd->GetNewDoub <<
206 } <<
207 <<
208 if (command == fAUpdateCmd) { <<
209 fDetector->UpdateGeometry(); <<
210 } <<
211 201
212 HistoManager* histo = HistoManager::GetPoint 202 HistoManager* histo = HistoManager::GetPointer();
213 if (command == fAaccCmd1) { << 203 if( command == fAaccCmd1 )
214 histo->SetEdepAndRMS(0, fAaccCmd1->GetNew3 << 204 { histo->SetEdepAndRMS(0,fAaccCmd1->GetNew3VectorValue(newValue));}
215 } << 205
216 << 206 if( command == fAaccCmd2 )
217 if (command == fAaccCmd2) { << 207 { histo->SetEdepAndRMS(1,fAaccCmd2->GetNew3VectorValue(newValue));}
218 histo->SetEdepAndRMS(1, fAaccCmd2->GetNew3 << 208
219 } << 209 if( command == fAaccCmd3 )
220 << 210 { histo->SetEdepAndRMS(2,fAaccCmd3->GetNew3VectorValue(newValue));}
221 if (command == fAaccCmd3) { << 211
222 histo->SetEdepAndRMS(2, fAaccCmd3->GetNew3 << 212
223 } <<
224 } 213 }
225 214
226 //....oooOO0OOooo........oooOO0OOooo........oo 215 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 216
227 217