Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 /// \file electromagnetic/TestEm3/src/Detector << 23 // $Id: DetectorMessenger.cc,v 1.10 2005/11/17 08:16:23 gcosmo Exp $ 27 /// \brief Implementation of the DetectorMesse << 24 // GEANT4 tag $Name: geant4-08-00 $ 28 // << 29 // 25 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 26 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 28 33 #include "DetectorMessenger.hh" 29 #include "DetectorMessenger.hh" 34 30 35 #include "DetectorConstruction.hh" << 31 #include <sstream> 36 32 37 #include "G4UIcmdWithADoubleAndUnit.hh" << 33 #include "DetectorConstruction.hh" 38 #include "G4UIcmdWithAnInteger.hh" << 39 #include "G4UIcmdWithoutParameter.hh" << 40 #include "G4UIcommand.hh" << 41 #include "G4UIdirectory.hh" 34 #include "G4UIdirectory.hh" >> 35 #include "G4UIcommand.hh" 42 #include "G4UIparameter.hh" 36 #include "G4UIparameter.hh" 43 << 37 #include "G4UIcmdWithAnInteger.hh" 44 #include <sstream> << 38 #include "G4UIcmdWithADoubleAndUnit.hh" >> 39 #include "G4UIcmdWithoutParameter.hh" 45 40 46 //....oooOO0OOooo........oooOO0OOooo........oo 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 42 48 DetectorMessenger::DetectorMessenger(DetectorC << 43 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) 49 { << 44 :Detector(Det) 50 fTestemDir = new G4UIdirectory("/testem/"); << 45 { 51 fTestemDir->SetGuidance("UI commands specifi << 46 testemDir = new G4UIdirectory("/testem/"); 52 << 47 testemDir->SetGuidance("UI commands specific to this example"); 53 fDetDir = new G4UIdirectory("/testem/det/"); << 48 54 fDetDir->SetGuidance("detector construction << 49 detDir = new G4UIdirectory("/testem/det/"); 55 << 50 detDir->SetGuidance("detector construction commands"); 56 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" << 51 57 fSizeYZCmd->SetGuidance("Set tranverse size << 52 SizeYZCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setSizeYZ",this); 58 fSizeYZCmd->SetParameterName("Size", false); << 53 SizeYZCmd->SetGuidance("Set tranverse size of the calorimeter"); 59 fSizeYZCmd->SetRange("Size>0."); << 54 SizeYZCmd->SetParameterName("Size",false); 60 fSizeYZCmd->SetUnitCategory("Length"); << 55 SizeYZCmd->SetRange("Size>0."); 61 fSizeYZCmd->AvailableForStates(G4State_PreIn << 56 SizeYZCmd->SetUnitCategory("Length"); 62 fSizeYZCmd->SetToBeBroadcasted(false); << 57 SizeYZCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 63 << 58 64 fNbLayersCmd = new G4UIcmdWithAnInteger("/te << 59 NbLayersCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfLayers",this); 65 fNbLayersCmd->SetGuidance("Set number of lay << 60 NbLayersCmd->SetGuidance("Set number of layers."); 66 fNbLayersCmd->SetParameterName("NbLayers", f << 61 NbLayersCmd->SetParameterName("NbLayers",false); 67 fNbLayersCmd->SetRange("NbLayers>0"); << 62 NbLayersCmd->SetRange("NbLayers>0"); 68 fNbLayersCmd->AvailableForStates(G4State_Pre << 63 NbLayersCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 69 fNbLayersCmd->SetToBeBroadcasted(false); << 64 70 << 65 NbAbsorCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfAbsor",this); 71 fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes << 66 NbAbsorCmd->SetGuidance("Set number of Absorbers."); 72 fNbAbsorCmd->SetGuidance("Set number of Abso << 67 NbAbsorCmd->SetParameterName("NbAbsor",false); 73 fNbAbsorCmd->SetParameterName("NbAbsor", fal << 68 NbAbsorCmd->SetRange("NbAbsor>0"); 74 fNbAbsorCmd->SetRange("NbAbsor>0"); << 69 NbAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 75 fNbAbsorCmd->AvailableForStates(G4State_PreI << 70 76 fNbAbsorCmd->SetToBeBroadcasted(false); << 71 AbsorCmd = new G4UIcommand("/testem/det/setAbsor",this); 77 << 72 AbsorCmd->SetGuidance("Set the absor nb, the material, the thickness."); 78 fAbsorCmd = new G4UIcommand("/testem/det/set << 73 AbsorCmd->SetGuidance(" absor number : from 1 to NbOfAbsor"); 79 fAbsorCmd->SetGuidance("Set the absor nb, th << 74 AbsorCmd->SetGuidance(" material name"); 80 fAbsorCmd->SetGuidance(" absor number : fro << 75 AbsorCmd->SetGuidance(" thickness (with unit) : t>0."); 81 fAbsorCmd->SetGuidance(" material name"); << 82 fAbsorCmd->SetGuidance(" thickness (with un << 83 // 76 // 84 G4UIparameter* AbsNbPrm = new G4UIparameter( << 77 G4UIparameter* AbsNbPrm = new G4UIparameter("AbsorNb",'i',false); 85 AbsNbPrm->SetGuidance("absor number : from 1 78 AbsNbPrm->SetGuidance("absor number : from 1 to NbOfAbsor"); 86 AbsNbPrm->SetParameterRange("AbsorNb>0"); 79 AbsNbPrm->SetParameterRange("AbsorNb>0"); 87 fAbsorCmd->SetParameter(AbsNbPrm); << 80 AbsorCmd->SetParameter(AbsNbPrm); 88 // 81 // 89 G4UIparameter* MatPrm = new G4UIparameter("m << 82 G4UIparameter* MatPrm = new G4UIparameter("material",'s',false); 90 MatPrm->SetGuidance("material name"); 83 MatPrm->SetGuidance("material name"); 91 fAbsorCmd->SetParameter(MatPrm); << 84 AbsorCmd->SetParameter(MatPrm); 92 // << 85 // 93 G4UIparameter* ThickPrm = new G4UIparameter( << 86 G4UIparameter* ThickPrm = new G4UIparameter("thickness",'d',false); 94 ThickPrm->SetGuidance("thickness of absorber 87 ThickPrm->SetGuidance("thickness of absorber"); 95 ThickPrm->SetParameterRange("thickness>0."); 88 ThickPrm->SetParameterRange("thickness>0."); 96 fAbsorCmd->SetParameter(ThickPrm); << 89 AbsorCmd->SetParameter(ThickPrm); 97 // 90 // 98 G4UIparameter* unitPrm = new G4UIparameter(" << 91 G4UIparameter* unitPrm = new G4UIparameter("unit",'s',false); 99 unitPrm->SetGuidance("unit of thickness"); 92 unitPrm->SetGuidance("unit of thickness"); 100 G4String unitList = G4UIcommand::UnitsList(G 93 G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("mm")); 101 unitPrm->SetParameterCandidates(unitList); 94 unitPrm->SetParameterCandidates(unitList); 102 fAbsorCmd->SetParameter(unitPrm); << 95 AbsorCmd->SetParameter(unitPrm); 103 // 96 // 104 fAbsorCmd->AvailableForStates(G4State_PreIni << 97 AbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 105 fAbsorCmd->SetToBeBroadcasted(false); << 98 >> 99 MagFieldCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setField",this); >> 100 MagFieldCmd->SetGuidance("Define magnetic field."); >> 101 MagFieldCmd->SetGuidance("Magnetic field will be in Z direction."); >> 102 MagFieldCmd->SetParameterName("Bz",false); >> 103 MagFieldCmd->SetUnitCategory("Magnetic flux density"); >> 104 MagFieldCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 105 >> 106 UpdateCmd = new G4UIcmdWithoutParameter("/testem/det/update",this); >> 107 UpdateCmd->SetGuidance("Update calorimeter geometry."); >> 108 UpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" "); >> 109 UpdateCmd->SetGuidance("if you changed geometrical value(s)."); >> 110 UpdateCmd->AvailableForStates(G4State_Idle); 106 } 111 } 107 112 108 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 109 114 110 DetectorMessenger::~DetectorMessenger() 115 DetectorMessenger::~DetectorMessenger() 111 { 116 { 112 delete fSizeYZCmd; << 117 delete SizeYZCmd; 113 delete fNbLayersCmd; << 118 delete NbLayersCmd; 114 delete fNbAbsorCmd; << 119 delete NbAbsorCmd; 115 delete fAbsorCmd; << 120 delete AbsorCmd; 116 delete fDetDir; << 121 delete MagFieldCmd; 117 delete fTestemDir; << 122 delete UpdateCmd; >> 123 delete detDir; >> 124 delete testemDir; 118 } 125 } 119 126 120 //....oooOO0OOooo........oooOO0OOooo........oo 127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 121 128 122 void DetectorMessenger::SetNewValue(G4UIcomman << 129 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue) 123 { 130 { 124 if (command == fSizeYZCmd) { << 131 if( command == SizeYZCmd ) 125 fDetector->SetCalorSizeYZ(fSizeYZCmd->GetN << 132 { Detector->SetCalorSizeYZ(SizeYZCmd->GetNewDoubleValue(newValue));} 126 } << 133 127 << 134 if( command == NbLayersCmd ) 128 if (command == fNbLayersCmd) { << 135 { Detector->SetNbOfLayers(NbLayersCmd->GetNewIntValue(newValue));} 129 fDetector->SetNbOfLayers(fNbLayersCmd->Get << 136 130 } << 137 if( command == NbAbsorCmd ) 131 << 138 { Detector->SetNbOfAbsor(NbAbsorCmd->GetNewIntValue(newValue));} 132 if (command == fNbAbsorCmd) { << 139 133 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe << 140 if (command == AbsorCmd) 134 } << 141 { 135 << 142 G4int num; G4double tick; 136 if (command == fAbsorCmd) { << 143 G4String unt, mat; 137 G4int num; << 144 std::istringstream is(newValue); 138 G4double tick; << 145 is >> num >> mat >> tick >> unt; 139 G4String unt, mat; << 146 G4String material=mat; 140 std::istringstream is(newValue); << 147 tick *= G4UIcommand::ValueOf(unt); 141 is >> num >> mat >> tick >> unt; << 148 Detector->SetAbsorMaterial (num,material); 142 G4String material = mat; << 149 Detector->SetAbsorThickness(num,tick); 143 tick *= G4UIcommand::ValueOf(unt); << 150 } 144 fDetector->SetAbsorMaterial(num, material) << 151 145 fDetector->SetAbsorThickness(num, tick); << 152 if( command == MagFieldCmd ) 146 } << 153 { Detector->SetMagField(MagFieldCmd->GetNewDoubleValue(newValue));} >> 154 >> 155 if( command == UpdateCmd ) >> 156 { Detector->UpdateGeometry();} 147 } 157 } 148 158 149 //....oooOO0OOooo........oooOO0OOooo........oo 159 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 150 160