# -*- coding: windows-1251 -*-

# Calculates delays from FIX 35=D,F,G to FIX 35=8 150=0,4,6,5 replies and builds latency probability distribution
# uses tshark output with FIX packets dissections 
# python 2.x and 3.x compatible on Windows and Linux 

#==========================================
from ctypes import *
import time
import sys
import os
import socket
import random
import math
import string
import re
import array


#===========================================
def simplecount(filename):
    lines= 0
    for line in open(filename):
        lines += 1
    return lines


ord_file_name =sys.argv[1]
out_fname=ord_file_name + ".log"
stat_fname=ord_file_name + ".stats"
out_file=open(out_fname,"w+")
stat_file=open(stat_fname,"w+")


LCount=simplecount(ord_file_name)
ord_file = open(ord_file_name)
record= ord_file.readlines()
ord_file.close()
st=[]
b=[]
dd=0
ddd=0
bg=20
distr=[]

lv=1000  # number of steps in distribution .stats file
step=5 # bin size in microseconds

while ddd<lv+step:
 distr.append(0)
 ddd=ddd+1
ii=1
cnt=0
tottime=0
delay=0

while ii< LCount-120:
    if record[ii].split(chr(9))[3] == "D" or record[ii].split(chr(9))[3] == "F" or record[ii].split(chr(9))[3] == "G" :
        iii=ii
        while iii<ii+110:
           iii=iii+1
           if record[iii].split(chr(9))[4] == record[ii].split(chr(9))[4] and record[ii].split(chr(9))[5]==record[iii].split(chr(9))[9] and (record[iii].split(chr(9))[3]=="9" or (record[iii].split(chr(9))[3]=="8" and (record[iii].split(chr(9))[6]=="0" or record[iii].split(chr(9))[6]=="6" or record[iii].split(chr(9))[6]=="8" or record[iii].split(chr(9))[6]=="4" or record[iii].split(chr(9))[6]=="5"))):
              delay=int(1000000/step*(float(record[iii].split(chr(9))[2]) - float(record[ii].split(chr(9))[2])) )*step
              cnt=cnt+1
              st.append(delay)
              tottime=tottime+delay
              if delay<lv*step+bg+step:
                  distr[int((delay-bg)/step)]=1+distr[int((delay-bg)/step)]
              else:
                  distr[lv+1]=distr[lv+1]+1

              out_file.write(str(delay)+chr(9)+str(record[ii].split(chr(9))[0])+chr(9)+str(iii-ii) +  chr(9)+'35='+str(record[ii].split(chr(9))[3])+ chr(9)+ str(record[ii].split(chr(9))[5]) + chr(9)+ str(record[ii].split(chr(9))[1]) +chr(9)+ str(record[iii].split(chr(9))[7]) + chr(9) + str(record[iii].split(chr(9))[8])+chr(9)+str(record[iii].split(chr(9))[9])+chr(9)+str(record[iii].split(chr(9))[2])+chr(9)+str(record[ii].split(chr(9))[2])+chr(9)+str(record[iii].split(chr(9))[10]))
              iii=ii+200				  

        
    ii=ii+1  
	
b=sorted(st)
print (" ")
print ("count="+chr(9)+str(cnt))
print ("avg="+chr(9)+str(int(tottime/cnt)))
print ("peak="+chr(9)+str(distr.index(max(distr))*step+bg))
print ("50%  "+chr(9)+str(b[int(0.5*cnt)]))
print ("80%  "+chr(9)+str(b[int(0.8*cnt)]))
print ("90%  "+chr(9)+str(b[int(0.90*cnt)]))
print ("99%  "+chr (9)+str(b[int(0.99*cnt)]))
print ("99.9%"+chr(9)+str(b[int(0.999*cnt)]))

stat_file.write( "count="+chr(9)+ str(cnt)+chr(10))
stat_file.write( "avg="+chr(9)+str(int(tottime/cnt))+chr(10))
stat_file.write("peak probability at" + chr(9)+str(distr.index(max(distr))*step+bg)+" microseconds"+chr(10))
stat_file.write( "50%   "+ chr(9)+str( b[int(0.5*cnt)])+chr(10))
stat_file.write( "80%   "+ chr(9)+str(b[int(0.8*cnt)])+chr(10))
stat_file.write( "90%   "+ chr(9)+str( b[int(0.90*cnt)])+chr(10))
stat_file.write( "99%   "+ chr (9)+ str(b[int(0.99*cnt)])+chr(10))
stat_file.write( "99.9% "+ chr(9)+str(b[int(0.999*cnt)])+chr(10))


stat_file.write(chr(10))
stat_file.write(chr(10))
stat_file.write("delay"+chr(9)+"FIX Ack delay probability"+chr(10))

ddd=0
while ddd<lv+2:
  stat_file.write(str(ddd*step+bg)+chr(9)+str(int(distr[ddd]*100000/cnt)/1000.0) + chr(10))
  ddd=ddd+1