qsedatadeglitcher.cpp

Go to the documentation of this file.

#include "qsedatadeglitcher.h"
#include "qsedataset.h"
#include "qsescan.h"

QseDataDeglitcher::QseDataDeglitcher
(QseDataSet* inData, QseDataSet* outData, const QString& name, QObject* parent)
  : QseEXAFSOperator(inData, outData, name, parent)
{
}

void QseDataDeglitcher::exec()
{
  m_OutputData -> clear();

  int ns = m_InputData -> scanCount();

  for (int s = 0; s < ns; s++) {
    QseScan *osc = m_OutputData -> appendScan();

    deglitchScan(s, osc);
  }

  m_OutputData -> emitOperationCompleted();
}

void QseDataDeglitcher::deglitchScan(int n, QseScan* osc)
{
//   QseScan *sc = m_InputData->scan(n);

//   if (sc) {
//     int nc = sc -> columnCount();
//     int energycol = -1;

//     for (int c = 0; c < nc; c++) {
//       QseColumn *col = sc -> column(c);

//       if (matchesEnergyColumn(col)) {
//      energycol = c;
//       }
//     }

//     if (energycol < 0) {
//       emit message(QString("No energy column found for scan %1").arg(n));
//     } else {
//       osc -> mergeHeaders(sc->header());
//       osc -> setScanNumber(sc->scanNumber());

//       QseColumn *energy = osc -> appendColumn("Energy", sc -> column(energycol));
//       QseColumn *normal = osc -> appendColumn("Norm");
//       QseColumn *detect = osc -> appendColumn("Detector");
//       QseColumn *seconds = osc -> appendColumn("Seconds");

//       energy -> multiplication(energyScalingRatio());

//       int nnorm = 0, ndet = 0, nsecs = 0;
      
//       for (int c = 0; c < nc; c++) {
//      QseColumn *col = sc -> column(c);

//      if (matchesNormalizationColumn(col)) {
//        nnorm += 1;
//        normal -> addition(col);
//      } else if (matchesDetectorColumn(col)) {
//        ndet += 1;
//        detect -> addition(col);
//      } else if (matchesSecondsColumn(col)) {
//        nsecs += 1;
//        seconds -> addition(col);
//      }
//       }
      
//       if (nnorm > 0) {
//      QseColumn *ratio = osc -> appendColumn("Ratio", detect);
        
//      ratio -> division(normal);
        
//      emit message(QString("Scan %1 had %2 detector columns, was normalized by %3 ion chambers").arg(n).arg(ndet).arg(nnorm));
//       } else if (nsecs > 0) {
//      QseColumn *ratio = osc -> appendColumn("Ratio", detect);
        
//      ratio -> division(seconds);
//      emit message(QString("Scan %1 had %2 detector columns, was normalized by counting time").arg(n).arg(ndet));
//       } else {
//      QseColumn *ratio = osc -> appendColumn("Ratio", detect);
        
//      emit message(QString("Scan %1 had %2 detector columns, no normalization performed").arg(n).arg(ndet));
//       }
//     }
//   }
}