###### John's Result Evaluation #######
log.info("Evaluation of QMix results.")
z = 6.45 - 4.21j
y = np.reciprocal(z)
#iLO_fr_cir = M.iLO_from_circuit_calc(y)
#iLO = M.masking(iLO_fr_cir,M.mask_photon_steps)
#iLO = np.average(iLO[1])
iLO = M.current_LO_from_Embedding_Circuit(
M.total_admittance([y.real, y.imag], M.ySIS_masked_positive),
M.pumping_Levels_Volt_masked_positive)
#TODO put the following three equaitons in an own function
vLO = M.vLO_from_circuit_calc(iLO=iLO,
yEmb=y,
unpumpedExpanded=M.Unpumped.binedDataExpanded,
iKKExpanded=M.Unpumped.iKKExpanded,
vrangeEvaluated=M.Pumped.binedIVData[0])
alp = M.pumping_Levels_calc(vLO)
pumprec = M.pumped_from_unpumped_calc(
alphas=alp, unpumpedExpanded=M.Unpumped.binedDataExpanded)
title = newfig('Comparison_with_Johns_Embedding_Impedance')
plot(M.Pumped.binedIVData, label='Binned Data')
plot(pumprec, label='QMix Fitting Result')
plot(M.pumped_from_embedding_circuit, label='Own Fitting Result')
pltsettings(directory + title,
xlabel=lbl['mV'],
ylabel=lbl['uA'],
