It is currently Mon Oct 07, 2024 7:59 pm

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 4 posts ] 
Author Message
PostPosted: Mon Apr 14, 2014 7:00 pm 
Offline

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
 Profile  
 
PostPosted: Mon Apr 14, 2014 11:24 pm 
Offline

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
 Profile  
 
PostPosted: Tue Apr 15, 2014 2:07 pm 
Offline

Joined: Thu Apr 05, 2012 10:40 pm
Posts: 14
Perfect! Thank you.


Top
 Profile  
 
PostPosted: Wed Apr 16, 2014 12:08 am 
Offline

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
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 4 posts ] 

All times are UTC - 8 hours [ DST ]


Who is online

Users browsing this forum: No registered users and 1 guest


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

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group