Hello,

I want to get flux of a primary/secondary particle at different (R,z) or (x,y,z) w.r.to Galactic centre how can I plot the flux distribution as a function of (R,z) or (x,y,z). I can see pion intensity skymaps which gives intensity as a function of (l,b,E)using ds9 software.But I want to know the flux at different (R,z) point in the Galaxy.

Please help me to configure the galdef file accordingly.

Hoping for your kind response

Paramita

The nuclei file in the output gives you the CR flux at the grid points in the plane. The units of this file is MeV/(cm^2 sr s). If you want z dependence as well you will need to specify output_gcr_full in the galdef file and use the nuclei_full file.

_________________
Gudlaugur Johannesson, GALPROP developer

Is the suggested procedure possible on WebRun, or needs one to download galprop and run it locally? Thanks.

If you go into advanced mode you can set the output_gcr_full option.

_________________
Gudlaugur Johannesson, GALPROP developer

Hello, above it was stated that the nuclei file gives the CR flux at the grid points in the Galactic plane in units of MeV/(cm2 s sr). So, from webrunning galprop, I expected to obtain the same quantity in the same units by multiplying the output spectrum, given in (1 / cm2 s sr MeV), times E^2 (given in MeV^2). However, DS9 reads values in the E-9 - E-10 range from the nuclei file, whereas by webrunning galprop (assuming typical values fror the parameters to be input) I get values in the E+2 - E-7 MeV/(cm2 s sr) range. If both procedures return CR fluxes in MeV/(cm2 s sr), why don't I obtain consistent results? Can anybody please help me? Thanks.

The units in the nuclei files are E^2 times flux as I stated above. I am not sure what you mean by webrunning galprop, that is just a regular install of galprop run on our cluster. Without further information I cannot really help you more. Note that the 3D nuclei_full file has the coordinates x,y,z,E,isotope so by default ds9 will show the x-y planes and you have to select different z,E, and isotope coordinates with the additional window that pops up.

_________________
Gudlaugur Johannesson, GALPROP developer

Thank you for your prompt reply. Reading the DS9 Users Manual, I think the pop-up window you are referring to is the "Bin --> Binning Parameters" dialog box. That would allow me to choose any column among x,y,z,E,A. Unfortunately, the "Bin" dialog box is disabled. Using the nuclei file, or the "primary electron (or proton) source function" (extracted from the nuclei_full file), the Bin dialog box still remains disabled. So I am not able to read the 5 columns contained in the nuclei_full_file. Can you figure out where I'm going wrong? Thank you for your help!

I suggest you read up on opening the fits file in your favorite script language (I use python and the fits module from astropy). It allows much more flexibility than ds9 that was really designed to view skymaps, not galprop nuclei files.

_________________
Gudlaugur Johannesson, GALPROP developer