I am simulating a single Gaussian proton point source and I would like to normalize the flux to the total injected proton luminosity (i.e. integrated over the entire galaxy.) In the explanatory supplement this is supposed to be possible using proton_norm_type=2 or 3 (where 1 is the standard normalization at the solar position.) When I look at the resulting pion_decay_mapcube I get the same results using proton_norm_type=1,2 or 3. The result is proportional to proton_norm_flux.

This section of the galdef looks like this:

proton_norm_type= 1 or 2 or 3
proton_norm_Ekin = 1.00e+3 proton kinetic energy for normalization (MeV)
proton_norm_flux = 1e-8 5.00e-9 to renorm nuclei/flux of protons at norm energy (cm^-2 sr^-1 s^-1 MeV^-1)

electron_norm_flux = 0 flux of electrons at normalization energy (cm^-2 sr^-1 s^-1 MeV^-1)
source_norm = 1.0 absolute normalization for proton CR source function (only if electron_norm_flux=proton_norm_flux=0)
electron_source_norm = .0 absolute normalization for electron CR source function (only if electron_norm_flux=proton_norm_flux=0)


Any ideas? I am using galprop r2423.

-Eric

Good point!
As written in the Explanatory Supplement, the parameters proton_norm_type and electron_norm_type were added at r603 but not used. It was never in fact implemented. It would be quite easy to do, thanks for the reminder.

Meanwhile it can be done using current GALPROP output by just integrating the proton injection over the galaxy to get the required scaling factor.
The source functions are output : files primary_protons_source_function_ , primary_electrons_source_function_ etc

NB for the units of those files - what is written in the Explanatory Supplement should be cm-2 sr-1 s-2 MeV-1 (not cm-3 sr-1 s-1 MeV-1).
since there is a c/4pi factor for the flux units - see description of GALPROP units in Sec 2.2, and the units of the source term q in equation (1) must match those, just with extra s-1. I will correct it in the Explanatory Supplement. For details see also the routine cr_luminosity.cc.


Added later: in fact it is easy since GALPROP prints out the total CR luminosity of each species during the run. Example : galdef_54_reltest01
=====================================================
CR proton luminosity=3.06412e+41 erg s^-1
=====================================================
=====================================================
CR primary electron luminosity=1.38217e+39 erg s^-1
=====================================================

so you can just use this to normalize to your required value.

_________________
Andy Strong, MPE

Perfect! Thank you.

this was useful since it enabled me to correct an error in the documentation for future releases!
see the latest r2439.

_________________
Andy Strong, MPE