def dunlop(palint_object, ax,
component='m', norm_factor=1,
plt_idx=0,
out='show', folder=None, name='output.pdf',
plt_opt={},
**options):
components = {'x': 1, 'y': 2, 'z': 3, 'm': 4}
idx = components[component]
colors = helper.get_colors()
marker = helper.get_marker()
ax.plot(palint_object.th[:, 0], palint_object.th[:, idx] / norm_factor, '.-',
marker=marker[plt_idx],
color=colors[0], **plt_opt)
ax.plot(palint_object.sum[:, 0], palint_object.sum[:, idx] / norm_factor, '.-',
marker=marker[plt_idx],
color=colors[1], **plt_opt)
ax.plot(palint_object.ptrm[:, 0], palint_object.ptrm[:, idx] / norm_factor, '.-',
marker=marker[plt_idx],
color=colors[2], **plt_opt)
# todo write dunlop_add_tr()
# ax.plot(palint_object.tr[:, 0], palint_object.tr[:, idx] / norm_factor,
# color='k', marker=marker[plt_idx], alpha=0.8, ls='')
dunlop_std(palint_object=palint_object, ax=ax, component=component, norm_factor=norm_factor, plt_idx=plt_idx,
plt_opt=plt_opt)
return ax
def arai(palint_object, ax,
component='m', norm=None, norm_factor=[1, 1],
plt_idx=0,
t_min=20, t_max=700,
out='show', folder=None, name='output.pdf',
line=True, check=True,
plt_opt={}, **options):
label = options.get('label', '')
idx = palint_object.components[component]
markers = helper.get_marker()
colors = ['r', 'b', 'g'] # helper.get_colors()
xy = np.array(
[[i[idx], j[idx]] for i in palint_object.ptrm for j in palint_object.th if i[0] == j[0]]) / norm_factor
th = xy[:, 1]
ptrm = xy[:, 0]
if not 'ls' in plt_opt or 'linestyle' in plt_opt:
plt_opt.update({'ls': ''})
ax.plot(np.fabs(ptrm), np.fabs(th),
marker='o', color=colors[plt_idx], label=label, **plt_opt)
area = options.get('area')
temps = options.get('temps')
if area:
slopes, sigmas, y_intercept, x_intercept = palint_object.calculate_slope(b_min=t_min, t_max=t_max)
y_fit = (y_intercept[idx - 1] + slopes[idx - 1] * ptrm * norm_factor[0]) / norm_factor[1]
ax.fill_between(np.fabs(ptrm), np.fabs(th), y_fit, alpha=0.2, color=colors[plt_idx])
if line:
add_arai_line(palint_object=palint_object, ax=ax,
t_min=t_min, t_max=t_max,
plt_idx=plt_idx,
norm_factor=norm_factor,
component=component,
plt_opt=plt_opt)
if check:
add_ck_check(palint_object=palint_object, ax=ax,
plt_idx=plt_idx,
norm_factor=norm_factor,
component=component,
plt_opt=plt_opt)
add_ac_check(palint_object=palint_object, ax=ax,
plt_idx=plt_idx,
norm_factor=norm_factor,
component=component,
plt_opt=plt_opt)
if temps:
add_arai_temps(palint_object=palint_object, ax=ax,
plt_idx=plt_idx,
norm_factor=norm_factor,
component=component,
plt_opt=plt_opt)
return ax
def arai_stdev(palint_object, ax,
component='m', norm=None, norm_factor=[1, 1],
plt_idx=0,
t_min=20, t_max=700,
line=True, check=True,
plt_opt={}, **options):
idx = palint_object.components[component] - 1
markers = helper.get_marker()
# colors = helper.get_colors()
colors = ['r', 'b', 'g'] # helper.get_colors()
y, x = palint_object._get_th_ptrm_data(t_min=t_min, t_max=t_max, m=True)
x /= norm_factor[0]
y /= norm_factor[1]
y_stdev, x_stdev = palint_object._get_th_ptrm_stdev_data(t_min=t_min, t_max=t_max)
ax.fill_between(x[:, idx], y[:, idx] + y_stdev[:, idx], y[:, idx] - y_stdev[:, idx], color=colors[plt_idx],
alpha=0.2)
return ax
