%%% Script to look at a set of regularized P(r) transforms %%%
%%% Set the filename of the *.param file to be read below  %%%
%%%
%%% Author: Jan Lipfert
%%% Date: 11/1/2007
clf;clc;clear;
hold on; box on;

%%% --- Set the filename of the *.param file here --%%%
par = load('filename.param');  

%%% --- For comparison, enter the Rg obtained from Guinier analysis --- %%%
expRg = 13.7; % Rg from Guinier analysis in Angstroem


subplot(3,2,1)
hold on; box on;
title('Discrepancy (Ideal 0.7)', 'fontsize', 16)
plot(par(:,1), par(:,2))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);

subplot(3,2,2)
hold on; box on;
title('Oscillations (Ideal 1.1)', 'fontsize', 16)
plot(par(:,1), par(:,3))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);

subplot(3,2,3)
hold on; box on;
title('Stability (Ideal 0.0)', 'fontsize', 16)
plot(par(:,1), par(:,4))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);

subplot(3,2,4)
hold on; box on;
title('Systematic deviations (Ideal 1.0)', 'fontsize', 16)
plot(par(:,1), par(:,5))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);

subplot(3,2,5)
hold on; box on;
title('Positivity (Ideal 1.0)', 'fontsize', 16)
plot(par(:,1), par(:,6))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);

subplot(3,2,6)
hold on; box on;
title('Real space Rg', 'fontsize', 16)
plot(par(:,1), par(:,8))
set(gca, 'fontsize', 16)
xlabel('Dmax', 'fontsize', 16);
line([par(1,1) par(size(par,1),1)],[ expRg expRg])