NWB: Insert recreated epoch data into LNB on NWB_ExportAllData

Epoch data is recreated for each sweep and inserted with the additional
LNB column "PostProcessed" set to 1.

Author: MichaelHuth

Packages/MIES/MIES_Epochs.ipf

@@ -1741,7 +1741,7 @@ End
 /// @param sweepDFR        single sweep folder, e.g. for measurement with a device this wold be DFREF sweepDFR = GetSingleSweepFolder(deviceDFR, sweepNo)
 /// @param sweepNo         sweep number
 /// @returns recreated 4D epoch wave
-static Function/WAVE EP_RecreateEpochsFromLoadedData(WAVE numericalValues, WAVE/T textualValues, DFREF sweepDFR, variable sweepNo)
+Function/WAVE EP_RecreateEpochsFromLoadedData(WAVE numericalValues, WAVE/T textualValues, DFREF sweepDFR, variable sweepNo)
 
 	STRUCT DataConfigurationResult s
 	variable channelNr, plannedTime, acquiredTime, adSize, firstUnacquiredIndex

Packages/MIES/MIES_MiesUtilities_Logbook.ipf

@@ -1938,3 +1938,42 @@ Function ParseLogbookMode(string modeText)
 			break
 	endswitch
 End
+
+Function InsertRecreatedEpochsIntoLBN(WAVE numericalValues, WAVE/T textualValues, string device, variable sweepNo)
+
+	string epochList
+	variable channelNumber, channelType
+	variable rowIndex, colIndexEpochs, headstage
+
+	DFREF  deviceDFR = GetDeviceDataPath(device)
+	DFREF  sweepDFR  = GetSingleSweepFolder(deviceDFR, sweepNo)
+	WAVE/Z recEpochs = EP_RecreateEpochsFromLoadedData(numericalValues, textualValues, sweepDFR, sweepNo)
+	if(!WaveExists(recEpochs))
+		print "Could not recreate Epochs."
+		return NaN
+	endif
+
+	[rowIndex, colIndexEpochs] = ED_InsertPostProcessRow(textualValues, device, sweepNo, DATA_ACQUISITION_MODE, EPOCHS_ENTRY_KEY)
+
+	for(channelType = 0; channelType < XOP_CHANNEL_TYPE_COUNT; channelType += 1)
+		for(channelNumber = 0; channelNumber < NUM_DA_TTL_CHANNELS; channelNumber += 1)
+			// Currently only implemented for DAC channel type
+			if(channelType != XOP_CHANNEL_TYPE_DAC)
+				continue
+			endif
+
+			epochList = EP_EpochWaveToStr(recEpochs, channelNumber, channelType)
+			if(IsEmpty(epochList))
+				continue
+			endif
+
+			headstage = AFH_GetHeadstageFromChannelNumberAndType(device, channelType, channelNumber)
+			if(!IsValidHeadstage(headstage))
+				printf "Could not resolve headstage for device %s, channel %d, type %d but got epoch info for it.\r", device, channelNumber, channelType
+				continue
+			endif
+
+			textualValues[rowIndex][colIndexEpochs][headstage] = epochList
+		endfor
+	endfor
+End

Packages/MIES/MIES_NeuroDataWithoutBorders.ipf

@@ -467,7 +467,6 @@ static Function NWB_AddDeviceSpecificData(STRUCT NWBAsyncParameters &s, variable
 
 	AddModificationTimeEntry(s.locationID, s.nwbVersion)
 	NWB_AddDevice(s)
-	NWB_WriteLabnotebooksAndComments(s)
 	NWB_WriteTestpulseData(s, writeStoredTestPulses)
 End
 
@@ -724,6 +723,8 @@ Function NWB_ExportAllData(variable nwbVersion, [string device, string overrideF
 			WAVE s.DAQDataWave   = TextSweepToWaveRef(sweepWave)
 			WAVE s.DAQConfigWave = configWave
 
+			InsertRecreatedEpochsIntoLBN(s.numericalValues, s.textualValues, s.device, s.sweep)
+
 			NWB_AppendSweepLowLevel(s)
 			stimsetList += AB_GetStimsetFromPanel(device, sweep)
 
@@ -732,6 +733,8 @@ Function NWB_ExportAllData(variable nwbVersion, [string device, string overrideF
 		endfor
 		LOG_AddEntry(PACKAGE_MIES, "export", keys = {"size [MiB]"}, values = {num2str(NWB_GetExportedFileSize(device))})
 
+		NWB_WriteLabnotebooksAndComments(s)
+
 		NWB_AppendStimset(nwbVersion, s.locationID, stimsetList, compressionMode)
 
 		if(writeIgorHistory)