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 "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........oooOO0OOooo........oooOO0OOooo...... 47 48 DetectorMessenger::DetectorMessenger(DetectorConstruction* Det) : fDetector(Det) 49 { 50 fTestemDir = new G4UIdirectory("/testhadr/"); 51 fTestemDir->SetGuidance("UI commands specific to this example"); 52 53 fDetDir = new G4UIdirectory("/testhadr/det/"); 54 fDetDir->SetGuidance("detector construction commands"); 55 56 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/det/setSizeYZ", this); 57 fSizeYZCmd->SetGuidance("Set tranverse size of the calorimeter"); 58 fSizeYZCmd->SetParameterName("Size", false); 59 fSizeYZCmd->SetRange("Size>0."); 60 fSizeYZCmd->SetUnitCategory("Length"); 61 fSizeYZCmd->AvailableForStates(G4State_PreInit); 62 fSizeYZCmd->SetToBeBroadcasted(false); 63 64 fNbLayersCmd = new G4UIcmdWithAnInteger("/testhadr/det/setNbOfLayers", this); 65 fNbLayersCmd->SetGuidance("Set number of layers."); 66 fNbLayersCmd->SetParameterName("NbLayers", false); 67 fNbLayersCmd->SetRange("NbLayers>0"); 68 fNbLayersCmd->AvailableForStates(G4State_PreInit); 69 fNbLayersCmd->SetToBeBroadcasted(false); 70 71 fNbAbsorCmd = new G4UIcmdWithAnInteger("/testhadr/det/setNbOfAbsor", this); 72 fNbAbsorCmd->SetGuidance("Set number of Absorbers."); 73 fNbAbsorCmd->SetParameterName("NbAbsor", false); 74 fNbAbsorCmd->SetRange("NbAbsor>0"); 75 fNbAbsorCmd->AvailableForStates(G4State_PreInit); 76 fNbAbsorCmd->SetToBeBroadcasted(false); 77 78 fAbsorCmd = new G4UIcommand("/testhadr/det/setAbsor", this); 79 fAbsorCmd->SetGuidance("Set the absor nb, the material, the thickness."); 80 fAbsorCmd->SetGuidance(" absor number : from 1 to NbOfAbsor"); 81 fAbsorCmd->SetGuidance(" material name"); 82 fAbsorCmd->SetGuidance(" thickness (with unit) : t>0."); 83 // 84 G4UIparameter* AbsNbPrm = new G4UIparameter("AbsorNb", 'i', false); 85 AbsNbPrm->SetGuidance("absor number : from 1 to NbOfAbsor"); 86 AbsNbPrm->SetParameterRange("AbsorNb>0"); 87 fAbsorCmd->SetParameter(AbsNbPrm); 88 // 89 G4UIparameter* MatPrm = new G4UIparameter("material", 's', false); 90 MatPrm->SetGuidance("material name"); 91 fAbsorCmd->SetParameter(MatPrm); 92 // 93 G4UIparameter* ThickPrm = new G4UIparameter("thickness", 'd', false); 94 ThickPrm->SetGuidance("thickness of absorber"); 95 ThickPrm->SetParameterRange("thickness>0."); 96 fAbsorCmd->SetParameter(ThickPrm); 97 // 98 G4UIparameter* unitPrm = new G4UIparameter("unit", 's', false); 99 unitPrm->SetGuidance("unit of thickness"); 100 G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("mm")); 101 unitPrm->SetParameterCandidates(unitList); 102 fAbsorCmd->SetParameter(unitPrm); 103 // 104 fAbsorCmd->AvailableForStates(G4State_PreInit); 105 fAbsorCmd->SetToBeBroadcasted(false); 106 107 fIsotopeCmd = new G4UIcommand("/testhadr/det/setIsotopeMat", this); 108 fIsotopeCmd->SetGuidance("Build and select a material with single isotope"); 109 fIsotopeCmd->SetGuidance(" symbol of isotope, Z, A, density of material"); 110 // 111 G4UIparameter* symbPrm = new G4UIparameter("isotope", 's', false); 112 symbPrm->SetGuidance("isotope symbol"); 113 fIsotopeCmd->SetParameter(symbPrm); 114 // 115 G4UIparameter* ZPrm = new G4UIparameter("Z", 'i', false); 116 ZPrm->SetGuidance("Z"); 117 ZPrm->SetParameterRange("Z>0"); 118 fIsotopeCmd->SetParameter(ZPrm); 119 // 120 G4UIparameter* APrm = new G4UIparameter("A", 'i', false); 121 APrm->SetGuidance("A"); 122 APrm->SetParameterRange("A>0"); 123 fIsotopeCmd->SetParameter(APrm); 124 // 125 G4UIparameter* densityPrm = new G4UIparameter("density", 'd', false); 126 densityPrm->SetGuidance("density of material"); 127 densityPrm->SetParameterRange("density>0."); 128 fIsotopeCmd->SetParameter(densityPrm); 129 // 130 G4UIparameter* unitPrm1 = new G4UIparameter("unit", 's', false); 131 unitPrm1->SetGuidance("unit of density"); 132 G4String unitList1 = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("g/cm3")); 133 unitPrm1->SetParameterCandidates(unitList1); 134 fIsotopeCmd->SetParameter(unitPrm1); 135 // 136 fIsotopeCmd->AvailableForStates(G4State_PreInit, G4State_Idle); 137 } 138 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 140 141 DetectorMessenger::~DetectorMessenger() 142 { 143 delete fSizeYZCmd; 144 delete fNbLayersCmd; 145 delete fNbAbsorCmd; 146 delete fAbsorCmd; 147 delete fIsotopeCmd; 148 delete fDetDir; 149 delete fTestemDir; 150 } 151 152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 153 154 void DetectorMessenger::SetNewValue(G4UIcommand* command, G4String newValue) 155 { 156 if (command == fSizeYZCmd) { 157 fDetector->SetCalorSizeYZ(fSizeYZCmd->GetNewDoubleValue(newValue)); 158 } 159 160 if (command == fNbLayersCmd) { 161 fDetector->SetNbOfLayers(fNbLayersCmd->GetNewIntValue(newValue)); 162 } 163 164 if (command == fNbAbsorCmd) { 165 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNewIntValue(newValue)); 166 } 167 168 if (command == fAbsorCmd) { 169 G4int num; 170 G4double tick; 171 G4String unt, mat; 172 std::istringstream is(newValue); 173 is >> num >> mat >> tick >> unt; 174 G4String material = mat; 175 tick *= G4UIcommand::ValueOf(unt); 176 fDetector->SetAbsorMaterial(num, material); 177 fDetector->SetAbsorThickness(num, tick); 178 } 179 180 if (command == fIsotopeCmd) { 181 G4int Z; 182 G4int A; 183 G4double dens; 184 G4String name, unt; 185 std::istringstream is(newValue); 186 is >> name >> Z >> A >> dens >> unt; 187 dens *= G4UIcommand::ValueOf(unt); 188 fDetector->MaterialWithSingleIsotope(name, name, dens, Z, A); 189 } 190 } 191 192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 193