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/Hadr00/src/DetectorMessenger.cc 27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class 28 // 28 // 29 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo << 30 // $Id: DetectorMessenger.cc 77210 2013-11-22 01:58:38Z adotti $ 31 //....oooOO0OOooo........oooOO0OOooo........oo << 31 // >> 32 // >> 33 ///////////////////////////////////////////////////////////////////////// >> 34 // >> 35 // DetectorMessenger >> 36 // >> 37 // Created: 20.06.08 V.Ivanchenko >> 38 // >> 39 // Modified: >> 40 // >> 41 //////////////////////////////////////////////////////////////////////// >> 42 // 32 43 33 #include "DetectorMessenger.hh" 44 #include "DetectorMessenger.hh" 34 45 35 #include "DetectorConstruction.hh" 46 #include "DetectorConstruction.hh" 36 << 47 #include "G4UIdirectory.hh" 37 #include "G4UIcmdWithADoubleAndUnit.hh" << 48 #include "G4UIcmdWithABool.hh" >> 49 #include "G4UIcmdWithAString.hh" 38 #include "G4UIcmdWithAnInteger.hh" 50 #include "G4UIcmdWithAnInteger.hh" >> 51 #include "G4UIcmdWith3Vector.hh" >> 52 #include "G4UIcmdWithADoubleAndUnit.hh" 39 #include "G4UIcmdWithoutParameter.hh" 53 #include "G4UIcmdWithoutParameter.hh" 40 #include "G4UIcommand.hh" << 54 #include "HistoManager.hh" 41 #include "G4UIdirectory.hh" << 42 #include "G4UIparameter.hh" << 43 << 44 #include <sstream> << 45 55 46 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 57 48 DetectorMessenger::DetectorMessenger(DetectorC << 58 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) >> 59 :G4UImessenger(), fDetector(Det) 49 { 60 { 50 fTestemDir = new G4UIdirectory("/testhadr/") << 61 ftestDir = new G4UIdirectory("/testhadr/"); 51 fTestemDir->SetGuidance("UI commands specifi << 62 ftestDir->SetGuidance(" Hadronic Extended Example."); 52 63 53 fDetDir = new G4UIdirectory("/testhadr/det/" << 64 fmatCmd = new G4UIcmdWithAString("/testhadr/TargetMat",this); 54 fDetDir->SetGuidance("detector construction << 65 fmatCmd->SetGuidance("Select Material for the target"); >> 66 fmatCmd->SetParameterName("tMaterial",false); >> 67 fmatCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 68 fmatCmd->SetToBeBroadcasted(false); >> 69 >> 70 fmat1Cmd = new G4UIcmdWithAString("/testhadr/WorldMat",this); >> 71 fmat1Cmd->SetGuidance("Select Material for world"); >> 72 fmat1Cmd->SetParameterName("wMaterial",false); >> 73 fmat1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 74 fmat1Cmd->SetToBeBroadcasted(false); >> 75 >> 76 frCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetRadius",this); >> 77 frCmd->SetGuidance("Set radius of the target"); >> 78 frCmd->SetParameterName("radius",false); >> 79 frCmd->SetUnitCategory("Length"); >> 80 frCmd->SetRange("radius>0"); >> 81 frCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 82 frCmd->SetToBeBroadcasted(false); >> 83 >> 84 flCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetLength",this); >> 85 flCmd->SetGuidance("Set length of the target"); >> 86 flCmd->SetParameterName("length",false); >> 87 flCmd->SetUnitCategory("Length"); >> 88 flCmd->SetRange("length>0"); >> 89 flCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 90 flCmd->SetToBeBroadcasted(false); 55 91 56 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" << 57 fSizeYZCmd->SetGuidance("Set tranverse size << 58 fSizeYZCmd->SetParameterName("Size", false); << 59 fSizeYZCmd->SetRange("Size>0."); << 60 fSizeYZCmd->SetUnitCategory("Length"); << 61 fSizeYZCmd->AvailableForStates(G4State_PreIn << 62 fSizeYZCmd->SetToBeBroadcasted(false); << 63 << 64 fNbLayersCmd = new G4UIcmdWithAnInteger("/te << 65 fNbLayersCmd->SetGuidance("Set number of lay << 66 fNbLayersCmd->SetParameterName("NbLayers", f << 67 fNbLayersCmd->SetRange("NbLayers>0"); << 68 fNbLayersCmd->AvailableForStates(G4State_Pre << 69 fNbLayersCmd->SetToBeBroadcasted(false); << 70 << 71 fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes << 72 fNbAbsorCmd->SetGuidance("Set number of Abso << 73 fNbAbsorCmd->SetParameterName("NbAbsor", fal << 74 fNbAbsorCmd->SetRange("NbAbsor>0"); << 75 fNbAbsorCmd->AvailableForStates(G4State_PreI << 76 fNbAbsorCmd->SetToBeBroadcasted(false); << 77 << 78 fAbsorCmd = new G4UIcommand("/testhadr/det/s << 79 fAbsorCmd->SetGuidance("Set the absor nb, th << 80 fAbsorCmd->SetGuidance(" absor number : fro << 81 fAbsorCmd->SetGuidance(" material name"); << 82 fAbsorCmd->SetGuidance(" thickness (with un << 83 // << 84 G4UIparameter* AbsNbPrm = new G4UIparameter( << 85 AbsNbPrm->SetGuidance("absor number : from 1 << 86 AbsNbPrm->SetParameterRange("AbsorNb>0"); << 87 fAbsorCmd->SetParameter(AbsNbPrm); << 88 // << 89 G4UIparameter* MatPrm = new G4UIparameter("m << 90 MatPrm->SetGuidance("material name"); << 91 fAbsorCmd->SetParameter(MatPrm); << 92 // << 93 G4UIparameter* ThickPrm = new G4UIparameter( << 94 ThickPrm->SetGuidance("thickness of absorber << 95 ThickPrm->SetParameterRange("thickness>0."); << 96 fAbsorCmd->SetParameter(ThickPrm); << 97 // << 98 G4UIparameter* unitPrm = new G4UIparameter(" << 99 unitPrm->SetGuidance("unit of thickness"); << 100 G4String unitList = G4UIcommand::UnitsList(G << 101 unitPrm->SetParameterCandidates(unitList); << 102 fAbsorCmd->SetParameter(unitPrm); << 103 // << 104 fAbsorCmd->AvailableForStates(G4State_PreIni << 105 fAbsorCmd->SetToBeBroadcasted(false); << 106 << 107 fIsotopeCmd = new G4UIcommand("/testhadr/det << 108 fIsotopeCmd->SetGuidance("Build and select a << 109 fIsotopeCmd->SetGuidance(" symbol of isotop << 110 // << 111 G4UIparameter* symbPrm = new G4UIparameter(" << 112 symbPrm->SetGuidance("isotope symbol"); << 113 fIsotopeCmd->SetParameter(symbPrm); << 114 // << 115 G4UIparameter* ZPrm = new G4UIparameter("Z", << 116 ZPrm->SetGuidance("Z"); << 117 ZPrm->SetParameterRange("Z>0"); << 118 fIsotopeCmd->SetParameter(ZPrm); << 119 // << 120 G4UIparameter* APrm = new G4UIparameter("A", << 121 APrm->SetGuidance("A"); << 122 APrm->SetParameterRange("A>0"); << 123 fIsotopeCmd->SetParameter(APrm); << 124 // << 125 G4UIparameter* densityPrm = new G4UIparamete << 126 densityPrm->SetGuidance("density of material << 127 densityPrm->SetParameterRange("density>0."); << 128 fIsotopeCmd->SetParameter(densityPrm); << 129 // << 130 G4UIparameter* unitPrm1 = new G4UIparameter( << 131 unitPrm1->SetGuidance("unit of density"); << 132 G4String unitList1 = G4UIcommand::UnitsList( << 133 unitPrm1->SetParameterCandidates(unitList1); << 134 fIsotopeCmd->SetParameter(unitPrm1); << 135 // << 136 fIsotopeCmd->AvailableForStates(G4State_PreI << 137 } 92 } 138 93 139 //....oooOO0OOooo........oooOO0OOooo........oo 94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 140 95 141 DetectorMessenger::~DetectorMessenger() 96 DetectorMessenger::~DetectorMessenger() 142 { 97 { 143 delete fSizeYZCmd; << 98 delete fmatCmd; 144 delete fNbLayersCmd; << 99 delete fmat1Cmd; 145 delete fNbAbsorCmd; << 100 delete frCmd; 146 delete fAbsorCmd; << 101 delete flCmd; 147 delete fIsotopeCmd; << 102 delete ftestDir; 148 delete fDetDir; << 149 delete fTestemDir; << 150 } 103 } 151 104 152 //....oooOO0OOooo........oooOO0OOooo........oo 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 153 106 154 void DetectorMessenger::SetNewValue(G4UIcomman 107 void DetectorMessenger::SetNewValue(G4UIcommand* command, G4String newValue) 155 { 108 { 156 if (command == fSizeYZCmd) { << 109 if( command == fmatCmd ) { 157 fDetector->SetCalorSizeYZ(fSizeYZCmd->GetN << 110 fDetector->SetTargetMaterial(newValue); 158 } << 111 } else if( command == fmat1Cmd ) { 159 << 112 fDetector->SetWorldMaterial(newValue); 160 if (command == fNbLayersCmd) { << 113 } else if( command == frCmd ) { 161 fDetector->SetNbOfLayers(fNbLayersCmd->Get << 114 fDetector->SetTargetRadius(frCmd->GetNewDoubleValue(newValue)); 162 } << 115 } else if( command == flCmd ) { 163 << 116 fDetector->SetTargetLength(flCmd->GetNewDoubleValue(newValue)); 164 if (command == fNbAbsorCmd) { << 165 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe << 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, << 189 } 117 } 190 } 118 } 191 119 192 //....oooOO0OOooo........oooOO0OOooo........oo 120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 121 193 122