Exemplo n.º 1
0
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;
  }