Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-11-25 02:30:05

0001 import sys
0002 if len(sys.argv) !=3 :
0003   print("usage:" + sys.argv[0] + "   bla bla")
0004   exit
0005 # go
0006 mm=0.
0007 n=0
0008 ave=0
0009 tot=0
0010 ln=0
0011 ml=0
0012 file1 = open(sys.argv[1])
0013 file2 = open(sys.argv[2])
0014 for l1 in file1.readlines() :
0015     ln+=1
0016     l2 = file2.readline()
0017     a1 =l1.split('(')
0018     a2 = l2.split('(')
0019     if len(a1)<2 :
0020         continue
0021     p1 = a1[1].split(',')
0022     p1[2] = p1[2].rstrip(') \n')
0023     p2 = a2[1].split(',')
0024     p2[2] = p2[2].rstrip(') \n')
0025     f1 = [float(i) for i in p1]
0026     f2 = [float(i) for i in p2]
0027     diff = []
0028     tot+=1
0029     for i in range(len(f1)):
0030         diff.append(f1[i]-f2[i])
0031     for x in diff:
0032         if (abs(x)>mm):
0033             mm=x
0034             ml=ln
0035         if (x!=0) :
0036             n+=1
0037             ave+=abs(x)
0038 print(tot, mm, n, ave/n)
0039 #
0040 ml = ml-5
0041 file1.seek(0)
0042 file2.seek(0)
0043 l1 = file1.readlines()
0044 l2 = file2.readlines()
0045 for i in range(10) :
0046     print(l1[ml+i].rstrip('\n'))
0047     print(l2[ml+i].rstrip('\n'))