TotalExpTime.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// SVN: $Id: TotalExpTime.hxx 4695 2008-04-09 15:18:42Z EckhardSutorius $
//-----------------------------------------------------------------------------
#ifndef TOTALEXPTIME_H
#define TOTALEXPTIME_H

#include "Constants.h"
#include "DataMethod.hxx"
#include "FitsFile.h"
//-----------------------------------------------------------------------------
template<typename DataType>
class TotalExpTime : public DataMethod<DataType>
{
public:
  TotalExpTime(const FitsFile& aFile)
    : mFile(aFile)
  {
    DataMethod<DataType>::methodName = "TotalExpTime";
  }

  virtual ~TotalExpTime() { }

  void setSource(IntMap& tbls, IntMap& cols, StringMap units)
  {
    mExtnumColno = cols[0];
    if (mExtnumColno < 0)
    {
      throw DataMethodException("TotalExpTime method: cannot source");
    }
  }

  void setTarget(IntMap& cols, StringMap units)
  {
    mTotalExpTimeColno = cols[0];
  }

  void doit(TableData<DataType>& data, int row1, int row2)
  {
    double extNum, texptime, exptime;
    double njitter, nustep;
    bool noData = false;

    try
    {
      // get data from primary header
      mFile.movabsHdu(1);
      mFile.readKey("NJITTER", njitter);
      mFile.readKey("NUSTEP", nustep);
      mFile.readKey("EXP_TIME", exptime);
    }
    catch (FitsIOException)
    {
      noData = true;
    }

    for (int rowNo = row1; rowNo < row2; ++rowNo)
    {
      if (mFile.getFileName().find(PreSuffixNames::stackSuffix,
                                   mFile.getFileName().rfind("/"))
          != std::string::npos)
      {
        // move to the extension
        data.value(mExtnumColno, rowNo, extNum);
        mFile.movabsHdu(static_cast<int>(extNum));
        try
        { // if keyword is available get data from its value
          mFile.readKey(FitsFile::TEXPTIME, texptime);
        }
        catch (FitsIOException)
        {
          texptime = (!noData) ? njitter * nustep * exptime
                               : DefaultValues::realdef;
        }
      }
      else
      {
        texptime = DefaultValues::realdef;
      }
      data.assign(mTotalExpTimeColno, rowNo, texptime);
    }
  }

  // No data from raw data included
  void doit2(TableData<DataType>& rawData, TableData<DataType>& data, int row1, int row2) { }

private:
  FitsFile mFile;
  int mTotalExpTimeColno;
  int mExtnumColno;
};
//-----------------------------------------------------------------------------
#endif
//-----------------------------------------------------------------------------
// Change log:
//
// 10-Mar-2008, ETWS: First version.

---