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