Update eloss authored by Jeppe Schultz Nielsen's avatar Jeppe Schultz Nielsen
......@@ -30,7 +30,7 @@ int main(void) {
Any questions ? Then read on...
# Overview
ausa/eloss/ allows you to calculate energy losses for practically any ion-target combination using several different tabulations of stopping powers and ranges. The basic structure of the program is shown the figure below. Tabulations of stopping power and range may be generated for any ion-target combination using one of three sources: SRIM, GEANT or ICRU. The tabulations are fed to the `StoppingPowerCalculator` and `RangeCalculator` which calculate the stopping power and range for an arbitrary energy, typically, by cubic spline interpolation. Finally, the `EnergLossCalculator` calculates the energy loss for an arbitrary energy (initial or final) and thickness, either by integrating the stopping powers or by inverting the range function. For complete documentation of the code see http://docs.kern.phys.au.dk/.
ausa/eloss/ allows you to calculate energy losses for practically any ion-target combination using several different tabulations of stopping powers and ranges. The basic structure of the program is shown the figure below. Tabulations of stopping power and range may be generated for any ion-target combination using one of three sources: SRIM, GEANT or ICRU. The tabulations are fed to the `StoppingPowerCalculator` and `RangeCalculator` which calculate the stopping power and range for an arbitrary energy, typically, by cubic spline interpolation. Finally, the `EnergLossCalculator` calculates the energy loss for an arbitrary energy (initial or final) and thickness, either by integrating the stopping powers or by inverting the range function.
![eloss](uploads/c433f98b106be1898bc88267aac406f5/eloss.jpg)
# Tutorial
......@@ -113,7 +113,7 @@ double de = DEsp.getTotalEnergyLoss( 3000., 100E-6 );
```
The result should be close to 100 keV, the exact value depending on the tabulation you are using.
More examples can be found [here](uploads/66bc7fb78ed597a789249a0e2870cc03/snippet7.cpp).
For a complete list of available get methods see [EnergyLossCalculator.h](http://docs.kern.phys.au.dk/classAUSA_1_1EnergyLoss_1_1EnergyLossCalculator.html).
For a complete list of available get methods see [EnergyLossCalculator.h](https://gitlab.au.dk/ausa/ausalib/-/blob/master/include/ausa/eloss/EnergyLossCalculator.h).
## 6. Default calculator
If you are in a rush and do not care too much about how the energy loss is being determined, you can skip steps 2-4 and instead use the [default energy-loss calculator](blob/master/include/ausa/eloss/Default.h) which determines energy loss by numerical integration of GEANT stopping powers. To determine, say, the energy loss of 1 MeV protons in 100 nm aluminium, you would write
......
......