Page MenuHomeHEPForge
Files F11221464

get_chisq.f90
No OneTemporary
Actions

get_chisq.f90
View Options

program getchisq
!gfortran chi2-CL-stuff_fromOliver/chi2fromCL.f get_chisq.f90 -o get_chisq.exe
implicit none
integer :: ifault
double precision :: ppchi2
double precision :: chisq,clsb
clsb=1.0D-6
if(abs(1.0D0-2.0D0*clsb) .le. 0.000002D0) then ! 0.499999 <clsb< 0.500001
chisq=0.0D0!ppchi2 not valid: see comment at the beginning of function ppchi2
elseif(abs(1.0D0-2.0D0*clsb) .ge. 0.999998D0)then !|clsb| <=1.0D-6
chisq=1.0D5!ppchi2 not valid: see comment at the beginning of function ppchi2
elseif(clsb.le. -1.0D-3) then
chisq= -1.0D0
elseif(clsb .le. 0.5D0) then
chisq=ppchi2(1.0D0-2.0D0*clsb,1.0D0,ifault)
else
chisq=ppchi2(1.0D0-2.0D0*(1.0D0-clsb),1.0D0,ifault)
endif
write(*,*)'clsb=',clsb
write(*,*)'chisq=',chisq
end program getchisq

File Metadata

Mime Type
text/plain
Expires
Wed, May 14, 10:28 AM (1 d, 20 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5090861
Default Alt Text
get_chisq.f90 (820 B)

Event Timeline