ISAJET/code/mbset.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE MBSET
   3 C
   4 C          SET PARAMETERS FOR GENERATING MINBIAS EVENTS OR BEAM JETS,
   5 C          ALLOWING DIFFERENT PARAMETERS FOR TWO CASES.
   6 C
   7 #include "isajet/itapes.inc"
   8 #include "isajet/mbpar.inc"
   9 #include "isajet/mbgen.inc"
  10 #include "isajet/primar.inc"
  11 #include "isajet/totals.inc"
  12 #include "isajet/keys.inc"
  13 C
  14 C
  15 C          DN/DY INCREASES WITH LOG(S). INCLUDED IN SPLITTING FUNCTION
  16 C          BECAUSE AVERAGE MULTIPLICITY COMES FROM SINGLE CHAIN GRAPH.
  17       XGEN0(1)=.9
  18       XGEN0(2)=1.+0.35*ALOG(ECM/60.)
  19 C
  20 C          POMWT ARE (RELATIVE) PROBABILITIES FOR N CUT POMERONS.
  21 C          PDIFFR IS DIFFRACTIVE PROBABILITY.
  22 C          SIGQT0 IS MEAN PT.
  23       IF(KEYS(4)) THEN
  24         PDIFFR=.15
  25         SIGQT0=.35
  26         PSUM=0.
  27         DO 100 I=1,LIMPOM
  28         POMWT(I)=(1.+4.*I**2)*EXP(-1.8*I)
  29         PSUM=PSUM+POMWT(I)
  30 100     CONTINUE
  31       ELSE
  32         PDIFFR=0.
  33         SIGQT0=.45
  34         PSUM=0.
  35         DO 110 I=1,LIMPOM
  36         POMWT(I)=(1.+4.*I**2)*EXP(-1.8*I)
  37         PSUM=PSUM+POMWT(I)
  38 110     CONTINUE
  39         POMWT(1)=.1*POMWT(1)
  40         POMWT(2)=.2*POMWT(2)
  41         POMWT(3)=.5*POMWT(3)
  42       ENDIF
  43 C
  44 C          RENORMALIZE POMWT.
  45       PSUM=1./PSUM
  46       DO 200 I=1,LIMPOM
  47       POMWT(I)=PSUM*POMWT(I)
  48 200   CONTINUE
  49       PSUM=0.
  50       DO 210 I=MNPOM,MXPOM
  51       PSUM=PSUM+POMWT(I)
  52 210   CONTINUE
  53 C
  54 C          POMGEN IS USED TO SELECT NUMBER OF POMERONS.
  55       PGEN=0.
  56       PSUM=1./PSUM
  57       DO 300 I=1,LIMPOM
  58       POMGEN(I)=0.
  59 300   CONTINUE
  60       DO 310 I=MNPOM,MXPOM
  61       PGEN=PGEN+PSUM*POMWT(I)
  62       POMGEN(I)=PGEN
  63 310   CONTINUE
  64       POMGEN(MXPOM)=1.
  65 C
  66 C          SET /TOTALS/ FOR MINBIAS EVENTS USING LOG**2(S) FIT TO
  67 C          TOTAL CROSS SECTION.
  68       IF(KEYS(4)) THEN
  69         SIGTOT=25.65*(1.+.0102*ALOG(SCM/1.76)**2)
  70         SIGTOT=PSUM*SIGTOT
  71         NKINPT=NEVENT
  72         SUMWT=SIGTOT*NKINPT
  73       ENDIF
  74 C
  75       RETURN
  76       END