Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 /// \file DetectorMessenger.cc 27 /// \brief Implementation of the DetectorMesse 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "DetectorMessenger.hh" 34 35 #include "DetectorConstruction.hh" 36 37 #include "G4UIcmdWithADoubleAndUnit.hh" 38 #include "G4UIcmdWithAnInteger.hh" 39 #include "G4UIcmdWithoutParameter.hh" 40 #include "G4UIcommand.hh" 41 #include "G4UIdirectory.hh" 42 #include "G4UIparameter.hh" 43 44 #include <sstream> 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 48 DetectorMessenger::DetectorMessenger(DetectorC 49 { 50 fTestemDir = new G4UIdirectory("/testhadr/") 51 fTestemDir->SetGuidance(" detector control." 52 53 fDetDir = new G4UIdirectory("/testhadr/det/" 54 fDetDir->SetGuidance("detector construction 55 56 fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes 57 fNbAbsorCmd->SetGuidance("Set number of Abso 58 fNbAbsorCmd->SetParameterName("NbAbsor", fal 59 fNbAbsorCmd->SetRange("NbAbsor>0"); 60 fNbAbsorCmd->AvailableForStates(G4State_PreI 61 fNbAbsorCmd->SetToBeBroadcasted(false); 62 63 fAbsorCmd = new G4UIcommand("/testhadr/det/s 64 fAbsorCmd->SetGuidance("Set the absor nb, th 65 fAbsorCmd->SetGuidance(" absor number : fro 66 fAbsorCmd->SetGuidance(" material name"); 67 fAbsorCmd->SetGuidance(" thickness (with un 68 // 69 G4UIparameter* AbsNbPrm = new G4UIparameter( 70 AbsNbPrm->SetGuidance("absor number : from 1 71 AbsNbPrm->SetParameterRange("AbsorNb>0"); 72 fAbsorCmd->SetParameter(AbsNbPrm); 73 // 74 G4UIparameter* MatPrm = new G4UIparameter("m 75 MatPrm->SetGuidance("material name"); 76 fAbsorCmd->SetParameter(MatPrm); 77 // 78 G4UIparameter* ThickPrm = new G4UIparameter( 79 ThickPrm->SetGuidance("thickness of absorber 80 ThickPrm->SetParameterRange("thickness>0."); 81 fAbsorCmd->SetParameter(ThickPrm); 82 // 83 G4UIparameter* unitPrm = new G4UIparameter(" 84 unitPrm->SetGuidance("unit of thickness"); 85 G4String unitList = G4UIcommand::UnitsList(G 86 unitPrm->SetParameterCandidates(unitList); 87 fAbsorCmd->SetParameter(unitPrm); 88 // 89 fAbsorCmd->AvailableForStates(G4State_PreIni 90 fAbsorCmd->SetToBeBroadcasted(false); 91 92 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" 93 fSizeYZCmd->SetGuidance("Set sizeYZ of the a 94 fSizeYZCmd->SetParameterName("SizeYZ", false 95 fSizeYZCmd->SetRange("SizeYZ>0."); 96 fSizeYZCmd->SetUnitCategory("Length"); 97 fSizeYZCmd->AvailableForStates(G4State_PreIn 98 fSizeYZCmd->SetToBeBroadcasted(false); 99 100 fIsotopeCmd = new G4UIcommand("/testhadr/det 101 fIsotopeCmd->SetGuidance("Build and select a 102 fIsotopeCmd->SetGuidance(" symbol of isotop 103 // 104 G4UIparameter* symbPrm = new G4UIparameter(" 105 symbPrm->SetGuidance("isotope symbol"); 106 fIsotopeCmd->SetParameter(symbPrm); 107 // 108 G4UIparameter* ZPrm = new G4UIparameter("Z", 109 ZPrm->SetGuidance("Z"); 110 ZPrm->SetParameterRange("Z>0"); 111 fIsotopeCmd->SetParameter(ZPrm); 112 // 113 G4UIparameter* APrm = new G4UIparameter("A", 114 APrm->SetGuidance("A"); 115 APrm->SetParameterRange("A>0"); 116 fIsotopeCmd->SetParameter(APrm); 117 // 118 G4UIparameter* densityPrm = new G4UIparamete 119 densityPrm->SetGuidance("density of material 120 densityPrm->SetParameterRange("density>0."); 121 fIsotopeCmd->SetParameter(densityPrm); 122 // 123 G4UIparameter* unitPrm1 = new G4UIparameter( 124 unitPrm1->SetGuidance("unit of density"); 125 G4String unitList1 = G4UIcommand::UnitsList( 126 unitPrm1->SetParameterCandidates(unitList1); 127 fIsotopeCmd->SetParameter(unitPrm1); 128 // 129 fIsotopeCmd->AvailableForStates(G4State_PreI 130 } 131 132 //....oooOO0OOooo........oooOO0OOooo........oo 133 134 DetectorMessenger::~DetectorMessenger() 135 { 136 delete fNbAbsorCmd; 137 delete fAbsorCmd; 138 delete fSizeYZCmd; 139 delete fIsotopeCmd; 140 delete fDetDir; 141 delete fTestemDir; 142 } 143 144 //....oooOO0OOooo........oooOO0OOooo........oo 145 146 void DetectorMessenger::SetNewValue(G4UIcomman 147 { 148 if (command == fNbAbsorCmd) { 149 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe 150 } 151 152 if (command == fAbsorCmd) { 153 G4int num; 154 G4double tick; 155 G4String unt, mat; 156 std::istringstream is(newValue); 157 is >> num >> mat >> tick >> unt; 158 G4String material = mat; 159 tick *= G4UIcommand::ValueOf(unt); 160 fDetector->SetAbsorMaterial(num, material) 161 fDetector->SetAbsorThickness(num, tick); 162 } 163 164 if (command == fSizeYZCmd) { 165 fDetector->SetAbsorSizeYZ(fSizeYZCmd->GetN 166 } 167 168 if (command == fIsotopeCmd) { 169 G4int Z; 170 G4int A; 171 G4double dens; 172 G4String name, unt; 173 std::istringstream is(newValue); 174 is >> name >> Z >> A >> dens >> unt; 175 dens *= G4UIcommand::ValueOf(unt); 176 fDetector->MaterialWithSingleIsotope(name, 177 } 178 } 179 180 //....oooOO0OOooo........oooOO0OOooo........oo 181