Page MenuHomeHEPForge
Files F8309735

test_pc_chisq.f90
No OneTemporary
Actions

test_pc_chisq.f90
View Options

program test_pc_chisq
use pc_chisq
implicit none
integer, allocatable :: indices_best(:), indices_tr(:)
double precision :: pccsq, prob
integer :: imin, imax, ncomb,ii,jj
Smin=2.0d0
pdf=3
call setup_mutables
! Look at 1 implemented analysis first.
imin = lbound(mutables(1)%mu,dim=1)
imax = ubound(mutables(1)%mu,dim=1)
call findpeaks(mutables(1), Smin, peaks)
! Model input -> Higgses. Here, we set them explicitly:
allocate(Higgses(2))
! Higgses(1)%m=120.0D0
Higgses(1)%dm=2.0D0
Higgses(1)%mu=0.907D0
! Higgses(2)%m=128.0D0
Higgses(2)%dm=2.0D0
Higgses(2)%mu=0.907D0
open(21, file='pc_csq_results.dat')
write(21,*) '#--------------------------------------------------------------------------'
write(21,*) '# mh1 dmh1 mu1 mh2 dmh2 mu2 chisq prob ncomb indices'
write(21,*) '#--------------------------------------------------------------------------'
if(allocated(peaks)) then
! pc_chisq(peaks(1), Higgses,pdf, indices_best, Nindices_in, indices_in)
! This finds the best combination automatically.
do ii=1,101
Higgses(1)%m=120.0D0+(ii-1)*0.1
do jj=1,101
Higgses(2)%m=120.0D0+(jj-1)*0.1
call calc_pc_chisq(pccsq, peaks(1), Higgses,pdf, indices_best)
prob = 0.5D0*(1 - gammp(dble(1)/2,pccsq/2))
call get_ncomb(ncomb, indices_tr, indices_best)
write(21,*) Higgses(1)%m, Higgses(1)%dm, Higgses(1)%mu, &
& Higgses(2)%m, Higgses(2)%dm, Higgses(2)%mu, &
& pccsq, prob, ncomb, indices_tr
! write(*,*) 'Peak 1 was fitted with', ncomb,' Higgs boson(s) with indices:',indices_tr
! write(*,*) 'with a chi^2 value of', pccsq
! write(*,*) 'the corresponding probability is',prob
deallocate(indices_best, indices_tr)
enddo
enddo
endif
close(21)
end program test_pc_chisq

File Metadata

Mime Type
text/plain
Expires
Sat, Dec 21, 4:21 PM (21 h, 7 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023427
Default Alt Text
test_pc_chisq.f90 (1 KB)

Event Timeline