*
* $Id: zbottom,v 1.1.1.1 1996/01/11 14:15:07 mclareni Exp $
*
* $Log: zbottom,v $
* Revision 1.1.1.1  1996/01/11 14:15:07  mclareni
* Cojets
*
*
* This directory was created from cojets.car patch zbottom
#ifndef CERNLIB_COJCONV
#define CERNLIB_COJCONV
#endif
+ADD,COJCONV,MAIN,38
C-
C-
C- SET UP FOR PRODUCTION OF BOTTOM QUARKS VIA Z DECAY
C-
      ITDKY=1
      ITLIS=6
      NLEAVE=10
      IEVTYP=2
      ECM=630.
      NFLAV=5
      ALAMB=0.2
      IEICHT=1
C- INPUT PARAMETERS FOR Z PRODUCTION AND DECAY
      IDECBO=11
C-
C- FORCE ELECTRON SEMI-LEPTONIC DECAY FOR BOTTOM PARTICLES
      CALL FORCSL('BT','E  ',2)
C-
C- STOP DECAYS OF PI0 AND K0S
      CALL STOPD(110)
      CALL STOPD(20)
C-
C-
+ADD,COJCONV,BOOKH,9
C-
C-
C- SET UP FOR PRODUCTION OF BOTTOM QUARKS VIA Z DECAY
C-
      COMMON//B(10000)
      CALL HLIMIT(10000)
      CALL HOUTPU(ITLIS)
      CALL HERMES(ITLIS)
C
C-- BOOK HISTOGRAMS
      CALL HBOOK1(11,'ETL Z-BOT.$',100,0.,100.)
C-
C-
+ADD,COJCONV,FILLH,15
C-
C-
C- SET UP FOR PRODUCTION OF BOTTOM QUARKS VIA Z DECAY
C-
      COMMON/CUTS/ETLC
      DATA ICALL/0/
C
      IF(ICALL.GT.0) GO TO 10
      ICALL=1
      ETLC=0.
10    CONTINUE
C
C- REJECT EVENTS NOT CONTAINING HEAVY QUARKS WITH ENOUGH ENERGY
      IFL=0
      DO 20 IQ=1,NQUA
      IF(PARQUA(IQ,4).LT.ETLC) GO TO 20
      IL=ABS(INT(PARQUA(IQ,6)))
      IF(IL.EQ.5) IFL=IL
20    CONTINUE
      IF(IFL.EQ.0) RETURN
C
C- BOOK ELECTRONS FROM BOTTOM QUARKS
C
      IF(NPART.EQ.0) RETURN
      DO 1 IPART=1,NPART
C- SELECT ELECTRON/POSITRON
      ITYPE=INT(PARHAD(IPART,6))
      IF(ITYPE.NE.7.AND.ITYPE.NE.8) GO TO 1
C- SELECT BOTTOM QUARK DECAYS ONLY
      IFQRK=ABS(PARQUA(INT(PARHAD(IPART,7)),6))
      IF(IFQRK.NE.5) GO TO 1
C
      ETL=SQRT(PARHAD(IPART,1)**2+PARHAD(IPART,2)**2)
      IF(ETL.LT.ETLC) GO TO 1
      CALL HFILL(11,ETL,0.,1.)
1     CONTINUE
C-
C-
+ADD,COJCONV,OUTH,10
C-
C-
C- SET UP FOR PRODUCTION OF BOTTOM QUARKS VIA Z DECAY
C-
      COMMON/CUTS/ETLC
C
      IF(NEVNTS.EQ.0) RETURN
      FACT=SIGTOT/FLOAT(NEVNTS)
      CALL HOPERA(11,'+',11,11,FACT,0.)
C
      WRITE(ITLIS,5) ETLC
5     FORMAT(/////1X,50(1H=)/1X,'ELECTRONS FROM B QUARKS'
     1,5X,6HETLC =,E12.3,' GEV'/1X,50(1H=)/////)
C
      CALL HISTDO
C-
C-
