###### 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'],