Package net.sf.mzmine.util

Examples of net.sf.mzmine.util.Range


    // <PEAK>
    if (qName.equals(PeakListElementName_2_0.PEAK.getElementName())) {

      DataPoint[] mzPeaks = new DataPoint[numOfMZpeaks];
      Range peakRTRange = null, peakMZRange = null, peakIntensityRange = null;
      RawDataFile dataFile = dataFilesIDMap.get(peakColumnID);

      if (dataFile == null)
        throw new SAXException("Error in project: data file "
            + peakColumnID + " not found");

      for (int i = 0; i < numOfMZpeaks; i++) {

        Scan sc = dataFile.getScan(scanNumbers[i]);
        double retentionTime = sc.getRetentionTime();

        double mz = masses[i];
        double intensity = intensities[i];

        if ((peakRTRange == null) || (peakIntensityRange == null)) {
          peakRTRange = new Range(retentionTime);
          peakIntensityRange = new Range(intensity);
        } else {
          peakRTRange.extendRange(retentionTime);
          peakIntensityRange.extendRange(intensity);
        }
        if (mz > 0.0) {
          mzPeaks[i] = new SimpleDataPoint(mz, intensity);
          if (peakMZRange == null)
            peakMZRange = new Range(mz);
          else
            peakMZRange.extendRange(mz);
        }
      }
View Full Code Here


  }

  @Nonnull
  @Override
  public Range getRawDataPointsRTRange() {
    return new Range(retentionTime, retentionTime);
  }
View Full Code Here

  }

  @Nonnull
  @Override
  public Range getRawDataPointsMZRange() {
    return new Range(ionMass, ionMass);
  }
View Full Code Here

  }

  @Nonnull
  @Override
  public Range getRawDataPointsIntensityRange() {
    return new Range(0, 0);
  }
View Full Code Here

        String name = FameData.FAME_NAMES[i];

        // Check for ion qualifier
        int qualifier = FameData.QUALIFIER_IONS[i];

        DataPoint[] p = s.getDataPointsByMass(new Range(qualifier,
            qualifier));

        // Confirm that the qualifier ion exists
        if (p.length != 1)
          continue;
View Full Code Here

    if (data.size() > n) {
      Collections.sort(data, new AlignStructMol());
      double min = data.get(0).RT;
      double max = data.get(data.size() - 1).RT;

      Range rtRange = new Range(min, ((max - min) / 2) + min);

      int cont = 0, bucle = 0;
      while (cont < n / 2 && bucle < 1000) {
        int index = (int) (data.size() * Math.random());
        if (!data.get(index).ransacMaybeInLiers
            && rtRange.contains(data.get(index).RT)) {
          data.get(index).ransacMaybeInLiers = true;
          cont++;

        }

        bucle++;
      }
      if (bucle >= 1000) {
        getN(data, (n / 2) - cont);
      }

      bucle = 0;
      rtRange = new Range(((max - min) / 2) + min, max);

      while (cont < n && bucle < 1000) {

        int index = (int) (data.size() * Math.random());
        if (!data.get(index).ransacMaybeInLiers
            && rtRange.contains(data.get(index).RT)) {
          data.get(index).ransacMaybeInLiers = true;
          cont++;
        }
        bucle++;
      }
View Full Code Here

    Map<Integer, Map<Range, List<MassCandidate>>> massRanges = new HashMap<Integer, Map<Range, List<MassCandidate>>>();

    for (List<MassCandidate> masses : massCandidates.values()) {
      for (MassCandidate m : masses) {
        int mass = m.getIonMass();
        Range range = getRIWindow(m);

        // Create a new HashMap entry with a simple Range Comparator if
        // this is a new mass
        if (!massRanges.containsKey(mass)) {
          massRanges.put(mass,
              new TreeMap<Range, List<MassCandidate>>(
                  new RangeComparator()));

          totalScans += 2;
        }

        // Otherwise, search through all ranges
        else {
          Map<Range, List<MassCandidate>> map = massRanges.get(mass);
          Iterator<Map.Entry<Range, List<MassCandidate>>> iter = map
              .entrySet().iterator();

          boolean matched = false;

          while (iter.hasNext()) {
            Map.Entry<Range, List<MassCandidate>> e = iter.next();

            // If the range overlaps, remove the old key, extend it,
            // and re-add the list with the new key
            if (rangeOverlaps(e.getKey(), range)) {
              List<MassCandidate> tmpList = e.getValue();
              tmpList.add(m);

              Range matchedRange = e.getKey();
              map.remove(e.getKey());

              matchedRange.extendRange(range);
              map.put(matchedRange, tmpList);

              matched = true;
              break;
            }
          }

          // If no match is found, add another range entry
          if (!matched) {
            List<MassCandidate> tmpList = new ArrayList<MassCandidate>();
            tmpList.add(m);

            massRanges.get(mass).put(range, tmpList);
          }
        }
      }

      processedScans++;
    }

    // Combine any possible overlapping ranges
    for (Integer mass : massRanges.keySet()) {
      Map<Range, List<MassCandidate>> map = massRanges.get(mass);
      Range[] keys = map.keySet().toArray(new Range[map.size()]);

      for (int i = 0; i < keys.length - 1;) {
        if (rangeOverlaps(keys[i], keys[i + 1])) {
          List<MassCandidate> tmpList = map.get(keys[i]);
          tmpList.addAll(map.get(keys[i + 1]));

          Range tmpRange = new Range(keys[i]);
          tmpRange.extendRange(keys[i + 1]);

          map.remove(keys[i]);
          map.remove(keys[i + 1]);
          map.put(tmpRange, tmpList);
View Full Code Here

   * @param rt
   *            retention time
   * @return time window about given retention time
   */
  private Range getTimeWindow(double rt) {
    return new Range(rt - timeWindow, rt + timeWindow);
  }
View Full Code Here

   */
  private Range getRIWindow(MassCandidate m) {
    CombinedRegression fit = m.getFit();
    double rt = m.getOriginalRetentionTime();

    return new Range(fit.getY(rt - timeWindow), fit.getY(rt + timeWindow));
  }
View Full Code Here

    // Tine window in minutes
    double window = 0.5;

    TICVisualizerModule.showNewTICVisualizerWindow(
        resultsList.getRawDataFiles(), null, null, 1,
        PlotType.BASEPEAK, new Range(rt - window, rt + window),
        new Range(0, 600));
  }
View Full Code Here

TOP

Related Classes of net.sf.mzmine.util.Range

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.