ISAJET/code/reset.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE RESET
   3 C          RESET ALL USER DEFINED VARIABLES
   4 #if defined(CERNLIB_IMPNONE)
   5       IMPLICIT NONE
   6 #endif
   7 #include "isajet/itapes.inc"
   8 #include "isajet/qcdpar.inc"
   9 #include "isajet/dylim.inc"
  10 #include "isajet/frgpar.inc"
  11 #include "isajet/hcon.inc"
  12 #include "isajet/jetlim.inc"
  13 #include "isajet/jetpar.inc"
  14 #include "isajet/nodcay.inc"
  15 #include "isajet/primar.inc"
  16 #include "isajet/qlmass.inc"
  17 #include "isajet/q1q2.inc"
  18 #include "isajet/seed.inc"
  19 #include "isajet/sspar.inc"
  20 #include "isajet/tcpar.inc"
  21 #include "isajet/types.inc"
  22 #include "isajet/wcon.inc"
  23 #include "isajet/force.inc"
  24 #include "isajet/mbgen.inc"
  25 #include "isajet/isloop.inc"
  26 #include "isajet/limevl.inc"
  27 #include "isajet/xmssm.inc"
  28 #include "isajet/eepar.inc"
  29 #include "isajet/mglims.inc"
  30 C
  31       INTEGER I,I1,I2,I3,J1,INDEX,J,K
  32       REAL UNDEF,AMASS
  33       CHARACTER*8 BLANK
  34       DATA BLANK/'        '/
  35       DATA UNDEF/-1.E9/
  36 C          RESET DYLIM
  37       DO 110 I=1,12
  38       BLIM1(I)=UNDEF
  39       SETLMQ(I)=.TRUE.
  40 110   CONTINUE
  41 C          RESET FRGPAR
  42       PUD=.43
  43       PBARY=.10
  44       SIGQT=.35
  45       PEND=.14
  46       XGEN(1)=.96
  47       XGEN(2)=3.
  48       XGEN(3)=0.
  49       XGEN(4)=.8
  50       XGEN(5)=.5
  51       XGEN(6)=.5
  52       XGEN(7)=.5
  53       XGEN(8)=.5
  54       DO 111 K=1,9
  55 111   XGENSS(K)=.5
  56       PSPIN1(1)=.5
  57       PSPIN1(2)=.5
  58       PSPIN1(3)=.5
  59       PSPIN1(4)=.75
  60       PSPIN1(5)=.75
  61       PSPIN1(6)=.75
  62       PSPIN1(7)=.75
  63       PSPIN1(8)=.75
  64       PMIXX1(1)=.25
  65       PMIXX1(2)=.25
  66       PMIXX1(3)=.5
  67       PMIXX1(4)=0.
  68       PMIXX1(5)=.5
  69       PMIXX1(6)=1.
  70       PMIXX2(1)=.5
  71       PMIXX2(2)=.5
  72       PMIXX2(3)=1.
  73       PMIXX2(4)=0.
  74       PMIXX2(5)=0.
  75       PMIXX2(6)=1.
  76 C          RESET ISLOOP
  77       NEVOLV=1
  78       NFRGMN=1
  79 C          RESET JETLIM
  80       DO 120 I=1,12*MXLIM
  81         BLIMS(I)=UNDEF
  82         SETLMJ(I)=.TRUE.
  83 120   CONTINUE
  84 C          RESET NODCAY
  85       NODCAY=.FALSE.
  86       NOETA=.FALSE.
  87       NOPI0=.FALSE.
  88       NONUNU=.FALSE.
  89       NOEVOL=.FALSE.
  90       NOHADR=.FALSE.
  91       NOGRAV=.FALSE.
  92 C          RESET PRIMAR
  93       IDIN(1)=1120
  94       IDIN(2)=1120
  95       NTRIES=1000
  96       NSIGMA=20
  97 C          RESET QCDPAR
  98       ALAM=.2
  99       ALAM2=ALAM**2
 100       CUTJET=6.
 101       ISTRUC=6
 102 C          RESET QLMASS
 103       AMLEP(6)=175.
 104       AMLEP(7)=-1.
 105       AMLEP(8)=-1.
 106       DO 125 I=1,9
 107       CALL FLAVOR(80+I,I1,I2,I3,J1,INDEX)
 108 125   AMLEP(INDEX)=0.
 109       CALL FLAVOR(29,I1,I2,I3,J1,INDEX)
 110       AMLEP(INDEX)=100.
 111       CALL FLAVOR(30,I1,I2,I3,J1,INDEX)
 112       AMLEP(INDEX)=0.
 113       CALL FLAVOR(39,I1,I2,I3,J1,INDEX)
 114       AMLEP(INDEX)=100.
 115       CALL FLAVOR(40,I1,I2,I3,J1,INDEX)
 116       AMLEP(INDEX)=100.
 117       DO 126 I=1,6
 118       CALL FLAVOR(20+I,I1,I2,I3,J1,INDEX)
 119       AMLEP(INDEX)=100.+AMASS(I)
 120       CALL FLAVOR(30+I,I1,I2,I3,J1,INDEX)
 121       AMLEP(INDEX)=100.+AMASS(I+10)
 122 126   CONTINUE
 123 C          RESET Q1Q2
 124       DO 130 I=1,MXGOQ
 125       DO 130 J=1,MXGOJ
 126 130   GOQ(I,J)=.TRUE.
 127       DO 131 I=1,MXGOJ
 128 131   GOALL(I)=.TRUE.
 129       GODY(1)=.TRUE.
 130       GODY(2)=.FALSE.
 131       GODY(3)=.FALSE.
 132       GODY(4)=.TRUE.
 133       DO 132 I=1,2
 134       ALLWW(I)=.TRUE.
 135       DO 132 J=1,25
 136 132   GOWW(J,I)=.TRUE.
 137       DO 133 I=1,3
 138       DO 133 J=1,25
 139 133   GOWMOD(J,I)=.TRUE.
 140 C          RESET TCPAR
 141       TCMRHO=1000.
 142       TCGRHO=100.
 143 C          RESET TYPES
 144       DO 140 I=1,NTYP
 145 140   LOC(I)=0
 146       DO 141 I=1,MXTYPE
 147       NJTTYP(I)=0
 148       JETYP(1,I)='ALL     '
 149       DO 141 K=2,30
 150 141   JETYP(K,I)=BLANK
 151       JWTYP=4
 152       DO 142 I=1,2
 153       NWWTYP(I)=0
 154       WWTYP(1,I)='ALL     '
 155       DO 142 K=2,4
 156 142   WWTYP(K,I)=BLANK
 157       DO 143 I=1,3
 158       NWMODE(I)=0
 159       WMODES(1,I)='ALL     '
 160       DO 143 K=2,30
 161 143   WMODES(K,I)=BLANK
 162 C          RESET WCON
 163       SIN2W=.232
 164       WMASS(2)=80.2
 165       WMASS(3)=WMASS(2)
 166       WMASS(4)=91.19
 167       CALL FLAVOR(80,I1,I2,I3,J,INDEX)
 168       AMLEP(INDEX)=WMASS(2)
 169       CALL FLAVOR(90,I1,I2,I3,J,INDEX)
 170       AMLEP(INDEX)=WMASS(4)
 171       CUTOFF=.200
 172       CUTPOW=1.0
 173       WFUDGE=1.75
 174 C          RESET MBGEN
 175       MNPOM=1
 176       MXPOM=LIMPOM
 177 C          RESET FORCE
 178       NFORCE=0
 179 C
 180 C          RESET QCD EVOLUTION CUTS
 181       USELIM=.FALSE.
 182       CONCUT=1.0
 183 C
 184 C          RESET SSPAR
 185       AMGVSS=1.E20
 186 C
 187 C          RESET XMSSM
 188       GOMSSM=.FALSE.
 189       GOSUG=.FALSE.
 190       GOAMSB=.FALSE.
 191       XM1SS=1.E20
 192       XM2SS=1.E20
 193       XMAJNR=1.E20
 194       XMGVTO=1.E20
 195 C
 196 C          RESET HCON
 197       IHTYPE=0
 198 C
 199 C          RESET EEPAR
 200       PLEP=0.
 201       PLEM=0.
 202 C
 203 C          RESET MGLIMS
 204       EHMGMN=-1.E9
 205       EHMGMX=-1.E9
 206       YHMGMN=-1.E9
 207       YHMGMX=-1.E9
 208       DO 150 I=1,MXLIM
 209         DO 151 J=1,MXLIM
 210           AMIJMN(I,J)=-1.E9
 211           AMIJMX(I,J)=-1.E9
 212 151     CONTINUE
 213 150   CONTINUE
 214 C
 215       RETURN
 216       END