Page MenuHomeHEPForge
Files F11221686

example_functions.F90
No OneTemporary
Actions

example_functions.F90
View Options

! This file is part of HiggsBounds
!******************************************************
program example_functions
#define NHIGGS 3
#define NHIGGSPLUS 1
#ifdef NAGf90Fortran
use F90_UNIX_ENV, only : iargc,getarg
use F90_UNIX_IO, only : flush
#endif
!******************************************************
implicit none
integer :: nH,nHplus,HBresult,chan,ncombined,i,j,Mhtemp,x
double precision :: obsratio,mass
double precision :: SMGammaTotal(NHIGGS),SMGamma_hgg(NHIGGS),Gamma_hgg(NHIGGS),SMGamma_h(NHIGGS), &
& SMBR_Hgg(NHIGGS),FGMGammaTotal(NHIGGS)
double precision :: Mh(NHIGGS),GammaTotal(NHIGGS), &
& BR_hjinvisible(NHIGGS),BR_hjhihi(NHIGGS,NHIGGS)
double precision :: CS_lep_hjZ_ratio(NHIGGS), &
& CS_lep_bbhj_ratio(NHIGGS),CS_lep_tautauhj_ratio(NHIGGS), &
& CS_lep_hjhi_ratio_nHbynH(NHIGGS,NHIGGS), &
& CS_gg_hj_ratio(NHIGGS),CS_bb_hj_ratio(NHIGGS), &
& CS_bg_hjb_ratio(NHIGGS), &
& CS_ud_hjWp_ratio(NHIGGS),CS_cs_hjWp_ratio(NHIGGS), &
& CS_ud_hjWm_ratio(NHIGGS),CS_cs_hjWm_ratio(NHIGGS), &
& CS_gg_hjZ_ratio(NHIGGS), &
& CS_dd_hjZ_ratio(NHIGGS),CS_uu_hjZ_ratio(NHIGGS), &
& CS_ss_hjZ_ratio(NHIGGS),CS_cc_hjZ_ratio(NHIGGS), &
& CS_bb_hjZ_ratio(NHIGGS), &
& CS_tev_vbf_ratio(NHIGGS),CS_tev_tthj_ratio(NHIGGS), &
& CS_lhc7_vbf_ratio(NHIGGS),CS_lhc7_tthj_ratio(NHIGGS), &
& BR_hjss(NHIGGS),BR_hjcc(NHIGGS), &
& BR_hjbb(NHIGGS),BR_hjmumu(NHIGGS),BR_hjtautau(NHIGGS), &
& BR_hjWW(NHIGGS),BR_hjZZ(NHIGGS), &
& BR_hjZga(NHIGGS), BR_hjgaga(NHIGGS),BR_hjgg(NHIGGS)
integer ::CP_value
character(LEN=5) :: wchinput
character(len=100) :: filename
wchinput='SLHA'
nH=NHIGGS
nHplus=NHIGGSPLUS
call initialize_HiggsBounds(nH,nHplus,'LandH')
!do i=2,1300
do i=1,1
Mh=dble(i-1)/10.0D0
!write(*,*)'*******hello******'
!write(*,*)'hello Mh(1)=',Mh(1)
CP_value=1
CS_lep_hjZ_ratio = 0.0D0
CS_lep_bbhj_ratio=0.0D0
CS_lep_tautauhj_ratio=0.0D0
CS_lep_hjhi_ratio_nHbynH=0.0D0
CS_gg_hj_ratio=0.0D0
CS_bb_hj_ratio=0.0D0
CS_bg_hjb_ratio=0.0D0
CS_ud_hjWp_ratio=0.0D0
CS_cs_hjWp_ratio=0.0D0
CS_ud_hjWm_ratio=0.0D0
CS_cs_hjWm_ratio=0.0D0
CS_gg_hjZ_ratio=0.0D0
CS_dd_hjZ_ratio=0.0D0
CS_uu_hjZ_ratio=0.0D0
CS_ss_hjZ_ratio=0.0D0
CS_cc_hjZ_ratio=0.0D0
CS_bb_hjZ_ratio=0.0D0
CS_tev_vbf_ratio=0.0D0
CS_tev_tthj_ratio=0.0D0
CS_lhc7_vbf_ratio=0.0D0
CS_lhc7_tthj_ratio=0.0D0
BR_hjss=0.0D0
BR_hjcc=0.0D0
BR_hjbb = 0.0D0
BR_hjmumu=0.0D0
BR_hjtautau=0.0D0
BR_hjWW=0.0D0
BR_hjZZ=0.0D0
BR_hjZga=0.0D0
BR_hjgaga=0.0D0
BR_hjgg=0.0D0
BR_hjinvisible=0.0D0
BR_hjhihi=0.0D0
Mh(2)=Mh(1)*2.0D0+2.0D0
BR_hjbb(1) = 1.0D0
BR_hjhihi(2,1)=1.0D0
CS_lep_hjhi_ratio_nHbynH(2,1)=0.1D0
CS_lep_hjhi_ratio_nHbynH(1,2)=0.1D0
do j=1,nH
Mhtemp=Mh(j)
SMGammaTotal=SMGamma_h(Mhtemp)
enddo
!write(*,*)'hello about to give input to HB' ; call flush(6)
select case(wchinput)
case('part')
call HiggsBounds_neutral_input_part(Mh,SMGammaTotal,CP_value, &
& CS_lep_hjZ_ratio, &
& CS_lep_bbhj_ratio,CS_lep_tautauhj_ratio, &
& CS_lep_hjhi_ratio_nHbynH, &
& CS_gg_hj_ratio,CS_bb_hj_ratio, &
& CS_bg_hjb_ratio, &
& CS_ud_hjWp_ratio,CS_cs_hjWp_ratio, &
& CS_ud_hjWm_ratio,CS_cs_hjWm_ratio, &
& CS_gg_hjZ_ratio, &
& CS_dd_hjZ_ratio,CS_uu_hjZ_ratio, &
& CS_ss_hjZ_ratio,CS_cc_hjZ_ratio, &
& CS_bb_hjZ_ratio, &
& CS_tev_vbf_ratio,CS_tev_tthj_ratio, &
& CS_lhc7_vbf_ratio,CS_lhc7_tthj_ratio, &
& BR_hjss,BR_hjcc, &
& BR_hjbb,BR_hjmumu,BR_hjtautau, &
& BR_hjWW,BR_hjZZ,BR_hjZga, BR_hjgaga,BR_hjgg, &
& BR_hjinvisible,BR_hjhihi )
case('SLHA')
filename='../example_data/example_slha_in.txt'
call HiggsBounds_input_SLHA(filename)
case default
stop'problem here in example_chisq'
end select
!write(*,*)'hello about to call run_HiggsBounds' ; call flush(6)
call run_HiggsBounds( HBresult,chan, &
& obsratio, ncombined )
write(*,*)'HBresult,chan,' &
& //'obsratio, ncombined:',HBresult,chan, &
& obsratio, ncombined
call flush(19)
enddo
call finish_HiggsBounds
end program example_functions

File Metadata

Mime Type
text/plain
Expires
Wed, May 14, 10:46 AM (1 d, 8 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5067243
Default Alt Text
example_functions.F90 (5 KB)

Event Timeline