Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
FITS_tutorial8.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_FITS
3/// \notebook
4/// Open a FITS file with columns containing variable-length arrays.
5///
6/// ABOUT THE DATA: To illustrate such a case we use a Redistribution
7/// Matrix File (RMF) conform with the HEASARC OGIP specifications
8/// OGIP definition
9/// https://heasarc.gsfc.nasa.gov/docs/heasarc/ofwg/docs/spectra/ogip_92_007/ogip_92_007.html
10/// RMF defintion
11/// https://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/docs/memos/cal_gen_92_002/cal_gen_92_002.html#tth_sEc3.1
12///
13/// Variable-length arrays are used to return the non-null row elements of a
14/// sparse matrix containing the probability of detecting a photon of a given
15/// energy at a given detector channel.
16///
17/// The test data file, rmf.fits, is taken from the test data of the sherpa X-ray
18/// analysis tools
19/// https://cxc.harvard.edu/sherpa/
20/// the original file is available under the repository
21/// https://github.com/sherpa/sherpa-test-data
22///
23/// \macro_code
24/// \macro_output
25///
26/// \author Cosimo Nigro
27
28void FITS_tutorial8()
29{
30 // let us read the first header data unit, opening the file with a FITS viewer
31 // we see that the following cells contain variable-length arrays with values
32 // row 215 column "F_CHAN": (1, 68)
33 // row 215 column "N_CHAN": (66, 130)
34 // row 215 column "MATRIX" : (5.8425176E-6, 7.290097E-6, 8.188037E-6, 9.157882E-6, 1.018355E-5, ..)
35 TString dir = gROOT->GetTutorialDir();
36 TFITSHDU* hdu = new TFITSHDU(dir + "/fitsio/rmf.fits", 1);
37 int rownum = 214; // FITS tables are indexed starting from 1
38 TString colname1 = "F_CHAN";
39 TString colname2 = "N_CHAN";
40 TString colname3 = "MATRIX";
41
42 printf("reading in row %d, column %s \n", rownum+1, colname1.Data());
43 TArrayD *arr1 = hdu->GetTabVarLengthVectorCell(rownum, colname1);
44 printf("(%f, %f) \n", arr1->At(0), arr1->At(1));
45
46 printf("reading in row %d, column %s \n", rownum+1, colname2.Data());
47 TArrayD *arr2 = hdu->GetTabVarLengthVectorCell(rownum, colname2);
48 printf("(%f, %f) \n", arr2->At(0), arr2->At(1));
49
50 // printing only the first 5 values in the array
51 printf("reading in row %d, column %s \n", rownum+1, colname3.Data());
52 TArrayD *arr3 = hdu->GetTabVarLengthVectorCell(rownum, colname3);
53 printf("(%e, %e, %e, %e, %e, ...) \n", arr3->At(0), arr3->At(1), arr3->At(2), arr3->At(3), arr3->At(4));
54
55}
#define gROOT
Definition TROOT.h:406
Array of doubles (64 bits per element).
Definition TArrayD.h:27
Double_t At(Int_t i) const
Definition TArrayD.h:79
FITS file interface class.
Definition TFITS.h:35
TArrayD * GetTabVarLengthVectorCell(Int_t rownum, Int_t colnum)
Get the variable-length array contained in a cell given by (row>=0, column name)
Definition TFITS.cxx:1583
Basic string class.
Definition TString.h:139
const char * Data() const
Definition TString.h:376