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 // $Id: DetectorMessenger.cc 77255 2013-11-22 10:09:14Z gcosmo $ 29 // 30 // 30 ////////////////////////////////////////////// 31 ///////////////////////////////////////////////////////////////////////// 31 // 32 // 32 // DetectorMessenger 33 // DetectorMessenger 33 // 34 // 34 // Created: 31.01.03 V.Ivanchenko 35 // Created: 31.01.03 V.Ivanchenko 35 // 36 // 36 // Modified: 37 // Modified: 37 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko) 38 // 16.11.2006 Add beamCmd (V.Ivanchenko) 39 // 16.11.2006 Add beamCmd (V.Ivanchenko) 39 // 40 // 40 ////////////////////////////////////////////// 41 //////////////////////////////////////////////////////////////////////// 41 // 42 // 42 43 43 #include "DetectorMessenger.hh" 44 #include "DetectorMessenger.hh" 44 45 45 #include "DetectorConstruction.hh" 46 #include "DetectorConstruction.hh" 46 #include "HistoManager.hh" << 47 #include "G4UIdirectory.hh" 47 << 48 #include "G4HadronicParameters.hh" << 49 #include "G4UIcmdWith3Vector.hh" << 50 #include "G4UIcmdWithABool.hh" 48 #include "G4UIcmdWithABool.hh" 51 #include "G4UIcmdWithADoubleAndUnit.hh" << 52 #include "G4UIcmdWithAString.hh" 49 #include "G4UIcmdWithAString.hh" 53 #include "G4UIcmdWithAnInteger.hh" 50 #include "G4UIcmdWithAnInteger.hh" >> 51 #include "G4UIcmdWith3Vector.hh" >> 52 #include "G4UIcmdWithADoubleAndUnit.hh" 54 #include "G4UIcmdWithoutParameter.hh" 53 #include "G4UIcmdWithoutParameter.hh" 55 #include "G4UIdirectory.hh" << 54 #include "HistoManager.hh" 56 55 57 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 58 57 59 DetectorMessenger::DetectorMessenger(DetectorC << 58 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det) 60 : G4UImessenger(), << 59 : G4UImessenger(), fDetector(Det), 61 fDetector(Det), << 60 fTestDir(0), 62 fTestDir(0), << 61 fMatCmd(0), 63 fMatCmd(0), << 62 fMat1Cmd(0), 64 fMat1Cmd(0), << 63 fRCmd(0), 65 fRCmd(0), << 64 fLCmd(0), 66 fLCmd(0), << 65 fEdepCmd(0), 67 fEdepCmd(0), << 66 fBinCmd(0), 68 fBinCmd(0), << 67 fNOfAbsCmd(0), 69 fNOfAbsCmd(0), << 68 fVerbCmd(0), 70 fVerbCmd(0), << 69 fBeamCmd(0) 71 fBeamCmd(0) << 72 { 70 { 73 fTestDir = new G4UIdirectory("/testhadr/"); 71 fTestDir = new G4UIdirectory("/testhadr/"); 74 fTestDir->SetGuidance(" Hadronic Extended Ex 72 fTestDir->SetGuidance(" Hadronic Extended Example."); 75 73 76 fMatCmd = new G4UIcmdWithAString("/testhadr/ << 74 fMatCmd = new G4UIcmdWithAString("/testhadr/TargetMat",this); 77 fMatCmd->SetGuidance("Select Material for th 75 fMatCmd->SetGuidance("Select Material for the target"); 78 fMatCmd->SetParameterName("tMaterial", false << 76 fMatCmd->SetParameterName("tMaterial",false); 79 fMatCmd->AvailableForStates(G4State_PreInit, << 77 fMatCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 80 78 81 fMat1Cmd = new G4UIcmdWithAString("/testhadr << 79 fMat1Cmd = new G4UIcmdWithAString("/testhadr/WorldMat",this); 82 fMat1Cmd->SetGuidance("Select Material for w 80 fMat1Cmd->SetGuidance("Select Material for world"); 83 fMat1Cmd->SetParameterName("wMaterial", fals << 81 fMat1Cmd->SetParameterName("wMaterial",false); 84 fMat1Cmd->AvailableForStates(G4State_PreInit << 82 fMat1Cmd->AvailableForStates(G4State_PreInit,G4State_Idle); 85 83 86 fRCmd = new G4UIcmdWithADoubleAndUnit("/test << 84 fRCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetRadius",this); 87 fRCmd->SetGuidance("Set radius of the target 85 fRCmd->SetGuidance("Set radius of the target"); 88 fRCmd->SetParameterName("radius", false); << 86 fRCmd->SetParameterName("radius",false); 89 fRCmd->SetUnitCategory("Length"); 87 fRCmd->SetUnitCategory("Length"); 90 fRCmd->SetRange("radius>0"); 88 fRCmd->SetRange("radius>0"); 91 fRCmd->AvailableForStates(G4State_PreInit); << 89 fRCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 92 90 93 fLCmd = new G4UIcmdWithADoubleAndUnit("/test << 91 fLCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/TargetLength",this); 94 fLCmd->SetGuidance("Set length of the target 92 fLCmd->SetGuidance("Set length of the target"); 95 fLCmd->SetParameterName("length", false); << 93 fLCmd->SetParameterName("length",false); 96 fLCmd->SetUnitCategory("Length"); 94 fLCmd->SetUnitCategory("Length"); 97 fLCmd->SetRange("length>0"); 95 fLCmd->SetRange("length>0"); 98 fLCmd->AvailableForStates(G4State_PreInit); << 96 fLCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 99 97 100 fBinCmd = new G4UIcmdWithAnInteger("/testhad << 98 fBinCmd = new G4UIcmdWithAnInteger("/testhadr/NumberOfBinsE",this); 101 fBinCmd->SetGuidance("Set number of bins for 99 fBinCmd->SetGuidance("Set number of bins for Energy"); 102 fBinCmd->SetParameterName("NEbins", false); << 100 fBinCmd->SetParameterName("NEbins",false); 103 fBinCmd->AvailableForStates(G4State_PreInit) << 101 fBinCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 104 102 105 fNOfAbsCmd = new G4UIcmdWithAnInteger("/test << 103 fNOfAbsCmd = new G4UIcmdWithAnInteger("/testhadr/NumberDivZ",this); 106 fNOfAbsCmd->SetGuidance("Set number of slice 104 fNOfAbsCmd->SetGuidance("Set number of slices"); 107 fNOfAbsCmd->SetParameterName("NZ", false); << 105 fNOfAbsCmd->SetParameterName("NZ",false); 108 fNOfAbsCmd->AvailableForStates(G4State_PreIn << 106 fNOfAbsCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 109 107 110 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/t << 108 fEdepCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/MaxEdep",this); 111 fEdepCmd->SetGuidance("Set max energy in his 109 fEdepCmd->SetGuidance("Set max energy in histogram"); 112 fEdepCmd->SetParameterName("edep", false); << 110 fEdepCmd->SetParameterName("edep",false); 113 fEdepCmd->SetUnitCategory("Energy"); 111 fEdepCmd->SetUnitCategory("Energy"); 114 fEdepCmd->SetRange("edep>0"); 112 fEdepCmd->SetRange("edep>0"); 115 fEdepCmd->AvailableForStates(G4State_PreInit << 113 fEdepCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 116 114 117 fBeamCmd = new G4UIcmdWithABool("/testhadr/D << 115 fBeamCmd = new G4UIcmdWithABool("/testhadr/DefaultBeamPosition",this); 118 fBeamCmd->SetGuidance("Use default beam posi << 116 fBeamCmd->SetGuidance("show inelastic and elastic cross sections"); 119 fBeamCmd->AvailableForStates(G4State_PreInit << 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) { << 148 fDetector->SetWorldMaterial(newValue); 176 fDetector->SetWorldMaterial(newValue); << 149 } else if( command == fRCmd ) { 177 } << 150 fDetector->SetTargetRadius(fRCmd->GetNewDoubleValue(newValue)); 178 else if (command == fRCmd) { << 151 } else if( command == fLCmd ) { 179 h->SetTargetRadius(fRCmd->GetNewDoubleValu << 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