Geant4 Cross Reference |
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 // 26 // This is the *BASIC* version of IORT, a Geant4-based application 27 // 28 // Main Authors: G.Russo(a,b), C.Casarino*(c), G.C. Candiano(c), G.A.P. Cirrone(d), F.Romano(d) 29 // Contributor Authors: S.Guatelli(e) 30 // Past Authors: G.Arnetta(c), S.E.Mazzaglia(d) 31 // 32 // (a) Fondazione Istituto San Raffaele G.Giglio, Cefalù, Italy 33 // (b) IBFM-CNR , Segrate (Milano), Italy 34 // (c) LATO (Laboratorio di Tecnologie Oncologiche), Cefalù, Italy 35 // (d) Laboratori Nazionali del Sud of the INFN, Catania, Italy 36 // (e) University of Wollongong, Australia 37 // 38 // *Corresponding author, email to carlo.casarino@polooncologicocefalu.it 39 ////////////////////////////////////////////////////////////////////////////////////////////// 40 41 #include "IORTPhysicsListMessenger.hh" 42 43 #include "IORTPhysicsList.hh" 44 #include "G4UIdirectory.hh" 45 #include "G4UIcmdWithADoubleAndUnit.hh" 46 #include "G4UIcmdWithAString.hh" 47 48 ///////////////////////////////////////////////////////////////////////////// 49 IORTPhysicsListMessenger::IORTPhysicsListMessenger(IORTPhysicsList* pPhys) 50 :pPhysicsList(pPhys) 51 { 52 physDir = new G4UIdirectory("/Physics/"); 53 physDir->SetGuidance("Commands to activate physics models and set cuts"); 54 55 gammaCutCmd = new G4UIcmdWithADoubleAndUnit("/Physics/setGCut",this); 56 gammaCutCmd->SetGuidance("Set gamma cut."); 57 gammaCutCmd->SetParameterName("Gcut",false); 58 gammaCutCmd->SetUnitCategory("Length"); 59 gammaCutCmd->SetRange("Gcut>0.0"); 60 gammaCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 61 62 electCutCmd = new G4UIcmdWithADoubleAndUnit("/Physics/setECut",this); 63 electCutCmd->SetGuidance("Set electron cut."); 64 electCutCmd->SetParameterName("Ecut",false); 65 electCutCmd->SetUnitCategory("Length"); 66 electCutCmd->SetRange("Ecut>0.0"); 67 electCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 68 69 protoCutCmd = new G4UIcmdWithADoubleAndUnit("/Physics/setPCut",this); 70 protoCutCmd->SetGuidance("Set positron cut."); 71 protoCutCmd->SetParameterName("Pcut",false); 72 protoCutCmd->SetUnitCategory("Length"); 73 protoCutCmd->SetRange("Pcut>0.0"); 74 protoCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 75 76 allCutCmd = new G4UIcmdWithADoubleAndUnit("/Physics/setCuts",this); 77 allCutCmd->SetGuidance("Set cut for all."); 78 allCutCmd->SetParameterName("cut",false); 79 allCutCmd->SetUnitCategory("Length"); 80 allCutCmd->SetRange("cut>0.0"); 81 allCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 82 83 allDetectorCmd = new G4UIcmdWithADoubleAndUnit("/Physics/setDetectorCuts",this); 84 allDetectorCmd->SetGuidance("Set cut for all. into Detector"); 85 allDetectorCmd->SetParameterName("cut",false); 86 allDetectorCmd->SetUnitCategory("Length"); 87 allDetectorCmd->SetRange("cut>0.0"); 88 allDetectorCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 89 90 pListCmd = new G4UIcmdWithAString("/Physics/addPhysics",this); 91 pListCmd->SetGuidance("Add physics list."); 92 pListCmd->SetParameterName("PList",false); 93 pListCmd->AvailableForStates(G4State_PreInit, G4State_Idle); 94 } 95 96 ///////////////////////////////////////////////////////////////////////////// 97 IORTPhysicsListMessenger::~IORTPhysicsListMessenger() 98 { 99 delete gammaCutCmd; 100 delete electCutCmd; 101 delete protoCutCmd; 102 delete allCutCmd; 103 delete allDetectorCmd; 104 delete pListCmd; 105 delete physDir; 106 } 107 108 ///////////////////////////////////////////////////////////////////////////// 109 void IORTPhysicsListMessenger::SetNewValue(G4UIcommand* command, 110 G4String newValue) 111 { 112 if( command == gammaCutCmd ) 113 { pPhysicsList->SetCutForGamma(gammaCutCmd->GetNewDoubleValue(newValue));} 114 115 else if( command == electCutCmd ) 116 { pPhysicsList->SetCutForElectron(electCutCmd->GetNewDoubleValue(newValue));} 117 118 else if( command == protoCutCmd ) 119 { pPhysicsList->SetCutForPositron(protoCutCmd->GetNewDoubleValue(newValue));} 120 121 else if( command == allCutCmd ) 122 { 123 G4double cut = allCutCmd->GetNewDoubleValue(newValue); 124 pPhysicsList->SetCutForGamma(cut); 125 pPhysicsList->SetCutForElectron(cut); 126 pPhysicsList->SetCutForPositron(cut); 127 } 128 else if( command == allDetectorCmd) 129 { 130 G4double cut = allDetectorCmd -> GetNewDoubleValue(newValue); 131 pPhysicsList -> SetDetectorCut(cut); 132 } 133 else if( command == pListCmd ) 134 { pPhysicsList->AddPhysicsList(newValue);} 135 } 136