本文整理汇总了C#中DataReader.GetFramesAndScanIntensitiesForAGivenMz方法的典型用法代码示例。如果您正苦于以下问题:C# DataReader.GetFramesAndScanIntensitiesForAGivenMz方法的具体用法?C# DataReader.GetFramesAndScanIntensitiesForAGivenMz怎么用?C# DataReader.GetFramesAndScanIntensitiesForAGivenMz使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类DataReader
的用法示例。
在下文中一共展示了DataReader.GetFramesAndScanIntensitiesForAGivenMz方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C#代码示例。
示例1: PeakShapeScore
/// <summary>
/// The Peak shape score. Evaluating how "good" the peak looks. A good
/// peak shape score indicates that the peak is not a result of noise
/// or instrument errors. Mostly the feature intensity along is sufficient
/// to exclude noise but a good shape score helps evaluating the experiment
/// and thus the reliability of the data analysis result.
/// </summary>
/// <param name="reader">
/// The reader.
/// </param>
/// <param name="massToleranceInPpm">
/// The mass Tolerance In Ppm.
/// </param>
/// <param name="driftTimeToleranceInMs">
/// The drift Time Tolerance In Scans.
/// </param>
/// <param name="imsPeak">
/// The imsPeak.
/// </param>
/// <param name="voltageGroup">
/// The voltage group.
/// </param>
/// <param name="targetMz">
/// The Target MZ.
/// </param>
/// <param name="globalMaxIntensities">
/// The global Max Intensities.
/// </param>
/// <param name="numberOfScans">
/// The number Of Scans.
/// </param>
/// <returns>
/// The <see cref="double"/>.
/// </returns>
public static double PeakShapeScore(StandardImsPeak imsPeak, DataReader reader, double massToleranceInPpm, double driftTimeToleranceInMs, VoltageGroup voltageGroup, double globalMaxIntensities, int numberOfScans)
{
int scanRep = imsPeak.PeakApex.DriftTimeCenterInScanNumber;
double toleranceInMz = massToleranceInPpm / 1e6 * imsPeak.PeakApex.MzCenterInDalton;
int scanWidth = (int)Math.Ceiling(driftTimeToleranceInMs / 1000 / voltageGroup.AverageTofWidthInSeconds);
int scanWindowSize = scanWidth * 2 + 1;
int scanNumberMin = scanRep - scanWidth;
int scanNumberMax = scanRep + scanWidth;
if ((scanNumberMin < 0) || (scanNumberMax > numberOfScans - 1))
{
return 0;
}
int[][] intensityWindow = reader.GetFramesAndScanIntensitiesForAGivenMz(
voltageGroup.FirstFrameNumber,
voltageGroup.LastFrameNumber,
DataReader.FrameType.MS1,
scanNumberMin,
scanNumberMax,
imsPeak.PeakApex.MzCenterInDalton,
toleranceInMz);
// Average the intensity window across frames
int frames = intensityWindow.GetLength(0);
double[] averagedPeak = new double[scanWindowSize];
double highestPeak = 0;
for (int i = 0; i < scanWindowSize; i++)
{
for (int j = 0; j < frames; j++)
{
averagedPeak[i] += intensityWindow[j][i];
}
averagedPeak[i] /= frames;
highestPeak = (averagedPeak[i] > highestPeak) ? averagedPeak[i] : highestPeak;
}
// For peaks with peak width lower than 3, return a peak score of 0
// TODO get the intensity threshold here from the noise level instead.
if (scanWindowSize >= 3)
{
if (averagedPeak[scanRep - scanNumberMin] < globalMaxIntensities * 0.0001
|| averagedPeak[scanRep - scanNumberMin - 1] < globalMaxIntensities * 0.0001
|| averagedPeak[scanRep - scanNumberMin + 1] < globalMaxIntensities * 0.0001)
{
return 0;
}
}
// Perform a statistical normality test
double normalityScore = NormalityTest.PeakNormalityTest(averagedPeak, NormalityTest.JaqueBeraTest, 100, globalMaxIntensities);
return normalityScore;
}