SANCinterfaceF.F

00001       subroutine 
00002      $  flagset(iqedx,iewx,ibornx,gfschemex,ifggx,ncx,fcx,tlmu2x)
00003       implicit none!
00004 #include "s2n_declare.h"
00005       integer iqedx,iewx,ibornx,gfschemex,ifggx
00006       real*8 ncx,fcx,tlmu2x
00007 *-- Flag settings 
00008       iqed     = iqedx
00009       iew      = iewx
00010       iborn    = ibornx
00011       gfscheme = gfschemex
00012       ifgg     = ifggx
00013 * Input constants (masses, couplings)
00014 #include "leshw_input.h"
00015       call s2n_init()
00016       call ffini()
00017       nc=ncx
00018       fc=fcx
00019       tlmu2=tlmu2x
00020       end
00021 
00022       subroutine paraget(mtax,conhcx,pix)
00023       implicit none!
00024       real*8 mtax,conhcx,pix
00025 #include "s2n_declare.h"
00026       mtax=mta
00027       conhcx=conhc
00028       pix= pi
00029       end
00030 
00031       subroutine PrintConsts (pmg)
00032       implicit none!
00033       integer pmg,NOUT
00034 #include "s2n_declare.h"
00035 
00036       if ((pmg.eq.0).or.(pmg.eq.1)) then
00037          print *
00038          print *, "flag iqed     = ",iqed
00039          print *, "flag iew      = ",iew
00040          print *, "flag qcd      = ",iqcd
00041          print *, "flag gfscheme = ",gfscheme
00042          print *, "thmu2(GeV^2)  = ",thmu2
00043          print *, "tlmu2(GeV^2)  = ",tlmu2
00044          print *, "mup(GeV)      = ",mup
00045          print *, "mta(GeV)      = ",mta
00046        endif
00047       if ((pmg.eq.0).or.(pmg.eq.2)) then
00048          print *
00049          print *, "alphai    = ",alphai
00050          print *, "alpha     = ",alpha
00051          print *, "GF(GeV-^2)= ",gf
00052          print *, "g         = ",g
00053          print *, "e         = ",e
00054          print *, "alphas    = ",alpha
00055          print *, "gs        = ",gs
00056          print *, "cwm       = ",ctw
00057          print *, "swm2      = ",stw2
00058          print *, "mh(GeV)   = ",mh
00059          print *, "mw(GeV)   = ",mw
00060          print *, "mz(GeV)   = ",mz
00061          print *, "wz(GeV)   = ",wz
00062          print *, "ww(GeV)   = ",ww
00063          print *, "wh(GeV)   = ",wh
00064          print *, "wtp(GeV)  = ",wtp
00065       endif
00066       
00067       NOUT=11
00068       OPEN(NOUT   ,file='var.dump')
00069          write(NOUT,*) "TAUOLA interace: SANC input parameters:"
00070       if ((pmg.eq.0).or.(pmg.eq.1)) then
00071          write(NOUT,*) " "
00072          write(NOUT,*) "flag iqed     = ",iqed
00073          write(NOUT,*) "flag iew      = ",iew
00074          write(NOUT,*) "flag qcd      = ",iqcd
00075          write(NOUT,*) "flag gfscheme = ",gfscheme
00076          write(NOUT,*) "thmu2(GeV^2)  = ",thmu2
00077          write(NOUT,*) "tlmu2(GeV^2)  = ",tlmu2
00078          write(NOUT,*) "mup(GeV)      = ",mup
00079          write(NOUT,*) "mta(GeV)      = ",mta
00080        endif
00081       if ((pmg.eq.0).or.(pmg.eq.2)) then
00082          write(NOUT,*) " "
00083          write(NOUT,*) "alphai    = ",alphai
00084          write(NOUT,*) "alpha     = ",alpha
00085          write(NOUT,*) "GF(GeV-^2)= ",gf
00086          write(NOUT,*) "g         = ",g
00087          write(NOUT,*) "e         = ",e
00088          write(NOUT,*) "alphas    = ",alpha
00089          write(NOUT,*) "gs        = ",gs
00090          write(NOUT,*) "cwm       = ",ctw
00091          write(NOUT,*) "swm2      = ",stw2
00092          write(NOUT,*) "mh(GeV)   = ",mh
00093          write(NOUT,*) "mw(GeV)   = ",mw
00094          write(NOUT,*) "mz(GeV)   = ",mz
00095          write(NOUT,*) "wz(GeV)   = ",wz
00096          write(NOUT,*) "ww(GeV)   = ",ww
00097          write(NOUT,*) "wh(GeV)   = ",wh
00098          write(NOUT,*) "wtp(GeV)  = ",wtp
00099       endif
00100  
00101       CLOSE(NOUT)
00102 
00103       return
00104       end