Skip to content
Snippets Groups Projects
Commit 9c1c3d60 authored by Hans Rabus's avatar Hans Rabus
Browse files

Upload New File

parent 8f5dde7a
No related branches found
No related tags found
1 merge request!1Add new directory
## How to get the code to run
### Step 1: Download and compile the codes
1. Download all *.f files into the same directory. <br />
2. Open a terminal and navigate to this directory. <br />
3. Compile the codes FilterExtrAp.f and PTra_c3h8_20170608HIT.f using a Fortran 90/95 compatible compiler. <br />
Example using gfortran and producing an executable PTra_HIT: <br />
$ *gfortran **PTra_c3h8_20170608HIT.f** > ~/bin/PTra_HIT* <br />
Example using gfortran and producing an executable PTra_HIT: <br />
$ *gfortran **FilterExtrAp.f** > ~/bin/FilterExtrAp*
### Step 2: Filter events containing charged particles from the Geant4 output files
Run the executable produced from FilterExtrAp.f best with a steering text file. <br />
Example with the executable produced above and the sample steering file **cox.dat**: <br />
$ *~/bin/FilterExtrAp < cox.dat* <br />
The first two lines are the numbers of files per case (PMMA thickness and primary energy) and of subdirectories (each corresponding to a combination of PMMA thickness and primary energy). <br />
The third line is a flag that should be set to 0 if the results from the Geant4 simulations have been produced in different runs and are merged during the filtering. (In this case files named F_ExtractionScoringPlane_t*.dat resdiding in the output directories must be deleted before running FilterExtrAp.) <br />
The remaining even lines are the paths to the Geant4 simulation output files and the ensuing odd lines are the corresponding paths where the output files are written.
### Step 3: Run the simulations
1. Download the data file **nanodosimeter3_eff_HIT_1.2mbar-C3H8_drift_V4c_part-1-2_ds2v4508.dat** which contains the detection efficiency map. <br />
2. Create a text file containing the info on the simulation parameters. <br />
An example is provided here as the file **HITEIN_240513.txt**. It may be necessary to change the last line of the section "DETECTION EFFICIENCY" containing the path to the efficiency map.
3. Create a steering text file (to avoid manual input). <br />
An example is provided here as the file **bco.dat**. (It will be necessary to change the first line to reflect the correct path to the text file containing the info on the simulation parameters.)
4. Open a terminal and navigate to the directory where the data files (produced by FilterExrAp) are situated. <br />
5. Start the simulation.<br />
Example with the steering file bco.dat: <br />
$ *~/bin/PTra_HIT < bco.dat* <br />
The output files of the simulation are written to the directory in which the filtered datafiles reside.
## Contact
Hans Rabus - hans.rabus@ptb.de <br />
## License
GNU GENERAL PUBLIC LICENSE. See `LICENSE.txt` for more information.
## Citation
If you use the code for your research, please consider citing:
Miriam Schwarze, Gerhard Hilgers, and Hans Rabus. "Nanodosimetric investigation of the track structure of therapeutic carbon ion radiation part 2: Detailed simulation." 10.1088/2057-1976/ad9152.<br />
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment