private MSScan GetMSScan(int argScanNo)
{
MSScan msScan = new MSScan(argScanNo);
// GlypID.Readers.clsRawData rawData = new GlypID.Readers.clsRawData() ;
GlypID.HornTransform.clsHornTransform massTransform = new GlypID.HornTransform.clsHornTransform();
GlypID.Peaks.clsPeakProcessor peakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.HornTransform.clsHornTransformParameters transform_parameters = new GlypID.HornTransform.clsHornTransformParameters();
GlypID.Peaks.clsPeakProcessorParameters peak_parameters = new GlypID.Peaks.clsPeakProcessorParameters();
// Declarations
float[] _mzs = null;
float[] _intensities = null;
GlypID.Peaks.clsPeak[] peaks;
GlypID.HornTransform.clsHornTransformResults[] _transformResults;
// Settings
peak_parameters.SignalToNoiseThreshold = 3.0;
peak_parameters.PeakBackgroundRatio = 5.0;
peakProcessor.SetOptions(peak_parameters);
peakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
transform_parameters.PeptideMinBackgroundRatio = 5.0;
transform_parameters.MaxCharge = 10;
massTransform.TransformParameters = transform_parameters;
// Load
Raw = new GlypID.Readers.clsRawData(_fullFilePath, GlypID.Readers.FileType.FINNIGAN);
if (Raw.GetMSLevel(argScanNo) == 1)
{
// Get spectra
Raw.GetSpectrum(argScanNo, ref _mzs, ref _intensities);
// Get peaks
peaks = new GlypID.Peaks.clsPeak[1];
peakProcessor.DiscoverPeaks(ref _mzs, ref _intensities, ref peaks, Convert.ToSingle(transform_parameters.MinMZ), Convert.ToSingle(transform_parameters.MaxMZ), true);
// Deisotope
double min_background_intensity = peakProcessor.GetBackgroundIntensity(ref _intensities);
double min_peptide_intensity = min_background_intensity * transform_parameters.PeptideMinBackgroundRatio;
_transformResults = new GlypID.HornTransform.clsHornTransformResults[1];
massTransform.PerformTransform(Convert.ToSingle(min_background_intensity), Convert.ToSingle(min_peptide_intensity), ref _mzs, ref _intensities, ref peaks, ref _transformResults);
Array.Clear(_transformResults, 0, _transformResults.Length);
Array.Clear(peaks, 0, peaks.Length);
}
return msScan;
}
/// <summary>
/// Function to normalize peaks into bins of set width
/// </summary>
/// <param name="?"></param>
public void BinNormalize(ref GlypID.Peaks.clsPeak[] peaks, float min_mz, float max_mz, double bin_size)
{
int num_bins = (int)((max_mz - min_mz) / bin_size)+1;
GlypID.Peaks.clsPeak[] baseSpectrumPeaks = new GlypID.Peaks.clsPeak[num_bins];
for (int i = 0; i < baseSpectrumPeaks.Length; i++)
{
baseSpectrumPeaks[i] = new GlypID.Peaks.clsPeak();
baseSpectrumPeaks[i].mdbl_mz = min_mz + i * bin_size;
baseSpectrumPeaks[i].mdbl_intensity = 0;
}
for (int i = 0; i < peaks.Length; i++)
{
int bin = (int) ((peaks[i].mdbl_mz - min_mz) / bin_size);
baseSpectrumPeaks[i].mdbl_intensity += peaks[i].mdbl_intensity;
}
Array.Clear(peaks, 0, peaks.Length);
peaks = baseSpectrumPeaks;
}
private static MSScan GetScanFromFile(int argScanNo, double argSingleToNoise, double argPeakBackground, double argPeptideBackground, short argMaxCharge, GlypID.Readers.clsRawData Raw)
{
float[] _cidMzs = null;
float[] _cidIntensities = null;
GlypID.Peaks.clsPeak[] _cidPeaks = new GlypID.Peaks.clsPeak[1];
GlypID.Peaks.clsPeak[] _parentPeaks = new GlypID.Peaks.clsPeak[1];
GlypID.HornTransform.clsHornTransform mobjTransform = new GlypID.HornTransform.clsHornTransform();
GlypID.HornTransform.clsHornTransformParameters mobjTransformParameters = new GlypID.HornTransform.clsHornTransformParameters() ;
GlypID.HornTransform.clsHornTransformResults[] _transformResult;
GlypID.Peaks.clsPeakProcessor cidPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.Peaks.clsPeakProcessorParameters cidPeakParameters = new GlypID.Peaks.clsPeakProcessorParameters();
GlypID.Peaks.clsPeakProcessor parentPeakProcessor = new GlypID.Peaks.clsPeakProcessor() ;
GlypID.Peaks.clsPeakProcessorParameters parentPeakParameters = new GlypID.Peaks.clsPeakProcessorParameters();
//Start Read Scan
MSScan scan = new MSScan(argScanNo);
Raw.GetSpectrum(argScanNo, ref _cidMzs, ref _cidIntensities);
scan.MsLevel = Raw.GetMSLevel(Convert.ToInt32(argScanNo));
double min_peptide_intensity = 0;
scan.Time = Raw.GetScanTime(scan.ScanNo);
scan.ScanHeader = Raw.GetScanDescription(scan.ScanNo);
if (scan.MsLevel != 1)
{
float[] _parentRawMzs = null;
float[] _parentRawIntensitys = null;
string Header = Raw.GetScanDescription(argScanNo);
cidPeakProcessor.ProfileType = GlypID.enmProfileType.CENTROIDED;
if (Header.Substring(Header.IndexOf("+") + 1).Trim().StartsWith("p"))
{
cidPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
}
// cidPeakProcessor.DiscoverPeaks(ref _cidMzs, ref _cidIntensities, ref _cidPeaks,
// Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), false);
for (int chNum = 0; chNum < _cidMzs.Length; chNum++)
{
scan.MSPeaks.Add(new MSPeak(
Convert.ToSingle(_cidMzs[chNum]),
Convert.ToSingle(_cidIntensities[chNum])));
}
//for (int chNum = 0; chNum < _cidMzs.Length; chNum++)
//{
// scan.MSPeaks.Add(new MSPeak(
// Convert.ToSingle(_cidMzs[chNum]),
// Convert.ToSingle(_cidIntensities[chNum])));
//}
// Get parent information
scan.ParentScanNo = Raw.GetParentScan(scan.ScanNo);
Raw.GetSpectrum(scan.ParentScanNo, ref _parentRawMzs, ref _parentRawIntensitys);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _parentRawMzs, ref _parentRawIntensitys, ref _parentPeaks, Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
float _parentBackgroundIntensity = (float)parentPeakProcessor.GetBackgroundIntensity(ref _parentRawIntensitys);
_transformResult = new GlypID.HornTransform.clsHornTransformResults[1];
bool found = false;
if (Raw.IsFTScan(scan.ParentScanNo))
{
// High resolution data
found = mobjTransform.FindPrecursorTransform(Convert.ToSingle(_parentBackgroundIntensity), Convert.ToSingle(min_peptide_intensity), ref _parentRawMzs, ref _parentRawIntensitys, ref _parentPeaks, Convert.ToSingle(scan.ParentMZ), ref _transformResult);
}
if (!found)//de-isotope fail
{
// Low resolution data or bad high res spectra
short cs = Raw.GetMonoChargeFromHeader(scan.ScanNo);
double monoMZ = Raw.GetMonoMzFromHeader(scan.ScanNo);
List<float> ParentMzs = new List<float>(_parentRawMzs);
int CloseIdx = MassUtility.GetClosestMassIdx(ParentMzs, Convert.ToSingle(monoMZ));
if (cs > 0)
{
short[] charges = new short[1];
charges[0] = cs;
mobjTransform.AllocateValuesToTransform(Convert.ToSingle(scan.ParentMZ), Convert.ToInt32(_parentRawIntensitys[CloseIdx]), ref charges, ref _transformResult);
}
else
{
// instrument has no charge just store 2 and 3.
short[] charges = new short[2];
charges[0] = 2;
charges[1] = 3;
mobjTransform.AllocateValuesToTransform(Convert.ToSingle(scan.ParentMZ), Convert.ToInt32(_parentRawIntensitys[CloseIdx]), ref charges, ref _transformResult);
}
}
if (_transformResult[0].mint_peak_index == -1) //De-isotope parent scan
{
//Get parent info
MSScan _parentScan = GetScanFromFile(scan.ParentScanNo, argSingleToNoise, argPeakBackground, argPeptideBackground, argMaxCharge,Raw);
float[] _MSMzs = null;
float[] _MSIntensities = null;
Raw.GetSpectrum(scan.ParentScanNo, ref _MSMzs, ref _MSIntensities);
// Now find peaks
parentPeakParameters.SignalToNoiseThreshold = 0;
parentPeakParameters.PeakBackgroundRatio = 0.01;
parentPeakProcessor.SetOptions(parentPeakParameters);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _MSMzs, ref _MSIntensities, ref _cidPeaks,
Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
//Look for charge and mono.
float[] monoandcharge = FindChargeAndMono(_cidPeaks, Convert.ToSingle(Raw.GetParentMz(scan.ScanNo)), scan.ScanNo,Raw);
//scan.ParentMonoMW = _parentScan.MSPeaks[ClosedIdx].MonoMass;
//scan.ParentAVGMonoMW = _parentScan.MSPeaks[ClosedIdx].;
scan.ParentMZ = monoandcharge[0];
if (monoandcharge[1] == 0.0f)
{
scan.ParentCharge = Convert.ToInt32(Raw.GetMonoChargeFromHeader(scan.ParentScanNo));
}
else
{
scan.ParentCharge = Convert.ToInt32(monoandcharge[1]);
}
scan.ParentMonoMW = (monoandcharge[0] - Atoms.ProtonMass) * monoandcharge[1];
}
else
{
scan.ParentMonoMW = (float)_transformResult[0].mdbl_mono_mw;
scan.ParentAVGMonoMW = (float)_transformResult[0].mdbl_average_mw;
scan.ParentMZ = (float)_transformResult[0].mdbl_mz;
scan.ParentCharge = (int)_transformResult[0].mshort_cs;
}
scan.IsCIDScan = Raw.IsCIDScan(argScanNo);
scan.IsFTScan = Raw.IsFTScan(argScanNo);
Array.Clear(_transformResult, 0, _transformResult.Length);
Array.Clear(_cidPeaks, 0, _cidPeaks.Length);
Array.Clear(_cidMzs, 0, _cidMzs.Length);
Array.Clear(_cidIntensities, 0, _cidIntensities.Length);
Array.Clear(_parentRawMzs, 0, _parentRawMzs.Length);
Array.Clear(_parentRawIntensitys, 0, _parentRawIntensitys.Length);
}
else //MS Scan
{
scan.ParentMZ = 0.0f;
double mdbl_current_background_intensity = 0;
// Now find peaks
parentPeakParameters.SignalToNoiseThreshold = argSingleToNoise;
parentPeakParameters.PeakBackgroundRatio = argPeakBackground;
parentPeakProcessor.SetOptions(parentPeakParameters);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _cidMzs, ref _cidIntensities, ref _cidPeaks,
Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
mdbl_current_background_intensity = parentPeakProcessor.GetBackgroundIntensity(ref _cidIntensities);
// Settings
min_peptide_intensity = mdbl_current_background_intensity * mobjTransformParameters.PeptideMinBackgroundRatio;
if (mobjTransformParameters.UseAbsolutePeptideIntensity)
{
if (min_peptide_intensity < mobjTransformParameters.AbsolutePeptideIntensity)
min_peptide_intensity = mobjTransformParameters.AbsolutePeptideIntensity;
}
mobjTransformParameters.PeptideMinBackgroundRatio = argPeptideBackground;
mobjTransformParameters.MaxCharge = argMaxCharge;
mobjTransform.TransformParameters = mobjTransformParameters;
// Now perform deisotoping
_transformResult = new GlypID.HornTransform.clsHornTransformResults[1];
mobjTransform.PerformTransform(Convert.ToSingle(mdbl_current_background_intensity), Convert.ToSingle(min_peptide_intensity), ref _cidMzs, ref _cidIntensities, ref _cidPeaks, ref _transformResult);
// for getting results
for (int chNum = 0; chNum < _transformResult.Length; chNum++)
{
double sumintensity = 0.0;
double mostIntenseIntensity = 0.0;
for (int i = 0; i < _transformResult[chNum].marr_isotope_peak_indices.Length; i++)
{
sumintensity = sumintensity + _cidPeaks[_transformResult[chNum].marr_isotope_peak_indices[i]].mdbl_intensity;
if (Math.Abs(_transformResult[chNum].mdbl_most_intense_mw -
(_cidPeaks[_transformResult[chNum].marr_isotope_peak_indices[i]].mdbl_mz * _transformResult[chNum].mshort_cs - Atoms.ProtonMass * _transformResult[chNum].mshort_cs))
< 1.0 / _transformResult[chNum].mshort_cs)
{
mostIntenseIntensity = _cidPeaks[_transformResult[chNum].mint_peak_index].mdbl_intensity;
}
}
scan.MSPeaks.Add(new MSPeak(
Convert.ToSingle(_transformResult[chNum].mdbl_mono_mw),
_transformResult[chNum].mint_mono_intensity,
_transformResult[chNum].mshort_cs,
Convert.ToSingle(_transformResult[chNum].mdbl_mz),
Convert.ToSingle(_transformResult[chNum].mdbl_fit),
Convert.ToSingle(_transformResult[chNum].mdbl_most_intense_mw),
mostIntenseIntensity,
sumintensity
));
}
Array.Clear(_transformResult, 0, _transformResult.Length);
Array.Clear(_cidPeaks, 0, _cidPeaks.Length);
Array.Clear(_cidMzs, 0, _cidMzs.Length);
Array.Clear(_cidIntensities, 0, _cidIntensities.Length);
}
return scan;
}
/// <summary>
/// No HCD info the return will be null
/// </summary>
/// <param name="argReader"></param>
/// <param name="argScanNo"></param>
/// <returns></returns>
private static HCDInfo GetHCDInfo(GlypID.Readers.clsRawData argReader, int argScanNo)
{
if(!argReader.IsHCDScan(argScanNo))
{
return null;
}
HCDInfo HCDinfo = null;
double _hcd_score;
GlypID.enmGlycanType _gType;
int ParentScan = argReader.GetParentScan(argScanNo);
// Scorers and transforms
GlypID.HCDScoring.clsHCDScoring HCDScoring = new GlypID.HCDScoring.clsHCDScoring();
GlypID.HornTransform.clsHornTransform Transform = new GlypID.HornTransform.clsHornTransform();
// mzs , intensities
float[] hcd_mzs = null;
float[] hcd_intensities = null;
float[] parent_mzs = null;
float[] parent_intensities = null;
// Peaks
GlypID.Peaks.clsPeak[] parent_peaks;
GlypID.Peaks.clsPeak[] hcd_peaks;
// Peak Processors
GlypID.Peaks.clsPeakProcessor hcdPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.Peaks.clsPeakProcessor parentPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
// Results
GlypID.HCDScoring.clsHCDScoringScanResults[] hcd_scoring_results;
GlypID.HornTransform.clsHornTransformResults[] transform_results;
// Params
GlypID.Scoring.clsScoringParameters scoring_parameters = new GlypID.Scoring.clsScoringParameters();
GlypID.HornTransform.clsHornTransformParameters transform_parameters = new GlypID.HornTransform.clsHornTransformParameters();
scoring_parameters.MinNumPeaksToConsider = 2;
scoring_parameters.PPMTolerance = 10;
scoring_parameters.UsePPM = true;
// Init
Transform.TransformParameters = transform_parameters;
// Loading parent
int parent_scan = argReader.GetParentScan(argScanNo);
double parent_mz = argReader.GetParentMz(argScanNo);
int scan_level = argReader.GetMSLevel(argScanNo);
int parent_level = argReader.GetMSLevel(parent_scan);
argReader.GetSpectrum(parent_scan, ref parent_mzs, ref parent_intensities);
// Parent processing
parent_peaks = new GlypID.Peaks.clsPeak[1];
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref parent_mzs, ref parent_intensities, ref parent_peaks,
Convert.ToSingle(transform_parameters.MinMZ), Convert.ToSingle(transform_parameters.MaxMZ), true);
double bkg_intensity = parentPeakProcessor.GetBackgroundIntensity(ref parent_intensities);
double min_peptide_intensity = bkg_intensity * transform_parameters.PeptideMinBackgroundRatio;
transform_results = new GlypID.HornTransform.clsHornTransformResults[1];
bool found = Transform.FindPrecursorTransform(Convert.ToSingle(bkg_intensity), Convert.ToSingle(min_peptide_intensity), ref parent_mzs, ref parent_intensities, ref parent_peaks, Convert.ToSingle(parent_mz), ref transform_results);
if (!found && (argReader.GetMonoChargeFromHeader(ParentScan) > 0))
{
found = true;
double mono_mz = argReader.GetMonoMzFromHeader(ParentScan);
if (mono_mz == 0)
mono_mz = parent_mz;
short[] charges = new short[1];
charges[0] = argReader.GetMonoChargeFromHeader(ParentScan);
Transform.AllocateValuesToTransform(Convert.ToSingle(mono_mz), 0, ref charges, ref transform_results); // Change abundance value from 0 to parent_intensity if you wish
}
if (found && transform_results.Length == 1)
{
// Score HCD scan first
argReader.GetSpectrum(argScanNo, ref hcd_mzs, ref hcd_intensities);
double hcd_background_intensity = GlypID.Utils.GetAverage(ref hcd_intensities, ref hcd_mzs, Convert.ToSingle(scoring_parameters.MinHCDMz), Convert.ToSingle(scoring_parameters.MaxHCDMz));
hcdPeakProcessor.SetPeakIntensityThreshold(hcd_background_intensity);
hcd_peaks = new GlypID.Peaks.clsPeak[1];
//Check Header
string Header = argReader.GetScanDescription(argScanNo);
hcdPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
if (Header.Substring(Header.IndexOf("+") + 1).Trim().StartsWith("c"))
{
hcdPeakProcessor.ProfileType = GlypID.enmProfileType.CENTROIDED;
}
hcdPeakProcessor.DiscoverPeaks(ref hcd_mzs, ref hcd_intensities, ref hcd_peaks, Convert.ToSingle
(scoring_parameters.MinHCDMz), Convert.ToSingle(scoring_parameters.MaxHCDMz), false);
hcdPeakProcessor.InitializeUnprocessedData();
hcd_scoring_results = new GlypID.HCDScoring.clsHCDScoringScanResults[1];
HCDScoring.ScoringParameters = scoring_parameters;
_hcd_score = HCDScoring.ScoreHCDSpectra(ref hcd_peaks, ref hcd_mzs, ref hcd_intensities, ref transform_results, ref hcd_scoring_results);
_gType = (GlypID.enmGlycanType)hcd_scoring_results[0].menm_glycan_type;
enumGlycanType GType; //Convert from GlypID.enumGlycanType to MassLib.enumGlycanType;
if (_gType == GlypID.enmGlycanType.CA)
{
GType = enumGlycanType.CA;
}
else if (_gType == GlypID.enmGlycanType.CS)
{
GType = enumGlycanType.CS;
}
else if (_gType == GlypID.enmGlycanType.HM)
{
GType = enumGlycanType.HM;
}
else if (_gType == GlypID.enmGlycanType.HY)
{
GType = enumGlycanType.HY;
}
else
{
GType = enumGlycanType.NA;
}
HCDinfo = new HCDInfo(argScanNo, GType, _hcd_score);
}
return HCDinfo;
}
/// <summary>
/// Given a cluster index, return the representative of that cluster. Currently, the spectrum with maximum SNR in a given range is returned
/// </summary>
/// <param name="clusterIndex"></param>
/// <param name="repPeaks"></param>
/// <param name="minmz"></param>
/// <param name="maxmz"></param>
/// <param name="return_original_id">true or false. If true, return the original index(among all spectra), else returns the index within the cluster</param>
/// <returns></returns>
public int GetRepresentativePeaksFromCluster(int clusterIndex, ref GlypID.Peaks.clsPeak[] repPeaks, double minmz, double maxmz, bool return_original_id)
{
double max_snr = 0;
int spectra_with_max_snr = -1;
int orig_index = -1;
string clustername = _ClusterNames[clusterIndex];
string[] spectra = new string[1];
if (clustername.Contains('-'))
spectra = clustername.Split('-');
else
spectra[0] = clustername;
if (spectra.Length > 0)
{
foreach (string s in spectra)
{
List<GlypID.Peaks.clsPeak> this_peaks = new List<GlypID.Peaks.clsPeak>();
string[] parts = s.Split('_');
int spectra_index = Convert.ToInt32(parts[2]);
this_peaks = _AllOriginalSpectra[spectra_index];
double noise_level = CalculateAverage(ref this_peaks, 0, minmz, maxmz);
double signal_level = CalculateAverage(ref this_peaks, noise_level, minmz, maxmz);
double snr = signal_level / noise_level;
if (snr > max_snr)
{
max_snr = snr;
spectra_with_max_snr = spectra_index;
orig_index = Convert.ToInt32(parts[1]);
}
}
if (spectra_with_max_snr >= 0)
{
repPeaks = new GlypID.Peaks.clsPeak[_AllOriginalSpectra[spectra_with_max_snr].Count];
for (int i = 0; i < _AllOriginalSpectra[spectra_with_max_snr].Count; i++)
{
repPeaks[i] = _AllOriginalSpectra[spectra_with_max_snr][i];
}
}
}
if (return_original_id)
return orig_index;
else
return spectra_with_max_snr;
}
/// <summary>
/// Function to calculate representative fragmentation spectra for each multi record.
/// </summary>
/// <param name="_glycoMap"></param>
public void CalculateRepresentativeFragmentationSpectraThroughSum(ref Classes.MapRecord _glycoMap)
{
for (int i = 0; i < _glycoMap._AllMLNRecords.Count; i++)
{
MultiAlignRecord m = _glycoMap._AllMLNRecords[i];
int num_records = m._AssociatedUMCRecords.Count;
int num_cid_peaks = (int)((_TransformParameters.MaxMZ - _TransformParameters.MinMZ) / _binSize) + 1;
int num_hcd_peaks = (int) ((_ScoringParameters.MaxHCDMz - _ScoringParameters.MinHCDMz) / _binSize) + 1;
int num_etd_peaks = (int)((_TransformParameters.MaxMZ - _TransformParameters.MinMZ) / _binSize) + 1;
GlypID.Peaks.clsPeak[] summedCIDPeaks = new GlypID.Peaks.clsPeak[num_cid_peaks];
GlypID.Peaks.clsPeak[] summedHCDPeaks = new GlypID.Peaks.clsPeak[num_hcd_peaks];
GlypID.Peaks.clsPeak[] summedETDPeaks = new GlypID.Peaks.clsPeak[num_etd_peaks];
for (int j = 0; j < num_cid_peaks; j++)
{
summedCIDPeaks[j] = new GlypID.Peaks.clsPeak();
summedCIDPeaks[j].mdbl_mz = _TransformParameters.MinMZ + j * _binSize;
summedCIDPeaks[j].mdbl_intensity = 0;
}
for (int j = 0; j < num_etd_peaks; j++)
{
summedETDPeaks[j] = new GlypID.Peaks.clsPeak();
summedETDPeaks[j].mdbl_mz = _TransformParameters.MinMZ + j * _binSize;
summedETDPeaks[j].mdbl_intensity = 0;
}
for (int j = 0; j < num_hcd_peaks; j++)
{
summedHCDPeaks[j] = new GlypID.Peaks.clsPeak();
summedHCDPeaks[j].mdbl_mz = _TransformParameters.MinMZ + j * _binSize;
summedHCDPeaks[j].mdbl_intensity = 0;
}
for (int j = 0; j < num_records; j++)
{
if (m._AssociatedUMCRecords[j].CIDPeaks[0] != null)
{
GlypID.Peaks.clsPeak[] thisCIDPeaks = new GlypID.Peaks.clsPeak[m._AssociatedUMCRecords[j].CIDPeaks.Length];
Array.Copy(m._AssociatedUMCRecords[j].CIDPeaks, thisCIDPeaks, m._AssociatedUMCRecords[j].CIDPeaks.Length);
BinNormalize(ref thisCIDPeaks, (float)_TransformParameters.MinMZ, (float)_TransformParameters.MaxMZ, _binSize);
AddToPeaks(ref summedCIDPeaks, ref thisCIDPeaks);
}
if (m._AssociatedUMCRecords[j].HCDPeaks[0] !=null)
{
GlypID.Peaks.clsPeak[] thisHCDPeaks = new GlypID.Peaks.clsPeak[m._AssociatedUMCRecords[j].HCDPeaks.Length];
Array.Copy(m._AssociatedUMCRecords[j].HCDPeaks, thisHCDPeaks, m._AssociatedUMCRecords[j].HCDPeaks.Length);
BinNormalize(ref thisHCDPeaks, (float)_ScoringParameters.MinHCDMz, (float)_ScoringParameters.MaxHCDMz, _binSize);
AddToPeaks(ref summedHCDPeaks, ref thisHCDPeaks);
}
if (m._AssociatedUMCRecords[j].ETDPeaks[0] !=null)
{
GlypID.Peaks.clsPeak[] thisETDPeaks = new GlypID.Peaks.clsPeak[m._AssociatedUMCRecords[j].ETDPeaks.Length];
Array.Copy(m._AssociatedUMCRecords[j].ETDPeaks, thisETDPeaks, m._AssociatedUMCRecords[j].ETDPeaks.Length);
BinNormalize(ref thisETDPeaks, (float)_TransformParameters.MinMZ, (float)_TransformParameters.MaxMZ, _binSize);
AddToPeaks(ref summedETDPeaks, ref thisETDPeaks);
}
}
m._RepresentativeCIDPeaks = summedCIDPeaks;
m._RepresentativeETDPeaks = summedETDPeaks;
m._RepresentativeHCDPeaks = summedHCDPeaks;
}
}
public HCDInfo(GlypID.Readers.clsRawData argReader, int argScanNo)
{
int ParentScan = argReader.GetParentScan(argScanNo);
// Scorers and transforms
GlypID.HCDScoring.clsHCDScoring HCDScoring = new GlypID.HCDScoring.clsHCDScoring();
GlypID.HornTransform.clsHornTransform Transform = new GlypID.HornTransform.clsHornTransform();
// mzs , intensities
float[] hcd_mzs = null;
float[] hcd_intensities = null;
float[] parent_mzs = null;
float[] parent_intensities = null;
// Peaks
GlypID.Peaks.clsPeak[] parent_peaks;
GlypID.Peaks.clsPeak[] hcd_peaks;
// Peak Processors
GlypID.Peaks.clsPeakProcessor hcdPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.Peaks.clsPeakProcessor parentPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
// Results
GlypID.HCDScoring.clsHCDScoringScanResults[] hcd_scoring_results;
GlypID.HornTransform.clsHornTransformResults[] transform_results;
// Params
GlypID.Scoring.clsScoringParameters scoring_parameters = new GlypID.Scoring.clsScoringParameters();
GlypID.HornTransform.clsHornTransformParameters transform_parameters = new GlypID.HornTransform.clsHornTransformParameters();
scoring_parameters.MinNumPeaksToConsider = 2;
scoring_parameters.PPMTolerance = 10;
scoring_parameters.UsePPM = true;
// Init
Transform.TransformParameters = transform_parameters;
// Loading parent
int parent_scan = argReader.GetParentScan(argScanNo);
double parent_mz = argReader.GetParentMz(argScanNo);
int scan_level = argReader.GetMSLevel(argScanNo);
int parent_level = argReader.GetMSLevel(parent_scan);
argReader.GetSpectrum(parent_scan, ref parent_mzs, ref parent_intensities);
// Parent processing
parent_peaks = new GlypID.Peaks.clsPeak[1];
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref parent_mzs, ref parent_intensities, ref parent_peaks,
Convert.ToSingle(transform_parameters.MinMZ), Convert.ToSingle(transform_parameters.MaxMZ), true);
double bkg_intensity = parentPeakProcessor.GetBackgroundIntensity(ref parent_intensities);
double min_peptide_intensity = bkg_intensity * transform_parameters.PeptideMinBackgroundRatio;
transform_results = new GlypID.HornTransform.clsHornTransformResults[1];
bool found = Transform.FindPrecursorTransform(Convert.ToSingle(bkg_intensity), Convert.ToSingle(min_peptide_intensity), ref parent_mzs, ref parent_intensities, ref parent_peaks, Convert.ToSingle(parent_mz), ref transform_results);
if (!found && (argReader.GetMonoChargeFromHeader(ParentScan) > 0))
{
found = true;
double mono_mz = argReader.GetMonoMzFromHeader(ParentScan);
if (mono_mz == 0)
{
mono_mz = parent_mz;
}
short[] charges = new short[1];
charges[0] = argReader.GetMonoChargeFromHeader(ParentScan);
Transform.AllocateValuesToTransform(Convert.ToSingle(mono_mz), 0, ref charges, ref transform_results); // Change abundance value from 0 to parent_intensity if you wish
}
if (found && transform_results.Length == 1)
{
// Score HCD scan first
argReader.GetSpectrum(argScanNo, ref hcd_mzs, ref hcd_intensities);
double hcd_background_intensity = GlypID.Utils.GetAverage(ref hcd_intensities, ref hcd_mzs, Convert.ToSingle(scoring_parameters.MinHCDMz), Convert.ToSingle(scoring_parameters.MaxHCDMz));
hcdPeakProcessor.SetPeakIntensityThreshold(hcd_background_intensity);
hcd_peaks = new GlypID.Peaks.clsPeak[1];
//Check Header
string Header = argReader.GetScanDescription(argScanNo);
hcdPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
if (Header.Substring(Header.IndexOf("+") + 1).Trim().StartsWith("c"))
{
hcdPeakProcessor.ProfileType = GlypID.enmProfileType.CENTROIDED;
}
hcdPeakProcessor.DiscoverPeaks(ref hcd_mzs, ref hcd_intensities, ref hcd_peaks, Convert.ToSingle
(scoring_parameters.MinHCDMz), Convert.ToSingle(scoring_parameters.MaxHCDMz), false);
hcdPeakProcessor.InitializeUnprocessedData();
hcd_scoring_results = new GlypID.HCDScoring.clsHCDScoringScanResults[1];
HCDScoring.ScoringParameters = scoring_parameters;
_hcd_score = HCDScoring.ScoreHCDSpectra(ref hcd_peaks, ref hcd_mzs, ref hcd_intensities, ref transform_results, ref hcd_scoring_results);
_gType = (GlypID.enmGlycanType)hcd_scoring_results[0].menm_glycan_type;
}
}
private static MSScan GetScanFromFile(int argScanNo, double argSingleToNoise, double argPeakBackground, double argPeptideBackground, short argMaxCharge, GlypID.Readers.clsRawData Raw)
{
float[] _cidMzs = null;
float[] _cidIntensities = null;
GlypID.Peaks.clsPeak[] _cidPeaks = new GlypID.Peaks.clsPeak[1];
GlypID.Peaks.clsPeak[] _parentPeaks = new GlypID.Peaks.clsPeak[1];
GlypID.HornTransform.clsHornTransform mobjTransform = new GlypID.HornTransform.clsHornTransform();
GlypID.HornTransform.clsHornTransformParameters mobjTransformParameters = new GlypID.HornTransform.clsHornTransformParameters();
GlypID.HornTransform.clsHornTransformResults[] _transformResult;
GlypID.Peaks.clsPeakProcessor cidPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.Peaks.clsPeakProcessorParameters cidPeakParameters = new GlypID.Peaks.clsPeakProcessorParameters();
GlypID.Peaks.clsPeakProcessor parentPeakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.Peaks.clsPeakProcessorParameters parentPeakParameters = new GlypID.Peaks.clsPeakProcessorParameters();
//Start Read Scan
MSScan scan = new MSScan(argScanNo);
Raw.GetSpectrum(argScanNo, ref _cidMzs, ref _cidIntensities);
scan.MsLevel = Raw.GetMSLevel(Convert.ToInt32(argScanNo));
double min_peptide_intensity = 0;
scan.Time = Raw.GetScanTime(scan.ScanNo);
scan.ScanHeader = Raw.GetScanDescription(scan.ScanNo);
if (scan.MsLevel != 1)
{
float[] _parentRawMzs = null;
float[] _parentRawIntensitys = null;
string Header = Raw.GetScanDescription(argScanNo);
cidPeakProcessor.ProfileType = GlypID.enmProfileType.CENTROIDED;
if (Header.Substring(Header.IndexOf("+") + 1).Trim().StartsWith("p"))
{
cidPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
}
// cidPeakProcessor.DiscoverPeaks(ref _cidMzs, ref _cidIntensities, ref _cidPeaks,
// Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), false);
for (int chNum = 0; chNum < _cidMzs.Length; chNum++)
{
scan.MSPeaks.Add(new MSPeak(
Convert.ToSingle(_cidMzs[chNum]),
Convert.ToSingle(_cidIntensities[chNum])));
}
//for (int chNum = 0; chNum < _cidMzs.Length; chNum++)
//{
// scan.MSPeaks.Add(new MSPeak(
// Convert.ToSingle(_cidMzs[chNum]),
// Convert.ToSingle(_cidIntensities[chNum])));
//}
// Get parent information
scan.ParentScanNo = Raw.GetParentScan(scan.ScanNo);
Raw.GetSpectrum(scan.ParentScanNo, ref _parentRawMzs, ref _parentRawIntensitys);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _parentRawMzs, ref _parentRawIntensitys, ref _parentPeaks, Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
float _parentBackgroundIntensity = (float)parentPeakProcessor.GetBackgroundIntensity(ref _parentRawIntensitys);
_transformResult = new GlypID.HornTransform.clsHornTransformResults[1];
bool found = false;
if (Raw.IsFTScan(scan.ParentScanNo))
{
// High resolution data
found = mobjTransform.FindPrecursorTransform(Convert.ToSingle(_parentBackgroundIntensity), Convert.ToSingle(min_peptide_intensity), ref _parentRawMzs, ref _parentRawIntensitys, ref _parentPeaks, Convert.ToSingle(scan.ParentMZ), ref _transformResult);
}
if (!found)//de-isotope fail
{
// Low resolution data or bad high res spectra
short cs = Raw.GetMonoChargeFromHeader(scan.ScanNo);
double monoMZ = Raw.GetMonoMzFromHeader(scan.ScanNo);
List <float> ParentMzs = new List <float>(_parentRawMzs);
int CloseIdx = MassUtility.GetClosestMassIdx(ParentMzs, Convert.ToSingle(monoMZ));
if (cs > 0)
{
short[] charges = new short[1];
charges[0] = cs;
mobjTransform.AllocateValuesToTransform(Convert.ToSingle(scan.ParentMZ), Convert.ToInt32(_parentRawIntensitys[CloseIdx]), ref charges, ref _transformResult);
}
else
{
// instrument has no charge just store 2 and 3.
short[] charges = new short[2];
charges[0] = 2;
charges[1] = 3;
mobjTransform.AllocateValuesToTransform(Convert.ToSingle(scan.ParentMZ), Convert.ToInt32(_parentRawIntensitys[CloseIdx]), ref charges, ref _transformResult);
}
}
if (_transformResult[0].mint_peak_index == -1) //De-isotope parent scan
{
//Get parent info
MSScan _parentScan = GetScanFromFile(scan.ParentScanNo, argSingleToNoise, argPeakBackground, argPeptideBackground, argMaxCharge, Raw);
float[] _MSMzs = null;
float[] _MSIntensities = null;
Raw.GetSpectrum(scan.ParentScanNo, ref _MSMzs, ref _MSIntensities);
// Now find peaks
parentPeakParameters.SignalToNoiseThreshold = 0;
parentPeakParameters.PeakBackgroundRatio = 0.01;
parentPeakProcessor.SetOptions(parentPeakParameters);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _MSMzs, ref _MSIntensities, ref _cidPeaks,
Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
//Look for charge and mono.
float[] monoandcharge = FindChargeAndMono(_cidPeaks, Convert.ToSingle(Raw.GetParentMz(scan.ScanNo)), scan.ScanNo, Raw);
//scan.ParentMonoMW = _parentScan.MSPeaks[ClosedIdx].MonoMass;
//scan.ParentAVGMonoMW = _parentScan.MSPeaks[ClosedIdx].;
scan.ParentMZ = monoandcharge[0];
if (monoandcharge[1] == 0.0f)
{
scan.ParentCharge = Convert.ToInt32(Raw.GetMonoChargeFromHeader(scan.ParentScanNo));
}
else
{
scan.ParentCharge = Convert.ToInt32(monoandcharge[1]);
}
scan.ParentMonoMW = (monoandcharge[0] - Atoms.ProtonMass) * monoandcharge[1];
}
else
{
scan.ParentMonoMW = (float)_transformResult[0].mdbl_mono_mw;
scan.ParentAVGMonoMW = (float)_transformResult[0].mdbl_average_mw;
scan.ParentMZ = (float)_transformResult[0].mdbl_mz;
scan.ParentCharge = (int)_transformResult[0].mshort_cs;
}
scan.IsCIDScan = Raw.IsCIDScan(argScanNo);
scan.IsFTScan = Raw.IsFTScan(argScanNo);
Array.Clear(_transformResult, 0, _transformResult.Length);
Array.Clear(_cidPeaks, 0, _cidPeaks.Length);
Array.Clear(_cidMzs, 0, _cidMzs.Length);
Array.Clear(_cidIntensities, 0, _cidIntensities.Length);
Array.Clear(_parentRawMzs, 0, _parentRawMzs.Length);
Array.Clear(_parentRawIntensitys, 0, _parentRawIntensitys.Length);
}
else //MS Scan
{
scan.ParentMZ = 0.0f;
double mdbl_current_background_intensity = 0;
// Now find peaks
parentPeakParameters.SignalToNoiseThreshold = argSingleToNoise;
parentPeakParameters.PeakBackgroundRatio = argPeakBackground;
parentPeakProcessor.SetOptions(parentPeakParameters);
parentPeakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
parentPeakProcessor.DiscoverPeaks(ref _cidMzs, ref _cidIntensities, ref _cidPeaks,
Convert.ToSingle(mobjTransformParameters.MinMZ), Convert.ToSingle(mobjTransformParameters.MaxMZ), true);
mdbl_current_background_intensity = parentPeakProcessor.GetBackgroundIntensity(ref _cidIntensities);
// Settings
min_peptide_intensity = mdbl_current_background_intensity * mobjTransformParameters.PeptideMinBackgroundRatio;
if (mobjTransformParameters.UseAbsolutePeptideIntensity)
{
if (min_peptide_intensity < mobjTransformParameters.AbsolutePeptideIntensity)
{
min_peptide_intensity = mobjTransformParameters.AbsolutePeptideIntensity;
}
}
mobjTransformParameters.PeptideMinBackgroundRatio = argPeptideBackground;
mobjTransformParameters.MaxCharge = argMaxCharge;
mobjTransform.TransformParameters = mobjTransformParameters;
// Now perform deisotoping
_transformResult = new GlypID.HornTransform.clsHornTransformResults[1];
mobjTransform.PerformTransform(Convert.ToSingle(mdbl_current_background_intensity), Convert.ToSingle(min_peptide_intensity), ref _cidMzs, ref _cidIntensities, ref _cidPeaks, ref _transformResult);
// for getting results
for (int chNum = 0; chNum < _transformResult.Length; chNum++)
{
double sumintensity = 0.0;
double mostIntenseIntensity = 0.0;
for (int i = 0; i < _transformResult[chNum].marr_isotope_peak_indices.Length; i++)
{
sumintensity = sumintensity + _cidPeaks[_transformResult[chNum].marr_isotope_peak_indices[i]].mdbl_intensity;
if (Math.Abs(_transformResult[chNum].mdbl_most_intense_mw -
(_cidPeaks[_transformResult[chNum].marr_isotope_peak_indices[i]].mdbl_mz * _transformResult[chNum].mshort_cs - Atoms.ProtonMass * _transformResult[chNum].mshort_cs))
< 1.0 / _transformResult[chNum].mshort_cs)
{
mostIntenseIntensity = _cidPeaks[_transformResult[chNum].mint_peak_index].mdbl_intensity;
}
}
scan.MSPeaks.Add(new MSPeak(
Convert.ToSingle(_transformResult[chNum].mdbl_mono_mw),
_transformResult[chNum].mint_mono_intensity,
_transformResult[chNum].mshort_cs,
Convert.ToSingle(_transformResult[chNum].mdbl_mz),
Convert.ToSingle(_transformResult[chNum].mdbl_fit),
Convert.ToSingle(_transformResult[chNum].mdbl_most_intense_mw),
mostIntenseIntensity,
sumintensity
));
}
Array.Clear(_transformResult, 0, _transformResult.Length);
Array.Clear(_cidPeaks, 0, _cidPeaks.Length);
Array.Clear(_cidMzs, 0, _cidMzs.Length);
Array.Clear(_cidIntensities, 0, _cidIntensities.Length);
}
return(scan);
}
private MSScan GetMSScan(int argScanNo)
{
MSScan msScan = new MSScan(argScanNo);
// GlypID.Readers.clsRawData rawData = new GlypID.Readers.clsRawData() ;
GlypID.HornTransform.clsHornTransform massTransform = new GlypID.HornTransform.clsHornTransform();
GlypID.Peaks.clsPeakProcessor peakProcessor = new GlypID.Peaks.clsPeakProcessor();
GlypID.HornTransform.clsHornTransformParameters transform_parameters = new GlypID.HornTransform.clsHornTransformParameters();
GlypID.Peaks.clsPeakProcessorParameters peak_parameters = new GlypID.Peaks.clsPeakProcessorParameters();
// Declarations
float[] _mzs = null;
float[] _intensities = null;
GlypID.Peaks.clsPeak[] peaks;
GlypID.HornTransform.clsHornTransformResults[] _transformResults;
// Settings
peak_parameters.SignalToNoiseThreshold = 3.0;
peak_parameters.PeakBackgroundRatio = 5.0;
peakProcessor.SetOptions(peak_parameters);
peakProcessor.ProfileType = GlypID.enmProfileType.PROFILE;
transform_parameters.PeptideMinBackgroundRatio = 5.0;
transform_parameters.MaxCharge = 10;
massTransform.TransformParameters = transform_parameters;
// Load
Raw = new GlypID.Readers.clsRawData(_fullFilePath, GlypID.Readers.FileType.FINNIGAN);
if (Raw.GetMSLevel(argScanNo) == 1)
{
// Get spectra
Raw.GetSpectrum(argScanNo, ref _mzs, ref _intensities);
// Get peaks
peaks = new GlypID.Peaks.clsPeak[1];
peakProcessor.DiscoverPeaks(ref _mzs, ref _intensities, ref peaks, Convert.ToSingle(transform_parameters.MinMZ), Convert.ToSingle(transform_parameters.MaxMZ), true);
// Deisotope
double min_background_intensity = peakProcessor.GetBackgroundIntensity(ref _intensities);
double min_peptide_intensity = min_background_intensity * transform_parameters.PeptideMinBackgroundRatio;
_transformResults = new GlypID.HornTransform.clsHornTransformResults[1];
massTransform.PerformTransform(Convert.ToSingle(min_background_intensity), Convert.ToSingle(min_peptide_intensity), ref _mzs, ref _intensities, ref peaks, ref _transformResults);
Array.Clear(_transformResults, 0, _transformResults.Length);
Array.Clear(peaks, 0, peaks.Length);
}
return msScan;
}
/// <summary>
/// Function to cluster CID spectra in each record.
/// </summary>
/// <param name="_glycoMap"></param>
/// <param name="_params"></param>
public void ClusterRecordsOnCID(ref Classes.MapRecord _glycoMap, ref Classes.Params _params)
{
Classes.MapRecord _tempMap = new MapRecord();
_tempMap._AssociatedDatasetNames = _glycoMap._AssociatedDatasetNames;
_tempMap._IsCID = _glycoMap._IsCID;
_tempMap._IsETD = _glycoMap._IsETD;
_tempMap._IsHCD = _glycoMap._IsHCD;
int Id = 0;
FragEvents tempE = new FragEvents();
for (int i = 0; i < _glycoMap._AllMLNRecords.Count; i++)
{
MultiAlignRecord m = new MultiAlignRecord();
m = _glycoMap._AllMLNRecords[i];
int num_records = m._AssociatedUMCRecords.Count;
short mincs = m.MinChargeStateObserved();
short maxcs = m.MaxChargeStateObserved();
for (short thiscs = mincs; thiscs <= maxcs; thiscs++)
{
int spectra_num = 0;
string cum_spectra_names = null;
List<string> cluster_names = new List<string>();
List<int> orphan_ids = new List<int>();
for (int j = 0; j < num_records; j++)
{
if (m._AssociatedUMCRecords[j]._AssociatedFragEvents.Count == 0)
{
// Not fragmented in this UMC but keep track of this but
// // This happens which means there was no fragmentation event then (or) HCD score was bad
if (m._AssociatedUMCRecords[j].Abundance > 0)
orphan_ids.Add(m._AssociatedUMCRecords[j].DatasetID);
}
for (int k = 0; k < m._AssociatedUMCRecords[j]._AssociatedFragEvents.Count; k++)
{
if ((m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].CIDPeaks[0] != null) &&
(m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].TransformResult.mshort_cs == thiscs))
{
string spectra_name = "Spectra_" + m._AssociatedUMCRecords[j].DatasetID + "_" + m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].ID + "_" +
spectra_num + "_" + m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].CIDScan;
GlypID.Peaks.clsPeak[] thisCIDPeaks = new GlypID.Peaks.clsPeak[m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].CIDPeaks.Length];
Array.Copy(m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].CIDPeaks, thisCIDPeaks, m._AssociatedUMCRecords[j]._AssociatedFragEvents[k].CIDPeaks.Length);
_CIDSpectralUtilities.AddPeaksToList(ref thisCIDPeaks, spectra_name);
spectra_num++;
if (cum_spectra_names != null)
cum_spectra_names = cum_spectra_names + "-" + spectra_name;
else
cum_spectra_names = spectra_name;
}
}
}
int num_clusters = _CIDSpectralUtilities.ClusterSpectraInList();
_CIDSpectralUtilities.GetClusterNames(ref cluster_names);
if (num_clusters > 1)
{
// Indicates glycoforms
bool debug = true;
debug = true;
}
for (int c = 0; c < num_clusters; c++)
{
m._ClusterNames.Add(cluster_names[c]);
}
_CIDSpectralUtilities.Clear();
}
_tempMap.AddRecord(m);
}
// Restore them.
_glycoMap.ClearRecords();
for (int i = 0; i < _tempMap._AllMLNRecords.Count; i++)
{
_glycoMap.AddRecord(_tempMap._AllMLNRecords[i]);
}
/*else
{
/* MultiAlignRecord _tempM = new MultiAlignRecord(m);
_tempM.ID = Id;
_tempMap.AddRecord(_tempM);
Id++; */
/* if (cum_spectra_names != null)
{
// Choose the spectra with the greatest SNR
num_clusters = 1;
cluster_names.Add(cum_spectra_names);
_CIDSpectralUtilities.AssignClusters(ref cluster_names);
}
}*/
/* for (int k = 0; k < num_clusters; k++)
{
// ----- Get a representative spectrum index for each cluster ---//
MultiAlignRecord _tempM = new MultiAlignRecord(m);
_tempM.ID = Id;
_tempM._AssociatedUMCRecords.Clear();
GlypID.Peaks.clsPeak[] repCIDPeaks = new GlypID.Peaks.clsPeak[0];
int repOrigIndex = _CIDSpectralUtilities.GetRepresentativePeaksFromCluster(k, ref repCIDPeaks, _params.ScoringParams.MinCIDMz, _params.ScoringParams.MaxCIDMz, true);
if (repCIDPeaks.Length > 1)
{
_tempM._RepresentativeCIDPeaks = repCIDPeaks;
_tempM._RepresentativeDatasetID_CID = repOrigIndex;
}
// Attach UMCs corresponding to that cluster.
List<int> allOrigIDs = new List<int>();
_CIDSpectralUtilities.GetOriginalIDFromCluster(k, ref allOrigIDs);
for (int j = 0; j < num_records; j++)
{
int id = m._AssociatedUMCRecords[j].DatasetID;
if (allOrigIDs.Exists(element => element == id))
{
UMCRecord _tempUMC = new UMCRecord();
_tempUMC = m._AssociatedUMCRecords[j];
_tempM._AssociatedUMCRecords.Add(_tempUMC);
}
else if (orphan_ids.Exists(element => element == id)) // This takes care of non fragmentation but still has stuff present
{
UMCRecord _tempUMC = new UMCRecord();
_tempUMC = m._AssociatedUMCRecords[j];
_tempM._AssociatedUMCRecords.Add(_tempUMC);
}
}
_tempMap.AddRecord(_tempM);
Id++;
}*/
}
