#include <qxrdhistogramdialog.h>

Inheritance diagram for QxrdHistogramDialog:

Collaboration diagram for QxrdHistogramDialog:

Public Slots
void	histogramSelectionChanged (QRectF rect)

void	updateHistogramNeeded ()

Public Member Functions
	QxrdHistogramDialog (QxrdHistogramDialogSettingsWPtr settings, QxrdExperimentWPtr expt, QWidget *parent)

virtual	~QxrdHistogramDialog ()

void	onProcessedImageAvailable (QcepDoubleImageDataPtr image, QcepMaskDataPtr overflow)

Private Member Functions
void	recalculateHistogram ()

Private Attributes
QxrdExperimentWPtr	m_Experiment

QxrdHistogramDialogSettingsWPtr	m_HistogramDialogSettings

QcepDoubleImageDataPtr	m_Image

Detailed Description

Definition at line 13 of file qxrdhistogramdialog.h.

Constructor & Destructor Documentation

QxrdHistogramDialog::QxrdHistogramDialog	(	QxrdHistogramDialogSettingsWPtr	settings,
		QxrdExperimentWPtr	expt,
		QWidget *	parent
	)

explicit

Definition at line 10 of file qxrdhistogramdialog.cpp.

References DEBUG_CONSTRUCTORS, m_HistogramDialogSettings, and qcepDebug().

                                                           :
   QDockWidget(parent),
   m_Experiment(expt),
   m_HistogramDialogSettings(settings)
 {
   if (qcepDebug(DEBUG_CONSTRUCTORS)) {
     printf("QxrdHistogramDialog::QxrdHistogramDialog(%p)\n", this);
   }
 
   setupUi(this);
 
   QxrdHistogramDialogSettingsPtr set(m_HistogramDialogSettings);
 
   if (set) {
     m_HistogramPlot->init(set->histogramPlotSettings());
   }
 }

Here is the call graph for this function:

QxrdHistogramDialog::~QxrdHistogramDialog ( )

virtual

Definition at line 30 of file qxrdhistogramdialog.cpp.

References DEBUG_CONSTRUCTORS, and qcepDebug().

 {
   if (qcepDebug(DEBUG_CONSTRUCTORS)) {
     printf("QxrdHistogramDialog::~QxrdHistogramDialog(%p)\n", this);
   }
 }

Here is the call graph for this function:

Member Function Documentation

void QxrdHistogramDialog::histogramSelectionChanged ( QRectF rect )

slot

Definition at line 44 of file qxrdhistogramdialog.cpp.

References m_HistogramDialogSettings, and recalculateHistogram().

Referenced by QxrdWindow::initialize().

 {
   QxrdHistogramDialogSettingsPtr set(m_HistogramDialogSettings);
 
   if (set) {
     set->set_HistogramRect(rect);
 
     recalculateHistogram();
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QxrdHistogramDialog::onProcessedImageAvailable	(	QcepDoubleImageDataPtr	image,
		QcepMaskDataPtr	overflow
	)

Definition at line 37 of file qxrdhistogramdialog.cpp.

References m_Image, and recalculateHistogram().

 {
   m_Image = image;
 
   recalculateHistogram();
 }

Here is the call graph for this function:

void QxrdHistogramDialog::recalculateHistogram ( )

private

Definition at line 60 of file qxrdhistogramdialog.cpp.

References DEBUG_HISTOGRAM, m_Experiment, m_HistogramDialogSettings, m_Image, and qcepDebug().

Referenced by histogramSelectionChanged(), onProcessedImageAvailable(), and updateHistogramNeeded().

 {
   QxrdHistogramDialogSettingsPtr set(m_HistogramDialogSettings);
 
   if (set && m_Image) {
     QxrdExperimentPtr expt(m_Experiment);
 
     if (expt) {
       QTime tic;
       tic.start();
 
       QRectF rect = set->get_HistogramRect();
 
       int nsum = m_Image->get_SummedExposures();
 
       if (nsum < 1) {
         nsum = 1;
       }
 
       QxrdAcquisitionPtr acq(expt->acquisition());
 
       double satlev = 60000;
 
       if (acq) {
         satlev = acq->get_OverflowLevel();
       }
 
       double min = 0, max = satlev*nsum;
 
       const int nbins=1000;
 
       QcepDoubleVector x0(nbins), h0(nbins), h1(nbins);
 
       for (int i=0; i<nbins; i++) {
 //        double x = min+i*(max-min)/1000.0;
         x0[i] = i*100.0/(double)nbins;
         h0[i] = 0;
         h1[i] = 0;
       }
 
       int width = m_Image->get_Width();
       int height= m_Image->get_Height();
 
       double *data = m_Image->data();
 
       for (int i=0; i<width; i++) {
         for (int j=0; j<height; j++) {
           double v = *data++;
 
           int n;
 
           if (v<min) {
             n=0;
           } else if (v>max) {
             n=nbins-1;
           } else {
             n=(nbins-1)*(v-min)/(max-min);
           }
 
           if (n<0) {
             n=0;
           }
 
           if (n>=nbins) {
             n=nbins-1;
           }
 
           h0[n]++;
 
           if (rect.contains(i,j)) {
             h1[n]++;
           }
         }
       }
 
       m_HistogramPlot->detachItems();
 
       QwtPlotPiecewiseCurve *pc0 = new QwtPlotPiecewiseCurve(m_HistogramPlot, "Entire Image");
 
       pc0->setSamples(x0, h0);
       pc0->setPen(QPen(Qt::red));
 
       pc0->attach(m_HistogramPlot);
 
       QwtPlotPiecewiseCurve *pc1 = new QwtPlotPiecewiseCurve(m_HistogramPlot,
                                                              tr("[%1,%2]-[%3,%4]")
                                                              .arg(rect.left()).arg(rect.bottom())
                                                              .arg(rect.right()).arg(rect.top()));
 
       pc1->setSamples(x0, h1);
       pc1->setPen(QPen(Qt::darkRed));
 
       pc1->attach(m_HistogramPlot);
 
       m_HistogramPlot->replot();
 
       if (qcepDebug(DEBUG_HISTOGRAM)) {
         expt -> printMessage(tr("Histogram of data took %1 msec").arg(tic.elapsed()));
       }
     }
   }
 }