Colormaps = [plt.get_cmap(Map) for Map in ['Reds', 'Greens', 'Blues', 'Purples', 'Greys', 'Oranges']] Fig, Axes = plt.subplots(len(DataInfo['NoiseFrequency']), sharex=True, figsize=(6, 12)) for FKey in FreqList: Freq = FreqList.index(FKey) FigTitle = 'Power spectral density' AxTitle = FKey + ' Hz' YLabel = 'PSD [V$^{2}$/Hz]' XLabel = 'Frequency [Hz]' # AmpFList = list(SoundIntensity[FKey].keys()); AmpFList.sort() Intensities = [80, 70, 60, 50, 40, 0] AmpFList = ControlSoundBoard.dBToAmpF(Intensities, FileName) for AFKey in AmpFList: AmpFList[AFKey] = [str(_) for _ in AmpFList[AFKey]] for AKey in AmpFList[FKey]: AmpF = AmpFList[FKey].index(AKey) if AKey == '0.0': AKey = '0' Colors = [Colormaps[Freq](255-(AmpF*255//len(AmpFList)))] LineLabel = str(round(SoundIntensity[FKey][AKey])) + ' dB' for IAmpF in Intensity[FKey]: IdB = Intensity[FKey][IAmpF]['dB'] if SoundIntensity[FKey][AKey] == IdB: # print(IAmpF)
SoundPostPauseDur = 0.039 SoundPulseNo = 2 SoundStimBlockNo = 1 SoundPauseBetweenStimBlocksDur = 1 Intensities = [80] NoiseFrequency = [[8000, 10000]] SBOutAmpF = 1 LaserPrePauseDur = 0.004 LaserPulseDur = 0.01 LaserPostPauseDur = 0.036 LaserPulseNo = 2 LaserStimBlockNo = 1 LaserPauseBetweenStimBlocksDur = 1 SoundAmpF = ControlSoundBoard.dBToAmpF(Intensities, CalibrationFile) SoundPulse = ControlSoundBoard.GenNoise(Rate, SoundPulseDur) SoundPulseFiltered = ControlSoundBoard.BandpassFilterSound(SoundPulse, Rate, NoiseFrequency) SoundUnit = ControlSoundBoard.ApplySoundAmpF(SoundPulseFiltered, Rate, SoundAmpF, NoiseFrequency, SBOutAmpF, SoundPrePauseDur, SoundPostPauseDur) SoundTTLUnit = ControlSoundBoard.GenTTL(Rate, SoundPulseDur, TTLAmpF, SoundTTLVal, SoundBoard, SBOutAmpF, SoundPrePauseDur, SoundPostPauseDur) LaserUnit = ControlSoundBoard.GenTTL(Rate, LaserPulseDur, TTLAmpF, LaserTTLVal, SoundBoard, SBOutAmpF, LaserPrePauseDur, LaserPostPauseDur) SoundTTLLaserUnit = [LaserUnit[_]+SoundTTLUnit[_] \ for _ in range(len(SoundTTLUnit))]