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