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 // Hadrontherapy advanced example for Geant4 << 26 // HadrontherapyPhysicsListMessenger.cc 27 // See more at: https://twiki.cern.ch/twiki/bi << 27 // See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy 28 28 29 #include "HadrontherapyPhysicsListMessenger.hh 29 #include "HadrontherapyPhysicsListMessenger.hh" >> 30 30 #include "HadrontherapyPhysicsList.hh" 31 #include "HadrontherapyPhysicsList.hh" 31 #include "G4UIdirectory.hh" 32 #include "G4UIdirectory.hh" 32 #include "G4UIcmdWithADoubleAndUnit.hh" 33 #include "G4UIcmdWithADoubleAndUnit.hh" 33 #include "G4UIcmdWithAString.hh" 34 #include "G4UIcmdWithAString.hh" 34 35 35 ////////////////////////////////////////////// 36 ///////////////////////////////////////////////////////////////////////////// 36 HadrontherapyPhysicsListMessenger::Hadronthera 37 HadrontherapyPhysicsListMessenger::HadrontherapyPhysicsListMessenger(HadrontherapyPhysicsList* pPhys) 37 :pPhysicsList(pPhys) 38 :pPhysicsList(pPhys) 38 { 39 { 39 physDir = new G4UIdirectory("/Physics/"); << 40 physDir = new G4UIdirectory("/physic/"); 40 physDir->SetGuidance("Commands to activate p 41 physDir->SetGuidance("Commands to activate physics models and set cuts"); >> 42 >> 43 gammaCutCmd = new G4UIcmdWithADoubleAndUnit("/physic/setGCut",this); >> 44 gammaCutCmd->SetGuidance("Set gamma cut."); >> 45 gammaCutCmd->SetParameterName("Gcut",false); >> 46 gammaCutCmd->SetUnitCategory("Length"); >> 47 gammaCutCmd->SetRange("Gcut>0.0"); >> 48 gammaCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 49 >> 50 electCutCmd = new G4UIcmdWithADoubleAndUnit("/physic/setECut",this); >> 51 electCutCmd->SetGuidance("Set electron cut."); >> 52 electCutCmd->SetParameterName("Ecut",false); >> 53 electCutCmd->SetUnitCategory("Length"); >> 54 electCutCmd->SetRange("Ecut>0.0"); >> 55 electCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 41 56 42 pListCmd = new G4UIcmdWithAString("/Physics/ << 57 protoCutCmd = new G4UIcmdWithADoubleAndUnit("/physic/setPCut",this); >> 58 protoCutCmd->SetGuidance("Set positron cut."); >> 59 protoCutCmd->SetParameterName("Pcut",false); >> 60 protoCutCmd->SetUnitCategory("Length"); >> 61 protoCutCmd->SetRange("Pcut>0.0"); >> 62 protoCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 63 >> 64 allCutCmd = new G4UIcmdWithADoubleAndUnit("/physic/setCuts",this); >> 65 allCutCmd->SetGuidance("Set cut for all."); >> 66 allCutCmd->SetParameterName("cut",false); >> 67 allCutCmd->SetUnitCategory("Length"); >> 68 allCutCmd->SetRange("cut>0.0"); >> 69 allCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); >> 70 >> 71 pListCmd = new G4UIcmdWithAString("/physic/addPhysics",this); 43 pListCmd->SetGuidance("Add physics list."); 72 pListCmd->SetGuidance("Add physics list."); 44 pListCmd->SetParameterName("PList",false); 73 pListCmd->SetParameterName("PList",false); 45 pListCmd->AvailableForStates(G4State_PreInit << 74 pListCmd->AvailableForStates(G4State_PreInit); >> 75 >> 76 packageListCmd = new G4UIcmdWithAString("/physic/addPackage",this); >> 77 packageListCmd->SetGuidance("Add physics package."); >> 78 packageListCmd->SetParameterName("package",false); >> 79 packageListCmd->AvailableForStates(G4State_PreInit); 46 } 80 } 47 81 48 ////////////////////////////////////////////// 82 ///////////////////////////////////////////////////////////////////////////// 49 HadrontherapyPhysicsListMessenger::~Hadronther 83 HadrontherapyPhysicsListMessenger::~HadrontherapyPhysicsListMessenger() 50 { 84 { 51 delete physDir; << 85 delete gammaCutCmd; >> 86 delete electCutCmd; >> 87 delete protoCutCmd; >> 88 delete allCutCmd; >> 89 delete pListCmd; >> 90 delete physDir; >> 91 delete packageListCmd; 52 } 92 } 53 93 54 ////////////////////////////////////////////// 94 ///////////////////////////////////////////////////////////////////////////// 55 void HadrontherapyPhysicsListMessenger::SetNew 95 void HadrontherapyPhysicsListMessenger::SetNewValue(G4UIcommand* command, 56 G4St 96 G4String newValue) 57 { << 97 { 58 if( command == pListCmd ) << 98 if( command == gammaCutCmd ) >> 99 { pPhysicsList->SetCutForGamma(gammaCutCmd->GetNewDoubleValue(newValue));} >> 100 >> 101 if( command == electCutCmd ) >> 102 { pPhysicsList->SetCutForElectron(electCutCmd->GetNewDoubleValue(newValue));} >> 103 >> 104 if( command == protoCutCmd ) >> 105 { pPhysicsList->SetCutForPositron(protoCutCmd->GetNewDoubleValue(newValue));} >> 106 >> 107 if( command == allCutCmd ) >> 108 { >> 109 G4double cut = allCutCmd->GetNewDoubleValue(newValue); >> 110 pPhysicsList->SetCutForGamma(cut); >> 111 pPhysicsList->SetCutForElectron(cut); >> 112 pPhysicsList->SetCutForPositron(cut); >> 113 } >> 114 >> 115 if( command == pListCmd ) 59 { pPhysicsList->AddPhysicsList(newValue);} 116 { pPhysicsList->AddPhysicsList(newValue);} 60 } << 61 117 >> 118 >> 119 if( command == packageListCmd ) >> 120 { pPhysicsList->AddPackage(newValue);} >> 121 >> 122 >> 123 } 62 124