quantumaharonovbohmplotwidget.cpp

Go to the documentation of this file.

/**
 * @file quantumaharonovbohmplotwidget.cpp
 * @brief Implements the quantumaharonovbohmplotwidget visualization component used by QPE.
 * @ingroup qpe_graphics
 *
 * Source: @c graphics/renderer/quantumaharonovbohmplotwidget.cpp. This file is part of Quantum Phenomena Explorer (QPE),
 * maintained by Officina Turini. See LICENSE_IT.txt / LICENSE.txt for the
 * QPE freeware source-available terms and THIRD_PARTY_NOTICES.txt for
 * third-party components.
 */
 
#include "quantumaharonovbohmplotwidget.h"
#include <QPainter>
#include <QPaintEvent>
#include <QPainterPath>
#include <algorithm>
#include <cmath>
 
namespace {
/**
 * @brief Performs heat.
 * @param v The v value used by the operation.
 */
QColor heat(double v)
{
    v = std::max(0.0, std::min(1.0, v));
    if (v < 0.5) {
        const double u = v / 0.5;
        return QColor(static_cast<int>(8 + 20*u), static_cast<int>(12 + 90*u), static_cast<int>(30 + 180*u));
    }
    const double u = (v - 0.5) / 0.5;
    return QColor(static_cast<int>(28 + 227*u), static_cast<int>(102 + 130*u), static_cast<int>(210 - 80*u));
}
}
 
QuantumAharonovBohmPlotWidget::QuantumAharonovBohmPlotWidget(QWidget *parent):QWidget(parent){}
void QuantumAharonovBohmPlotWidget::setData(const qpe::QuantumAharonovBohmHistory2D &h){history_=&h;update();}
void QuantumAharonovBohmPlotWidget::clear(){history_=nullptr;update();}
QSize QuantumAharonovBohmPlotWidget::minimumSizeHint() const{return QSize(650,430);}
 
void QuantumAharonovBohmPlotWidget::paintEvent(QPaintEvent *event)
{
    Q_UNUSED(event);
    QPainter p(this); p.setRenderHint(QPainter::Antialiasing,true); p.fillRect(rect(),QColor(10,14,20));
    QFont f=font();f.setBold(true);p.setFont(f);p.setPen(QColor(225,235,244));
    p.drawText(QRect(12,7,width()-24,24),Qt::AlignLeft|Qt::AlignVCenter,QStringLiteral("Aharonov-Bohm — detector interference vs enclosed flux"));
    if(!history_||!history_->configured()||history_->empty()){
        p.setFont(font());p.setPen(QColor(150,166,180));
        p.drawText(QRect(20,50,width()-40,height()-70),Qt::AlignCenter|Qt::TextWordWrap,
                   QStringLiteral("Load Aharonov-Bohm Effect / Magnetic Flux and press Play"));return;
    }
    const auto &samples=history_->samples(); const auto &cfg=history_->config();
    const int left=64,right=18,top=40,gap=34; const int traceH=130;
    QRect map(left,top,std::max(10,width()-left-right),std::max(100,height()-top-gap-traceH-28));
    QRect trace(left,map.bottom()+gap,map.width(),std::max(90,height()-(map.bottom()+gap)-28));
    p.setPen(QColor(55,72,88));p.setBrush(QColor(5,9,15));p.drawRect(map);p.drawRect(trace);
    double vmax=1e-15,cmax=1e-15;
    for(const auto&s:samples){for(float v:s.detectorProfileY)vmax=std::max(vmax,(double)v);cmax=std::max(cmax,s.detectorCenterSignal);}
    const std::size_t ny=history_->gridHeight();
    for(std::size_t si=0;si<samples.size();++si){
        const auto&s=samples[si];
        int x0=map.left()+int((double(si)/samples.size())*map.width());
        int x1=map.left()+int((double(si+1)/samples.size())*map.width());
        int ww=std::max(1,x1-x0+1);
        for(std::size_t yi=0;yi<ny&&yi<s.detectorProfileY.size();++yi){
            int y1=map.bottom()-int((double(yi)/ny)*map.height());
            int y0=map.bottom()-int((double(yi+1)/ny)*map.height());
            int hh=std::max(1,y1-y0+1);
            double a=std::sqrt(std::max(0.0,(double)s.detectorProfileY[yi])/vmax);
            p.fillRect(QRect(x0,y0,ww,hh),heat(a));
        }
    }
    const double tmax=std::max(cfg.analysisEndTime,samples.back().time);
    auto tx=[&](double t){return trace.left()+t/std::max(1e-12,tmax)*trace.width();};
    auto cy=[&](double v){return trace.bottom()-v/cmax*trace.height();};
    QPainterPath path;for(std::size_t i=0;i<samples.size();++i){QPointF pt(tx(samples[i].time),cy(samples[i].detectorCenterSignal));if(i==0)path.moveTo(pt);else path.lineTo(pt);} 
    p.setPen(QPen(QColor(90,225,255),2.0));p.drawPath(path);
    p.setPen(QColor(180,198,212));p.setFont(font());
    p.drawText(QRect(4,map.top(),left-10,20),Qt::AlignRight|Qt::AlignVCenter,QStringLiteral("+y"));
    p.drawText(QRect(4,map.bottom()-20,left-10,20),Qt::AlignRight|Qt::AlignVCenter,QStringLiteral("-y"));
    p.drawText(QRect(map.left(),map.bottom()+4,map.width(),20),Qt::AlignCenter,QStringLiteral("time  —  detector density profile rho(xd,y,t)"));
    p.drawText(QRect(trace.left(),trace.bottom()+4,trace.width(),20),Qt::AlignCenter,QStringLiteral("time  —  central fringe intensity"));
    const auto &last=samples.back(); const double deg=360.0*cfg.fluxQuanta;
    p.setPen(QColor(232,238,244));
    p.drawText(QRect(trace.left()+8,trace.top()+5,trace.width()-16,22),Qt::AlignLeft|Qt::AlignVCenter,
               QStringLiteral("Phi/Phi0=%1   AB phase=%2 deg   center=%3   line=%4   norm=%5")
               .arg(cfg.fluxQuanta,0,'f',3).arg(deg,0,'f',1).arg(last.detectorCenterSignal,0,'g',5)
               .arg(last.detectorSignal,0,'g',5).arg(last.norm,0,'f',6));
}