void InstrumentWindowPickTab::updatePlot(int detid)
{
if (m_instrWindow->blocked())
{
m_plot->clearCurve();
return;
}
if (m_plotPanel->isCollapsed()) return;
InstrumentActor* instrActor = m_instrWindow->getInstrumentActor();
Mantid::API::MatrixWorkspace_const_sptr ws = instrActor->getWorkspace();
if (detid >= 0)
{
if (m_one->isChecked() || m_peak->isChecked())
{// plot spectrum of a single detector
plotSingle(detid);
}
else if (m_tube->isChecked())
{// plot integrals
plotTube(detid);
}
}
else
{
m_plot->clearCurve();
}
m_plot->recalcAxisDivs();
m_plot->replot();
}
void InstrumentWindowPickTab::getBinMinMaxIndex(size_t wi,size_t& imin, size_t& imax)
{
InstrumentActor* instrActor = m_instrWindow->getInstrumentActor();
Mantid::API::MatrixWorkspace_const_sptr ws = instrActor->getWorkspace();
const Mantid::MantidVec& x = ws->readX(wi);
if (instrActor->wholeRange())
{
imin = 0;
imax = x.size() - 1;
}
else
{
Mantid::MantidVec::const_iterator x_begin = std::lower_bound(x.begin(),x.end(),instrActor->minBinValue());
Mantid::MantidVec::const_iterator x_end = std::lower_bound(x.begin(),x.end(),instrActor->maxBinValue());
imin = static_cast<size_t>(x_begin - x.begin());
imax = static_cast<size_t>(x_end - x.begin()) - 1;
}
}
/**
* Plot data for a detector.
* @param detid :: ID of the detector to be plotted.
*/
void InstrumentWindowPickTab::plotSingle(int detid)
{
m_plot->clearLabels();
InstrumentActor* instrActor = m_instrWindow->getInstrumentActor();
Mantid::API::MatrixWorkspace_const_sptr ws = instrActor->getWorkspace();
size_t wi;
try {
wi = instrActor->getWorkspaceIndex(detid);
} catch (Mantid::Kernel::Exception::NotFoundError) {
return; // Detector doesn't have a workspace index relating to it
}
// get the data
const Mantid::MantidVec& x = ws->readX(wi);
const Mantid::MantidVec& y = ws->readY(wi);
// find min and max for x
size_t imin,imax;
getBinMinMaxIndex(wi,imin,imax);
Mantid::MantidVec::const_iterator y_begin = y.begin() + imin;
Mantid::MantidVec::const_iterator y_end = y.begin() + imax;
m_plot->setXScale(x[imin],x[imax]);
// fins min and max for y
Mantid::MantidVec::const_iterator min_it = std::min_element(y_begin,y_end);
Mantid::MantidVec::const_iterator max_it = std::max_element(y_begin,y_end);
// set the data
m_plot->setData(&x[0],&y[0],static_cast<int>(y.size()));
m_plot->setYScale(*min_it,*max_it);
// find any markers
ProjectionSurface* surface = mInstrumentDisplay->getSurface();
if (surface)
{
QList<PeakMarker2D*> markers = surface->getMarkersWithID(detid);
foreach(PeakMarker2D* marker,markers)
{
m_plot->addLabel(new PeakLabel(marker));
//std::cerr << marker->getLabel().toStdString() << std::endl;
}
}
/**
* Update the info window with information for a selected detector.
* @param detid :: ID of the selected detector.
*/
void InstrumentWindowPickTab::updateSelectionInfo(int detid)
{
if (m_freezePlot)
{ // freeze the plot for one update
m_freezePlot = false;
return;
}
if (m_instrWindow->blocked())
{
m_selectionInfoDisplay->clear();
return;
}
if (detid >= 0)
{
InstrumentActor* instrActor = m_instrWindow->getInstrumentActor();
Mantid::Geometry::IDetector_const_sptr det = instrActor->getInstrument()->getDetector(detid);
QString text = "Selected detector: " + QString::fromStdString(det->getName()) + "\n";
text += "Detector ID: " + QString::number(detid) + '\n';
QString wsIndex;
try {
wsIndex = QString::number(instrActor->getWorkspaceIndex(detid));
updatePlot(detid); // Update the plot if the detector links to some data
} catch (Mantid::Kernel::Exception::NotFoundError &) {
// Detector doesn't have a workspace index relating to it
wsIndex = "None";
m_plot->clearCurve(); // Clear the plot window
m_plot->replot();
}
text += "Workspace index: " + wsIndex + '\n';
Mantid::Kernel::V3D pos = det->getPos();
text += "xyz: " + QString::number(pos.X()) + "," + QString::number(pos.Y()) + "," + QString::number(pos.Z()) + '\n';
double r,t,p;
pos.getSpherical(r,t,p);
text += "rtp: " + QString::number(r) + "," + QString::number(t) + "," + QString::number(p) + '\n';
Mantid::Geometry::ICompAssembly_const_sptr parent = boost::dynamic_pointer_cast<const Mantid::Geometry::ICompAssembly>(det->getParent());
if (parent)
{
QString textpath;
while (parent)
{
textpath="/"+QString::fromStdString(parent->getName())+textpath;
parent=boost::dynamic_pointer_cast<const Mantid::Geometry::ICompAssembly>(parent->getParent());
}
text += "Component path:" +textpath+"/"+ QString::fromStdString(det->getName()) +'\n';
}
const double integrated = instrActor->getIntegratedCounts(detid);
const QString counts = integrated == -1.0 ? "N/A" : QString::number(integrated);
text += "Counts: " + counts + '\n';
QString xUnits;
if (m_selectionType > SingleDetectorSelection && !m_plotSum)
{
switch(m_tubeXUnits)
{
case DETECTOR_ID:
xUnits = "Detector ID";
break;
case LENGTH:
xUnits = "Length";
break;
case PHI:
xUnits = "Phi";
break;
default:
xUnits = "Detector ID";
}
}
else
{
xUnits = QString::fromStdString(instrActor->getWorkspace()->getAxis(0)->unit()->caption());
//xUnits = "Time of flight";
}
text += "X units: " + xUnits + '\n';
m_selectionInfoDisplay->setText(text);
}
else
{
m_selectionInfoDisplay->clear();
m_plot->clearCurve(); // Clear the plot window
m_plot->replot();
}
}