function [ results ] = plot_field_integrals_raw( y_offset_scan_Q, Bmap, relative_field_rolloff )

if(y_offset_scan_Q==0) % x offsets

    x0s = 0.999*linspace(-1,1,11) * 1e-3;
    fis = [];
    for x0 = x0s
    fi = field_integral_raw( x0, 0, Bmap, relative_field_rolloff );
    fis = [fis; fi(end,2:end)];
    end
    figure; hold on; plot(1e3*x0s, fis(:,3)); plot(1e3*x0s, fis(:,4)); plot(1e3*x0s, 1e-3 + 0*x0s); plot(1e3*x0s, -1e-3 + 0*x0s); hold off; xlabel('x offset (mm)'); ylabel('1st integral (T m)'); legend('By', 'Bx', '+limit', '-limit'); enhance_plot('arial',20,2,8);
    figure; hold on; plot(1e3*x0s, fis(:,1)); plot(1e3*x0s, fis(:,2)); plot(1e3*x0s, 1.5e-3 + 0*x0s); plot(1e3*x0s, -1.5e-3 + 0*x0s); hold off; xlabel('x offset (mm)'); ylabel('2nd integral (T m^2)'); legend('By', 'Bx', '+limit', '-limit'); enhance_plot('arial',20,2,8);

else % y offsets

    x0s = 0.999*linspace(-0.2,0.2,5) * 1e-3;
    fis = [];
    for x0 = x0s
    fi = field_integral_raw( 0, x0, Bmap, 0*relative_field_rolloff );
    fis = [fis; fi(end,2:end)];
    end
    figure; hold on; plot(1e3*x0s, fis(:,3)); plot(1e3*x0s, fis(:,4)); plot(1e3*x0s, 1e-3 + 0*x0s); plot(1e3*x0s, -1e-3 + 0*x0s); hold off; xlabel('y offset (mm)'); ylabel('1st integral (T m)'); legend('By', 'Bx', '+limit', '-limit'); enhance_plot('arial',20,2,8);
    figure; hold on; plot(1e3*x0s, fis(:,1)); plot(1e3*x0s, fis(:,2)); plot(1e3*x0s, 1.5e-3 + 0*x0s); plot(1e3*x0s, -1.5e-3 + 0*x0s); hold off; xlabel('y offset (mm)'); ylabel('2nd integral (T m^2)'); legend('By', 'Bx', '+limit', '-limit'); enhance_plot('arial',20,2,8);

end

results = [];

end