c	
c	DELTAS.FOR (delta**2 test/Secondary events)- April/May 2001
c
	implicit double precision (a-h),(o-z)
      dimension ta(71),ac(71)
	dimension tc(71),cc(71)   
	open(unit=1,file='deltas15.dat',status='new')
	open(unit=2,file='cross15.dat',status='old')
	open(unit=3,file='autos15.dat',status='old')
	read(2,*)(tc(i),cc(i),i=1,71)
	read(3,*)(ta(i),ac(i),i=1,71)
	do teta=400.D0,440.D0,1.D0
	delta=0.D0
	adat=0.D0
	do 10 i=1,71
	do 20 j=1,71
	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.0231D0
	end do
	stop
100	format(E10.3,4x,E10.3)
	end