Page MenuHomeHEPForge
Files F11222119

calculate_ratios.f90
No OneTemporary
Actions

calculate_ratios.f90
View Options

program calculate_ratios
! Simple program for transforming upper limits on Cross Section x BR
! into upper limits normalised to the SM.
! SM cross section and BRs are taken from one table, observed and expected
! upper limits on Cross Section x BR are given in the other table.
! Compiled with
! gfortran calculate_ratios.f90 -o calculate_ratios.exe
implicit none
character(LEN=150) :: description
double precision :: Mh_in1, Mh_in2, exp, obs, xsec, BR, expSMnorm, obsSMnorm
integer :: i,skip,lines
skip=5
lines=44
open(10,file='../../docs/Expt/10433_files/CDF_ggH-WW_7.1fb_10433.txt')
open(11,file='../../docs/Expt/10433_files/CDF_ggH-WW_7.1fb_10433xsec.txt')
open(14,file='../../docs/Expt/10433_files/CDF_ggH-WW_7.1fb_10433_SMnormalized.txt')
do i=1,lines
write(*,*) i
if(i.le.skip) then
read(10,'(a)') description
read(11,*)
write(14,*) description
else
read(10,*)Mh_in1, obs, exp
read(11,*)Mh_in2, xsec, BR
if(abs(Mh_in1-Mh_in2).gt.1.0D-10) stop'Higgs masses do not match'
obsSMnorm = obs/(xsec*BR)*1.0D-03
expSMnorm = exp/(xsec*BR)*1.0D-03
write(14,'(1D16.3,2D26.5)')Mh_in1, obsSMnorm, expSMnorm
endif
enddo
close(10)
close(11)
close(14)
end program calculate_ratios

File Metadata

Mime Type
text/plain
Expires
Wed, May 14, 11:26 AM (14 h, 36 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5068972
Default Alt Text
calculate_ratios.f90 (1 KB)

Event Timeline