JosephBateman: /* Monoenergetic Proton Pencil Beam */

2024-11-06T12:27:42Z

Monoenergetic Proton Pencil Beam

← Older revision		Revision as of 12:27, 6 November 2024
Line 16:		Line 16:

	The image shows the water box divided into slices using class <code>G4PVReplica</code>.		The image shows the water box divided into slices using class <code>G4PVReplica</code>.
	Protons are in blue, photons are in green.		Protons are in blue, photons are in green.

	== Running the Tutorial ==		== Running the Tutorial ==

SimonJolly at 17:49, 5 October 2016

2016-10-05T17:49:07Z

← Older revision		Revision as of 17:49, 5 October 2016
Line 39:		Line 39:

	<pre>		<pre>
	[username@plus1 ProtonPBFolder]$ source /unix/pbt/software/scripts/pbt.sh		[username@plus1 ProtonPBFolder]$ source /unix/pbt/software/scripts/pbt.sh
	</pre>		</pre>

SimonJolly at 14:38, 12 May 2016

2016-05-12T14:38:33Z

← Older revision		Revision as of 14:38, 12 May 2016
Line 478:		Line 478:
	.q		.q
	</pre>		</pre>

			=== Plotting the dose distribution ===

	Folder <code>RootScripts</code> contains several ROOT scripts which plot dose deposition in data and simulation.		Folder <code>RootScripts</code> contains several ROOT scripts which plot dose deposition in data and simulation.

SimonJolly at 18:00, 1 April 2016

2016-04-01T18:00:20Z

← Older revision		Revision as of 18:00, 1 April 2016
Line 1:		Line 1:
	== ~~Introduction~~ ==		== Monoenergetic Proton Pencil Beam ==

	This example shows the dose distribution in water along the incident proton beam.		This example shows the dose distribution in water along the incident proton beam.

SimonJolly at 18:51, 31 March 2016

2016-03-31T18:51:27Z

← Older revision		Revision as of 18:51, 31 March 2016
Line 837:		Line 837:
	== Files ==		== Files ==

	[[List of ~~monoenergetic proton pencil beam~~ files ~~with brief description~~]]		[[Software/Geant4/Tutorials/Basic/Proton Pencil Beam Files\|List of files for Mononenergetic Proton Pencil Beam simulation]].

SimonJolly at 18:49, 31 March 2016

2016-03-31T18:49:26Z

← Older revision		Revision as of 18:49, 31 March 2016
Line 39:		Line 39:

	<pre>		<pre>
	[username@plus1 ProtonPBFolder]$ source /unix/pbt/software/~~dev/bin~~/pbt~~-dev~~.sh		[username@plus1 ProtonPBFolder]$ source /unix/pbt/software/scripts/pbt.sh
	</pre>		</pre>

Line 45:		Line 45:

	<pre>		<pre>
	[username@plus1 ProtonPBFolder]$ cp -r /unix/pbt/tutorials/basic/ProtonPB .		[username@plus1 ProtonPBFolder]$ cp -rv /unix/pbt/tutorials/basic/ProtonPB .

	[username@plus1 ProtonPBFolder]$ mv ProtonPB ProtonPB_source		[username@plus1 ProtonPBFolder]$ mv ProtonPB ProtonPB_source

SimonJolly at 18:48, 31 March 2016

2016-03-31T18:48:25Z

← Older revision		Revision as of 18:48, 31 March 2016
Line 18:		Line 18:
	Protons are in blue, photons are in green.		Protons are in blue, photons are in green.

	== Running the Tutorial ==		== Running the Tutorial ==

			All the source files for the tutorial can be found in http://www.hep.ucl.ac.uk/pbt/wikiData/code/ProtonPB/.

			The corresponding data files are available from http://www.hep.ucl.ac.uk/pbt/wikiData/data/ProtonPB/.

	=== Create folder ProtonPBFolder ===		=== Create folder ProtonPBFolder ===

SimonJolly at 18:45, 31 March 2016

2016-03-31T18:45:55Z

← Older revision		Revision as of 18:45, 31 March 2016
Line 463:		Line 463:
	==== The energy deposition along the beam in the absorber ====		==== The energy deposition along the beam in the absorber ====

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/BraggPeak_PB1~~.png		[[File:ProtonPB_BraggPeak_PB1.png]]

	==== The energy deposition along the beam zoomed around the peak ====		==== The energy deposition along the beam zoomed around the peak ====

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/BraggPeak_PBzoom1~~.png		[[File:ProtonPB_BraggPeak_PBzoom1.png]]

	You can close your ROOT session by typing		You can close your ROOT session by typing
Line 494:		Line 494:
	This will create <code>Simulation.root</code> file with the following plot:		This will create <code>Simulation.root</code> file with the following plot:

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/Simulation1~~.png		[[File:ProtonPB_Simulation1.png]]

	You can also plot the file <code>PlotDose.txt</code> using MATLAB.		You can also plot the file <code>PlotDose.txt</code> using MATLAB.
Line 513:		Line 513:
	This plot will be created with added axis labels and a legend:		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		[[File:ProtonPB_matlab1proton.png]]

	You can also plot the data in <code>DoseLongitudinalMesh.txt</code> and <code>DoseLateralMesh.txt</code> which were created using commands <code>/score/</code> in the macro <code>proton.mac</code>.		You can also plot the data in <code>DoseLongitudinalMesh.txt</code> and <code>DoseLateralMesh.txt</code> which were created using commands <code>/score/</code> in the macro <code>proton.mac</code>.
Line 532:		Line 532:
	This will create <code>LateralDose_Mesh.root</code> file with the following plot:		This will create <code>LateralDose_Mesh.root</code> file with the following plot:

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/SimulationLateralMesh~~.png		[[File:ProtonPB_SimulationLateralMesh.png]]

	=== Run with different settings ===		=== Run with different settings ===
Line 713:		Line 713:
	</pre>		</pre>

	If you modify files in directory ProtonPB_source (for example you change the value of MaxLayer) you need to ~~comppile~~ the code.		If you modify files in directory ProtonPB_source (for example you change the value of MaxLayer) you need to compile the code.
	In directory PhotonPB_build do:		In directory PhotonPB_build do:

Line 744:		Line 744:
	This will create your first visualisation:		This will create your first visualisation:

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_00_6000e.eps~~		[[File:ProtonPB_g4_00_6000e.gif\|500px]]

	The second file g4_01.prim will not open automatically.		The second file g4_01.prim will not open automatically.
Line 756:		Line 756:
	In the opened window press OK and this will create the image:		In the opened window press OK and this will create the image:

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_02_6000e.eps~~		[[File:ProtonPB_g4_02_6000e.gif\|500px]]

	You can modify the .prim files in DAWN.		You can modify the .prim files in DAWN.
Line 762:		Line 762:
	1) (polar angle, azimuthal angle) = (0,90) will create this image:		1) (polar angle, azimuthal angle) = (0,90) will create this image:

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_01_6000e.eps~~		[[File:ProtonPB_g4_01_6000e.gif\|500px]]

	2) (polar angle, azimuthal angle) = (90,0) will create this image:		2) (polar angle, azimuthal angle) = (90,0) will create this image:

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_04_6000e.eps~~		[[File:ProtonPB_g4_04_6000e.gif\|500px]]

	Now, in <code>proton.mac</code> macro uncomment <code>/score/drawProjection waterMeshlateral doseDeposit</code> and comment <code>/score/drawProjection waterMeshlongitudinal doseDeposit</code>.		Now, in <code>proton.mac</code> macro uncomment <code>/score/drawProjection waterMeshlateral doseDeposit</code> and comment <code>/score/drawProjection waterMeshlongitudinal doseDeposit</code>.
Line 772:		Line 772:
	The image with the lateral dose projections will look like that:		The image with the lateral dose projections will look like that:

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_03_6000e.eps~~		[[File:ProtonPB_g4_03_6000e.gif\|500px]]

	== Data from The Clatterbridge Cancer Centre ==		== Data from The Clatterbridge Cancer Centre ==

	~~This is~~ [http://www.hep.ucl.ac.uk/pbt/~~RadiotherapyWorkbook~~/~~skins/common/images~~/Clatterbridge/ClatterbridgeBraggPeak.txt data] from the Clatterbridge Cancer Center.		The file [http://www.hep.ucl.ac.uk/pbt/wikiData/data/Clatterbridge/ClatterbridgeBraggPeak.txt <code>ClatterbridgeBraggPeak.txt</code>] contains Bragg Peak data from the Clatterbridge Cancer Center.
	The first column represents the values in mm in depth, the second column represents the normalized values of deposited dose and normalized at the peak.		The first column represents the values in mm in depth, the second column represents the normalized values of deposited dose and normalized at the peak.
	You can use script <code>PlotData.C</code> from folder <code>RootScripts</code> to plot the data.		You can use script <code>PlotData.C</code> from folder <code>RootScripts</code> to plot the data: these are available in the directory http://www.hep.ucl.ac.uk/pbt/wikiData/code/ProtonPB/.
	Copy files <code>ClatterbridgeData.txt</code> and <code>PlotData.C</code> to your current ProtonPB_build directory and run the script:		Copy files <code>ClatterbridgeData.txt</code> and <code>PlotData.C</code> to your current ProtonPB_build directory and run the script:

Line 793:		Line 793:
	This will create <code>ClatterbridgeData.root</code> file with the following plot:		This will create <code>ClatterbridgeData.root</code> file with the following plot:

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/Clatterbridge/Clatterbridge~~.png		[[File:ClatterbridgeBraggPeak.png]]

	== Comparison between data and simulation ==		== Comparison between data and simulation ==

	The scripts <code>PlotDataAndSim.C</code> and <code>PlotDataAndSimMesh.C</code> in folder <code>RootScripts</code> compare data with simulation.		The scripts <code>PlotDataAndSim.C</code> and <code>PlotDataAndSimMesh.C</code> in folder <code>RootScripts</code> compare data with simulation.
	You can use script <code>PlotDataAndSim.C</code> to compare data (ClatterbridgeData.txt) and simulation (<code>PlotDose.txt</code>).		You can use script <code>PlotDataAndSim.C</code> to compare data (<code>ClatterbridgeData.txt</code>) and simulation (<code>PlotDose.txt</code>).
	Both text files must be in the folder where you run the script.		Both text files must be in the folder where you run the script.

Line 811:		Line 811:
	This will create <code>BraggPeakComparison.root</code> file with the following plot:		This will create <code>BraggPeakComparison.root</code> file with the following plot:

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/DataSimulation~~.png		[[File:ProtonPB_DataSimulation.png]]

	You can also compare data (ClatterbridgeData.txt) with simulation done with scoring mesh (<code>DoseLongitudinalMesh.txt</code>).		You can also compare data (ClatterbridgeData.txt) with simulation done with scoring mesh (<code>DoseLongitudinalMesh.txt</code>).
Line 829:		Line 829:
	This creates <code>BraggPeakComparison_Mesh.root</code> file with the following plot:		This creates <code>BraggPeakComparison_Mesh.root</code> file with the following plot:

	~~http~~:~~//www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/DataSimulationMesh~~.png		[[File:ProtonPB_DataSimulationMesh.png]]

	== Files ==		== Files ==

	[[List of monoenergetic proton pencil beam files with brief description]]		[[List of monoenergetic proton pencil beam files with brief description]]

SimonJolly at 18:25, 31 March 2016

2016-03-31T18:25:31Z

Show changes

SimonJolly at 18:17, 31 March 2016

2016-03-31T18:17:30Z

← Older revision		Revision as of 18:17, 31 March 2016
Line 13:		Line 13:
	There is an option to chose among several <code>EM</code> and the <code>QGSP_BIC_EMY</code> physics lists.		There is an option to chose among several <code>EM</code> and the <code>QGSP_BIC_EMY</code> physics lists.

	~~http~~:~~//www~~.~~hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/g4_00_6000e.png~~		[[File:ProtonPB_g4_00_6000e.gif\|500px]]

	The image shows the water box divided into slices using class <code>G4PVReplica</code>.		The image shows the water box divided into slices using class <code>G4PVReplica</code>.

Software/Geant4/Tutorials/Basic/Monoenergetic Proton Pencil Beam - Revision history