Page MenuHomeHEPForge
Files F11221615

SM_effC.F
No OneTemporary
Actions

SM_effC.F
View Options

! This file is part of HiggsBounds
! -OS
!
! gfortran -L../HiggsBounds_KW/ -I../HiggsBounds_KW/ SM_list.F -o SM_list.exe -lHB
!******************************************************
program SM_list
!******************************************************
use usefulbits, only : analysislist
implicit none
integer nHzero,nHplus,HBresult,chan,ncombined,i
parameter (nHzero = 1)
parameter (nHplus = 0)
real*8 obsratio,mass
real*8 SMGammaTotal,SMGamma_hgg,Gamma_hgg,SMGamma_h,
& SMBR_Hgg,FGMGammaTotal
real*8 Mh,GammaTotal,g2hjss_s,g2hjss_p,g2hjcc_s,g2hjcc_p,
& g2hjbb_s,g2hjbb_p,g2hjtt_s,g2hjtt_p,
& g2hjmumu_s,g2hjmumu_p,g2hjtautau_s,g2hjtautau_p,
& g2hjWW,g2hjZZ,g2hjZga,g2hjgaga,g2hjgg,g2hjggZ,
& g2hjhiZ,BR_hjhihi,BR_hjinvisible
double precision R
c used by set_mass_uncertainties
double precision dmhneut(nHzero)
double precision dmhch(nHplus)
! Allocate list of analyses ID numbers to include
! and initialize with bogus (-1) numbers
allocate(analysislist(10))
do i=1,ubound(analysislist,1)
analysislist(i) = -1
enddo
analysislist(1) = 12045
R = 1.0d0
c Call init with new 'list ' setting. Note stupid trailing space.
c call initialize_HiggsBounds(nHzero,nHplus,'list ')
c call initialize_HiggsBounds(nHzero,nHplus,'onlyL')
call initialize_HiggsBounds(nHzero,nHplus,'LandH')
c call initialize_HiggsBounds(nHzero,nHplus,'onlyP')
open(21,file='SM-list-results.dat')
do i=800,1500
Mh=dble(i)/10
SMGammaTotal=SMGamma_h(Mh)
! SMGamma_h(Mh), SMBR_Hgg(Mh), SMBR_Hgg(Mh) are set to -1 if called
! with Mh out of range [0.8 GeV, 500 GeV]. The calculation is then bypassed.
if(.not. (SMGammaTotal .lt. 0)) then
! SM
GammaTotal=SMGammaTotal
c CP even
g2hjss_s=1d0
g2hjss_p=0d0
g2hjcc_s=1d0
g2hjcc_p=0d0
g2hjbb_s=1d0
g2hjbb_p=0d0
g2hjtt_s=1d0
g2hjtt_p=0d0
g2hjmumu_s=1d0
g2hjmumu_p=0d0
g2hjtautau_s=1d0
g2hjtautau_p=0d0
g2hjWW=R*1.d0
g2hjZZ=R*1.d0
g2hjZga=R*1.d0
g2hjgaga=R*1d0
g2hjgg=R*1.d0
g2hjggZ=R*1.d0
g2hjhiZ=0d0
BR_hjhihi=0d0
BR_hjinvisible=0d0
call HiggsBounds_neutral_input_effC(Mh,GammaTotal,
& g2hjss_s,g2hjss_p,g2hjcc_s,g2hjcc_p,g2hjbb_s,g2hjbb_p,
& g2hjtt_s,g2hjtt_p,
& g2hjmumu_s,g2hjmumu_p,g2hjtautau_s,g2hjtautau_p,
& g2hjWW,g2hjZZ,g2hjZga,g2hjgaga,g2hjgg,g2hjggZ,
& g2hjhiZ, BR_hjinvisible,BR_hjhihi)
dMhneut = (/2.D0/)
! dMhch = (/0.D0/)
call set_mass_uncertainties(dmhneut,dmhch)
call run_HiggsBounds(HBresult,chan,obsratio,ncombined)
write(21,*) Mh,HBresult,chan,obsratio
endif
enddo
close(21)
call finish_HiggsBounds
end

File Metadata

Mime Type
text/x-fortran
Expires
Wed, May 14, 10:37 AM (1 d, 8 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5070630
Default Alt Text
SM_effC.F (3 KB)

Event Timeline