# You have to provide a DPL94 results state file
res_folder = "resultsR1"
res_folder = "results_clustering"
res_state = os.path.join(res_folder, "DPL94", "results")
spin_inte = ":52@N"
# Calculate the offset data
theta_spin_dic, Domega_spin_dic, w_eff_spin_dic, dic_key_list = calc_rotating_frame_params(spin=cdp.myspin, spin_id=spin_inte, verbosity=01)
# Save the data in cdp to access it after execution of script.
cdp.myspin.theta_spin_dic = theta_spin_dic
w_eff = w_eff_spin_dic[param_key]
cdp.mydic[exp_type][frq][offset]['w_eff'].append(w_eff)
# Average resonance spin_lock_offset
#print Domega_spin_dic[param_key]
# Y val
cdp.mydic[exp_type][frq][offset]['R1_rho_R2eff_inter'].append(R1_rho_R2eff)
#if oi == 0: #print exp_type, frq, offset, point, theta, w_eff
####### PLOT ####