It is currently Mon Dec 06, 2021 10:24 pm

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 8 posts ] 
Author Message
PostPosted: Thu May 07, 2015 9:53 am 
Offline

Joined: Wed Mar 18, 2015 9:17 am
Posts: 2
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


Top
 Profile  
 
PostPosted: Fri Jun 02, 2017 3:31 am 
Offline
Site Admin

Joined: Fri Jul 18, 2008 3:04 pm
Posts: 60
Location: Stanford
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


Top
 Profile  
 
PostPosted: Wed Sep 13, 2017 8:00 am 
Offline

Joined: Wed Oct 27, 2010 3:19 am
Posts: 3
Is the suggested procedure possible on WebRun, or needs one to download galprop and run it locally? Thanks.


Top
 Profile  
 
PostPosted: Wed Sep 13, 2017 8:05 am 
Offline
Site Admin

Joined: Fri Jul 18, 2008 3:04 pm
Posts: 60
Location: Stanford
If you go into advanced mode you can set the output_gcr_full option.

_________________
Gudlaugur Johannesson, GALPROP developer


Top
 Profile  
 
PostPosted: Mon Oct 16, 2017 7:12 am 
Offline

Joined: Wed Oct 27, 2010 3:19 am
Posts: 3
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.


Top
 Profile  
 
PostPosted: Mon Oct 16, 2017 9:12 am 
Offline
Site Admin

Joined: Fri Jul 18, 2008 3:04 pm
Posts: 60
Location: Stanford
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


Top
 Profile  
 
PostPosted: Tue Oct 17, 2017 8:30 am 
Offline

Joined: Wed Oct 27, 2010 3:19 am
Posts: 3
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!


Top
 Profile  
 
PostPosted: Tue Oct 17, 2017 8:32 am 
Offline
Site Admin

Joined: Fri Jul 18, 2008 3:04 pm
Posts: 60
Location: Stanford
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


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