QPSK demod constellation problem

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/jpus7fn5fdze1.png?width=1178&format=png&auto=webp&s=c98fa280993e0aa716bc234e5103581eb0b4bdd7

thanks