[807] | 1 | // |
---|
| 2 | // ******************************************************************** |
---|
| 3 | // * License and Disclaimer * |
---|
| 4 | // * * |
---|
| 5 | // * The Geant4 software is copyright of the Copyright Holders of * |
---|
| 6 | // * the Geant4 Collaboration. It is provided under the terms and * |
---|
| 7 | // * conditions of the Geant4 Software License, included in the file * |
---|
| 8 | // * LICENSE and available at http://cern.ch/geant4/license . These * |
---|
| 9 | // * include a list of copyright holders. * |
---|
| 10 | // * * |
---|
| 11 | // * Neither the authors of this software system, nor their employing * |
---|
| 12 | // * institutes,nor the agencies providing financial support for this * |
---|
| 13 | // * work make any representation or warranty, express or implied, * |
---|
| 14 | // * regarding this software system or assume any liability for its * |
---|
| 15 | // * use. Please see the license in the file LICENSE and URL above * |
---|
| 16 | // * for the full disclaimer and the limitation of liability. * |
---|
| 17 | // * * |
---|
| 18 | // * This code implementation is the result of the scientific and * |
---|
| 19 | // * technical work of the GEANT4 collaboration. * |
---|
| 20 | // * By using, copying, modifying or distributing the software (or * |
---|
| 21 | // * any work based on the software) you agree to acknowledge its * |
---|
| 22 | // * use in resulting scientific publications, and indicate your * |
---|
| 23 | // * acceptance of all terms of the Geant4 Software license. * |
---|
| 24 | // ******************************************************************** |
---|
| 25 | // |
---|
[1230] | 26 | // HadrontherapyPhysicsListMessenger.cc |
---|
| 27 | // See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy |
---|
| 28 | |
---|
[807] | 29 | #include "HadrontherapyPhysicsListMessenger.hh" |
---|
[1230] | 30 | |
---|
[807] | 31 | #include "HadrontherapyPhysicsList.hh" |
---|
| 32 | #include "G4UIdirectory.hh" |
---|
| 33 | #include "G4UIcmdWithADoubleAndUnit.hh" |
---|
| 34 | #include "G4UIcmdWithAString.hh" |
---|
| 35 | |
---|
[1230] | 36 | ///////////////////////////////////////////////////////////////////////////// |
---|
| 37 | HadrontherapyPhysicsListMessenger::HadrontherapyPhysicsListMessenger(HadrontherapyPhysicsList* pPhys) |
---|
| 38 | :pPhysicsList(pPhys) |
---|
| 39 | { |
---|
| 40 | physDir = new G4UIdirectory("/physic/"); |
---|
| 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); |
---|
[807] | 49 | |
---|
[1230] | 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); |
---|
| 56 | |
---|
| 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 | |
---|
[1313] | 71 | allDetectorCmd = new G4UIcmdWithADoubleAndUnit("/physic/setDetectorCuts",this); |
---|
| 72 | allDetectorCmd->SetGuidance("Set cut for all. into Detector"); |
---|
| 73 | allDetectorCmd->SetParameterName("cut",false); |
---|
| 74 | allDetectorCmd->SetUnitCategory("Length"); |
---|
| 75 | allDetectorCmd->SetRange("cut>0.0"); |
---|
| 76 | allDetectorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); |
---|
| 77 | |
---|
[1230] | 78 | pListCmd = new G4UIcmdWithAString("/physic/addPhysics",this); |
---|
| 79 | pListCmd->SetGuidance("Add physics list."); |
---|
| 80 | pListCmd->SetParameterName("PList",false); |
---|
[1313] | 81 | pListCmd->AvailableForStates(G4State_PreInit, G4State_Idle); |
---|
[1230] | 82 | |
---|
| 83 | packageListCmd = new G4UIcmdWithAString("/physic/addPackage",this); |
---|
| 84 | packageListCmd->SetGuidance("Add physics package."); |
---|
| 85 | packageListCmd->SetParameterName("package",false); |
---|
[1313] | 86 | packageListCmd->AvailableForStates(G4State_PreInit, G4State_Idle); |
---|
[807] | 87 | } |
---|
| 88 | |
---|
[1230] | 89 | ///////////////////////////////////////////////////////////////////////////// |
---|
[807] | 90 | HadrontherapyPhysicsListMessenger::~HadrontherapyPhysicsListMessenger() |
---|
| 91 | { |
---|
[1230] | 92 | delete gammaCutCmd; |
---|
| 93 | delete electCutCmd; |
---|
| 94 | delete protoCutCmd; |
---|
| 95 | delete allCutCmd; |
---|
[1313] | 96 | delete allDetectorCmd; |
---|
[1230] | 97 | delete pListCmd; |
---|
| 98 | delete physDir; |
---|
| 99 | delete packageListCmd; |
---|
[807] | 100 | } |
---|
| 101 | |
---|
[1230] | 102 | ///////////////////////////////////////////////////////////////////////////// |
---|
| 103 | void HadrontherapyPhysicsListMessenger::SetNewValue(G4UIcommand* command, |
---|
| 104 | G4String newValue) |
---|
| 105 | { |
---|
| 106 | if( command == gammaCutCmd ) |
---|
| 107 | { pPhysicsList->SetCutForGamma(gammaCutCmd->GetNewDoubleValue(newValue));} |
---|
| 108 | |
---|
[1313] | 109 | else if( command == electCutCmd ) |
---|
[1230] | 110 | { pPhysicsList->SetCutForElectron(electCutCmd->GetNewDoubleValue(newValue));} |
---|
| 111 | |
---|
[1313] | 112 | else if( command == protoCutCmd ) |
---|
[1230] | 113 | { pPhysicsList->SetCutForPositron(protoCutCmd->GetNewDoubleValue(newValue));} |
---|
[807] | 114 | |
---|
[1313] | 115 | else if( command == allCutCmd ) |
---|
[1230] | 116 | { |
---|
| 117 | G4double cut = allCutCmd->GetNewDoubleValue(newValue); |
---|
| 118 | pPhysicsList->SetCutForGamma(cut); |
---|
| 119 | pPhysicsList->SetCutForElectron(cut); |
---|
| 120 | pPhysicsList->SetCutForPositron(cut); |
---|
| 121 | } |
---|
[1313] | 122 | else if( command == allDetectorCmd) |
---|
| 123 | { |
---|
| 124 | G4double cut = allDetectorCmd -> GetNewDoubleValue(newValue); |
---|
| 125 | pPhysicsList -> SetDetectorCut(cut); |
---|
| 126 | } |
---|
| 127 | else if( command == pListCmd ) |
---|
[1230] | 128 | { pPhysicsList->AddPhysicsList(newValue);} |
---|
[807] | 129 | |
---|
| 130 | |
---|
[1313] | 131 | else if( command == packageListCmd ) |
---|
[1230] | 132 | { pPhysicsList->AddPackage(newValue);} |
---|
[807] | 133 | |
---|
| 134 | |
---|
[1230] | 135 | } |
---|