In [1]:

    

import sys
sys.path.append('/mnt/labshare/Programs/python/HelperFunctions/')
from myfunctions import *

from os import listdir, path

import numpy as np


    

In [2]:

    

sample = 'CH010'
datadir = '/mnt/labshare/Experiments/JPM/CH010/'
listdir(datadir)


    

    
Out[2]:





['Flux Tuning of CH010.eps',
 'fluxat-55dBm',
 'CH010 Bias 2.eps',
 'flux2at-50dBm',
 'dispshift',
 'Flux tuning of CH010 - Bias 2.eps',
 '._.DS_Store',
 '.DS_Store',
 'CH010 Bias 1.eps']

In [3]:

    

dsdir = path.join(datadir, 'dispshift')
plotDispersiveShift(dsdir)


    

In [4]:

    

fdp = path.join(datadir, 'flux2at-50dBm')
plotFluxDependence(fdp);

fluxes = []
data = []

files = sorted([f for f in listdir(fdp)], key=lambda x: float(x.split('.s2p')[0]))
for f in files:
    fluxes.append(float(f.split('.s2p')[0]))
    fig, x, y = plotS21(path.join(fdp, f), show=False)
    data.append(10**(y/20))
    
layout, trace = HeatmapPlot(title=sample + ' Bias 2')

trace['z'] = data
trace['y'] = fluxes
trace['x'] = x

iplot(Figure(data=[trace], layout=layout))


    

In [39]:

    

fdp = path.join(datadir, 'fluxat-55dBm')
plotFluxDependence(fdp);

fluxes = []
data = []

files = sorted([f for f in listdir(fdp)], key=lambda x: float(x.split('.s2p')[0]))
for f in files:
    fluxes.append(float(f.split('.s2p')[0]))
    fig, x, y = plotS21(path.join(fdp, f), show=False)
    data.append(10**(y/20))
    
layout, trace = HeatmapPlot(title=sample)

trace['z'] = data


iplot(Figure(data=[trace], layout=layout))


    

    




---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-39-d7771bc93d3e> in <module>()
      1 fdp = path.join(datadir, 'fluxat-55dBm')
----> 2 plotFluxDependence(fdp);
      3 
      4 fluxes = []
      5 data = []

/mnt/labshare/Programs/python/HelperFunctions/myfunctions.py in plotFluxDependence(datapath, samplename, param)
    143             volt = float(path.splitext(f)[0])
    144             volts.append(volt)
--> 145             result = analyze(path.join(datapath,f), show=False)
    146             freqs.append(result[param])
    147     if param == "f0":

/mnt/labshare/Programs/python/HelperFunctions/myfunctions.py in analyze(filename, show, x, y)
     41     #     Calculate fit error
     42 
---> 43     perr = sqrt(diag(pcov))
     44 
     45     f0err = perr[1]*1e-9

/usr/lib/python3/dist-packages/numpy/lib/twodim_base.py in diag(v, k)
    287         return v.diagonal(k)
    288     else:
--> 289         raise ValueError("Input must be 1- or 2-d.")
    290 
    291 def diagflat(v, k=0):

ValueError: Input must be 1- or 2-d.

In [17]:

    

tdp = '/mnt/labshare/Experiments/JPM/Test Chips/T1/'
plotS21(path.join(tdp, 'T1 Spectrum'));
result = analyze(path.join(tdp, 'T1 Resonator'));


    

    












    












    




f0 = 4.704378+/- 0.00001 (0.0003%) GHz
kappa = 21147417.24560203+/- 186036.05489 (0.88%)
Q = 35.405+/- 0.012 (0.034%)






    

In [23]:

    

result['kappa']/(2*pi*1e6)


    

    
Out[23]:





3.3657159882643573

In [ ]: