| GALPROP https://galprop.stanford.edu/forum/ |
|
| CR from DM annihilation with a cuspy DM profile. https://galprop.stanford.edu/forum/viewtopic.php?f=7&t=79 |
Page 1 of 1 |
| Author: | jhhuh [ Fri Sep 25, 2009 4:38 am ] |
| Post subject: | CR from DM annihilation with a cuspy DM profile. |
I used galprop v.50p for my recent paper about CR from decaying DM. While I used it, I found some bug in the code averaging DM density. ( "DM_profile_av()" in "gen_DM_source.cc") The calculation with 3-d coordinate or non-cuspy profile is OK. But the calculation with cuspy profile like NFW in 2-d coordinate make unacceptably large CR fluxes. It maybe the problem in averaging DM near the origin (Jacobi factor become significant in this region.). By simply changing "int nuse" into "double ruse" and "nuse++" into "ruse+=rr", code produce appropriate result. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ double Galprop::DM_profile_av(double r,double z,double dr,double dz,double dzz) { double DM_profile_av_=0.0; int nuse=0; double ruse=0; for (double zz=z-dz/2.; zz<=z+dz/2.; zz+=dzz) for (double rr=r-dr/2.; rr<=r+dr/2.; rr+=dr/10.) { if (rr<0.) continue; DM_profile_av_+=rr*DM_profile(rr,0,zz); ruse+=rr; } return (DM_profile_av_/ruse); } ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|
| Page 1 of 1 | All times are UTC - 8 hours [ DST ] |
| Powered by phpBB® Forum Software © phpBB Group https://www.phpbb.com/ |
|