Computed tomography
From UCL HEP PBT Wiki
m |
m |
||
| Line 58: | Line 58: | ||
</pre> | </pre> | ||
| - | ; | + | ; Run macro run.mac |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| Line 68: | Line 64: | ||
</pre> | </pre> | ||
| - | * | + | == <span style="color:#000080"> How to analyze data </span> == |
| + | |||
| + | This is a [http://geant4.web.cern.ch/geant4/UserDocumentation/Doxygen/examples_doc/html/ExampleDICOM.html link] to the official DICOM explanation notes. We recommend to read it before proceeding with this tutorial. | ||
| + | |||
| + | === <span style="color:#000080"> Text files </span> === | ||
| + | |||
| + | The macro run.mac produces a text file '''dicom.out''' which includes the dose deposition in several voxels. | ||
| + | |||
| + | === <span style="color:#000080"> Run with different settings </span> === | ||
| + | |||
| + | ; Chose among several voxelization algorithms | ||
| + | |||
| + | * <span style="color:#ff0000"> G4RegularNavigation </span>: | ||
| + | |||
| + | '''G4RegularNavigation''' class is the default class for this tutorial. This algorithm skips frontiers between voxels when they have the same material i.e. "replacing group of voxels with a smaller number of larger voxels". | ||
| + | |||
| + | * <span style="color:#ff0000"> G4NestedParameterization </span>: | ||
| + | |||
| + | To run with this voxelization algorithm you need to set the variable '''DICOM_NESTED_PARAM''' | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| Line 74: | Line 88: | ||
</pre> | </pre> | ||
| - | You can check if | + | You can check if it is set correctly by typing in the command line |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| Line 80: | Line 94: | ||
</pre> | </pre> | ||
| - | Then, run the code: | + | Then, compile and run the code: |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| + | [username@plus1 DICOM_build]$ make | ||
| + | |||
[username@plus1 DICOM_build]$ ./DICOM run.mac | [username@plus1 DICOM_build]$ ./DICOM run.mac | ||
</pre> | </pre> | ||
| - | * | + | * <span style="color:#ff0000"> G4SmartVoxel / G4VoxelNavigation </span>: |
Using smart voxels required a huge amount of memory. At /home/username/DICOMFolder/DICOM_source/src/DicomRegularDetectorConstruction.cc set | Using smart voxels required a huge amount of memory. At /home/username/DICOMFolder/DICOM_source/src/DicomRegularDetectorConstruction.cc set | ||
| Line 102: | Line 118: | ||
</pre> | </pre> | ||
| - | + | ; Change the type, energy and position of the incident particle | |
| - | + | At /home/username/DICOMFolder/DICOM_source/src/DicomPrimaryGeneratorAction.cc change the following lines: | |
| - | + | <pre style="color: #800000; background-color: #dcdcdc"> | |
| + | G4ParticleDefinition* particle = particleTable->FindParticle(particleName="e-"); | ||
| - | + | fParticleGun->SetParticleEnergy(100.*MeV); | |
| + | |||
| + | fParticleGun->SetParticlePosition(G4ThreeVector(0.,0.,0.)); | ||
| + | </pre> | ||
| + | |||
| + | You can chose among several particles like "gamma", "e+", "alpha" and "He3". Then, compile and run the code | ||
| + | |||
| + | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| + | [username@plus1 DICOM_build]$ make | ||
| + | |||
| + | [username@plus1 DICOM_build]$ ./DICOM run.mac | ||
| + | </pre> | ||
=== <span style="color:#000080"> Visualisation </span> === | === <span style="color:#000080"> Visualisation </span> === | ||
