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 // passiveProtonBeamLineMessenger.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 "PassiveProtonBeamLineMessenger.hh" 29 #include "PassiveProtonBeamLineMessenger.hh" 30 #include "PassiveProtonBeamLine.hh" 30 #include "PassiveProtonBeamLine.hh" 31 #include "G4UIdirectory.hh" 31 #include "G4UIdirectory.hh" 32 #include "G4UIcmdWithADoubleAndUnit.hh" 32 #include "G4UIcmdWithADoubleAndUnit.hh" 33 #include "G4UIcmdWithAString.hh" 33 #include "G4UIcmdWithAString.hh" 34 #include "G4SystemOfUnits.hh" << 35 34 36 PassiveProtonBeamLineMessenger::PassivePro << 37 :passiveProton(beamLine) << 38 35 39 { << 36 PassiveProtonBeamLineMessenger::PassiveProtonBeamLineMessenger(PassiveProtonBeamLine* beamLine) 40 changeTheBeamLineDir = new G4UIdirectory(" << 37 :passiveProton(beamLine) 41 changeTheBeamLineDir -> SetGuidance("Comma << 42 38 >> 39 { >> 40 changeTheBeamLineDir = new G4UIdirectory("/ChangeBeamLine/"); >> 41 changeTheBeamLineDir -> SetGuidance("Command to change the transport beam line"); >> 42 43 changeTheBeamLineNameCmd = new G4UIcmdWith 43 changeTheBeamLineNameCmd = new G4UIcmdWithAString("/ChangeBeamLine/beamLineName",this); 44 changeTheBeamLineNameCmd -> SetGuidance("I << 44 changeTheBeamLineNameCmd -> SetGuidance("Insert the name of the beam line you want simulate"); 45 changeTheBeamLineNameCmd -> SetParameterNa << 45 changeTheBeamLineNameCmd -> SetParameterName("List",false); 46 changeTheBeamLineNameCmd -> AvailableForSt << 46 changeTheBeamLineNameCmd -> AvailableForStates(G4State_PreInit); 47 << 47 48 modulatorDir = new G4UIdirectory("/modulat << 48 modulatorDir = new G4UIdirectory("/modulator/"); 49 modulatorDir -> SetGuidance("Command to ro << 49 modulatorDir -> SetGuidance("Command to rotate the modulator wheel"); 50 << 50 51 beamLineDir = new G4UIdirectory("/beamLine << 51 beamLineDir = new G4UIdirectory("/beamLine/"); 52 beamLineDir -> SetGuidance("set specificat << 52 beamLineDir -> SetGuidance("set specification of range shifter"); 53 << 53 54 rangeShifterDir = new G4UIdirectory("/beam << 54 rangeShifterDir = new G4UIdirectory("/beamLine/RangeShifter/"); 55 rangeShifterDir -> SetGuidance("set specif << 55 rangeShifterDir -> SetGuidance("set specification of range shifter"); 56 << 56 57 firstScatteringFoilDir = new G4UIdirectory << 57 firstScatteringFoilDir = new G4UIdirectory("/beamLine/ScatteringFoil1/"); 58 firstScatteringFoilDir -> SetGuidance("set << 58 firstScatteringFoilDir -> SetGuidance("set specification of first scattering foil"); 59 << 59 60 secondScatteringFoilDir = new G4UIdirector << 60 secondScatteringFoilDir = new G4UIdirectory("/beamLine/ScatteringFoil2/"); 61 secondScatteringFoilDir -> SetGuidance("se << 61 secondScatteringFoilDir -> SetGuidance("set specification of second scattering foil"); 62 << 62 63 rangeStopperDir = new G4UIdirectory("/beam << 63 rangeStopperDir = new G4UIdirectory("/beamLine/Stopper/"); 64 rangeStopperDir -> SetGuidance("set specif << 64 rangeStopperDir -> SetGuidance("set specification of stopper"); 65 << 65 66 finalCollimatorDir = new G4UIdirectory("/b << 66 finalCollimatorDir = new G4UIdirectory("/beamLine/FinalCollimator/"); 67 finalCollimatorDir -> SetGuidance("set spe << 67 finalCollimatorDir -> SetGuidance("set specification of final collimator"); 68 << 68 69 modulatorAngleCmd = new G4UIcmdWithADouble << 69 modulatorAngleCmd = new G4UIcmdWithADoubleAndUnit("/modulator/angle",this); 70 modulatorAngleCmd -> SetGuidance("Set Modu << 70 modulatorAngleCmd -> SetGuidance("Set Modulator Angle"); 71 modulatorAngleCmd -> SetParameterName("Siz << 71 modulatorAngleCmd -> SetParameterName("Size",false); 72 modulatorAngleCmd -> SetRange("Size>=0."); << 72 modulatorAngleCmd -> SetRange("Size>=0."); 73 modulatorAngleCmd -> SetUnitCategory("Angl << 73 modulatorAngleCmd -> SetUnitCategory("Angle"); 74 modulatorAngleCmd -> AvailableForStates(G4 << 74 modulatorAngleCmd -> AvailableForStates(G4State_Idle); 75 << 75 76 rangeShifterMatCmd = new G4UIcmdWithAStrin << 76 rangeShifterMatCmd = new G4UIcmdWithAString("/beamLine/RangeShifter/RSMat",this); 77 rangeShifterMatCmd -> SetGuidance("Set mat << 77 rangeShifterMatCmd -> SetGuidance("Set material of range shifter"); 78 rangeShifterMatCmd -> SetParameterName("ch << 78 rangeShifterMatCmd -> SetParameterName("choice",false); 79 rangeShifterMatCmd -> AvailableForStates(G << 79 rangeShifterMatCmd -> AvailableForStates(G4State_Idle); 80 << 80 81 rangeShifterXSizeCmd = new G4UIcmdWithADou << 81 rangeShifterXSizeCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/RangeShifter/thickness",this); 82 rangeShifterXSizeCmd -> SetGuidance("Set h << 82 rangeShifterXSizeCmd -> SetGuidance("Set half of the thickness of range shifter along X axis"); 83 rangeShifterXSizeCmd -> SetParameterName(" << 83 rangeShifterXSizeCmd -> SetParameterName("Size",false); 84 rangeShifterXSizeCmd -> SetDefaultUnit("mm << 84 rangeShifterXSizeCmd -> SetDefaultUnit("mm"); 85 rangeShifterXSizeCmd -> SetUnitCandidates( << 85 rangeShifterXSizeCmd -> SetUnitCandidates("mm cm m"); 86 rangeShifterXSizeCmd -> AvailableForStates << 86 rangeShifterXSizeCmd -> AvailableForStates(G4State_Idle); 87 << 87 88 rangeShifterXPositionCmd = new G4UIcmdWith << 88 rangeShifterXPositionCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/RangeShifter/position",this); 89 rangeShifterXPositionCmd -> SetGuidance("S << 89 rangeShifterXPositionCmd -> SetGuidance("Set position of range shifter"); 90 rangeShifterXPositionCmd -> SetParameterNa << 90 rangeShifterXPositionCmd -> SetParameterName("Size",false); 91 rangeShifterXPositionCmd -> SetDefaultUnit << 91 rangeShifterXPositionCmd -> SetDefaultUnit("mm"); 92 rangeShifterXPositionCmd -> SetUnitCandida << 92 rangeShifterXPositionCmd -> SetUnitCandidates("mm cm m"); 93 rangeShifterXPositionCmd -> AvailableForSt << 93 rangeShifterXPositionCmd -> AvailableForStates(G4State_Idle); 94 << 94 95 firstScatteringFoilXSizeCmd = new G4UIcmdW << 95 firstScatteringFoilXSizeCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/ScatteringFoil1/thickness",this); 96 firstScatteringFoilXSizeCmd -> SetGuidance << 96 firstScatteringFoilXSizeCmd -> SetGuidance("Set half thickness of first scattering foil"); 97 firstScatteringFoilXSizeCmd -> SetParamete << 97 firstScatteringFoilXSizeCmd -> SetParameterName("Size",false); 98 firstScatteringFoilXSizeCmd -> SetDefaultU << 98 firstScatteringFoilXSizeCmd -> SetDefaultUnit("mm"); 99 firstScatteringFoilXSizeCmd -> SetUnitCand << 99 firstScatteringFoilXSizeCmd -> SetUnitCandidates("mm cm m"); 100 firstScatteringFoilXSizeCmd -> AvailableFo << 100 firstScatteringFoilXSizeCmd -> AvailableForStates(G4State_Idle); 101 << 101 102 secondScatteringFoilXSizeCmd = new G4UIcmd << 102 secondScatteringFoilXSizeCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/ScatteringFoil2/thickness",this); 103 secondScatteringFoilXSizeCmd -> SetGuidanc << 103 secondScatteringFoilXSizeCmd -> SetGuidance("Set half thickness of second scattering foil"); 104 secondScatteringFoilXSizeCmd -> SetParamet << 104 secondScatteringFoilXSizeCmd -> SetParameterName("Size",false); 105 secondScatteringFoilXSizeCmd -> SetDefault << 105 secondScatteringFoilXSizeCmd -> SetDefaultUnit("mm"); 106 secondScatteringFoilXSizeCmd -> SetUnitCan << 106 secondScatteringFoilXSizeCmd -> SetUnitCandidates("mm cm m"); 107 secondScatteringFoilXSizeCmd -> AvailableF << 107 secondScatteringFoilXSizeCmd -> AvailableForStates(G4State_Idle); 108 << 108 109 outerRadiusStopperCmd = new G4UIcmdWithADo << 109 outerRadiusStopperCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/Stopper/outRadius",this); 110 outerRadiusStopperCmd -> SetGuidance("Set << 110 outerRadiusStopperCmd -> SetGuidance("Set size of outer radius"); 111 outerRadiusStopperCmd -> SetParameterName( << 111 outerRadiusStopperCmd -> SetParameterName("Size",false); 112 outerRadiusStopperCmd -> SetDefaultUnit("m << 112 outerRadiusStopperCmd -> SetDefaultUnit("mm"); 113 outerRadiusStopperCmd -> SetUnitCandidates << 113 outerRadiusStopperCmd -> SetUnitCandidates("mm cm m"); 114 outerRadiusStopperCmd -> AvailableForState << 114 outerRadiusStopperCmd -> AvailableForStates(G4State_Idle); 115 << 115 116 innerRadiusFinalCollimatorCmd = new G4UIcm << 116 innerRadiusFinalCollimatorCmd = new G4UIcmdWithADoubleAndUnit("/beamLine/FinalCollimator/halfInnerRad",this); 117 innerRadiusFinalCollimatorCmd -> SetGuidan << 117 innerRadiusFinalCollimatorCmd -> SetGuidance("Set size of inner radius ( max 21.5 mm)"); 118 innerRadiusFinalCollimatorCmd -> SetParame << 118 innerRadiusFinalCollimatorCmd -> SetParameterName("Size",false); 119 innerRadiusFinalCollimatorCmd -> SetDefaul << 119 innerRadiusFinalCollimatorCmd -> SetDefaultUnit("mm"); 120 innerRadiusFinalCollimatorCmd -> SetUnitCa << 120 innerRadiusFinalCollimatorCmd -> SetUnitCandidates("mm cm m"); 121 innerRadiusFinalCollimatorCmd -> Available << 121 innerRadiusFinalCollimatorCmd -> AvailableForStates(G4State_Idle); 122 } 122 } 123 123 124 PassiveProtonBeamLineMessenger::~PassiveProton 124 PassiveProtonBeamLineMessenger::~PassiveProtonBeamLineMessenger() 125 { 125 { 126 delete innerRadiusFinalCollimatorCmd; << 126 delete innerRadiusFinalCollimatorCmd; 127 delete outerRadiusStopperCmd; << 127 delete outerRadiusStopperCmd; 128 delete secondScatteringFoilXSizeCmd; << 128 delete secondScatteringFoilXSizeCmd; 129 delete firstScatteringFoilXSizeCmd; << 129 delete firstScatteringFoilXSizeCmd; 130 delete rangeShifterXPositionCmd; << 130 delete rangeShifterXPositionCmd; 131 delete rangeShifterXSizeCmd; << 131 delete rangeShifterXSizeCmd; 132 delete rangeShifterMatCmd; << 132 delete rangeShifterMatCmd; 133 delete modulatorAngleCmd; << 133 delete modulatorAngleCmd; 134 delete finalCollimatorDir; << 134 delete finalCollimatorDir; 135 delete rangeStopperDir; << 135 delete rangeStopperDir; 136 delete secondScatteringFoilDir; << 136 delete secondScatteringFoilDir; 137 delete firstScatteringFoilDir; << 137 delete firstScatteringFoilDir; 138 delete rangeShifterDir; << 138 delete rangeShifterDir; 139 delete beamLineDir; << 139 delete beamLineDir; 140 delete modulatorDir; << 140 delete modulatorDir; 141 } 141 } 142 142 143 143 144 144 145 145 146 void PassiveProtonBeamLineMessenger::SetNewVal 146 void PassiveProtonBeamLineMessenger::SetNewValue(G4UIcommand* command,G4String newValue) 147 { 147 { 148 if( command == modulatorAngleCmd ) << 148 if( command == modulatorAngleCmd ) 149 { passiveProton -> SetModulatorAngle 149 { passiveProton -> SetModulatorAngle 150 (modulatorAngleCmd -> GetNewDoubleValue(newV << 150 (modulatorAngleCmd -> GetNewDoubleValue(newValue));} 151 151 152 else if( command == rangeShifterMatCmd ) << 152 else if( command == rangeShifterMatCmd ) 153 { passiveProton -> SetRSMaterial(newValue) 153 { passiveProton -> SetRSMaterial(newValue);} 154 154 155 else if( command == rangeShifterXSizeCmd ) << 155 else if( command == rangeShifterXSizeCmd ) 156 { passiveProton -> SetRangeShifterXSize 156 { passiveProton -> SetRangeShifterXSize 157 (rangeShifterXSizeCmd -> GetNewDoubleValue(n << 157 (rangeShifterXSizeCmd -> GetNewDoubleValue(newValue));} 158 158 159 else if( command == rangeShifterXPositionC << 159 else if( command == rangeShifterXPositionCmd ) 160 { passiveProton -> SetRangeShifterXPositio 160 { passiveProton -> SetRangeShifterXPosition 161 (rangeShifterXPositionCmd -> GetNewDoubleVal << 161 (rangeShifterXPositionCmd -> GetNewDoubleValue(newValue));} 162 162 163 else if( command == firstScatteringFoilXSi << 163 else if( command == firstScatteringFoilXSizeCmd ) 164 { passiveProton -> SetFirstScatteringFoilX 164 { passiveProton -> SetFirstScatteringFoilXSize 165 (firstScatteringFoilXSizeCmd -> GetNewDouble << 165 (firstScatteringFoilXSizeCmd -> GetNewDoubleValue(newValue));} 166 166 167 else if( command == secondScatteringFoilXS << 167 else if( command == secondScatteringFoilXSizeCmd ) 168 { passiveProton -> SetSecondScatteringFoil 168 { passiveProton -> SetSecondScatteringFoilXSize 169 (secondScatteringFoilXSizeCmd -> GetNewDoubl << 169 (secondScatteringFoilXSizeCmd -> GetNewDoubleValue(newValue));} 170 170 171 else if( command == outerRadiusStopperCmd << 171 else if( command == outerRadiusStopperCmd ) 172 { passiveProton -> SetOuterRadiusStopper( 172 { passiveProton -> SetOuterRadiusStopper( 173 outerRadiusStopperCmd -> GetNewDoubleVal << 173 outerRadiusStopperCmd -> GetNewDoubleValue(newValue));} 174 174 175 else if( command == innerRadiusFinalCollim << 175 else if( command == innerRadiusFinalCollimatorCmd ) 176 { passiveProton -> SetInnerRadiusFinalColl 176 { passiveProton -> SetInnerRadiusFinalCollimator 177 (innerRadiusFinalCollimatorCmd -> GetNewDoub << 177 (innerRadiusFinalCollimatorCmd -> GetNewDoubleValue(newValue));} 178 } 178 } 179 179 180 180