/** * if there is more than one peak in the histogram, reduce the histogram to only that peak, else * leave unaltered. * * @param hist * @param values * @param valueErrors * @return */ public static HistogramHolder reduceHistogramToFirstPeak( HistogramHolder hist, float[] values, float[] valueErrors) { int yPeakIdx = findFirstPeakIndex(hist); if (yPeakIdx == -1) { return hist; } int yMinPeakIdx = findFirstMinimaFollowingPeak(hist, yPeakIdx); if (yMinPeakIdx == -1) { return hist; } int n = yMinPeakIdx + 1; HistogramHolder tmp = new HistogramHolder(); tmp.setXHist(Arrays.copyOfRange(hist.getXHist(), 0, n)); tmp.setYHistFloat(Arrays.copyOfRange(hist.getYHistFloat(), 0, n)); tmp.setXErrors(Arrays.copyOfRange(hist.getXErrors(), 0, n)); tmp.setYErrors(Arrays.copyOfRange(hist.getYErrors(), 0, n)); return tmp; }
public static HistogramHolder createSimpleHistogram( final float xMin, final float xMax, int nBins, float[] values, float[] valueErrors) { if (values == null || valueErrors == null || values.length != valueErrors.length) { throw new IllegalArgumentException( "values and valueErrors cannot be null and must be the same length"); } float binWidth = calculateBinWidth(xMin, xMax, nBins); float[] xHist = new float[nBins]; int[] yHist = new int[nBins]; Histogram.createHistogram(values, nBins, xMin, xMax, xHist, yHist, binWidth); float[] yHistFloat = new float[yHist.length]; for (int i = 0; i < yHist.length; i++) { yHistFloat[i] = (float) yHist[i]; } float[] yErrors = new float[xHist.length]; float[] xErrors = new float[xHist.length]; calulateHistogramBinErrors(xHist, yHist, values, valueErrors, xErrors, yErrors); HistogramHolder histogram = new HistogramHolder(); histogram.setXHist(xHist); histogram.setYHist(yHist); histogram.setYHistFloat(yHistFloat); histogram.setYErrors(yErrors); histogram.setXErrors(xErrors); return histogram; }
public static HistogramHolder createSimpleHistogram(int binWidth, List<Integer> theValues) { if (theValues == null || theValues.isEmpty()) { throw new IllegalArgumentException( "values and valueErrors cannot be null and must be the same length"); } float[] values = new float[theValues.size()]; for (int i = 0; i < theValues.size(); ++i) { int v = theValues.get(i).intValue(); values[i] = v; } float[] valueErrors = Errors.populateYErrorsBySqrt(values); float[] minMax = MiscMath.calculateOuterRoundedMinAndMax(values); int nBins = (int) Math.ceil(((minMax[1] - minMax[0])) / binWidth); if (nBins < 0) { nBins *= -1; } float[] xHist = new float[nBins]; int[] yHist = new int[nBins]; Histogram.createHistogram(values, nBins, minMax[0], minMax[1], xHist, yHist, binWidth); float[] yHistFloat = new float[yHist.length]; for (int i = 0; i < yHist.length; i++) { yHistFloat[i] = (float) yHist[i]; } float[] yErrors = new float[xHist.length]; float[] xErrors = new float[xHist.length]; calulateHistogramBinErrors(xHist, yHist, values, valueErrors, xErrors, yErrors); HistogramHolder histogram = new HistogramHolder(); histogram.setXHist(xHist); histogram.setYHist(yHist); histogram.setYHistFloat(yHistFloat); histogram.setYErrors(yErrors); histogram.setXErrors(xErrors); return histogram; }
public static HistogramHolder createSimpleHistogram( float binWidth, float[] values, float[] valueErrors) { if (values == null || valueErrors == null || values.length != valueErrors.length) { throw new IllegalArgumentException( "values and valueErrors cannot be null and must be the same length"); } float[] minMax = MiscMath.calculateOuterRoundedMinAndMax(values); int nBins = (int) Math.ceil(((minMax[1] - minMax[0])) / binWidth); if (nBins < 0) { nBins *= -1; } float[] xHist = new float[nBins]; int[] yHist = new int[nBins]; Histogram.createHistogram(values, nBins, minMax[0], minMax[1], xHist, yHist, binWidth); float[] yHistFloat = new float[yHist.length]; for (int i = 0; i < yHist.length; i++) { yHistFloat[i] = (float) yHist[i]; } float[] yErrors = new float[xHist.length]; float[] xErrors = new float[xHist.length]; calulateHistogramBinErrors(xHist, yHist, values, valueErrors, xErrors, yErrors); HistogramHolder histogram = new HistogramHolder(); histogram.setXHist(xHist); histogram.setYHist(yHist); histogram.setYHistFloat(yHistFloat); histogram.setYErrors(yErrors); histogram.setXErrors(xErrors); return histogram; }
public static HistogramHolder calculateSturgesHistogramRemoveZeroTail( float[] values, float[] valueErrors) { if (values == null || valueErrors == null || values.length != valueErrors.length) { throw new IllegalArgumentException( "values and valueErrors cannot be null and must be the same length"); } int nIntervalsSturges = (int) Math.ceil(Math.log(values.length) / Math.log(2)); // int nItervalsRice = (int)(2*Math.pow(values.length, 0.3333)); int nBins = 25; if (values.length > 10000) { nBins = 40; } nBins = Math.max(nIntervalsSturges, nBins); float[] xHist = new float[nBins]; int[] yHist = new int[nBins]; float minx = MiscMath.findMin(values); float maxx = MiscMath.findMax(values); float binWidth = calculateBinWidth(minx, maxx, nBins); Histogram.createHistogram(values, nBins, minx, maxx, xHist, yHist, binWidth); float maxy = MiscMath.findMax(yHist); int minCountsLimit = (int) Math.max(5, 0.03f * maxy); int countsBelowMinAtTail = 0; int lastLowCountIdx = yHist.length - 1; for (int i = (yHist.length - 1); i > -1; i--) { if (yHist[i] < minCountsLimit) { countsBelowMinAtTail++; lastLowCountIdx = i; } else { break; } } if (countsBelowMinAtTail > 0) { maxx = xHist[lastLowCountIdx]; // keep nbins the same? binWidth = calculateBinWidth(minx, maxx, nBins); Histogram.createHistogram(values, nBins, minx, maxx, xHist, yHist, binWidth); if (countsBelowMinAtTail > (nBins >> 1)) { // one more round of trimming maxy = MiscMath.findMax(yHist); minCountsLimit = (int) Math.max(5, 0.03f * maxy); countsBelowMinAtTail = 0; lastLowCountIdx = yHist.length - 1; for (int i = (yHist.length - 1); i > -1; i--) { if (yHist[i] < minCountsLimit) { countsBelowMinAtTail++; lastLowCountIdx = i; } else { break; } } if (countsBelowMinAtTail > 0) { maxx = xHist[lastLowCountIdx]; binWidth = calculateBinWidth(minx, maxx, nBins); Histogram.createHistogram(values, nBins, minx, maxx, xHist, yHist, binWidth); } } } if (values.length > 100) { // if there are a large number of points, we'd like to increase the // resolution of the peak if needed int nLeftOfPeak = MiscMath.findYMaxIndex(yHist); int nIter = 0; while (nIter < 30 && nLeftOfPeak < 3 && (yHist[nLeftOfPeak] > 100)) { binWidth *= 0.8f; Histogram.createHistogram(values, nBins, minx, maxx, xHist, yHist, binWidth); nLeftOfPeak = MiscMath.findYMaxIndex(yHist); nIter++; } } float[] yHistFloat = new float[yHist.length]; for (int i = 0; i < yHist.length; i++) { yHistFloat[i] = (float) yHist[i]; } float[] yErrors = new float[xHist.length]; float[] xErrors = new float[xHist.length]; calulateHistogramBinErrors(xHist, yHist, values, valueErrors, xErrors, yErrors); HistogramHolder histogram = new HistogramHolder(); histogram.setXHist(xHist); histogram.setYHist(yHist); histogram.setYHistFloat(yHistFloat); histogram.setYErrors(yErrors); histogram.setXErrors(xErrors); return histogram; }