PacCoincidenceSpectrum.h

Go to the documentation of this file.

/***************************************************************************
 *   Copyright (C) 2008-2010 by Matthias Nagl                              *
 *   mnagl@uni-goettingen.de                                               *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

#ifndef PACCOINCIDENCEPARAMETERS_H
#define PACCOINCIDENCEPARAMETERS_H

#include <QObject>
#include <QVector>
#include <QPointer>
#include <QPair>
#include <QMutex>
#include <qwt_double_rect.h>
#include "PacParameterState.h"
#include <Minuit2/FunctionMinimum.h>
//#include "PacCoincidenceControlWidget.h"
#include "PacCoincidenceParameterDialog.h"
#include "pacCoincidenceDefinitions.h"

class QPushButton;
class QRadioButton;
class QDoubleSpinBox;
class QDomElement;
class QwtPlotMarker;
class QwtPlotCurve;
class PacPlot;
class PacCoincidenceContainer;
class PacAbstractFitFunction;
class PacBasicMeasurementInfo;

class PacCoincidenceSpectrum : public QObject
{
  Q_OBJECT
public:
  PacCoincidenceSpectrum(
      PacCoincidenceContainer *container,
      const PacBasicMeasurementInfo *basicinfo,
      QPair<quint8,quint8> &detectoridxs
  );
  ~PacCoincidenceSpectrum();

  QWidget * plotWidget();
  QWidget * controlWidget();
  QWidget * parameterDialogWidget();

  void setSpectrum(const QVector<quint64> &spectY, bool addendum);


  QVector<double> spectrumX() const;
  QVector<quint64> spectrumY() const;
  QVector<double> doubleSpectrumY() const;
  quint64 maxY() const;
  quint32 maxIdx() const;
  double angle() const;
  quint8 startNodeId() const;
  quint8 stopNodeId() const;
  qreal timeResolutionFWHM() const;


  bool ignoreThisSpectrum() const;
  Pac::ParameterDeterminationMethod determinationMethod() const;
  qreal t0FitCutterPosition() const;
  qreal backgroundFitStartPosition() const;
  qreal timeOffset() const;
  qreal backgroundLevel() const;


  void receiveFitResult(const ROOT::Minuit2::ParametricFunction *fitfunc, const ROOT::Minuit2::FunctionMinimum &result);

  void updateButtonBackground();
  void appendToDomNode(QDomElement &node) const;
  void initializeFromDomNode(const QDomElement &node);

public slots:

  void setIgnoreThisSpectrum(bool ignore);
  void setDeterminationMethod(Pac::ParameterDeterminationMethod m);
  void setT0FitCutterPosition(double percent);
  void setBackgroundFitStartPosition(double starttime, bool refit=true);
  void setTimeOffset(double timeoffset);
  void setBackgroundLevel(double counts);
  void t0ParametersUpdated();
  void backgroundParametersUpdated();


  void autoEstimateT0FitParameters();
  void autoEstimateBgFitParameters();


  void restartFitting();
  void sendT0FunctionToFitThread(bool force = false);
  void sendBackgroundFunctionToFitThread(bool force = false);


  void setButton(QPushButton *button);

private slots:
  void loadReceivedT0FitResults();
  void loadReceivedBgFitResults();
  void selectMarker(QwtDoublePoint p);
  void applyNewMarkerPosition(QwtDoublePoint p);
  void methodChanged(bool usefit);
  void updateLockButtons();
  void setMethodLock(bool locked);
  void setT0Lock(bool locked);
  void setBgLock(bool locked);

signals:
  void timeOffsetChanged(double newtimeoffset);
  void backgroundLevelChanged(double newcounts);
  void updateRTSpectrum();

private:
  void setT0FitCurve();
  void setBackgroundFitCurve();
  void plotSpectrum();
  void lockT0Widgets();
  void lockBackgroundWidgets();
  qreal autoEstimateBackground() const;
  qreal autoEstimateHalfLife() const;
  qreal autoEstimateExponentialAmplitude() const;


  PacCoincidenceContainer *container;
  const PacBasicMeasurementInfo *basicinfo;
  QPointer<QPushButton> button; 



  QPair<quint8,quint8> detectoridxs; 



  QVector<quint64> spectY;


  Pac::ParameterDeterminationMethod method;
  qreal t0fitcutterposition; 
  qreal bgfitcutterposition; 
  qreal toffset; 
  qreal bglevel;
  bool ignorethisspectrum; 

  PacParameterState t0fitparams;
  bool t0successfullyfitted;
  PacParameterState bgfitparams;
  bool bgsuccessfullyfitted;


  QPointer<QWidget> plotwidget;
  QPointer<QWidget> controlwidget;
  QPointer<PacCoincidenceParameterDialog> parameterdialog;

  PacPlot *plot;
  QwtPlotMarker *convfitcutter, *expfitcutter, *backgroundlevelmarker, *toffsetmarker;
  QwtPlotCurve *coincidencespectrum, *convfitcurve, *expfitcurve;

  QRadioButton *usefitmethod, *manuallysetparams;
  QDoubleSpinBox *tospinbox, *bgspinbox, *t0fitcutterpos, *bgfitstartpos;
  QToolButton *methodlock, *t0lock, *bglock;
  QGroupBox *methodgroupbox;


  mutable qreal detectorangle;
  mutable QVector<double> doublespecty;
  mutable quint64 max;
  mutable qint32 maxidx;
  mutable qreal bgautoestimate;
  mutable qreal hlautoestimate;
  mutable qreal expamplitudeautoestimate;
  const PacAbstractFitFunction *t0fitfunction;
  const PacAbstractFitFunction *bgfitfunction;
  ROOT::Minuit2::FunctionMinimum *receivedt0minimum, *receivedbgminimum;
  QTimer *processt0fit, *processbgfit, *startdeferredt0fit, *startdeferredbgfit;
  bool restartt0fitasap, restartbgfitasap;
  double commitbgfitstartposchangeasap;
  QMutex newparammutex;
  bool t0locked, bglocked; 
  enum CoincidenceMarkerSelection {
    Background,
    Offset,
    ConvolvedCutter,
    ExponentialCutter
  } selectedcoincidencemarker;
};

#endif // PACCOINCIDENCEPARAMETERS_H