GALPROP https://galprop.stanford.edu/forum/ 

Proton Norm Type 2,3 not working https://galprop.stanford.edu/forum/viewtopic.php?f=7&t=187 
Page 1 of 1 
Author:  erccarls [ Mon Apr 14, 2014 7:00 pm ] 
Post subject:  Proton Norm Type 2,3 not working 
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 = 1e8 5.00e9 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 
Author:  strong [ Mon Apr 14, 2014 11:24 pm ] 
Post subject:  Re: Proton Norm Type 2,3 not working 
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 cm2 sr1 s2 MeV1 (not cm3 sr1 s1 MeV1). 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 s1. 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. 
Author:  erccarls [ Tue Apr 15, 2014 2:07 pm ] 
Post subject:  Re: Proton Norm Type 2,3 not working 
Perfect! Thank you. 
Author:  strong [ Wed Apr 16, 2014 12:08 am ] 
Post subject:  Re: Proton Norm Type 2,3 not working 
this was useful since it enabled me to correct an error in the documentation for future releases! see the latest r2439. 
Page 1 of 1  All times are UTC  8 hours [ DST ] 
Powered by phpBB® Forum Software © phpBB Group https://www.phpbb.com/ 