function plotTodHistogramPanel(data)

    width = 100;
    height = 76;
    %outDir = "D:\Harvard\Manuscripts\Summit paper\Figures\Figure 1\";
    outDir = "D:\Harvard\Manuscripts\Summit paper\Figures\Fig 1 Supp\";

    numCad = length(data.ZTlastMvmnts);
    if numCad >1
        figure('Units','Pixels','Position',[0 0 (width*5)+200 (height*5)+200]);
        axesH = axes;
        set(axesH,'Units','pixels','Position',[50,50,width*5,height*5]);
        hold on;
        
        plotTimes = mod(data.ZTlastMvmnts,24);
        
        areaDay = 0:12;
        areaNight = 12:24;
        areaY = ones(1,length(areaNight));
        area(areaNight,areaY,'FaceColor','k','FaceAlpha',0.3,'HandleVisibility','off');   
        area(areaDay,areaY,'FaceColor','y','FaceAlpha',0.3,'HandleVisibility','off'); 
        h = histogram(plotTimes,48,'Normalization','probability','BinLimits',[0 24],'FaceColor',[1 0 0]);
        %'FaceColor',[0.93 0.13 0.14]
        
        %title(['Times of last movements',data.graphTitle],'FontSize',data.graph.titleFont);
        xlabel('Time (ZT hours)','FontSize',data.graph.labelFont);
        ylabel('Probability','FontSize',data.graph.labelFont);
        set(gca,'FontSize',data.graph.labelFont,'TickLength',[0 0],'xtick',0:4:24);
        
        xlim([0 24]);
        %ylim([0 200]); % countdensity normalization
        ylim([0 0.2]); %probability normalization
        
        legend(strcat("N = ", num2str(numCad)));
        print(strcat(outDir,data.filePre,'_tod_histogram2.png'),'-dpng');
        print(strcat(outDir,data.filePre,'_tod_histogram2.pdf'),'-dpdf');
        close();
    end
end