program fom_calc
c
      implicit none
c
      real*8 ac
      real*8 eff,th,efferr,sumeff,sumac2eff
      integer i
c
      sumeff = 0
      sumac2eff = 0
      open(unit=1,file='eff_ch2.dat',status='old')
      read(1,*) th,eff,efferr
      do i=1,24
      read(1,*) th,eff,efferr
c ch2
      ac = -0.01453+0.04803*th-0.005183*th*th
      ac = ac + 0.0001872*th*th*th-2.093e-06*(th)**4
      sumeff = sumeff + eff
      sumac2eff = sumac2eff + ac*ac*eff
      write(*,'(6g15.6)') th,eff,ac, ac*ac*eff,sumeff
      enddo
      write(*,*) ' eff = ',sumeff
      write(*,*) ' Ac^2 = ',sumac2eff/sumeff
      write(*,*) ' Ac  = ',sqrt(sumac2eff/sumeff)
c
      close(1)
      sumeff = 0
      sumac2eff = 0
c
      write(*,*) ' carbon'
      open(unit=1,file='eff_car.dat',status='old')
      read(1,*) th,eff,efferr
      do i=1,24
      read(1,*) th,eff,efferr
c carbon
      ac = 0.09623 -9.629e-05*th-0.0001892*th*th
      sumeff = sumeff + eff
      sumac2eff = sumac2eff + ac*ac*eff
      write(*,'(6g15.6)') th,eff,ac, ac*ac*eff,sumeff
      enddo
      write(*,*) ' eff = ',sumeff
      write(*,*) ' Ac^2 = ',sumac2eff/sumeff
      write(*,*) ' Ac  = ',sqrt(sumac2eff/sumeff)
c
      end