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