It is currently Wed Nov 20, 2019 7:04 pm

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 20 posts ]  Go to page 1, 2  Next
Author Message
PostPosted: Mon Sep 21, 2015 4:53 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

From the example topic in the webpage, I downloaded "nuclei_44_999726pub.gz", this contains a .fits file, how I can extract (in a mac) the antiproton information? I mean, Energy v/s Flux?

Best

Nicolas


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 5:38 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
That is a very old example, but anyway the format is explained in Section 3.4 of the GALPROP Explanatory Supplement, the latest version and examples are at
http://sourceforge.net/projects/galprop/


As explained there, in file nuclei_.... the nuclei are ordered on axis 4 (for 2D) and antiprotons have
NUCZxxx = -1
NUCAxxx = 1
Axis 1 corresponds to R so choose the bin at R=8.5 for local values (using CRVAL1 and CDELT2).
Axis 2 has just one entry, for z=0
Axis 3 is the flux as a function of energy, log scale, see the explanation in the Supplement.
The units are also explained in the Supplement.


To read the files you can use the HEASARC fv tool, available for mac.
There is also a standalone c++ converter for these files at http://sourceforge.net/projects/gcrconvert/.

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 5:47 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

Ok, there is a more easy way to obtain from GALPROP the antiproton flux v/s energy in a table format,, in order to directly use,,, handle .fits format is not easy to me.


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 5:51 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
all GALPROP output is in FITS only. But fv is very easy and can dump it as a table.

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 5:56 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

I have fv open, could you guide me in order to obtain the table, and on which example?

Best and thanks

Nicolás


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 6:16 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
download case reltest19 included in the above distribution, then I will describe that.

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 6:22 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
In the FITS folder I only see files with reltest02,,,,


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 6:25 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Ok, I found it


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 7:12 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

Recently I can open the fits with pyrites,,, but I don't know how to continue


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 7:35 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
I don't know pyrites. Anyway you should view the FITS header with it, and display the data as tables. The table is 4D (R,z,energy,nucleus) : can you see it?

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 7:43 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
YEs, I saw 5 columns, see the attached file


Attachments:
Screenshot 2015-09-21 11.42.49.png
Screenshot 2015-09-21 11.42.49.png [ 85.24 KiB | Viewed 10355 times ]
Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 8:06 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
look at the header to find the number N for antiprotons with NUCZ=-1 NUCA=1
then the spectrum is in array (8, 1, *, N)
where 8 is R=8 kpc roughly (since gridsize is CDELT1=1 kpc) (or interpolate between 8 and 9)
1 is z=0 i.e. in the Galactic plane
* are the 77 spectrum values (log10(E)=CRVAL3+ CDELT3 * (i-1), i=1-77

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 9:01 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

Ok, in which folder of Galprop_54 I can find the example that you said "reltest19" what I found is "galdef_54_reltest19" in the Galdef folder, but that file is not a fits file.

Best

Nicolás


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 9:23 am 
Offline

Joined: Wed May 17, 2006 7:02 am
Posts: 285
Location: MPE Garching
in the tar file there should be a folder FITS which contains the complete output files for this case.

_________________
Andy Strong, MPE


Top
 Profile  
 
PostPosted: Mon Sep 21, 2015 9:32 am 
Offline

Joined: Fri Jul 18, 2014 6:14 pm
Posts: 14
Hi,

Ok, N=2 , so what I want is in the array (8, 1, *, 2), and how I can extract this info in fv??

Best and Thanks

Nicolás


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 20 posts ]  Go to page 1, 2  Next

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:  
cron
Powered by phpBB® Forum Software © phpBB Group