Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download
Views: 143
Kernel: Python 3 (Anaconda)
import numpy as np import matplotlib.pyplot as plt from glob import glob
nPts = [] averages = [] stds = [] distances = [] cpm = [] ucpm = [] files = glob("*.csv") for dataFile in files: data = np.genfromtxt(dataFile) nPts.append(len(data)) avg = np.average(data) averages.append(avg) stds.append(np.std(data)) distances.append(dataFile) avgerages = np.array(avg) stds = np.array(stds) unAvgs = stds/(np.sqrt(nPts)) v = [] for dist in distances: v.append(float(dist.split('.')[0])) cpm = 60.0/np.array(averages) uncpm = (np.array(unAvgs)/np.array(averages))*np.array(cpm)
figure = plt.figure(figsize=(15,5)) ax = figure.add_axes([0.1, 0.1, 0.8, 0.8]) plt.errorbar(v, cpm, yerr=uncpm, fmt = '.') plt.title("Simplest er, 0.2 in x, 0.4 in y") ax.set_xlabel("Wait Time (s)") ax.set_ylabel("Counts (N)")
Text(0,0.5,'Counts (N)')
Image in a Jupyter notebook
r = np.sort(np.array(v)) def model(params,r): i0 = params[0] return (i0/r**2)+ 3.9172 params = np.array((6500000.0,)) rAxis = np.linspace(27,160,200)
figure = plt.figure(figsize=(15,5)) ax = figure.add_axes([0.1, 0.1, 0.8, 0.8]) plt.errorbar(v, cpm, yerr=uncpm, fmt = '.') plt.title("Simplest er, 0.2 in x, 0.4 in y") ax.set_xlabel("Distance (mm)") ax.set_ylabel("Intensity (counts/min)") ax.plot(rAxis, model(params, rAxis), 'r-')
[<matplotlib.lines.Line2D at 0x7fcbc1d88a90>]
Image in a Jupyter notebook
r = np.sort(r) print(r)