Confused Constellation on QPSK

I'm trying to QPSK demod a signal, i used a Frequency offsets are correction using Fourier transforms.
rrc_iq=apply_rrc_filter(iq_corrected, sample_rate, symbol_rate)

RRC filter with 0.25 rolloff
_, rrc = rrcosfilter(rrcos_length, 0.25, 1/(kbaud * resampling_factor), kbaud * resampling_factor * samples_per_symbol)

and finally
PLL
initial_bandwidth = 400.0 # Hz (fast acquisition)

damping_factor = 0.5 # Critically damped

iq_pll, phase_estimates = costas_loop1(rrc_iq, initial_bandwidth, damping_factor, sample_rate, symbol_rate)

But this is the result!
any suggestion? where I'm wrong?

https://preview.redd.it/fz4s7o4gfdze1.png?width=1178&format=png&auto=webp&s=7d56575ae3670d6a0fe32e4343d48dad110d6e6a

Thanks