DoubleReal ElutionPeakDetection::computeApexSNR(const MassTrace& tr) { DoubleReal snr(0.0); DoubleReal noise_level(computeMassTraceNoise(tr)); DoubleReal smoothed_apex_int(tr.getMaxIntensity(true)); if (noise_level > 0.0) { snr = smoothed_apex_int/noise_level; } // std::cout << "snr " << snr << " "; return snr; }
double ElutionPeakDetection::computeMassTraceSNR(const MassTrace& tr) { double snr(0.0); if (tr.getSize() > 0) { double noise_area = computeMassTraceNoise(tr) * tr.getTraceLength(); double signal_area = tr.computePeakArea(); snr = signal_area / noise_area; } // std::cout << "snr " << snr << " "; return snr; }