| 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 DetectorMessenger.cc 26 /// \file DetectorMessenger.cc
27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class
28 // 28 //
29 // 29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 32
33 #include "DetectorMessenger.hh" 33 #include "DetectorMessenger.hh"
34 34
35 #include "DetectorConstruction.hh" 35 #include "DetectorConstruction.hh"
36 <<
37 #include "G4UIcmdWithADoubleAndUnit.hh" <<
38 #include "G4UIcmdWithAString.hh" <<
39 #include "G4UIcmdWithoutParameter.hh" <<
40 #include "G4UIcommand.hh" <<
41 #include "G4UIdirectory.hh" 36 #include "G4UIdirectory.hh"
>> 37 #include "G4UIcommand.hh"
42 #include "G4UIparameter.hh" 38 #include "G4UIparameter.hh"
>> 39 #include "G4UIcmdWithAString.hh"
>> 40 #include "G4UIcmdWithADoubleAndUnit.hh"
>> 41 #include "G4UIcmdWithoutParameter.hh"
43 42
44 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45 44
46 DetectorMessenger::DetectorMessenger(DetectorC << 45 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det)
47 { << 46 :G4UImessenger(),
>> 47 fDetector(Det), fTesthadrDir(0), fDetDir(0), fMaterCmd1(0), fMaterCmd2(0),
>> 48 fSizeCmd1(0), fSizeCmd2(0), fSizeCmd3(0),
>> 49 fIsotopeCmd(0)
>> 50 {
48 fTesthadrDir = new G4UIdirectory("/testhadr/ 51 fTesthadrDir = new G4UIdirectory("/testhadr/");
49 fTesthadrDir->SetGuidance("commands specific 52 fTesthadrDir->SetGuidance("commands specific to this example");
50 53
51 G4bool broadcast = false; 54 G4bool broadcast = false;
52 fDetDir = new G4UIdirectory("/testhadr/det/" << 55 fDetDir = new G4UIdirectory("/testhadr/det/",broadcast);
53 fDetDir->SetGuidance("detector construction 56 fDetDir->SetGuidance("detector construction commands");
54 57
55 fMaterCmd1 = new G4UIcmdWithAString("/testha << 58 fMaterCmd1 = new G4UIcmdWithAString("/testhadr/det/setAbsorMat",this);
56 fMaterCmd1->SetGuidance("Select absorber mat 59 fMaterCmd1->SetGuidance("Select absorber material");
57 fMaterCmd1->SetParameterName("choice", false << 60 fMaterCmd1->SetParameterName("choice",false);
58 fMaterCmd1->AvailableForStates(G4State_PreIn << 61 fMaterCmd1->AvailableForStates(G4State_PreInit,G4State_Idle);
59 62
60 fMaterCmd2 = new G4UIcmdWithAString("/testha << 63 fMaterCmd2 = new G4UIcmdWithAString("/testhadr/det/setContMat",this);
61 fMaterCmd2->SetGuidance("Select container ma 64 fMaterCmd2->SetGuidance("Select container material");
62 fMaterCmd2->SetParameterName("choice", false << 65 fMaterCmd2->SetParameterName("choice",false);
63 fMaterCmd2->AvailableForStates(G4State_PreIn << 66 fMaterCmd2->AvailableForStates(G4State_PreInit,G4State_Idle);
64 67
65 fSizeCmd1 = new G4UIcmdWithADoubleAndUnit("/ << 68 fSizeCmd1 = new G4UIcmdWithADoubleAndUnit("/testhadr/det/setAbsorRadius",this);
66 fSizeCmd1->SetGuidance("Set absorber radius" 69 fSizeCmd1->SetGuidance("Set absorber radius");
67 fSizeCmd1->SetParameterName("Radius", false) << 70 fSizeCmd1->SetParameterName("Radius",false);
68 fSizeCmd1->SetRange("Radius>0."); 71 fSizeCmd1->SetRange("Radius>0.");
69 fSizeCmd1->SetUnitCategory("Length"); 72 fSizeCmd1->SetUnitCategory("Length");
70 fSizeCmd1->AvailableForStates(G4State_PreIni << 73 fSizeCmd1->AvailableForStates(G4State_PreInit,G4State_Idle);
71 74
72 fSizeCmd2 = new G4UIcmdWithADoubleAndUnit("/ << 75 fSizeCmd2 = new G4UIcmdWithADoubleAndUnit("/testhadr/det/setAbsorLength",this);
73 fSizeCmd2->SetGuidance("Set absorber length" 76 fSizeCmd2->SetGuidance("Set absorber length");
74 fSizeCmd2->SetParameterName("Length", false) << 77 fSizeCmd2->SetParameterName("Length",false);
75 fSizeCmd2->SetRange("Length>0."); 78 fSizeCmd2->SetRange("Length>0.");
76 fSizeCmd2->SetUnitCategory("Length"); 79 fSizeCmd2->SetUnitCategory("Length");
77 fSizeCmd2->AvailableForStates(G4State_PreIni << 80 fSizeCmd2->AvailableForStates(G4State_PreInit,G4State_Idle);
78 81
79 fSizeCmd3 = new G4UIcmdWithADoubleAndUnit("/ << 82 fSizeCmd3 = new G4UIcmdWithADoubleAndUnit("/testhadr/det/setContThick",this);
80 fSizeCmd3->SetGuidance("Set container thickn 83 fSizeCmd3->SetGuidance("Set container thickness");
81 fSizeCmd3->SetParameterName("Thick", false); << 84 fSizeCmd3->SetParameterName("Thick",false);
82 fSizeCmd3->SetRange("Thick>0."); 85 fSizeCmd3->SetRange("Thick>0.");
83 fSizeCmd3->SetUnitCategory("Length"); 86 fSizeCmd3->SetUnitCategory("Length");
84 fSizeCmd3->AvailableForStates(G4State_PreIni << 87 fSizeCmd3->AvailableForStates(G4State_PreInit,G4State_Idle);
85 << 88
86 fIsotopeCmd = new G4UIcommand("/testhadr/det << 89 fIsotopeCmd = new G4UIcommand("/testhadr/det/setIsotopeMat",this);
87 fIsotopeCmd->SetGuidance("Build and select a 90 fIsotopeCmd->SetGuidance("Build and select a material with single isotope");
88 fIsotopeCmd->SetGuidance(" symbol of isotop 91 fIsotopeCmd->SetGuidance(" symbol of isotope, Z, A, density of material");
89 // 92 //
90 G4UIparameter* symbPrm = new G4UIparameter(" << 93 G4UIparameter* symbPrm = new G4UIparameter("isotope",'s',false);
91 symbPrm->SetGuidance("isotope symbol"); 94 symbPrm->SetGuidance("isotope symbol");
92 fIsotopeCmd->SetParameter(symbPrm); 95 fIsotopeCmd->SetParameter(symbPrm);
93 // << 96 //
94 G4UIparameter* ZPrm = new G4UIparameter("Z", << 97 G4UIparameter* ZPrm = new G4UIparameter("Z",'i',false);
95 ZPrm->SetGuidance("Z"); 98 ZPrm->SetGuidance("Z");
96 ZPrm->SetParameterRange("Z>0"); 99 ZPrm->SetParameterRange("Z>0");
97 fIsotopeCmd->SetParameter(ZPrm); 100 fIsotopeCmd->SetParameter(ZPrm);
98 // << 101 //
99 G4UIparameter* APrm = new G4UIparameter("A", << 102 G4UIparameter* APrm = new G4UIparameter("A",'i',false);
100 APrm->SetGuidance("A"); 103 APrm->SetGuidance("A");
101 APrm->SetParameterRange("A>0"); 104 APrm->SetParameterRange("A>0");
102 fIsotopeCmd->SetParameter(APrm); << 105 fIsotopeCmd->SetParameter(APrm);
103 // << 106 //
104 G4UIparameter* densityPrm = new G4UIparamete << 107 G4UIparameter* densityPrm = new G4UIparameter("density",'d',false);
105 densityPrm->SetGuidance("density of material 108 densityPrm->SetGuidance("density of material");
106 densityPrm->SetParameterRange("density>0."); 109 densityPrm->SetParameterRange("density>0.");
107 fIsotopeCmd->SetParameter(densityPrm); 110 fIsotopeCmd->SetParameter(densityPrm);
108 // 111 //
109 G4UIparameter* unitPrm = new G4UIparameter(" << 112 G4UIparameter* unitPrm = new G4UIparameter("unit",'s',false);
110 unitPrm->SetGuidance("unit of density"); 113 unitPrm->SetGuidance("unit of density");
111 G4String unitList = G4UIcommand::UnitsList(G 114 G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("g/cm3"));
112 unitPrm->SetParameterCandidates(unitList); 115 unitPrm->SetParameterCandidates(unitList);
113 fIsotopeCmd->SetParameter(unitPrm); 116 fIsotopeCmd->SetParameter(unitPrm);
114 // 117 //
115 fIsotopeCmd->AvailableForStates(G4State_PreI << 118 fIsotopeCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
116 } 119 }
117 120
118 //....oooOO0OOooo........oooOO0OOooo........oo 121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
119 122
120 DetectorMessenger::~DetectorMessenger() 123 DetectorMessenger::~DetectorMessenger()
121 { 124 {
122 delete fMaterCmd1; 125 delete fMaterCmd1;
123 delete fMaterCmd2; 126 delete fMaterCmd2;
124 delete fSizeCmd1; 127 delete fSizeCmd1;
125 delete fSizeCmd2; 128 delete fSizeCmd2;
126 delete fSizeCmd3; 129 delete fSizeCmd3;
127 delete fIsotopeCmd; 130 delete fIsotopeCmd;
128 delete fDetDir; 131 delete fDetDir;
129 delete fTesthadrDir; 132 delete fTesthadrDir;
130 } 133 }
131 134
132 //....oooOO0OOooo........oooOO0OOooo........oo 135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
133 136
134 void DetectorMessenger::SetNewValue(G4UIcomman << 137 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
135 { << 138 {
136 if (command == fMaterCmd1) { << 139 if( command == fMaterCmd1 )
137 fDetector->SetAbsorMaterial(newValue); << 140 { fDetector->SetAbsorMaterial(newValue);}
138 } << 141 if( command == fMaterCmd2 )
139 if (command == fMaterCmd2) { << 142 { fDetector->SetContainMaterial(newValue);}
140 fDetector->SetContainMaterial(newValue); << 143
141 } << 144 if( command == fSizeCmd1 )
142 << 145 { fDetector->SetAbsorRadius(fSizeCmd1->GetNewDoubleValue(newValue));}
143 if (command == fSizeCmd1) { << 146
144 fDetector->SetAbsorRadius(fSizeCmd1->GetNe << 147 if( command == fSizeCmd2 )
145 } << 148 { fDetector->SetAbsorLength(fSizeCmd2->GetNewDoubleValue(newValue));}
146 << 149
147 if (command == fSizeCmd2) { << 150 if( command == fSizeCmd3 )
148 fDetector->SetAbsorLength(fSizeCmd2->GetNe << 151 { fDetector->SetContainThickness(fSizeCmd3->GetNewDoubleValue(newValue));}
149 } << 152
150 << 153 if (command == fIsotopeCmd)
151 if (command == fSizeCmd3) { << 154 {
152 fDetector->SetContainThickness(fSizeCmd3-> << 155 G4int Z; G4int A; G4double dens;
153 } << 156 G4String name, unt;
154 << 157 std::istringstream is(newValue);
155 if (command == fIsotopeCmd) { << 158 is >> name >> Z >> A >> dens >> unt;
156 G4int Z; << 159 dens *= G4UIcommand::ValueOf(unt);
157 G4int A; << 160 fDetector->MaterialWithSingleIsotope (name,name,dens,Z,A);
158 G4double dens; << 161 fDetector->SetAbsorMaterial(name);
159 G4String name, unt; << 162 }
160 std::istringstream is(newValue); <<
161 is >> name >> Z >> A >> dens >> unt; <<
162 dens *= G4UIcommand::ValueOf(unt); <<
163 fDetector->MaterialWithSingleIsotope(name, <<
164 fDetector->SetAbsorMaterial(name); <<
165 } <<
166 } 163 }
167 164
168 //....oooOO0OOooo........oooOO0OOooo........oo 165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169 166