#!/usr/bin/python
# compare_stars.py
# Jure Skvarc, May 2000
#
from string import *
import fileinput

class Star:
    pass

stars = {}
for a in fileinput.input():
    row = split(a)
    intensity = float(row[3])
    star_id = row[7]
    if stars.has_key(star_id):
        stars[star_id].mag.append(intensity)
        stars[star_id].catmag = float(row[9])
    else:
	stars[star_id] = Star()
	stars[star_id].mag = []
        stars[star_id].mag.append(intensity)

w = [0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 1000]
ne = [0, 0, 0, 0, 0, 0, 0]
for a in stars.keys():
    if len(stars[a].mag) > 1:
	e = abs(stars[a].mag[0] - stars[a].mag[1]) / stars[a].mag[0]
	print a, stars[a].mag[0], stars[a].mag[1], e, stars[a].catmag
	for i in range(0, len(w)):
	    if e < w[i]:
		ne[i] = ne[i] + 1


fp = open('errors.dat', 'a')
for i in range(0, len(ne)):
    fp.write('%.3f %4d\n' %  (w[i], ne[i]))
fp.write('\n')
fp.close()