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 hadronic/Hadr01/src/DetectorMessenge 26 /// \file hadronic/Hadr01/src/DetectorMessenger.cc 27 /// \brief Implementation of the DetectorMesse 27 /// \brief Implementation of the DetectorMessenger class 28 // 28 // 29 // 29 // 30 ////////////////////////////////////////////// 30 ///////////////////////////////////////////////////////////////////////// 31 // 31 // 32 // DetectorMessenger 32 // DetectorMessenger 33 // 33 // 34 // Created: 31.01.03 V.Ivanchenko 34 // Created: 31.01.03 V.Ivanchenko 35 // 35 // 36 // Modified: 36 // Modified: 37 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 37 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko) 38 // 16.11.2006 Add beamCmd (V.Ivanchenko) 38 // 16.11.2006 Add beamCmd (V.Ivanchenko) 39 // 39 // 40 ////////////////////////////////////////////// 40 //////////////////////////////////////////////////////////////////////// 41 // 41 // 42 42 43 #include "DetectorMessenger.hh" 43 #include "DetectorMessenger.hh" 44 44 45 #include "DetectorConstruction.hh" 45 #include "DetectorConstruction.hh" 46 #include "HistoManager.hh" << 46 #include "G4UIdirectory.hh" 47 << 48 #include "G4HadronicParameters.hh" << 49 #include "G4UIcmdWith3Vector.hh" << 50 #include "G4UIcmdWithABool.hh" 47 #include "G4UIcmdWithABool.hh" 51 #include "G4UIcmdWithADoubleAndUnit.hh" << 52 #include "G4UIcmdWithAString.hh" 48 #include "G4UIcmdWithAString.hh" 53 #include "G4UIcmdWithAnInteger.hh" 49 #include "G4UIcmdWithAnInteger.hh" >> 50 #include "G4UIcmdWith3Vector.hh" >> 51 #include "G4UIcmdWithADoubleAndUnit.hh" 54 #include "G4UIcmdWithoutParameter.hh" 52 #include "G4UIcmdWithoutParameter.hh" 55 #include "G4UIdirectory.hh" << 53 #include "HistoManager.hh" 56 54 57 //....oooOO0OOooo........oooOO0OOooo........oo 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 58 56 59 DetectorMessenger::DetectorMessenger(DetectorC << 57 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) 60 : G4UImessenger(), << 58 : G4UImessenger(), fDetector(Det), 61 fDetector(Det), << 59 fTestDir(0), 62 fTestDir(0), << 60 fMatCmd(0), 63 fMatCmd(0), << 61 fMat1Cmd(0), 64 fMat1Cmd(0), << 62 fRCmd(0), 65 fRCmd(0), << 63 fLCmd(0), 66 fLCmd(0), << 64 fEdepCmd(0), 67 fEdepCmd(0), << 65 fBinCmd(0), 68 fBinCmd(0), << 66 fNOfAbsCmd(0), 69 fNOfAbsCmd(0), << 67 fVerbCmd(0), 70 fVerbCmd(0), << 68 fBeamCmd(0) 71 fBeamCmd(0) << 72 { 69 { 73 fTestDir = new G4UIdirectory("/testhadr/"); 70 fTestDir = new G4UIdirectory("/testhadr/"); 74 fTestDir->SetGuidance(" Hadronic Extended Ex 71 fTestDir->SetGuidance(" Hadronic Extended Example."); 75 72 76 fMatCmd = new G4UIcmdWithAString("/testhadr/ << 73 fMatCmd = new G4UIcmdWithAString("/testhadr/TargetMat",this); 77 fMatCmd->SetGuidance("Select Material for th 74 fMatCmd->SetGuidance("Select Material for the target"); 78 fMatCmd->SetParameterName("tMaterial", false << 75 fMatCmd->SetParameterName("tMaterial",false); 79 fMatCmd->AvailableForStates(G4State_PreInit, << 76 fMatCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 80 77 81 fMat1Cmd = new G4UIcmdWithAString("/testhadr << 78 fMat1Cmd = new G4UIcmdWithAString("/testhadr/WorldMat",this); 82 fMat1Cmd->SetGuidance("Select Material for w 79 fMat1Cmd->SetGuidance("Select Material for world"); 83 fMat1Cmd->SetParameterName("wMaterial", fals << 80 fMat1Cmd->SetParameterName("wMaterial",false); 84 fMat1Cmd->AvailableForStates(G4State_PreInit << 81 fMat1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 85 82 86 fRCmd = new G4UIcmdWithADoubleAndUnit("/test << 83 fRCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetRadius",this); 87 fRCmd->SetGuidance("Set radius of the target 84 fRCmd->SetGuidance("Set radius of the target"); 88 fRCmd->SetParameterName("radius", false); << 85 fRCmd->SetParameterName("radius",false); 89 fRCmd->SetUnitCategory("Length"); 86 fRCmd->SetUnitCategory("Length"); 90 fRCmd->SetRange("radius>0"); 87 fRCmd->SetRange("radius>0"); 91 fRCmd->AvailableForStates(G4State_PreInit); << 88 fRCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 92 89 93 fLCmd = new G4UIcmdWithADoubleAndUnit("/test << 90 fLCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetLength",this); 94 fLCmd->SetGuidance("Set length of the target 91 fLCmd->SetGuidance("Set length of the target"); 95 fLCmd->SetParameterName("length", false); << 92 fLCmd->SetParameterName("length",false); 96 fLCmd->SetUnitCategory("Length"); 93 fLCmd->SetUnitCategory("Length"); 97 fLCmd->SetRange("length>0"); 94 fLCmd->SetRange("length>0"); 98 fLCmd->AvailableForStates(G4State_PreInit); 95 fLCmd->AvailableForStates(G4State_PreInit); 99 96 100 fBinCmd = new G4UIcmdWithAnInteger("/testhad << 97 fBinCmd = new G4UIcmdWithAnInteger("/testhadr/NumberOfBinsE",this); 101 fBinCmd->SetGuidance("Set number of bins for 98 fBinCmd->SetGuidance("Set number of bins for Energy"); 102 fBinCmd->SetParameterName("NEbins", false); << 99 fBinCmd->SetParameterName("NEbins",false); 103 fBinCmd->AvailableForStates(G4State_PreInit) 100 fBinCmd->AvailableForStates(G4State_PreInit); 104 101 105 fNOfAbsCmd = new G4UIcmdWithAnInteger("/test << 102 fNOfAbsCmd = new G4UIcmdWithAnInteger("/testhadr/NumberDivZ",this); 106 fNOfAbsCmd->SetGuidance("Set number of slice 103 fNOfAbsCmd->SetGuidance("Set number of slices"); 107 fNOfAbsCmd->SetParameterName("NZ", false); << 104 fNOfAbsCmd->SetParameterName("NZ",false); 108 fNOfAbsCmd->AvailableForStates(G4State_PreIn 105 fNOfAbsCmd->AvailableForStates(G4State_PreInit); 109 106 110 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/t << 107 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/MaxEdep",this); 111 fEdepCmd->SetGuidance("Set max energy in his 108 fEdepCmd->SetGuidance("Set max energy in histogram"); 112 fEdepCmd->SetParameterName("edep", false); << 109 fEdepCmd->SetParameterName("edep",false); 113 fEdepCmd->SetUnitCategory("Energy"); 110 fEdepCmd->SetUnitCategory("Energy"); 114 fEdepCmd->SetRange("edep>0"); 111 fEdepCmd->SetRange("edep>0"); 115 fEdepCmd->AvailableForStates(G4State_PreInit << 112 fEdepCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 116 113 117 fBeamCmd = new G4UIcmdWithABool("/testhadr/D << 114 fBeamCmd = new G4UIcmdWithABool("/testhadr/DefaultBeamPosition",this); 118 fBeamCmd->SetGuidance("Use default beam posi 115 fBeamCmd->SetGuidance("Use default beam position"); 119 fBeamCmd->AvailableForStates(G4State_PreInit << 116 fBeamCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 120 117 121 fVerbCmd = new G4UIcmdWithAnInteger("/testha << 118 fVerbCmd = new G4UIcmdWithAnInteger("/testhadr/Verbose",this); 122 fVerbCmd->SetGuidance("Set verbose for "); 119 fVerbCmd->SetGuidance("Set verbose for "); 123 fVerbCmd->SetParameterName("verb", false); << 120 fVerbCmd->SetParameterName("verb",false); 124 fVerbCmd->AvailableForStates(G4State_PreInit << 121 fVerbCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 125 << 126 fBertCmd = new G4UIcmdWithAnInteger("/testha << 127 fBertCmd->SetGuidance("Set printout of Berti << 128 fBertCmd->SetParameterName("bert", false); << 129 fBertCmd->AvailableForStates(G4State_PreInit << 130 << 131 fBCCmd = new G4UIcmdWithABool("/testhadr/BCP << 132 fBCCmd->SetGuidance("Enable b-, c- hadronic << 133 fBCCmd->SetParameterName("BCparticles", fals << 134 fBCCmd->AvailableForStates(G4State_PreInit); << 135 << 136 fHNCmd = new G4UIcmdWithABool("/testhadr/Hyp << 137 fHNCmd->SetGuidance("Enable hyper-nuclei had << 138 fHNCmd->SetParameterName("HyperNuclei", fals << 139 fHNCmd->AvailableForStates(G4State_PreInit); << 140 << 141 fNGPCmd = new G4UIcmdWithABool("/testhadr/Ne << 142 fNGPCmd->SetGuidance("Enable neutron general << 143 fNGPCmd->SetParameterName("NeutronGeneralPro << 144 fNGPCmd->AvailableForStates(G4State_PreInit) << 145 } 122 } 146 123 147 //....oooOO0OOooo........oooOO0OOooo........oo 124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 148 125 149 DetectorMessenger::~DetectorMessenger() 126 DetectorMessenger::~DetectorMessenger() 150 { 127 { 151 delete fMatCmd; 128 delete fMatCmd; 152 delete fMat1Cmd; 129 delete fMat1Cmd; 153 delete fRCmd; 130 delete fRCmd; 154 delete fLCmd; 131 delete fLCmd; 155 delete fBinCmd; 132 delete fBinCmd; 156 delete fNOfAbsCmd; 133 delete fNOfAbsCmd; 157 delete fTestDir; 134 delete fTestDir; 158 delete fBeamCmd; 135 delete fBeamCmd; 159 delete fVerbCmd; 136 delete fVerbCmd; 160 delete fBertCmd; << 161 delete fEdepCmd; 137 delete fEdepCmd; 162 delete fBCCmd; << 163 delete fHNCmd; << 164 delete fNGPCmd; << 165 } 138 } 166 139 167 //....oooOO0OOooo........oooOO0OOooo........oo 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 168 141 169 void DetectorMessenger::SetNewValue(G4UIcomman 142 void DetectorMessenger::SetNewValue(G4UIcommand* command, G4String newValue) 170 { 143 { 171 HistoManager* h = HistoManager::GetPointer() 144 HistoManager* h = HistoManager::GetPointer(); 172 if (command == fMatCmd) { << 145 if( command == fMatCmd ) { 173 fDetector->SetTargetMaterial(newValue); 146 fDetector->SetTargetMaterial(newValue); 174 } << 147 } else if( command == fMat1Cmd ) { 175 else if (command == fMat1Cmd) { << 176 fDetector->SetWorldMaterial(newValue); 148 fDetector->SetWorldMaterial(newValue); 177 } << 149 } else if( command == fRCmd ) { 178 else if (command == fRCmd) { << 150 fDetector->SetTargetRadius(fRCmd->GetNewDoubleValue(newValue)); 179 h->SetTargetRadius(fRCmd->GetNewDoubleValu << 151 } else if( command == fLCmd ) { 180 } << 181 else if (command == fLCmd) { << 182 h->SetTargetLength(fLCmd->GetNewDoubleValu 152 h->SetTargetLength(fLCmd->GetNewDoubleValue(newValue)); 183 } << 153 } else if( command == fNOfAbsCmd ) { 184 else if (command == fNOfAbsCmd) { << 185 h->SetNumberOfSlices(fNOfAbsCmd->GetNewInt 154 h->SetNumberOfSlices(fNOfAbsCmd->GetNewIntValue(newValue)); 186 } << 155 } else if( command == fBinCmd ) { 187 else if (command == fBinCmd) { << 188 h->SetNumberOfBinsE(fBinCmd->GetNewIntValu 156 h->SetNumberOfBinsE(fBinCmd->GetNewIntValue(newValue)); 189 } << 157 } else if( command == fVerbCmd ) { 190 else if (command == fVerbCmd) { << 191 h->SetVerbose(fVerbCmd->GetNewIntValue(new 158 h->SetVerbose(fVerbCmd->GetNewIntValue(newValue)); 192 } << 159 } else if (command == fBeamCmd) { 193 else if (command == fBertCmd) { << 194 h->SetPrintBertiniXS(fBertCmd->GetNewIntVa << 195 } << 196 else if (command == fBeamCmd) { << 197 h->SetDefaultBeamPositionFlag(fBeamCmd->Ge 160 h->SetDefaultBeamPositionFlag(fBeamCmd->GetNewBoolValue(newValue)); 198 } << 161 } else if (command == fEdepCmd) { 199 else if (command == fEdepCmd) { << 200 h->SetMaxEnergyDeposit(fEdepCmd->GetNewDou 162 h->SetMaxEnergyDeposit(fEdepCmd->GetNewDoubleValue(newValue)); 201 } 163 } 202 else if (command == fBCCmd) { << 203 G4HadronicParameters::Instance()->SetEnabl << 204 } << 205 else if (command == fHNCmd) { << 206 G4HadronicParameters::Instance()->SetEnabl << 207 } << 208 else if (command == fNGPCmd) { << 209 G4HadronicParameters::Instance()->SetEnabl << 210 fNGPCmd->GetNewBoolValue(newValue)); << 211 } << 212 } 164 } 213 165 214 //....oooOO0OOooo........oooOO0OOooo........oo 166 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 167 215 168