#!/bin/csh
###############################################
#   Submit to LSF with bsub < ZFITTER.JOB     #
###############################################
#BSUB -J ZFITTER
# Execute ~/.login and ~/.cshrc files
# # BSUB -L /usr/local/bin/csh
#  define the queue  (short=8nm, med=1nh or 8nh, long=1nd or 1nw)
#BSUB -q xu_8nm
# # BSUB -o zfitter.out
#
# Initialize
#source ~/.login
#source ~/.tcshrc
set echo
date
#
banner 'ZFITTER'
set echo  # echo commands.
echo $HOST
rm fun.F
rm fun
#
# === Fortran
#
cat <<EOF >fun.F
      PROGRAM FUN
C---------------------------------------------------------------------
C! Calculate Standard Model ee, mumu, tautau and qq cross-sections 
C! and asymmetries using ZFITTER.
C!
C! The ZFITTER writeup can be found in CERN-TH-6443/92.
C!
C! Author: I.Tomalin
C! Date: 5-Jun-1997
C---------------------------------------------------------------------
      IMPLICIT REAL*8(A-H,O-Z)
C
C Centre of mass energies.
      PARAMETER (NCMS=9)
      DIMENSION ECMS(NCMS)
      DATA ECMS/91.1862, 130.200, 136.200, 161.314, 172.086, 182.68,
     &          188.62,191.6,195.6/
C
C Standard Model inputs and their errors, taken from CERN-EP/99-15,
C except for Higgs mass which is taken from Dieu - private communication
C and DAL5H taken from p152 of Zfitter V6.10 manual.
      PARAMETER (ZMASS=91.1867,   EZMASS=0.0021)  !Z0 mass
      PARAMETER (TMASS=171.1,     ETMASS=4.9)     !Top quark mass
      PARAMETER (HMASS=127.0,     EHMASS=100.0)   !Higgs mass
      PARAMETER (ALEMINV=128.878, EALEMINV=0.090) !1/alpha QED at MZ
      DATA DAL5H /0.0280398/                      !Delta(alpha_k^5(M_Z)
      PARAMETER (ALFAS=0.119,     EALFAS=0.003)   !Strong coupling at MZ
C
      CHARACTER*9 FERMION(0:10)
      DATA FERMION /'nu nu', 'e+ e-', 'mu+ mu-', 'tau+ tau-',
     &      'u ubar', 'd dbar', 'c cbar', 's sbar', 't tbar', 'b bbar', 
     &      'hadrons' /
C
C Cuts on sqrt(s'/s), sqrt(s'/s)_pair and |cos(theta)!.
C
      PARAMETER (RSCUT=0.85, RSPCUT=RSCUT, CTHCUT=1.0)
C
C---------------------------------------------------------------------
C Initialize.
      ALEM = 1.0/ALEMINV
      CALL ZUINIT
C
      WRITE(6,10) RSCUT,RSPCUT,CTHCUT
   10 FORMAT(/,'f-fbar(gamma) cross-sections from ZFITTER v6.04 ',
     &/,' using sqrt(s''/s) cut of ',F5.3,' in normal events',
     &/,' using sqrt(s''/s) cut of ',F5.3,' in radiative pair events',
     &/,' and a |cos(theta)| cut of ',F5.3,
     &/,'==========================================================',
     &'==================')
C
C Define s' as invariant mass of propogator.
      CALL ZUFLAG('FINR',0)
C Switch off initial-final state photon interference.
      CALL ZUFLAG('INTF',0)
C Treat pair production from ISR photons as in old ZFITTER versions.
C (With the default option ISPP=1, the inclusive cross sections
C  depend very strongly on the argument SIPP of subroutine ZUCUTS).
C     CALL ZUFLAG('ISPP',-1)
C
      CALL ZUINFO(0)
C
C Do weak sector calculations.
      CALL ZUWEAK(ZMASS,TMASS,HMASS,DAL5H,ALFAS)
C
      WRITE(6,20) ZMASS,TMASS,HMASS,ALEMINV,DAL5H,ALFAS
   20 FORMAT(/,' Standard Model inputs ',
     &       /,' Z0 mass      = ',F9.5,' GeV',
     &       /,' t quark mass = ',F9.5,' GeV',
     &       /,' Higgs mass   = ',F9.5,' GeV',
     &       /,' 1/alpha_em   = ',F9.5,
     &       /,' DAL5H        = ',F9.5,
     &       /,' alpha_s      = ',F9.5)
C
C--- Loop over centre of mass energies.
      DO 80 IENE = 6,NCMS
C
C Calculate Standard Model cross-sections and asymmetries.
        WRITE(6,50) ECMS(IENE)
   50   FORMAT(/,' %%%%%%%%%%%%%%%% Results for centre of',
     &  ' mass energy = ',F7.3,' GeV %%%%%%%%%%%%%%%%',/,
     &  ' Fermion Type    Cross-Section (nb)     FB-Asymmetry')
C
        SPCUT  = (RSCUT*ECMS(IENE))**2
        SIPP   = (RSPCUT*ECMS(IENE))**2
        ANGMIN = ACOSD(CTHCUT)
        ANGMAX = 180.0D0 - ANGMIN
C
        DO 70 IFERM = 0,10
          CALL ZUCUTS(IFERM,1,180.0D0,0.0D0,SPCUT,ANGMIN,ANGMAX,SIPP)
C
          CALL ZUTHSM(IFERM,ECMS(IENE),ZMASS,TMASS,HMASS,DAL5H,ALFAS,
     &                SIGMA,ASYMM)
C
          WRITE(6,60) FERMION(IFERM),SIGMA,ASYMM
   60     FORMAT(3X,A9,6X,F10.6,10X,F10.5)
   70   CONTINUE
C
   80 CONTINUE
      END
EOF
#
# Compile
#
setenv BARDIN /afs/cern.ch/user/b/bardindy/public/ZF6_04
f77 -o fun \
    -trapuw -check_bounds -O -static -fpe -w1 -warn nouncalled -K \
    fun.F \
    $BARDIN/zfitr6_04.f \
    $BARDIN/dizet5_20.f \
    $BARDIN/acol6_04.f \
    $BARDIN/m2tcor5_11.f \
    $BARDIN/bcqcdl5_14.f \
    $BARDIN/bkqcdl5_14.f \
    $BARDIN/bhang4_62.f \
    $BARDIN/zf514_aux.f 
#
# f77 gives error return if duplicate subroutines exist.
if ($status > 1) goto EndJob
chmod u+x fun
#
# Execute
#
fun
#
# Conclude
EndJob:
date
pwd
time
#