/usr/share/psychtoolbox-3/PsychTests/MelanopsinFundamentalTest.m is in psychtoolbox-3-common 3.0.11.20140816.dfsg1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | % MelanopsinFundamentalTest
%
% Compare melanopsin fundamental generated by PTB's DefaultPhotoreceptors/FillInPhotoreceptors
% machinery to the tabulated function provided by Lucas.
%
% 11/11/13 ms Wrote it with DHB.
% 11/15/13 dhb Comment and clean up for insertion into PTB.
%% Clear
clear; close all;
%% Get the rod fundamental using PTB machinery
photoreceptors = DefaultPhotoreceptors('LivingHumanMelanopsin');
photoreceptors = FillInPhotoreceptors(photoreceptors);
%% Load in Lucas' melanopic efficiency function
load('T_melanopsin');
%% Make a plot
figure; clf;
set(gcf,'Position',[100 100 800 400]);
subplot(1,2,1); hold on
plot(SToWls(S_melanopsin), T_melanopsin', '-k', 'LineWidth', 2); hold on;
plot(SToWls(photoreceptors.nomogram.S),photoreceptors.energyFundamentals','r','LineWidth',1); hold on;
xlabel('Wavelength (nm)');
ylabel('Energy Sensitivity');
legend({'CIE' 'PTB'});
pbaspect([1 1 1]);
subplot(1,2,2); hold on
plot(SToWls(S_melanopsin), log10(T_melanopsin'), '-k', 'LineWidth', 2); hold on;
plot(SToWls(photoreceptors.nomogram.S),log10(photoreceptors.energyFundamentals'),'r','LineWidth',1); hold on;
xlabel('Wavelength (nm)');
ylabel('Log 10 Energy Sensitivity');
legend({'CIE' 'PTB'});
pbaspect([1 1 1]);
set(gcf, 'Color', [1 1 1]);
set(gcf, 'InvertHardCopy', 'off');
set(gcf, 'PaperPosition', [0 0 7 7]); %Position plot at left hand corner with width 5 and height 5.
set(gcf, 'PaperSize', [7 7]); %Set the paper to have width 5 and height 5.
|