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))]
