%%%% ACh and OXT speed scatterplot
clear all;clc
filename = ['F:\oxytocin\ACh_OXT004\ACh_OXT004_015_230411_150537',
            'F:\oxytocin\ACh_OXT004\ACh_OXT004_016_230412_145104',
            'F:\oxytocin\ACh_OXT004\ACh_OXT004_017_230412_151435',
            'F:\oxytocin\ACh_OXT005\ACh_OXT005_012_230411_153514',
            'F:\oxytocin\ACh_OXT005\ACh_OXT005_013_230412_154640',
            'F:\oxytocin\ACh_OXT005\ACh_OXT005_014_230412_160131',
            'F:\oxytocin\ACh_OXT006\ACh_OXT006_009_230411_155414',
            'F:\oxytocin\ACh_OXT006\ACh_OXT006_010_230412_163714',
            'F:\oxytocin\ACh_OXT006\ACh_OXT006_011_230412_164915',
            'F:\oxytocin\ACh_OXT007\ACh_OXT007_010_230411_161405',
            'F:\oxytocin\ACh_OXT007\ACh_OXT007_011_230412_170752',
            'F:\oxytocin\ACh_OXT007\ACh_OXT007_012_230412_172120'];

for ii = 1:size(filename,1)
    cd(filename(ii,:))
    basepath = pwd;
    basename = bz_BasenameFromBasepath(pwd);

    FluoData = Get_OxtAch_Fluos();
    spd = getSpd;
    figure;plot(FluoData.time,FluoData.Ach_dFF);
    hold on; plot(spd.time,spd.data/20); hold off;
    figure;plot(FluoData.time,FluoData.Oxt_dFF);  %%%% check the dFF time_0
    hold on; plot(spd.time,spd.data/20); hold off;

    

end
