Figure 2D: use dual_color recording file (plot every session and choose best one)
Figure 2EF: use single channel recording (ACh_popRate and Oxt_popRate during NREM sleep)


XXXXXXX_popRate.mat data structure

dFoverF:  100 Hz sampling rate and aligned already
lfp_fbtime: [Tstart  Tend]
popRate: 100 Hz sampling rate, not aligned.