Particle Energy Levels

import matplotlib.pyplot as plt
import numpy as np

np.random.seed(42)
particles = ['Electron', 'Muon', 'Tau', 'Proton', 'Neutron']
data = [
    np.random.gamma(2, 0.5, 180) * 100,
    np.random.gamma(3, 1, 180) * 100,
    np.random.gamma(4, 1.5, 180) * 100,
    np.random.gamma(5, 2, 180) * 100,
    np.random.gamma(5.5, 2.2, 180) * 100
]

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor('#0a0a0f')
fig.patch.set_facecolor('#0a0a0f')

colors = ['#5314E6', '#4927F5', '#276CF5', '#27D3F5', '#27F5B0']
bp = ax.boxplot(data, widths=0.5, patch_artist=True, showfliers=True, notch=True,
                flierprops=dict(marker='.', markerfacecolor='#F5D327', markersize=5, alpha=0.4),
                medianprops=dict(color='#F5D327', linewidth=2))

for patch, color in zip(bp['boxes'], colors):
    patch.set_facecolor(color)
    patch.set_alpha(0.8)
    patch.set_edgecolor('white')
for i, color in enumerate(colors):
    bp['whiskers'][i*2].set_color(color)
    bp['whiskers'][i*2+1].set_color(color)
    bp['caps'][i*2].set_color(color)
    bp['caps'][i*2+1].set_color(color)

ax.set_xticklabels(particles)
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color('#333333')
ax.spines['bottom'].set_color('#333333')
ax.yaxis.grid(True, color='#1a1a3f', linewidth=0.5, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors='#888888', labelsize=9, length=0, pad=8)
ax.set_ylabel('Energy (GeV)', fontsize=11, color='white', fontweight='500')
ax.set_title('Particle Energy Level Distribution', fontsize=14, color='white', fontweight='bold', pad=15)

plt.tight_layout()
plt.show()