View unanswered posts | View active topics
|
Page 1 of 1
|
[ 8 posts ] |
|
Author |
Message |
paramita
|
Posted: Thu May 07, 2015 9:53 am |
|
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 |
|
|
gudlaugu
|
Posted: Fri Jun 02, 2017 3:31 am |
|
|
Site Admin |
Joined: Fri Jul 18, 2008 3:04 pm Posts: 61 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 |
|
|
massimo.persic
|
Posted: Wed Sep 13, 2017 8:00 am |
|
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 |
|
|
gudlaugu
|
Posted: Wed Sep 13, 2017 8:05 am |
|
|
Site Admin |
Joined: Fri Jul 18, 2008 3:04 pm Posts: 61 Location: Stanford
|
If you go into advanced mode you can set the output_gcr_full option.
_________________ Gudlaugur Johannesson, GALPROP developer
|
|
Top |
|
|
massimo.persic
|
Posted: Mon Oct 16, 2017 7:12 am |
|
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 |
|
|
gudlaugu
|
Posted: Mon Oct 16, 2017 9:12 am |
|
|
Site Admin |
Joined: Fri Jul 18, 2008 3:04 pm Posts: 61 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 |
|
|
massimo.persic
|
Posted: Tue Oct 17, 2017 8:30 am |
|
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 |
|
|
gudlaugu
|
Posted: Tue Oct 17, 2017 8:32 am |
|
|
Site Admin |
Joined: Fri Jul 18, 2008 3:04 pm Posts: 61 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 |
|
|
|
Page 1 of 1
|
[ 8 posts ] |
|
Who is online |
Users browsing this forum: No registered users and 2 guests |
|
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
|