PtolemyGUI/Cleopatra/ExWindow.py

#!/usr/bin/python3

import os
import time
from PyQt6.QtWidgets import (
  QVBoxLayout, QWidget, QPushButton, QFileDialog
)
from PyQt6.QtCore import QUrl

from PyQt6.QtWebEngineWidgets import QWebEngineView
import plotly.graph_objects as go
import plotly.io as pio


from IAEANuclearData import IsotopeClass

class ExWindow(QWidget):
  def __init__(self):
    super().__init__()

    self.setWindowTitle("Ex Plot")
    # self.setGeometry(100, 100, 400, 800)
    self.resize(400, 800)

    self.ASym = ""
    self.maxEx = 0
    self.data = None
    self.Iso = IsotopeClass()

    self.save_button = QPushButton("Save Plot as Image")
    self.save_button.clicked.connect(self.save_plot)

    self.fig = go.Figure()
    self.html_file = None
    self.web_view = QWebEngineView()

    layout = QVBoxLayout(self)
    layout.addWidget(self.save_button)
    layout.addWidget(self.web_view)

    # self.plot_Ex_graph()

  def save_plot(self):
    file_path, file_type = QFileDialog.getSaveFileName(self, "Save Image", "", "PNG Files (*.png);;PDF Files (*.pdf);;All Files (*)")
    if file_path:

      if "PDF" in file_type:
        format = "pdf"
        extension = ".pdf"
      elif "PNG" in file_type:
        format = "png"
        extension = ".png"
      else:
        format = "png"
        extension = ".png"

      if not file_path.lower().endswith(extension):
        file_path += extension

      pio.write_image(self.fig, file_path, format=format)
      print(f"Plot saved to {file_path} as {format.upper()}")

  def GetEx(self, ASym :str, maxEx :float):
    self.ASym = ASym
    self.maxEx = maxEx
    self.data = self.Iso.GetExList(ASym, maxEx)    

    self.Iso.PrintIso(ASym)
    print(self.data)
    print("=============================")

  def plot_Ex_graph(self):

    if self.html_file and os.path.exists(self.html_file):
      os.remove(self.html_file)
    
    xShift = 0
    fontSize = 14
    plotHeight = 700
    plotWidth = 350
    yMin = -1
    
    A, Z = self.Iso.GetAZ(self.ASym)
    Sym = self.Iso.GetSymbol(A, Z)

    Sn = self.Iso.GetSn(self.ASym)
    Sp = self.Iso.GetSp(self.ASym)
    Sa = self.Iso.GetSa(self.ASym)

    ex=self.data['energy']/1000.
    jp=self.data['jp']
    
    # Create a Plotly figure
    self.fig.data = []

    self.fig.update_layout(plot_bgcolor='white', width=plotWidth, height = plotHeight, margin=dict(l=0, r=0, t=0, b=0))
    self.fig.update_layout(showlegend=False)
    self.fig.update_xaxes(showline=False,  visible= False, range=[-1, 2.5])
    self.fig.update_yaxes(showline=False,  visible= False, range=[yMin, self.maxEx+1])

    l=ex.last_valid_index()

    fontSizeMeV=fontSize/plotHeight*(self.maxEx+1-yMin)
    #print(fontSizeMeV)
    #adjust text label y-pos
    ypos = ex.copy()

    noOverlap = False
    loop = 0

    while noOverlap == False and loop < 2*l :
      #print("================= %d" % loop)
      for i in range(1, l+1) :
        diff = ypos[i] - ypos[i-1]
        #print("%2d | %.3f, %.3f | %.4f" % (i, ypos[i], ypos[i-1], diff))
        if diff < fontSizeMeV :
          ypos[i-1] += (diff - fontSizeMeV)/2
          ypos[i] += (fontSizeMeV - diff)/2
          if( ypos[i-1] < yMin + fontSizeMeV/2) :
            ypos[i-1] = yMin + fontSizeMeV/2
            ypos[i] = ypos[i-1] + fontSizeMeV
        #print("   | %.3f, %.3f" % (ypos[i], ypos[i-1]))

      #print(ypos)
      ###=======inspection
      count = 0
      for i in range(1, l+1) :
        diff = ypos[i] - ypos[i-1]
        if diff > fontSizeMeV :
          count = count +1

      if count == l :
        noOverlap = True

      loop += 1

    for i in range(0,l+1):
      self.fig.add_trace(go.Scatter(x=[xShift,1 + xShift], y=[ex[i],ex[i]],mode='lines',line=dict(color='black', width=1)))
      self.fig.add_trace(go.Scatter(x=[1.03 + xShift,1.1 + xShift, 1.19 + xShift], y=[ex[i],ypos[i],ypos[i]],mode='lines',line=dict(color='gray', width=1)))
      self.fig.add_annotation(x=1.2 + xShift, y=ypos[i], text=("%.3f, %s" % (ex[i], jp[i])), xanchor='left', font=dict(size=fontSize), showarrow=False)

    if( Sn < self.maxEx ):
      self.fig.add_trace(go.Scatter(x=[-0.6 + xShift,-0.1 + xShift], y=[Sn,Sn],mode='lines',line=dict(color='red', width=1)))
      self.fig.add_annotation(x=-0.6 + xShift, y=Sn, text=("Sn %.3f" % Sn), xanchor='left', yanchor='bottom', font=dict(size=fontSize, color='red'), showarrow=False)
    if( Sp < self.maxEx ):
      self.fig.add_trace(go.Scatter(x=[-0.6 + xShift,-0.1 + xShift], y=[Sp,Sp],mode='lines',line=dict(color='blue', width=1)))
      self.fig.add_annotation(x=-0.6 + xShift, y=Sp, text=("Sp %.3f" % Sp), xanchor='left', yanchor='bottom', font=dict(size=fontSize, color='blue'), showarrow=False)
    if( Sa < self.maxEx ):
      self.fig.add_trace(go.Scatter(x=[-0.6 + xShift,-0.1 + xShift], y=[Sa,Sa],mode='lines',line=dict(color='#9467bd', width=1)))
      self.fig.add_annotation(x=-0.6 + xShift, y=Sa, text=("Sa %.3f" % Sa), xanchor='left', yanchor='bottom', font=dict(size=fontSize, color='#9467bd'), showarrow=False)


    self.fig.add_annotation(x=0.5 + xShift, y=-0.6, text=("<sup>%s</sup>%s" % (A, Sym)), font=dict(size=2.5*fontSize), showarrow=False)

    # Save the plot as an HTML file in a temporary location
    timestamp = int(time.time() * 1000)  # Unique timestamp in milliseconds
    html_file = f"/tmp/Exwindow_{timestamp}.html"
    self.fig.write_html(html_file)
    self.html_file = html_file  # Store for cleanup
    self.web_view.setUrl(QUrl.fromLocalFile(html_file))

  def __del__(self):
    if self.html_file and os.path.exists(self.html_file):
       os.remove(self.html_file)