#!/usr/bin/env python

import optparse

import argweaver

from rasmus import intervals
from rasmus import stats
from rasmus import util


o = optparse.OptionParser()
conf, args = o.parse_args()


def iter_track(filename):
    for row in util.DelimReader(filename):
        yield int(row[0]), int(row[1]), float(row[2])

#=============================================================================

filenames = args

tracks = []
for filename in filenames:
    tracks.extend(iter_track(filename))
tracks.sort()

for start, end, group in intervals.iter_intersections(tracks):
    vals = util.cget(group, 2)
    util.print_row(start, end, stats.mean(vals),
                   stats.percentile(vals, .025), stats.percentile(vals, .975))
