fig, ax = plt.subplots(1, 1, figsize=(10, 7))
x = df.lz_ref
y = df.sqrtjR_ref
cond = (x > 500) & (y > 6)
x = x[cond]
y = y[cond]
ax.scatter(x, y, s=1, alpha=0.1)