It is currently Thu Mar 30, 2017 5:33 pm

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 1 post ] 
Author Message
PostPosted: Wed May 27, 2015 5:13 pm 
Offline

Joined: Thu Apr 05, 2012 10:40 pm
Posts: 14
Hi All,

I am trying to calculate propagation greens functions for antiprotons from dark matter using galprop. I have implemented a simple NFW profile and am using the gen_DM_Source.cc with the gen_DM_source routine containing the line:

particle.secondary_source_function.d2[ir][iz].s[ip]+=pow(DM_profile_av(galaxy.r[ir], galaxy.z[iz], galaxy.dr, galaxy.dz, dzz),2);

if we are in the energy bin of interest -- i.e. I inject antiprotons proportional to the NFW squared (in units GeV^2/cm^6) into a single energy bin and look at the resulting antiproton output, then repeat the galprop run for each energy bin.

I am trying to understand the resulting output units after propagation. The shape looks approximately correct if I compare to semi-analytic results for propagation functions (e.g. fig 9 of http://arxiv.org/pdf/1012.4515v4.pdf), but the normalization doesn't match. assuming a local rho_dm of 0.3 GeV/cm^3 I obtain a value at 1 GeV of around 10^21.

Any help would be much appreciated. Picture of resulting runs is also attached where one can see the Green functions for one of the models.


Attachments:
green_func.png
green_func.png [ 79.49 KiB | Viewed 2475 times ]
Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

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