Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/hadrontherapy/src/PassiveProtonBeamLineMessenger.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/advanced/hadrontherapy/src/PassiveProtonBeamLineMessenger.cc (Version 11.3.0) and /examples/advanced/hadrontherapy/src/PassiveProtonBeamLineMessenger.cc (Version 10.1.p1)


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