Plotting Zero-Crossing Roadmaps
Posted: Wed Dec 03, 2014 1:43 pm
I'm attempting to plot a simple "roadmap" of zero-crossings for an array of angles based on different rotation axes.
When I use a rotation axis that has the form, [1 0 x ], x/=0 I am only able to plot data for 45 degrees. When I plot [1 0 0] or [0 1 0] or [1 x 0], the code works fine, but for some reason on other axes, some spectra is "cut off".
I believe this may be a bug. This is the code that I am using.
When I use a rotation axis that has the form, [1 0 x ], x/=0 I am only able to plot data for 45 degrees. When I plot [1 0 0] or [0 1 0] or [1 x 0], the code works fine, but for some reason on other axes, some spectra is "cut off".
I believe this may be a bug. This is the code that I am using.
Code: Select all
% resonance roadmap from single crystal rotation
%================================================================
clear; clc; clf;
%Spin parameters
Sys.S = 5/2;
Sys.g = [2.0069];
Sys.lw = 3;
Sys.aFrame = 3;
Sys.D = [-768*3];
Sys.aF = [78*3 15*3];
% Experimental parameters
Exp.mwFreq = 9.401;
Exp.Range = [150 700];
Exp.CrystalSymmetry = 186;
Exp.nPoints = 1024;
% Generate orientations in a single rotation plane
rotN = [1 0 1]; % rotation axis
[phi,theta] = rotplane(rotN,[0 pi/2],91);
Exp.Orientations = [phi;theta];
% Simulate spectra
Opt.Output = 'separate'; % make sure spectra are not added up
Bres = resfields(Sys,Exp,Opt);
% plotting
plot(Bres,theta*180/pi,'b')
xlabel('magnetic field (mT)');
ylabel('theta (°)');