int main	(	int	argc,
		const char **	argv
	)

Definition at line 199 of file simpleFitsDumper.cc.

References CalculateBufferSize(), i, numRows, size, and writeValuesFromFits().

{
    //set the names of the file and table to be loaded
    string fileNameToLoad = "test.fits";
    string tableNameToLoad = "FACT-TIME_ETIENNE";
    //set the vector of columns to be dumped
    vector<string> columnsToDump;
    columnsToDump.push_back("Data0");
    columnsToDump.push_back("Data1");
    //set the name of the output text file
    string outputFile = "output.txt";

    //load the fits file
    CCfits::FITS* file = NULL;
    try
    {
        file = new CCfits::FITS(fileNameToLoad);
    }
    catch (CCfits::FitsException e)
    {
         cout << "Could not open FITS file " << fileNameToLoad << " reason: " << e.message() << endl;
         return -1;
    }
    //check if the selected table indeed exists in the loaded file. If so, load it. Otherwise display the existing tables
    CCfits::Table* table;
    const multimap< string, CCfits::ExtHDU * > extMap = file->extension();
    if (extMap.find(tableNameToLoad) == extMap.end())
    {
        cout << "Could not open table " << tableNameToLoad << ". Tables in file are: " << endl;
        for (std::multimap<string, CCfits::ExtHDU*>::const_iterator it=extMap.begin(); it != extMap.end(); it++)
            cout << it->first << " ";
        cout << endl;
        delete file;
        return -1;
    }
    else
        table = dynamic_cast<CCfits::Table*>(extMap.find(tableNameToLoad)->second);
    int numRows = table->rows();
    //check that the given column names are indeed part of that table
    map<string, CCfits::Column*> colMap = table->column();
    if (columnsToDump.size() != 0)
    {
        for (vector<string>::iterator it=columnsToDump.begin(); it!= columnsToDump.end(); it++)
        {
            if (colMap.find(*it) == colMap.end())
            {
                cout << "Config-given dump list contains invalid entry " << *it << endl;
                delete file;
                return -1;
            }
        }
    }
    //dump the requested columns
    table->makeThisCurrent();
    vector<int> offsets = CalculateBufferSize(colMap);
    int size = offsets[offsets.size()-1];
    offsets.pop_back();
    unsigned char* fitsBuffer = new unsigned char[size];

    ofstream targetFile(outputFile.c_str());
    int status = 0;

    for (int i=1;i<=table->rows(); i++)
    {
        fits_read_tblbytes(file->fitsPointer(), i, 1, size, fitsBuffer, &status);
        if (status)
        {
            cout << "An error occurred while reading fits row #" << i << " error code: " << status << endl;
            for (unsigned int j=0;j<offsets.size(); j++)
                cout << offsets[j] << " ";
            cout << endl;
        }
        writeValuesFromFits(offsets, targetFile, fitsBuffer, columnsToDump, colMap);
    }
    delete file;
    return 0;
}

Here is the call graph for this function: