Overlapping Distributions

import matplotlib.pyplot as plt
import numpy as np

COLORS = {'group_a': '#6366F1', 'group_b': '#10B981', 'background': '#FFFFFF', 'text': '#1E293B', 'text_muted': '#64748B', 'grid': '#F1F5F9'}

np.random.seed(42)
group_a = np.random.normal(45, 10, 500)
group_b = np.random.normal(55, 12, 500)

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor(COLORS['background'])
fig.patch.set_facecolor(COLORS['background'])

ax.hist(group_a, bins=25, color=COLORS['group_a'], edgecolor='white', linewidth=0.5, alpha=0.6, label='Group A')
ax.hist(group_b, bins=25, color=COLORS['group_b'], edgecolor='white', linewidth=0.5, alpha=0.6, label='Group B')

ax.axvline(np.mean(group_a), color=COLORS['group_a'], linewidth=2, linestyle='--')
ax.axvline(np.mean(group_b), color=COLORS['group_b'], linewidth=2, linestyle='--')

ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color(COLORS['grid'])
ax.spines['bottom'].set_color(COLORS['grid'])
ax.yaxis.grid(True, color=COLORS['grid'], linewidth=1, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors=COLORS['text_muted'], labelsize=9, length=0, pad=8)
ax.set_xlabel('Score', fontsize=10, color=COLORS['text'], labelpad=10)
ax.set_ylabel('Count', fontsize=10, color=COLORS['text'], labelpad=10)
ax.legend(loc='upper center', bbox_to_anchor=(0.5, -0.12), ncol=2, frameon=False, fontsize=10, labelcolor=COLORS['text_muted'])

plt.tight_layout()
plt.show()