Monoenergetic proton pencil beam
From UCL HEP PBT Wiki
m |
|||
| Line 7: | Line 7: | ||
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_00_6000e.png | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_00_6000e.png | ||
| - | The image shows the water box divided into slices using class '''G4PVReplica'''. | + | The image shows the water box divided into slices using class '''G4PVReplica'''. Protons are in blue, photons are in green. |
== <span style="color:#000080"> How to run the tutorial </span> == | == <span style="color:#000080"> How to run the tutorial </span> == | ||
| Line 53: | Line 53: | ||
</pre> | </pre> | ||
| - | ; Run macro proton.mac | + | ; Run macro proton.mac |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| Line 244: | Line 244: | ||
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/Simulation1.png | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/Simulation1.png | ||
| + | |||
| + | You can also plot the file '''PlotDose.txt''' using MATLAB. Similarly to the previous example first copy the text file to your computer. In the terminal at your computer write: | ||
| + | |||
| + | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| + | scp username@plus1.hep.ucl.ac.uk:/home/username/ProtonPBFolder/ProtonPB_build/PlotDose.txt . | ||
| + | </pre> | ||
| + | |||
| + | Then, open MATLAB and follow the procedure: | ||
| + | |||
| + | * Import the file: Chose 'HOME' tab and 'Import Data'. | ||
| + | * In the 'Import Data' window select the 'PlotDose.txt' file choosing the right path. | ||
| + | * In the opened window select the data points in the 'IMPORT' tab. If you like, you can change the name of the variables. For example, 'x' instead of 'VarName1' and 'Dose' instead of 'VarName2'. Then, press 'Import Selection'/'Import Data'. | ||
| + | * Close the Import Window and in the Command Window type plot(x,Dose). Press Enter. | ||
| + | |||
| + | This plot will be created with added axis labels and a legend: | ||
| + | |||
| + | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/matlab1proton.png | ||
You can also plot the data in '''DoseLongitudinalMesh.txt''' and '''DoseLateralMesh.txt''' which were created using commands '''/score/''' in the macro '''proton.mac'''. The file '''DoseLongitudinalMesh.txt''' will be used later to compare with data from the Clatterbridge Cancer Center. Now, use script '''PlotLateralDoseMesh.C''' to plot the lateral dose distribution. | You can also plot the data in '''DoseLongitudinalMesh.txt''' and '''DoseLateralMesh.txt''' which were created using commands '''/score/''' in the macro '''proton.mac'''. The file '''DoseLongitudinalMesh.txt''' will be used later to compare with data from the Clatterbridge Cancer Center. Now, use script '''PlotLateralDoseMesh.C''' to plot the lateral dose distribution. | ||
| Line 260: | Line 277: | ||
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/SimulationLateralMesh.png | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/SimulationLateralMesh.png | ||
| - | === <span style="color:#000080"> | + | === <span style="color:#000080"> Run with different settings </span> === |
You can change the physics process, incident proton energy and number of slices etc. by | You can change the physics process, incident proton energy and number of slices etc. by | ||
