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 electromagnetic/TestEm3/src/Detector 26 /// \file electromagnetic/TestEm3/src/DetectorMessenger.cc 27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class 28 // 28 // >> 29 // $Id: DetectorMessenger.cc 67268 2013-02-13 11:38:40Z 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 <sstream> 36 37 37 #include "G4UIcmdWithADoubleAndUnit.hh" << 38 #include "DetectorConstruction.hh" 38 #include "G4UIcmdWithAnInteger.hh" << 39 #include "G4UIcmdWithoutParameter.hh" << 40 #include "G4UIcommand.hh" << 41 #include "G4UIdirectory.hh" 39 #include "G4UIdirectory.hh" >> 40 #include "G4UIcommand.hh" 42 #include "G4UIparameter.hh" 41 #include "G4UIparameter.hh" 43 << 42 #include "G4UIcmdWithAnInteger.hh" 44 #include <sstream> << 43 #include "G4UIcmdWithADoubleAndUnit.hh" >> 44 #include "G4UIcmdWithoutParameter.hh" 45 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 47 48 DetectorMessenger::DetectorMessenger(DetectorC << 48 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) 49 { << 49 :G4UImessenger(),Detector(Det), >> 50 fTestemDir(0), >> 51 fDetDir(0), >> 52 fSizeYZCmd(0), >> 53 fNbLayersCmd(0), >> 54 fNbAbsorCmd(0), >> 55 fAbsorCmd(0), >> 56 fMagFieldCmd(0), >> 57 fUpdateCmd(0) >> 58 { 50 fTestemDir = new G4UIdirectory("/testem/"); 59 fTestemDir = new G4UIdirectory("/testem/"); 51 fTestemDir->SetGuidance("UI commands specifi 60 fTestemDir->SetGuidance("UI commands specific to this example"); 52 << 61 53 fDetDir = new G4UIdirectory("/testem/det/"); 62 fDetDir = new G4UIdirectory("/testem/det/"); 54 fDetDir->SetGuidance("detector construction 63 fDetDir->SetGuidance("detector construction commands"); 55 << 64 56 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" << 65 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setSizeYZ",this); 57 fSizeYZCmd->SetGuidance("Set tranverse size 66 fSizeYZCmd->SetGuidance("Set tranverse size of the calorimeter"); 58 fSizeYZCmd->SetParameterName("Size", false); << 67 fSizeYZCmd->SetParameterName("Size",false); 59 fSizeYZCmd->SetRange("Size>0."); 68 fSizeYZCmd->SetRange("Size>0."); 60 fSizeYZCmd->SetUnitCategory("Length"); 69 fSizeYZCmd->SetUnitCategory("Length"); 61 fSizeYZCmd->AvailableForStates(G4State_PreIn << 70 fSizeYZCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 62 fSizeYZCmd->SetToBeBroadcasted(false); << 71 63 << 72 fNbLayersCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfLayers",this); 64 fNbLayersCmd = new G4UIcmdWithAnInteger("/te << 65 fNbLayersCmd->SetGuidance("Set number of lay 73 fNbLayersCmd->SetGuidance("Set number of layers."); 66 fNbLayersCmd->SetParameterName("NbLayers", f << 74 fNbLayersCmd->SetParameterName("NbLayers",false); 67 fNbLayersCmd->SetRange("NbLayers>0"); 75 fNbLayersCmd->SetRange("NbLayers>0"); 68 fNbLayersCmd->AvailableForStates(G4State_Pre << 76 fNbLayersCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 69 fNbLayersCmd->SetToBeBroadcasted(false); << 77 70 << 78 fNbAbsorCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfAbsor",this); 71 fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes << 72 fNbAbsorCmd->SetGuidance("Set number of Abso 79 fNbAbsorCmd->SetGuidance("Set number of Absorbers."); 73 fNbAbsorCmd->SetParameterName("NbAbsor", fal << 80 fNbAbsorCmd->SetParameterName("NbAbsor",false); 74 fNbAbsorCmd->SetRange("NbAbsor>0"); 81 fNbAbsorCmd->SetRange("NbAbsor>0"); 75 fNbAbsorCmd->AvailableForStates(G4State_PreI << 82 fNbAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 76 fNbAbsorCmd->SetToBeBroadcasted(false); << 83 77 << 84 fAbsorCmd = new G4UIcommand("/testem/det/setAbsor",this); 78 fAbsorCmd = new G4UIcommand("/testem/det/set << 79 fAbsorCmd->SetGuidance("Set the absor nb, th 85 fAbsorCmd->SetGuidance("Set the absor nb, the material, the thickness."); 80 fAbsorCmd->SetGuidance(" absor number : fro 86 fAbsorCmd->SetGuidance(" absor number : from 1 to NbOfAbsor"); 81 fAbsorCmd->SetGuidance(" material name"); 87 fAbsorCmd->SetGuidance(" material name"); 82 fAbsorCmd->SetGuidance(" thickness (with un 88 fAbsorCmd->SetGuidance(" thickness (with unit) : t>0."); 83 // 89 // 84 G4UIparameter* AbsNbPrm = new G4UIparameter( << 90 G4UIparameter* AbsNbPrm = new G4UIparameter("AbsorNb",'i',false); 85 AbsNbPrm->SetGuidance("absor number : from 1 91 AbsNbPrm->SetGuidance("absor number : from 1 to NbOfAbsor"); 86 AbsNbPrm->SetParameterRange("AbsorNb>0"); 92 AbsNbPrm->SetParameterRange("AbsorNb>0"); 87 fAbsorCmd->SetParameter(AbsNbPrm); 93 fAbsorCmd->SetParameter(AbsNbPrm); 88 // 94 // 89 G4UIparameter* MatPrm = new G4UIparameter("m << 95 G4UIparameter* MatPrm = new G4UIparameter("material",'s',false); 90 MatPrm->SetGuidance("material name"); 96 MatPrm->SetGuidance("material name"); 91 fAbsorCmd->SetParameter(MatPrm); 97 fAbsorCmd->SetParameter(MatPrm); 92 // << 98 // 93 G4UIparameter* ThickPrm = new G4UIparameter( << 99 G4UIparameter* ThickPrm = new G4UIparameter("thickness",'d',false); 94 ThickPrm->SetGuidance("thickness of absorber 100 ThickPrm->SetGuidance("thickness of absorber"); 95 ThickPrm->SetParameterRange("thickness>0."); 101 ThickPrm->SetParameterRange("thickness>0."); 96 fAbsorCmd->SetParameter(ThickPrm); 102 fAbsorCmd->SetParameter(ThickPrm); 97 // 103 // 98 G4UIparameter* unitPrm = new G4UIparameter(" << 104 G4UIparameter* unitPrm = new G4UIparameter("unit",'s',false); 99 unitPrm->SetGuidance("unit of thickness"); 105 unitPrm->SetGuidance("unit of thickness"); 100 G4String unitList = G4UIcommand::UnitsList(G 106 G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("mm")); 101 unitPrm->SetParameterCandidates(unitList); 107 unitPrm->SetParameterCandidates(unitList); 102 fAbsorCmd->SetParameter(unitPrm); 108 fAbsorCmd->SetParameter(unitPrm); 103 // 109 // 104 fAbsorCmd->AvailableForStates(G4State_PreIni << 110 fAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 105 fAbsorCmd->SetToBeBroadcasted(false); << 111 >> 112 fMagFieldCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setField",this); >> 113 fMagFieldCmd->SetGuidance("Define magnetic field."); >> 114 fMagFieldCmd->SetGuidance("Magnetic field will be in Z direction."); >> 115 fMagFieldCmd->SetParameterName("Bz",false); >> 116 fMagFieldCmd->SetUnitCategory("Magnetic flux density"); >> 117 fMagFieldCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 118 >> 119 fUpdateCmd = new G4UIcmdWithoutParameter("/testem/det/update",this); >> 120 fUpdateCmd->SetGuidance("Update calorimeter geometry."); >> 121 fUpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" "); >> 122 fUpdateCmd->SetGuidance("if you changed geometrical value(s)."); >> 123 fUpdateCmd->AvailableForStates(G4State_Idle); 106 } 124 } 107 125 108 //....oooOO0OOooo........oooOO0OOooo........oo 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 109 127 110 DetectorMessenger::~DetectorMessenger() 128 DetectorMessenger::~DetectorMessenger() 111 { 129 { 112 delete fSizeYZCmd; 130 delete fSizeYZCmd; 113 delete fNbLayersCmd; 131 delete fNbLayersCmd; 114 delete fNbAbsorCmd; 132 delete fNbAbsorCmd; 115 delete fAbsorCmd; 133 delete fAbsorCmd; 116 delete fDetDir; << 134 delete fMagFieldCmd; >> 135 delete fUpdateCmd; >> 136 delete fDetDir; 117 delete fTestemDir; 137 delete fTestemDir; 118 } 138 } 119 139 120 //....oooOO0OOooo........oooOO0OOooo........oo 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 121 141 122 void DetectorMessenger::SetNewValue(G4UIcomman << 142 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue) 123 { 143 { 124 if (command == fSizeYZCmd) { << 144 if( command == fSizeYZCmd ) 125 fDetector->SetCalorSizeYZ(fSizeYZCmd->GetN << 145 { Detector->SetCalorSizeYZ(fSizeYZCmd->GetNewDoubleValue(newValue));} 126 } << 146 127 << 147 if( command == fNbLayersCmd ) 128 if (command == fNbLayersCmd) { << 148 { Detector->SetNbOfLayers(fNbLayersCmd->GetNewIntValue(newValue));} 129 fDetector->SetNbOfLayers(fNbLayersCmd->Get << 149 130 } << 150 if( command == fNbAbsorCmd ) 131 << 151 { Detector->SetNbOfAbsor(fNbAbsorCmd->GetNewIntValue(newValue));} 132 if (command == fNbAbsorCmd) { << 152 133 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe << 153 if (command == fAbsorCmd) 134 } << 154 { 135 << 155 G4int num; G4double tick; 136 if (command == fAbsorCmd) { << 156 G4String unt, mat; 137 G4int num; << 157 std::istringstream is(newValue); 138 G4double tick; << 158 is >> num >> mat >> tick >> unt; 139 G4String unt, mat; << 159 G4String material=mat; 140 std::istringstream is(newValue); << 160 tick *= G4UIcommand::ValueOf(unt); 141 is >> num >> mat >> tick >> unt; << 161 Detector->SetAbsorMaterial (num,material); 142 G4String material = mat; << 162 Detector->SetAbsorThickness(num,tick); 143 tick *= G4UIcommand::ValueOf(unt); << 163 } 144 fDetector->SetAbsorMaterial(num, material) << 164 145 fDetector->SetAbsorThickness(num, tick); << 165 if( command == fMagFieldCmd ) 146 } << 166 { Detector->SetMagField(fMagFieldCmd->GetNewDoubleValue(newValue));} >> 167 >> 168 if( command == fUpdateCmd ) >> 169 { Detector->UpdateGeometry();} 147 } 170 } 148 171 149 //....oooOO0OOooo........oooOO0OOooo........oo 172 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 150 173