c c DELTAP.FOR (delta**2 test/Main events)- April/May 2001 c implicit double precision (a-h),(o-z) dimension ta(51),ac(51) dimension tc(51),cc(51) open(unit=1,file='deltap15.dat',status='new') open(unit=2,file='crosp15.dat',status='old') open(unit=3,file='autop15.dat',status='old') read(2,*)(tc(i),cc(i),i=1,51) read(3,*)(ta(i),ac(i),i=1,51) do teta=400.D0,440.D0,1.D0 delta=0.D0 adat=0.D0 do 10 i=1,51 do 20 j=1,51 if(ta(j).eq.(tc(i)-teta))delta=delta+(cc(i)-ac(j))**2.D0 if(ta(j).eq.(tc(i)-teta))adat=adat+1.D0 20 continue 10 continue test=delta/adat write(1,100)teta,test/0.0249D0 end do stop 100 format(E10.3,4x,E10.3) end