Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/hadrontherapy/src/HadrontherapyDetectorMessenger.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/HadrontherapyDetectorMessenger.cc (Version 11.3.0) and /examples/advanced/hadrontherapy/src/HadrontherapyDetectorMessenger.cc (Version 10.4.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 //                                                 25 //
 26 // Hadrontherapy advanced example for Geant4       26 // Hadrontherapy advanced example for Geant4
 27 // See more at: https://twiki.cern.ch/twiki/bi     27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
 28                                                    28 
 29 #include "HadrontherapyDetectorMessenger.hh"       29 #include "HadrontherapyDetectorMessenger.hh"
 30 #include "HadrontherapyDetectorConstruction.hh     30 #include "HadrontherapyDetectorConstruction.hh"
 31 #include "G4UIdirectory.hh"                        31 #include "G4UIdirectory.hh"
 32 #include "G4UIcmdWith3VectorAndUnit.hh"            32 #include "G4UIcmdWith3VectorAndUnit.hh"
 33 #include "G4UIcmdWithoutParameter.hh"              33 #include "G4UIcmdWithoutParameter.hh"
 34 #include "G4UIcmdWithAString.hh"                   34 #include "G4UIcmdWithAString.hh"
 35 #include "G4SystemOfUnits.hh"                      35 #include "G4SystemOfUnits.hh"
 36 #include "G4UIcmdWithABool.hh"                     36 #include "G4UIcmdWithABool.hh"
 37                                                    37 
 38 //////////////////////////////////////////////     38 /////////////////////////////////////////////////////////////////////////////
 39 HadrontherapyDetectorMessenger::HadrontherapyD     39 HadrontherapyDetectorMessenger::HadrontherapyDetectorMessenger(HadrontherapyDetectorConstruction* detector)
 40 :hadrontherapyDetector(detector)               <<  40   :hadrontherapyDetector(detector)
 41 {                                                  41 {
 42     // Change Phantom size                         42     // Change Phantom size
 43     changeThePhantomDir = new G4UIdirectory("/     43     changeThePhantomDir = new G4UIdirectory("/changePhantom/");
 44     changeThePhantomDir -> SetGuidance("Comman     44     changeThePhantomDir -> SetGuidance("Command to change the Phantom Size/position");
 45     changeThePhantomSizeCmd = new G4UIcmdWith3     45     changeThePhantomSizeCmd = new G4UIcmdWith3VectorAndUnit("/changePhantom/size", this);
 46     changeThePhantomSizeCmd -> SetGuidance("In     46     changeThePhantomSizeCmd -> SetGuidance("Insert sizes X Y and Z"
 47                                            "\n <<  47                                      "\n   0 or negative values mean <<Don't change it!>>");
 48     changeThePhantomSizeCmd -> SetParameterNam <<  48     changeThePhantomSizeCmd -> SetParameterName("PhantomSizeAlongX", 
 49                                                <<  49             "PhantomSizeAlongY", 
 50                                                <<  50             "PhantomSizeAlongZ", false);
 51     changeThePhantomSizeCmd -> SetDefaultUnit(     51     changeThePhantomSizeCmd -> SetDefaultUnit("mm");
 52     changeThePhantomSizeCmd -> SetUnitCandidat <<  52     changeThePhantomSizeCmd -> SetUnitCandidates("nm um mm cm"); 
 53     changeThePhantomSizeCmd -> AvailableForSta     53     changeThePhantomSizeCmd -> AvailableForStates(G4State_Idle);
 54                                                <<  54 
 55                                                <<  55 
 56     // Change Phantom material                 <<  56     // Change Phantom material 
 57     changeThePhantomMaterialCmd = new G4UIcmdW     57     changeThePhantomMaterialCmd = new G4UIcmdWithAString("/changePhantom/material", this);
 58     changeThePhantomMaterialCmd -> SetGuidance <<  58     changeThePhantomMaterialCmd -> SetGuidance("Change the Phantom and the detector material"); 
 59     changeThePhantomMaterialCmd -> SetParamete     59     changeThePhantomMaterialCmd -> SetParameterName("PhantomMaterial", false);
 60     changeThePhantomMaterialCmd -> SetDefaultV     60     changeThePhantomMaterialCmd -> SetDefaultValue("G4_WATER");
 61     changeThePhantomMaterialCmd -> AvailableFo     61     changeThePhantomMaterialCmd -> AvailableForStates(G4State_Idle);
 62                                                <<  62 
 63     // Change Phantom position                     63     // Change Phantom position
 64     changeThePhantomPositionCmd = new G4UIcmdW     64     changeThePhantomPositionCmd = new G4UIcmdWith3VectorAndUnit("/changePhantom/position", this);
 65     changeThePhantomPositionCmd -> SetGuidance     65     changeThePhantomPositionCmd -> SetGuidance("Insert X Y and Z dimensions for the position of the center of the Phantom"
 66                                                <<  66             " respect to that of the \"World\""); 
 67     changeThePhantomPositionCmd -> SetParamete <<  67     changeThePhantomPositionCmd -> SetParameterName("PositionAlongX", 
 68                                                <<  68                 "PositionAlongY", 
 69                                                <<  69                 "PositionAlongZ", false);
 70     changeThePhantomPositionCmd -> SetDefaultU     70     changeThePhantomPositionCmd -> SetDefaultUnit("mm");
 71     changeThePhantomPositionCmd -> SetUnitCand <<  71     changeThePhantomPositionCmd -> SetUnitCandidates("um mm cm m"); 
 72     changeThePhantomPositionCmd -> AvailableFo     72     changeThePhantomPositionCmd -> AvailableForStates(G4State_Idle);
 73                                                <<  73 
 74                                                <<  74 
 75     updateCmd = new G4UIcmdWithoutParameter("/     75     updateCmd = new G4UIcmdWithoutParameter("/changePhantom/update",this);
 76     updateCmd->SetGuidance("Update Phantom/Det     76     updateCmd->SetGuidance("Update Phantom/Detector geometry.");
 77     updateCmd->SetGuidance("This command MUST      77     updateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
 78     updateCmd->SetGuidance("if you changed geo     78     updateCmd->SetGuidance("if you changed geometrical value(s).");
 79     updateCmd->AvailableForStates(G4State_Idle     79     updateCmd->AvailableForStates(G4State_Idle);
 80                                                <<  80 
 81     //  Change detector size                       81     //  Change detector size
 82     changeTheDetectorDir = new G4UIdirectory("     82     changeTheDetectorDir = new G4UIdirectory("/changeDetector/");
 83     changeTheDetectorDir -> SetGuidance("Comma     83     changeTheDetectorDir -> SetGuidance("Command to change the Detector's Size/position/Voxels");
 84                                                    84     
 85     changeTheDetectorSizeCmd = new G4UIcmdWith     85     changeTheDetectorSizeCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/size",this);
 86     changeTheDetectorSizeCmd -> SetGuidance("I     86     changeTheDetectorSizeCmd -> SetGuidance("Insert sizes for X Y and Z dimensions of the Detector"
 87                                             "\ <<  87               "\n   0 or negative values mean <<Don't change it>>");
 88     changeTheDetectorSizeCmd -> SetParameterNa     88     changeTheDetectorSizeCmd -> SetParameterName("DetectorSizeAlongX", "DetectorSizeAlongY", "DetectorSizeAlongZ", false);
 89     changeTheDetectorSizeCmd -> SetDefaultUnit     89     changeTheDetectorSizeCmd -> SetDefaultUnit("mm");
 90     changeTheDetectorSizeCmd -> SetUnitCandida <<  90     changeTheDetectorSizeCmd -> SetUnitCandidates("nm um mm cm"); 
 91     changeTheDetectorSizeCmd -> AvailableForSt     91     changeTheDetectorSizeCmd -> AvailableForStates(G4State_Idle);
 92                                                <<  92 
 93     //  Change the detector to phantom displac     93     //  Change the detector to phantom displacement
 94     changeTheDetectorToPhantomPositionCmd = ne     94     changeTheDetectorToPhantomPositionCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/displacement",this);
 95     changeTheDetectorToPhantomPositionCmd -> S     95     changeTheDetectorToPhantomPositionCmd -> SetGuidance("Insert X Y and Z displacements between Detector and Phantom"
 96                                                <<  96                                                    "\nNegative values mean <<Don't change it!>>"); 
 97     changeTheDetectorToPhantomPositionCmd -> S     97     changeTheDetectorToPhantomPositionCmd -> SetParameterName("DisplacementAlongX",
 98                                                <<  98                     "DisplacementAlongY", 
 99                                                <<  99                     "DisplacementAlongZ", false);
100     changeTheDetectorToPhantomPositionCmd -> S    100     changeTheDetectorToPhantomPositionCmd -> SetDefaultUnit("mm");
101     changeTheDetectorToPhantomPositionCmd -> S << 101     changeTheDetectorToPhantomPositionCmd -> SetUnitCandidates("nm um mm cm"); 
102     changeTheDetectorToPhantomPositionCmd -> A    102     changeTheDetectorToPhantomPositionCmd -> AvailableForStates(G4State_Idle);
103                                                   103     
104     // Change voxels by its size                  104     // Change voxels by its size
105     changeTheDetectorVoxelCmd = new G4UIcmdWit    105     changeTheDetectorVoxelCmd = new G4UIcmdWith3VectorAndUnit("/changeDetector/voxelSize",this);
106     changeTheDetectorVoxelCmd -> SetGuidance("    106     changeTheDetectorVoxelCmd -> SetGuidance("Insert Voxel sizes for X Y and Z dimensions"
107                                              " << 107                                  "\n   0 or negative values mean <<Don't change it!>>");
108     changeTheDetectorVoxelCmd -> SetParameterN    108     changeTheDetectorVoxelCmd -> SetParameterName("VoxelSizeAlongX", "VoxelSizeAlongY", "VoxelSizeAlongZ", false);
109     changeTheDetectorVoxelCmd -> SetDefaultUni    109     changeTheDetectorVoxelCmd -> SetDefaultUnit("mm");
110     changeTheDetectorVoxelCmd -> SetUnitCandid    110     changeTheDetectorVoxelCmd -> SetUnitCandidates("nm um mm cm");
111     changeTheDetectorVoxelCmd -> AvailableForS    111     changeTheDetectorVoxelCmd -> AvailableForStates(G4State_Idle);
                                                   >> 112 
                                                   >> 113 
112                                                   114     
113                                                   115     
114     changeTheSource = new G4UIdirectory("/chan << 116 
115     changeTheSource -> SetGuidance("Command to << 
116                                                << 
117                                                << 
118                                                << 
119                                                << 
120     ////////////////////////////////////////// << 
121                                                << 
122     VirtualLayer = new G4UIcmdWithABool("/chan << 
123     VirtualLayer -> SetParameterName("VirtualL << 
124     VirtualLayer -> SetDefaultValue("false");  << 
125     VirtualLayer -> SetGuidance("Set if you wa << 
126                                 "\n[usage]: /c << 
127     VirtualLayer -> AvailableForStates(G4State << 
128                                                << 
129                                                << 
130     VirtualLayerPosition = new G4UIcmdWith3Vec << 
131     VirtualLayerPosition -> SetGuidance("Inser << 
132                                         " resp << 
133     VirtualLayerPosition -> SetParameterName(" << 
134                                              " << 
135                                              " << 
136                                                << 
137     VirtualLayerPosition -> SetDefaultUnit("cm << 
138     VirtualLayerPosition -> SetUnitCandidates( << 
139     VirtualLayerPosition -> AvailableForStates << 
140                                                << 
141                                                << 
142                                                << 
143 }                                                 117 }
144                                                   118 
145 //////////////////////////////////////////////    119 /////////////////////////////////////////////////////////////////////////////
146 HadrontherapyDetectorMessenger::~Hadrontherapy    120 HadrontherapyDetectorMessenger::~HadrontherapyDetectorMessenger()
147 {                                                 121 {
148     delete changeThePhantomDir;                << 122     delete changeThePhantomDir; 
149     delete changeThePhantomSizeCmd;            << 123     delete changeThePhantomSizeCmd; 
150     delete changeThePhantomPositionCmd;        << 124     delete changeThePhantomPositionCmd; 
151     delete changeThePhantomMaterialCmd;        << 125     delete changeThePhantomMaterialCmd; 
152     delete updateCmd;                             126     delete updateCmd;
153     delete changeTheDetectorDir;               << 127     delete changeTheDetectorDir; 
154     delete changeTheDetectorSizeCmd;           << 128     delete changeTheDetectorSizeCmd; 
155     delete changeTheDetectorToPhantomPositionC << 129     delete changeTheDetectorToPhantomPositionCmd; 
156     delete changeTheDetectorVoxelCmd;             130     delete changeTheDetectorVoxelCmd;
157     delete VirtualLayer;                       << 131 
158     delete VirtualLayerPosition;               << 
159 }                                                 132 }
160                                                   133 
161 //////////////////////////////////////////////    134 /////////////////////////////////////////////////////////////////////////////
162 void HadrontherapyDetectorMessenger::SetNewVal    135 void HadrontherapyDetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
163 {                                                 136 {
164                                                << 137   
165     if( command == changeThePhantomSizeCmd)    << 138   if( command == changeThePhantomSizeCmd)
166     {                                          << 
167         G4ThreeVector size = changeThePhantomS << 
168         hadrontherapyDetector -> SetPhantomSiz << 
169     }                                          << 
170     else if (command == changeThePhantomPositi << 
171     {                                          << 
172                                                << 
173         G4ThreeVector size = changeThePhantomP << 
174         hadrontherapyDetector -> SetPhantomPos << 
175     }                                          << 
176     else if (command == changeThePhantomMateri << 
177     {                                             139     {
178         hadrontherapyDetector -> SetPhantomMat << 140   G4ThreeVector size = changeThePhantomSizeCmd -> GetNew3VectorValue(newValue);
                                                   >> 141   hadrontherapyDetector -> SetPhantomSize(size.getX(),size.getY(),size.getZ());
179     }                                             142     }
180     else if (command == changeTheDetectorSizeC << 143   else if (command == changeThePhantomPositionCmd )
181     {                                          << 144   {
182         G4ThreeVector size = changeTheDetector << 145    G4ThreeVector size = changeThePhantomPositionCmd -> GetNew3VectorValue(newValue);
                                                   >> 146          hadrontherapyDetector -> SetPhantomPosition(size);
                                                   >> 147   }
                                                   >> 148   else if (command == changeThePhantomMaterialCmd)
                                                   >> 149   {
                                                   >> 150       hadrontherapyDetector -> SetPhantomMaterial(newValue);
                                                   >> 151   }
                                                   >> 152   else if (command == changeTheDetectorSizeCmd)
                                                   >> 153   {
                                                   >> 154   G4ThreeVector size = changeTheDetectorSizeCmd  -> GetNew3VectorValue(newValue);
183         hadrontherapyDetector -> SetDetectorSi    155         hadrontherapyDetector -> SetDetectorSize(size.getX(),size.getY(),size.getZ());
184     }                                          << 156   }
185     else if (command == changeTheDetectorToPha << 157   else if (command == changeTheDetectorToPhantomPositionCmd)
186     {                                          << 158   {
187         G4ThreeVector size = changeTheDetector << 159   G4ThreeVector size = changeTheDetectorToPhantomPositionCmd-> GetNew3VectorValue(newValue);
188         hadrontherapyDetector -> SetDetectorTo    160         hadrontherapyDetector -> SetDetectorToPhantomPosition(size);
189                                                << 161   }
190     }                                          << 162   else if (command == changeTheDetectorVoxelCmd)
191     else if (command == changeTheDetectorVoxel << 163   {
192     {                                          << 164   G4ThreeVector size = changeTheDetectorVoxelCmd  -> GetNew3VectorValue(newValue);
193         G4ThreeVector size = changeTheDetector << 
194         hadrontherapyDetector -> SetVoxelSize(    165         hadrontherapyDetector -> SetVoxelSize(size.getX(),size.getY(),size.getZ());
195     }                                          << 166   }
196     else if (command == updateCmd)             << 167   else if (command == updateCmd)
197     {                                          << 168   {
198         hadrontherapyDetector -> UpdateGeometr << 169       hadrontherapyDetector -> UpdateGeometry();
199     }                                          << 170   }
200                                                << 171  
201     else if(command == VirtualLayer)           << 
202     {                                          << 
203         hadrontherapyDetector->VirtualLayer(Vi << 
204     }                                          << 
205     else if(command == VirtualLayerPosition)   << 
206     {                                          << 
207         G4ThreeVector size = VirtualLayerPosit << 
208                                                << 
209         hadrontherapyDetector -> SetVirtualLay << 
210     }                                          << 
211                                                << 
212 }                                                 172 }
213                                                   173