TotalExpTime.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// SVN: $Id: TotalExpTime.hxx 7109 2010-06-17 17:54:34Z NicholasCross $
//-----------------------------------------------------------------------------
#ifndef TOTALEXPTIME_H
#define TOTALEXPTIME_H

#include "Constants.h"
#include "DataMethod.hxx"
#include "FitsFile.h"
#include "Logger.h"
//-----------------------------------------------------------------------------
template<typename DataType>
class TotalExpTime : public DataMethod<DataType>
{
public:
  Logger log;
  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, noffset;
    bool noData = false;
    bool noOffs = false;

    try
    {
      // get data from primary header
      if (mFile.isTile()) {
        mFile.movabsHdu(2);
      } else {
        mFile.movabsHdu(1);
      }
      mFile.readKey(FitsFile::cEXP_TIME, exptime);
      mFile.readKey("NJITTER", njitter);
      mFile.readKey("NUSTEP", nustep);
    }
    catch (FitsIOException)
    {
      noData = true;
    }

    try
    {
      // get data from primary header
      if (mFile.isTile()) {
        mFile.movabsHdu(2);
      } else {
        mFile.movabsHdu(1);
      }
      mFile.readKey("NOFFSETS", noffset);
    }
    catch (FitsIOException)
    {
      noOffs = true;
    }
    // Calculate offset value. 1, if not tile.
    double offVal = (!noOffs && mFile.isTile()) ? noffset/3. : 1.;
    
    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 * offVal
                               : 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.