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

Update README.md

parent 9c1c3d60
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.
## How to use the GDL codes
The data analysis of the output files produced by the PTra simulations was performed using the GNU data language environment which can be generally directly installed in most Linux distributions.
1. Download all *.pro files into a directory in your default GDL path (encoded in the GDL constant !PATH) or into the directory mentioned in the next point. <br />
2. Open a terminal and navigate to the directory containing as subfolders the directories containing the data produced by the PTra simulations for the different cases (thickness of PMMA absorber and primary carbon ion energy). <br />
3. Start GDL and in the GDL prompt execute the routines Ptra4MS and Ptra4ICSD. <br />
GDL> Ptra4MS <br />
GDL> Ptra4ICSD <br />
## Contact
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment