def __init__(self, folder, filen_0, dim_1, dim_2, dim_3, label='_'):
        ''' Folder where to save the data files,
        filen_0 is the base name added to each file,
        dim_1,2,3 are the sweep objects used to grab data length and names,
        label is something that can be added to the filename
        '''
        gc.collect()
        self._folder = folder
        self.filen = []
        self.head = []
        self.filen.append(filen_0 + label + '_cI1I1' + '.mtx')
        self.filen.append(filen_0 + label + '_cQ1Q1' + '.mtx')
        self.filen.append(filen_0 + label + '_cI2I2' + '.mtx')
        self.filen.append(filen_0 + label + '_cQ2Q2' + '.mtx')
        self.filen.append(filen_0 + label + '_cI1Q1' + '.mtx')
        self.filen.append(filen_0 + label + '_cI2Q2' + '.mtx')
        self.filen.append(filen_0 + label + '_cI1I2' + '.mtx')
        self.filen.append(filen_0 + label + '_cQ1Q2' + '.mtx')
        self.filen.append(filen_0 + label + '_cI1Q2' + '.mtx')
        self.filen.append(filen_0 + label + '_cQ1I2' + '.mtx')
        self.filen.append(filen_0 + label + '_SqMag' + '.mtx')
        self.filen.append(filen_0 + label + '_SqPhs' + '.mtx')
        self.filen.append(filen_0 + label + '_SqMag2' + '.mtx')
        self.filen.append(filen_0 + label + '_cPow' + '.mtx')
        for i in range(10):
            self.head.append(make_header(dim_1, dim_2, dim_3, 'cCov'))

        self.head.append(make_header(dim_1, dim_2, dim_3, 'Sq-Mag'))
        self.head.append(make_header(dim_1, dim_2, dim_3, 'Sq-Phs'))
        self.head.append(make_header(dim_1, dim_2, dim_3, 'Sq-Mag2'))
        self.head.append(make_header(dim_1, dim_2, dim_3, 'cPower'))
        self.matrix3d = np.memmap(folder+'cov.dat', dtype=np.float32, mode='w+',
                                  shape=(14, dim_3.pt, dim_2.pt, dim_1.pt))
 def __init__(self, folder, filen_0, dim_1, dim_2, dim_3, label='_'):
     self._folder = folder
     self._filen_3 = filen_0 + label + '_mag' + '.mtx'
     self._filen_4 = filen_0 + label + '_phase' + '.mtx'
     self._head_3 = make_header(dim_1, dim_2, dim_3, 'S11 _mag')
     self._head_4 = make_header(dim_1, dim_2, dim_3, 'S11 _phase')
     self._matrix3d_3 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_4 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
 def __init__(self, folder, filen_0, dim_1, dim_2, dim_3):
     self._folder = folder
     self._filen_1 = filen_0 + '_real' + '.mtx'
     self._filen_2 = filen_0 + '_imag' + '.mtx'
     self._filen_3 = filen_0 + '_mag' + '.mtx'
     self._filen_4 = filen_0 + '_phase' + '.mtx'
     self._head_1 = make_header(dim_1, dim_2, dim_3, 'S11 _real')
     self._head_2 = make_header(dim_1, dim_2, dim_3, 'S11 _imag')
     self._head_3 = make_header(dim_1, dim_2, dim_3, 'S11 _mag')
     self._head_4 = make_header(dim_1, dim_2, dim_3, 'S11 _phase')
     self._matrix3d_1 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_2 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_3 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_4 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
예제 #4
0
 def prepare_data_save(self, folder, filen_0, dim1, dim_2, dim_3):
     if self.pt > 1:
         dim_1 = self
     else:
         dim_1 = dim1
     self._folder = folder
     self._filen_1 = filen_0 + '_r' + '.mtx'
     self._filen_2 = filen_0 + '_i' + '.mtx'
     self._filen_3 = filen_0 + '_m' + '.mtx'
     self._filen_4 = filen_0 + '_p' + '.mtx'
     self._head_1 = make_header(dim_1, dim_2, dim_3, 'S11 real')
     self._head_2 = make_header(dim_1, dim_2, dim_3, 'S11 imag')
     self._head_3 = make_header(dim_1, dim_2, dim_3, 'S11 mag')
     self._head_4 = make_header(dim_1, dim_2, dim_3, 'S11 phase')
     self._matrix3d_1 = np.memmap('VNAreal.mem', dtype=np.float32, mode='w+', shape=(dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_2 = np.memmap('VNAimag.mem', dtype=np.float32, mode='w+', shape=(dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_3 = np.memmap('VNAmag.mem', dtype=np.float32, mode='w+', shape=(dim_3.pt, dim_2.pt, dim_1.pt))
     self._matrix3d_4 = np.memmap('VNAphase.mem', dtype=np.float32, mode='w+', shape=(dim_3.pt, dim_2.pt, dim_1.pt))
 def __init__(self, folder, filen_0, dim_1, dim_2, dim_3,
              label='_', cname='Voltage x1k'):
     self._folder = folder
     self._filen_1 = filen_0 + label + '.mtx'
     self._head_1 = make_header(dim_1, dim_2, dim_3, cname)
     self._matrix3d_1 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
     self.UD = False
     if dim_1.UD is True:
         self._filen_2 = filen_0 + label + '_2' + '.mtx'
         self.UD = True
         self._matrix3d_2 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
예제 #6
0
            stop=sPar[4],
            pt=sPar[8],
            scale=1)
        for sPar in d.stepItems[n]]
'''
'''
d.n2[0].lin =  (d.n2[0].lin-xoff)/x1flux + 0.5
d.n2[0].start = d.n2[0].lin[0]
d.n2[0].stop =  d.n2[0].lin[-1]
d.n2[0].name = 'Flux/Flux0'
d.n2[1].lin =  (d.n2[1].lin-xoff)/x1flux + 0.5
d.n2[1].start = d.n2[1].lin[0]
d.n2[1].stop =  d.n2[1].lin[-1]
d.n2[1].name = 'Flux/Flux0'
'''
'''
M2 = np.zeros((MAT1.shape[0], d.n2[0].pt, d.n3.pt))
M4 = np.zeros((MAT1.shape[0], d.n2[1].pt, d.n3.pt))
M2 = MAT1[:, :d.n2[0].pt, :]
M3 = MAT1[:, d.n2[0].pt-1:, :]

header1 = make_header(d.n3, d.n2[0], d.n1, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '.mtx', M2, header=header1)
header2 = make_header(d.n3, d.n2[1], d.n1, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '2' + '.mtx', M3, header=header2)
'''

header1 = make_header(d.n1, d.n2, d.n3, meas_data=('Photon Flux'))
# header1 = make_header(d.n1, d.n2, d.n3, meas_data=('Pow [W]'))
savemtx('mtx_out//' + filein + '.mtx', MAT1, header=header1)
예제 #7
0
for i, vc.snr in enumerate(np.linspace(0, 8, num0)):
    vc.make_cvals()
    BLnM[i, :, :], BLnM2[i, :, :], MIneq[i, :, :], Mn1[i, :, :], Mn2[i, :, :] = NMatrix(vc, vc2)
vc.d1.name = 'Signal/Noise'
vc.d1.pt = num0
vc.d1.start = 0.0
vc.d1.stop = 8.0

# convert to effective Phi0
vc.d3.start = vc.d3.start*0.0113
vc.d3.stop = vc.d3.stop*0.0113

ext = 'med_'

savename3 = fname+ext+'n1.mtx'
headtxt3 = make_header(vc.d3, vc.d2, vc.d1, meas_data='n1')
savemtx(savename3, Mn1, headtxt3)

savename4 = fname+ext+'n2.mtx'
headtxt4 = make_header(vc.d3, vc.d2, vc.d1, meas_data='n2')
savemtx(savename4, Mn2, headtxt4)

savename5 = fname+ext+'IneQ.mtx'
headtxt5 = make_header(vc.d3, vc.d2, vc.d1, meas_data='Ineq')
savemtx(savename5, MIneq, headtxt5)

savename = fname+ext+'LogN.mtx'
headtxt = make_header(vc.d3, vc.d2, vc.d1, meas_data='Log-Negativity')
savemtx(savename, BLnM, headtxt)

savename2 = fname+ext+'ineqSq.mtx'
예제 #8
0
            stop=sPar[4],
            pt=sPar[8],
            scale=1)
        for sPar in d.stepItems[n]]
'''
'''
d.n2[0].lin =  (d.n2[0].lin-xoff)/x1flux + 0.5
d.n2[0].start = d.n2[0].lin[0]
d.n2[0].stop =  d.n2[0].lin[-1]
d.n2[0].name = 'Flux/Flux0'
d.n2[1].lin =  (d.n2[1].lin-xoff)/x1flux + 0.5
d.n2[1].start = d.n2[1].lin[0]
d.n2[1].stop =  d.n2[1].lin[-1]
d.n2[1].name = 'Flux/Flux0'
'''
'''
M2 = np.zeros((MAT1.shape[0], d.n2[0].pt, d.n3.pt))
M4 = np.zeros((MAT1.shape[0], d.n2[1].pt, d.n3.pt))
M2 = MAT1[:, :d.n2[0].pt, :]
M3 = MAT1[:, d.n2[0].pt-1:, :]

header1 = make_header(d.n3, d.n2[0], d.n1, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '.mtx', M2, header=header1)
header2 = make_header(d.n3, d.n2[1], d.n1, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '2' + '.mtx', M3, header=header2)
'''

# header1 = make_header(d.n1, d.n2, d.n3, meas_data=('Photon Flux'))
header1 = make_header(d.n1, d.n2, d.n3, meas_data=('Pow [W]'))
savemtx('mtx_out//' + filein + '.mtx', MAT1, header=header1)
'''
#Other Equipment
execfile('keithley2000.py')
vm = instrument3('GPIB0::29::INSTR')
vm.optimise()
sleep(0.1)
vm.testspeed()
'''

#for the VNA I want 4 data files for Real, Imag, MAG, Phase
filen_1 = filen_0 + '_real'  + '.mtx'
filen_2 = filen_0 + '_imag'  + '.mtx'
filen_3 = filen_0 + '_mag'   + '.mtx'
filen_4 = filen_0 + '_phase' + '.mtx'
head_1 = make_header(dim_1, dim_2, dim_3, 'S21 _real')
head_2 = make_header(dim_1, dim_2, dim_3, 'S21 _imag')
head_3 = make_header(dim_1, dim_2, dim_3, 'S21 _mag')
head_4 = make_header(dim_1, dim_2, dim_3, 'S21 _phase')
matrix3d_1 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
matrix3d_2 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
matrix3d_3 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
matrix3d_4 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))

ask_overwrite(folder+filen_1)
copy_file(thisfile, filen_0, folder) #backup this script

dim_1.lin = np.linspace(dim_1.start,dim_1.stop,dim_1.pt)
dim_2.lin = np.linspace(dim_2.start,dim_2.stop,dim_2.pt)
dim_3.lin = np.linspace(dim_3.start,dim_3.stop,dim_3.pt)
예제 #10
0
MAT1[1] = d.data[:, 5, :]
MAT1[2], MAT1[3] = get_MP(d, 6)
MAT1[4], MAT1[5] = get_MP(d, 8)
MAT1[6], MAT1[7] = get_MP(d, 10)
MAT1[8], MAT1[9] = get_MP(d, 12)

# MAT2[-1] = d.data[:, -1, :]

M2 = np.zeros((MAT1.shape[0], d.n2[0].pt, d.n3.pt))
M3 = np.zeros((MAT1.shape[0], d.n2[1].pt, d.n3.pt))
M3 = MAT1[:, :d.n2[0].pt, :]
M2 = MAT1[:, d.n2[0].pt-1:, :]

# meas specific to change mag field to flux
# simply comment this paragraph out
xoff = 140.5e-3  # 139.3e-3
x1flux = 479.6e-3
d.n2[0].lin =  (d.n2[0].lin-xoff)/x1flux + 0.5
d.n2[0].start = d.n2[0].lin[0]
d.n2[0].stop =  d.n2[0].lin[-1]
d.n2[0].name = 'Flux/Flux0'
d.n2[1].lin =  (d.n2[1].lin-xoff)/x1flux + 0.5
d.n2[1].start = d.n2[1].lin[0]
d.n2[1].stop =  d.n2[1].lin[-1]
d.n2[1].name = 'Flux/Flux0'

header2 = make_header(d.n3, d.n2[0], d.n1, meas_data=('a.u.'))
header1 = make_header(d.n3, d.n2[1], d.n1, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '.mtx', M3, header=header1)
savemtx('mtx_out//' + filein + '2' + '.mtx', M2, header=header2)
예제 #11
0
    MAT1[n*4+3] = np.angle(compx)

# store Power channels
MAT1[(numOfDigi)*4] = d.data[:, 5, :]
MAT1[(numOfDigi)*4+1] = d.data[:, 8, :]
MAT1[(numOfDigi)*4+2] = d.data[:, -1, :]  # Store measured Voltage

M2 = np.zeros((11, 1501, 70))
M3 = np.zeros((11, 1501, 70))

M3 = MAT1[:, :1501, :]
M2 = MAT1[:, 1500:, :]
d.n2.pt = (d.n2.pt-1)/2 +1

d.n2.stop = d.n2.stop*-1
d.n2.update_lin()
header1 = make_header(d.n1, d.n2, d.n2, meas_data=('a.u.'))
savemtx('mtx_out//' + filein + '.mtx', M3, header=header1)

#flip endpoints for the backsweep
d.n2.stop = d.n2.stop*-1
d.n2.start = d.n2.start*-1
d.n2.update_lin()
header1 = make_header(d.n1, d.n2, d.n2, meas_data=('a.u.'))
savemtx('mtx_out//' + filein +'2'+ '.mtx', M2, header=header1)

# for jj, filext in enumerate(d.channel[5:]):
#     MAT2[jj+2] = d.data[:, jj+5, :]
#     header1 = make_header(d.n1, d.n2, d.n2, meas_data=filext[1])
#     # savemtx(folder + filein + filext[0] + '.mtx', MAT1[jj+2], header=header1)
예제 #12
0
d.n1.start = d.n1.lin[0]
d.n1.stop = d.n1.lin[-1]
d.n1.name = 'Flux/Flux0'

MAT1 = np.zeros([8, d.shape[0], d.shape[1]])
MAT1[0] = d.D1mag
MAT1[1] = d.D2mag
MAT1[2] = d.D1phase
MAT1[3] = d.D2phase
MAT1[4] = d.D1Pow
MAT1[5] = d.D2Pow
MAT1[6] = d.D1Pow2
MAT1[7] = d.D2Pow2


header1 = make_header(d.n1, d.n2, d.n2, meas_data='est. Photon #')
savemtx('output//S1_480_V.mtx', MAT1, header=header1)

y = (d.n2.lin*d.n2.lin/50.0)  # position of the data

MAT2 = np.zeros([2, d.shape[0]*10, d.shape[1]])
y2, MAT2[0] = interp_y(y, d.D1Pow2, factor=10)
y2, MAT2[1] = interp_y(y, d.D2Pow2, factor=10)

d.dim_y2 = d.n2
d.dim_y2.start = y2[0]*1e3
d.dim_y2.stop = y2[-1]*1e3
d.dim_y2.pt = len(y2)
d.dim_y2.lin = y2*1e3
d.dim_y2.name = 'Pump power (mW)'
예제 #13
0
           start = -1,
           stop = 1,
           pt = 1001,
           scale = flux0)
freq = dim(name = 'Frequency (GHz)',
           start = 4,
           stop = 8,
           pt = 101,
           scale = 1e9)
dim_3 = handler(name = 'mag/phase',
           start = 0,
           stop = 10,
           pt = 11) #8 pts for S 4x2 values
dim_3._Z0 = 50

head1 = make_header(magnet, freq, dim_3, 'S11 S12 S21 S22 Z L')
dim_3.prepare_data_save(magnet, freq, dim_3)


for jj, f0 in enumerate(freq.lin):
    for ii, flux in enumerate(magnet.lin):
        # b = k = 2pi/wavelength; wavelength = velocity / frequency
        b = 2.0*pi*f0/v
        L = flux0 / (Ic*2.0*pi* abs(cos(pi*flux/flux0)))
        Ysq = (1.0/R + 1.0/(i*2.0*pi*f0*L +i*1e-90) + i*2.0*pi*f0*Cap)
        Zsq = 1.0/Ysq

        ABCD_Matrix = tline(70,b,0.01)*tline(50,b,0.3)*tline(10,b,900e-6)*sres(Zsq)*shunt(0.1)

        #record stuff into dim_3._SMat
        dim_3.record_SM(ABCD_Matrix,jj,ii)
예제 #14
0
 def prepare_data_save(self, folder, filen_0, dim_1, dim_2, dim_3, colour_name):
     '''Prepare an empty matrix which will be filled with measurement data'''
     self._folder = folder
     self._filen_1 = filen_0 + '.mtx'
     self._head_1 = make_header(dim_1, dim_2, dim_3, colour_name)
     self.matrix3d_1 = np.zeros((dim_3.pt, dim_2.pt, dim_1.pt))
d.shape = d.data.shape

xoff = 139.3e-3
x1flux = 479.6e-3
d.dim_1.lin = (d.dim_1.lin-xoff)/x1flux + 0.5
d.dim_1.start = d.dim_1.lin[0]
d.dim_1.stop = d.dim_1.lin[-1]

MAT1 = np.zeros([6, d.shape[0], d.shape[2]])
MAT1[0] = d.data[:, 3, :]
MAT1[1] = d.data[:, 4, :]
MAT1[2] = d.D1mag
MAT1[3] = d.D1ang
MAT1[4] = d.D1pow
MAT1[5] = d.D1lev

y = (d.dim_2.lin*d.dim_2.lin/50.0)  # position of the data
y2, MAT2 = interp_y(y, d.D1pow)
d.dim_y2 = d.dim_2
d.dim_y2.start = y2[0]
d.dim_y2.stop = y2[-1]
d.dim_y2.pt = len(y2)
d.dim_y2.lin = y2


header1 = make_header(d.dim_1, d.dim_2, d.dim_3, meas_data='(a.u)')
savemtx('S1_471_DCE_MAP2_volt.mtx', MAT1, header=header1)

header2 = make_header(d.dim_1, d.dim_y2, d.dim_3, meas_data='(a.u)')
savemtx('S1_471_DCE_MAP2_wats.mtx', MAT2, header=header2)