In [1]:
%matplotlib notebook

import matplotlib.pyplot as plt
import numpy as np
In [2]:
npz = np.load('ptopfit.npz')
dv = npz['dv']
dz = npz['dz']

def show(title, x, y):
    plt.figure(figsize=(9,8))
    plt.title(title)
    plt.scatter(x, y)
In [3]:
show('Original data', dv, dz)
In [4]:
show('Original data, rotated by $45^\circ$ around (1,1), default scale',
     ((dz-1)+(dv-1))/np.sqrt(2)+1, ((dz-1)-(dv-1))/np.sqrt(2)+1)
In [5]:
show('Original data, rotated by $45^\circ$ around (1,1), full scale for x-axis and y-axis',
     ((dz-1)+(dv-1))/np.sqrt(2)+1, ((dz-1)-(dv-1))/np.sqrt(2)+1)
plt.xlim(0.95, 1.05)
plt.ylim(0.99, 1.01);
In [6]:
show('Original data, rotated by $45^\circ$ around (1,1), identical scale for x-axis and y-axis',
     ((dz-1)+(dv-1))/np.sqrt(2)+1, ((dz-1)-(dv-1))/np.sqrt(2)+1)
plt.xlim(0.95, 1.05)
plt.ylim(0.95, 1.05);