/// <summary>
/// Test only method for creating a <see cref="XHunterLibrary"/> file
/// from another loaded <see cref="Library"/>. Should this move into test project?
/// </summary>
/// <param name="streamManager">Provides access to the file system</param>
/// <param name="path">Path to write to</param>
/// <param name="library">The loaded library to use as a data source</param>
/// <param name="lowIntensity">True to use 20 lowest intensity peaks for bad spectrum</param>
public static void Write(IStreamManager streamManager, string path, Library library, bool lowIntensity)
{
using (FileSaver fs = new FileSaver(path, streamManager))
using (Stream outStream = streamManager.CreateStream(fs.SafeName, FileMode.Create, true))
{
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int));
outStream.Write(BitConverter.GetBytes(library.SpectrumCount), 0, sizeof(int));
byte[] header = new byte[256 - 8];
const string headerText = @"HLF v=2 s=test.hlf d=2009.02.04";
Encoding.UTF8.GetBytes(headerText, 0, headerText.Length, header, 0);
outStream.Write(header, 0, header.Length);
SequenceMassCalc calc = new SequenceMassCalc(MassType.Monoisotopic);
byte[] seqBuffer = new byte[1024];
foreach (var key in library.Keys)
{
SpectrumPeaksInfo peaksInfo;
if (!library.TryLoadSpectrum(key, out peaksInfo))
{
continue;
}
// Fake X! Hunter filtering by choosing just the to 20 peaks
SpectrumPeaksInfo.MI[] peaks = peaksInfo.Peaks.ToArray();
// Sort by intensity
if (lowIntensity)
{
Array.Sort(peaks, (p1, p2) => Comparer.Default.Compare(p1.Intensity, p2.Intensity)); // ascending
}
else
{
Array.Sort(peaks, (p1, p2) => Comparer.Default.Compare(p2.Intensity, p1.Intensity)); // descending
}
float maxI = peaks.Length == 0 ? 0 : peaks[0].Intensity;
// Take 20 most intense peaks
SpectrumPeaksInfo.MI[] peaksFiltered = new SpectrumPeaksInfo.MI[Math.Min(20, peaks.Length)];
Array.Copy(peaks, peaksFiltered, peaksFiltered.Length);
// Resort by m/z (ineffient, but this is test code)
Array.Sort(peaksFiltered, (p1, p2) => Comparer.Default.Compare(p1.Mz, p2.Mz));
double totalI = 0;
byte[] peakBytes = new byte[(sizeof(byte) + sizeof(float)) * peaksFiltered.Length];
for (int i = 0; i < peaksFiltered.Length; i++)
{
var mi = peaksFiltered[i];
// Calculate the X! Hunter processed intensity value
float intensity = 100f * mi.Intensity / maxI;
totalI += intensity;
// Fill the peaks buffer
peakBytes[i] = (byte)(int)intensity;
Array.Copy(BitConverter.GetBytes((float)mi.Mz), 0, peakBytes, peaksFiltered.Length + i * 4, sizeof(float));
}
var sequence = key.Target.ToString();
// Only works for unmodified sequence
Debug.Assert(!key.IsModified);
double precursorMH = calc.GetPrecursorMass(sequence);
outStream.Write(BitConverter.GetBytes(precursorMH), 0, sizeof(double));
outStream.Write(BitConverter.GetBytes(key.Charge), 0, sizeof(int));
// Value rounded for consistent serialization round-tripping
float i2 = (float)Math.Round(Math.Sqrt(totalI), 4);
outStream.Write(BitConverter.GetBytes(i2), 0, sizeof(float));
outStream.Write(BitConverter.GetBytes(0.0001f), 0, sizeof(float));
outStream.Write(BitConverter.GetBytes(sequence.Length), 0, sizeof(int));
// Sequence
Encoding.UTF8.GetBytes(sequence, 0, sequence.Length, seqBuffer, 0);
outStream.Write(seqBuffer, 0, sequence.Length);
// Peaks
outStream.Write(BitConverter.GetBytes(peaksFiltered.Length), 0, sizeof(int));
outStream.Write(peakBytes, 0, peakBytes.Length);
// Modifications
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int));
// Homologs
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int));
}
streamManager.Finish(outStream);
fs.Commit();
}
}
/// <summary>
/// Test only method for creating a <see cref="BiblioSpecLibrary"/> file
/// from another loaded <see cref="Library"/>. Should this move into test project?
/// </summary>
/// <param name="streamManager">Provides access to the file system</param>
/// <param name="path">Path to write to</param>
/// <param name="library">The loaded library to use as a data source</param>
public static void Write(IStreamManager streamManager, string path, Library library)
{
using (FileSaver fs = new FileSaver(path, streamManager))
using (Stream outStream = streamManager.CreateStream(fs.SafeName, FileMode.Create, true))
{
outStream.Write(BitConverter.GetBytes(library.SpectrumCount), 0, sizeof(int)); // num_spectra
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // filtered
outStream.Write(BitConverter.GetBytes(1), 0, sizeof(int)); // version1
outStream.Write(BitConverter.GetBytes(1), 0, sizeof(int)); // version2
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // next_id
SequenceMassCalc calc = new SequenceMassCalc(MassType.Monoisotopic);
byte[] seqBuffer = new byte[1024];
int scanNum = 1;
foreach (var key in library.Keys)
{
SpectrumPeaksInfo peaksInfo;
if (!library.TryLoadSpectrum(key, out peaksInfo))
{
continue;
}
string sequence = key.Sequence;
// Only works for unmodified sequence
Debug.Assert(!key.IsModified);
double precursorMH = calc.GetPrecursorMass(sequence);
int charge = key.Charge;
float precursorMz = (float)SequenceMassCalc.GetMZ(precursorMH, charge);
outStream.Write(BitConverter.GetBytes(scanNum), 0, sizeof(int)); // scan_num
outStream.Write(BitConverter.GetBytes(2), 0, sizeof(int)); // scan_type
outStream.Write(BitConverter.GetBytes(precursorMz), 0, sizeof(float)); // pre_mz
outStream.Write(BitConverter.GetBytes(charge), 0, sizeof(int)); // scan_type
outStream.Write(BitConverter.GetBytes(0f), 0, sizeof(int)); // r_time
outStream.Write(BitConverter.GetBytes(peaksInfo.Peaks.Length), 0, sizeof(int)); // num_peaks
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // 32-bit peak_ptr
outStream.Write(BitConverter.GetBytes(sequence.Length), 0, sizeof(int)); // seq_len
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // annot
outStream.Write(BitConverter.GetBytes(scanNum), 0, sizeof(int)); // copies (bogus value for ranking)
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // lib_id
scanNum++;
// Sequence
int len = sequence.Length;
seqBuffer[len] = 0;
Encoding.UTF8.GetBytes(sequence, 0, len, seqBuffer, 0);
outStream.Write(seqBuffer, 0, len + 1);
// Modifications
const string zeros = "000000000000000000000000000000000000000000000000000"; // Not L10N
Encoding.UTF8.GetBytes(zeros.Substring(0, len), 0, len, seqBuffer, 0);
outStream.Write(seqBuffer, 0, len + 1);
// Peaks
foreach (var mi in peaksInfo.Peaks)
{
outStream.Write(BitConverter.GetBytes((float)mi.Mz), 0, sizeof(float));
outStream.Write(BitConverter.GetBytes(mi.Intensity), 0, sizeof(float));
}
}
streamManager.Finish(outStream);
fs.Commit();
}
}
/// <summary>
/// Test only method for creating a <see cref="BiblioSpecLibrary"/> file
/// from another loaded <see cref="Library"/>. Should this move into test project?
/// </summary>
/// <param name="streamManager">Provides access to the file system</param>
/// <param name="path">Path to write to</param>
/// <param name="library">The loaded library to use as a data source</param>
public static void Write(IStreamManager streamManager, string path, Library library)
{
using (FileSaver fs = new FileSaver(path, streamManager))
using (Stream outStream = streamManager.CreateStream(fs.SafeName, FileMode.Create, true))
{
outStream.Write(BitConverter.GetBytes(library.SpectrumCount), 0, sizeof(int)); // num_spectra
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // filtered
outStream.Write(BitConverter.GetBytes(1), 0, sizeof(int)); // version1
outStream.Write(BitConverter.GetBytes(1), 0, sizeof(int)); // version2
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // next_id
SequenceMassCalc calc = new SequenceMassCalc(MassType.Monoisotopic);
byte[] seqBuffer = new byte[1024];
int scanNum = 1;
foreach (var key in library.Keys)
{
SpectrumPeaksInfo peaksInfo;
if (!library.TryLoadSpectrum(key, out peaksInfo))
continue;
string sequence = key.Sequence;
// Only works for unmodified sequence
Debug.Assert(!key.IsModified);
double precursorMH = calc.GetPrecursorMass(sequence);
int charge = key.Charge;
float precursorMz = (float) SequenceMassCalc.GetMZ(precursorMH, charge);
outStream.Write(BitConverter.GetBytes(scanNum), 0, sizeof(int)); // scan_num
outStream.Write(BitConverter.GetBytes(2), 0, sizeof(int)); // scan_type
outStream.Write(BitConverter.GetBytes(precursorMz), 0, sizeof(float)); // pre_mz
outStream.Write(BitConverter.GetBytes(charge), 0, sizeof(int)); // scan_type
outStream.Write(BitConverter.GetBytes(0f), 0, sizeof(int)); // r_time
outStream.Write(BitConverter.GetBytes(peaksInfo.Peaks.Length), 0, sizeof(int)); // num_peaks
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // 32-bit peak_ptr
outStream.Write(BitConverter.GetBytes(sequence.Length), 0, sizeof(int)); // seq_len
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // annot
outStream.Write(BitConverter.GetBytes(scanNum), 0, sizeof(int)); // copies (bogus value for ranking)
outStream.Write(BitConverter.GetBytes(0), 0, sizeof(int)); // lib_id
scanNum++;
// Sequence
int len = sequence.Length;
seqBuffer[len] = 0;
Encoding.UTF8.GetBytes(sequence, 0, len, seqBuffer, 0);
outStream.Write(seqBuffer, 0, len + 1);
// Modifications
const string zeros = "000000000000000000000000000000000000000000000000000"; // Not L10N
Encoding.UTF8.GetBytes(zeros.Substring(0, len), 0, len, seqBuffer, 0);
outStream.Write(seqBuffer, 0, len + 1);
// Peaks
foreach (var mi in peaksInfo.Peaks)
{
outStream.Write(BitConverter.GetBytes((float)mi.Mz), 0, sizeof(float));
outStream.Write(BitConverter.GetBytes(mi.Intensity), 0, sizeof(float));
}
}
streamManager.Finish(outStream);
fs.Commit();
}
}