View unanswered posts | View active topics
|
Page 1 of 1
|
[ 4 posts ] |
|
Author |
Message |
erccarls
|
Posted: Mon Apr 14, 2014 7:00 pm |
|
Joined: Thu Apr 05, 2012 10:40 pm Posts: 14
|
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
|
|
Top |
|
|
strong
|
Posted: Mon Apr 14, 2014 11:24 pm |
|
Joined: Wed May 17, 2006 7:02 am Posts: 285 Location: MPE Garching
|
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
|
|
Top |
|
|
erccarls
|
Posted: Tue Apr 15, 2014 2:07 pm |
|
Joined: Thu Apr 05, 2012 10:40 pm Posts: 14
|
|
Top |
|
|
strong
|
Posted: Wed Apr 16, 2014 12:08 am |
|
Joined: Wed May 17, 2006 7:02 am Posts: 285 Location: MPE Garching
|
this was useful since it enabled me to correct an error in the documentation for future releases! see the latest r2439.
_________________ Andy Strong, MPE
|
|
Top |
|
|
|
Page 1 of 1
|
[ 4 posts ] |
|
Who is online |
Users browsing this forum: No registered users and 2 guests |
|
You cannot post new topics in this forum You cannot reply to topics in this forum You cannot edit your posts in this forum You cannot delete your posts in this forum You cannot post attachments in this forum
|