# Set the terminal for output; in this example, we'll use PNG
set terminal pngcairo enhanced font "Helvetica,16" size 1280,960
set output "makarna2.png"

# Define Fermi energy (set your value here)
Ef = -3.1840 # Replace 1.23 with your actual Fermi energy in eV


# Define range for the plot
set xrange [-0.02:3.2808]  # Adjust these values based on your high symmetry points
set yrange [-7:3]     # Adjust these based on your energy range

# Label axes
set xlabel "Wavevector"
set ylabel "Energy (eV)"

# Set style for the plot
set style line 1 lc rgb '#000000' lt 1 lw 2
set style line 2 lc rgb '#FF0000' lt 1 lw 1  # For Fermi level
set style line 3 lc rgb '#808080' lt 0 lw 1  # For vertical lines at symmetry points

# Remove legend
unset key

# Set grid for high symmetry points
set grid ytics lt 0 lw 1 lc rgb "#d3d3d3"
set grid xtics lt 0 lw 1 lc rgb "#d3d3d3"

# Define high symmetry points. Adjust these values according to your band structure:
set xtics ("Γ" 0.0000, "X" 0.5001, "M" 1.0002, "Γ" 1.7074, "R" 2.5736, "X" 3.2808)
set xtics rotate by -45 scale 1 offset 0,-0.5

# Draw vertical lines at high symmetry points
set arrow from 0.5001,-5 to 0.5001,10 nohead ls 3
set arrow from 1.0002,-5 to 1.0002,10 nohead ls 3
set arrow from 1.7074,-5 to 1.7074,10 nohead ls 3
set arrow from 2.5736,-5 to 2.5736,10 nohead ls 3

# Draw Fermi level dynamically based on Ef
set arrow from 0,Ef to 3.2808,Ef nohead ls 2

# Plot the band structure data
plot "gese.bands.dat.gnu" using 1:2 with lines ls 1

# Reset output to avoid locking the terminal
unset output