Table of Contents
Find out aluminium thickness from MSL RAD model, see paper J. Appel 2017
Mars Paper 2
Simulation preparation
* Lennart has prepared a PSF file on 12th of October
* Sasha has confirmed that the PSF file is OK
* Sasha has converted the PSF file to GDML
* Sasha has checked that the GDML file is OK
* Sasha has prepared extended GCR phantoms for ICRU and msl/rad detector B
* Sasha has created a 2mars directory in his phd git repository
* Based on the discussion with Jingnan, we are going to use the following energy range:
* Sasha has prepared a simulation script which takes: Z, number_of_threads, portion_factor and performs the needed run.
* We need for a single run the binary output! Also we need to benchmark if this is quicker or not.
* Based on the previous paper, we are going to use INCLXX for our calculations.
* INCLXX seams to have a bug when dealing with highest Z ions.
Simulation status
| Z | element name | status | machine | ETA |
|---|---|---|---|---|
| 1 | H | malanosa | 14.12. | |
| 2 | He | malanosa | 16.12. | |
| 3 | Li | optimus | 16.12. | |
| 4 | Be | tr1 | 17.12. | |
| 5 | B | tr1 | 18.12. | |
| 6 | C | tr1 | 18.12. | |
| 7 | N | tr1 | 19.12. | |
| 8 | O | optimus | 20.12. | |
| 9 | F | malanosa | 20.12. | |
| 10 | Ne | malanosa | 20.12. | |
| 11 | Natrium | terminator | 21.12. | |
| 12 | Magnesium | optimus | 23.12. | |
| 13 | Al | tr1 | 24.12. | |
| 14 | Si | tr1 | 25.12. | |
| 15 | P | tr1 | 26.12. | |
| 16 | S | tr1 | 26.12. | |
| 17 | Cl | tr2 | 27.12. | |
| 18 | Ar | tr2 | 31.12. | |
| 19 | K (potasium) | tr2 | 30.12. | |
| 20 | calcium | malanosa | 27.12. | |
| 21 | scandium | malanosa | 29.12. | |
| 22 | titanium | malanosa | 31.12. | |
| 23 | vanadium | optimus | 4.1. | |
| 24 | chromium | optimus | 4.1. | |
| 25 | manganese | optimus | 5.1. | |
| 26 | Ferum (Iron) | terminator | 1.1. | |
| 27 | cobalt | terminator | 2.1. | |
| 28 | nickel | terminator | 5.1. |
Mars SEP calculator paper
- 20 Mars atmospheres with approx 15Pa steps in pressure
- Only primary protons
- energy range between 0.1 MeV and 10 GeV, i.e. -1 to 4 decades
- 20 bins per decade ( total 100 bins)
- 5000 parts per decade (500.000 particles per run)
- 5 runs per atmosphere
- in total 100 runs
- Calculate results for all the GLEs
Goal: make a function that takes (pressure, diffspectrum) and calculates dose altitude profile for radb and for an icrusphere.
In paper: Talk about the atmospheres, talk about the convolution method, talk about the outlook: determining SEP spectra around mars