I quickly ran your script. A few more thoughts:
- The range is in kilometers and must therefore be multiplied by 1000 before being passed to Orekit. And I recommend you define a separate sigma value for the range
- But even after this fix, the estimator still does not converge
- When mixing different sort of measurements, it is very important to properly tune the weight and sigma of the different types of measurements
- You could also feed the Doppler data as an Orekit
RangeRate
measurement: RangeRate (ORbit Extrapolation KIT 11.1 API) - Your tracking data is a single pass whose duration is slightly more than 2 minutes. I am not an expert in orbit determination but this could prove challenging.