source: fact/tools/pyscripts/sandbox/dneise/compare_mapping/compare

#!/usr/bin/python -tti
#
# Dominik Neise
#
import sys
import numpy as np
from plotters import Plotter

def hard2chid(hardid):
    hardid = int(hardid)
    einer = hardid%10
    hardid -= einer
    zehner = hardid%100/10
    hardid -= zehner*10
    hunderter = hardid%1000/100
    hardid -= hunderter*100
    tausender = hardid%10000/1000

    chid = tausender*360 + hunderter*36 + zehner*9 + einer
    return chid

leftpath = sys.argv[1]
rightpath =sys.argv[2]

left = np.loadtxt(leftpath)
right = np.loadtxt(rightpath)

right = list(right)
for i,ri in enumerate(right):
    right[i] = list(ri)

for i in right:
    i.append( hard2chid(i[1]) )

right = sorted( right, key=lambda i: i[-1])

right = np.array(right)

pleft = Plotter('pleft', x=left[:,1])
pright= Plotter('pright',x=right[:,8])
pleft(left[:,2])
pright(-right[:,9])