*
* $Id: hbottom,v 1.1.1.1 1996/01/11 14:15:06 mclareni Exp $
*
* $Log: hbottom,v $
* Revision 1.1.1.1  1996/01/11 14:15:06  mclareni
* Cojets
*
*
* This directory was created from cojets.car patch hbottom
#ifndef CERNLIB_COJCONV
#define CERNLIB_COJCONV
#endif
+ADD,COJCONV,MAIN,38
C-
C-
C- SET UP FOR HADRONIC PRODUCTION OF BOTTOM QUARKS
C-
      ITDKY=1
      ITLIS=6
      NLEAVE=10
      IEVTYP=1
      ECM=630.
      NFLAV=5
      ALAMB=0.2
      IEICHT=1
      PTMGE=10.
C-
C- FUSION: FORCE BY CALLING FUSION
C- FLAVOR EXCITATION: FORCE BY CALLING FLAVEX
C- GLUON SPLITTING INTO HEAVY QUARK PAIRS: FORCE BY CALLING GSPLIT
C- HEAVY QUARK PRODUCTION: FORCE BY CALLING HEAVYQ
CCC   CALL HEAVYQ('BT')
CCC   CALL FUSION('BT')
      CALL FLAVEX('BT')
CCC   CALL GSPLIT('BT')
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 HADRONIC PRODUCTION OF BOTTOM QUARKS
      COMMON//B(10000)
      CALL HLIMIT(10000)
      CALL HOUTPU(ITLIS)
      CALL HERMES(ITLIS)
C-
C- BOOK HISTOGRAMS
      CALL HBOOK1(1,'PTHARD   $',100,0.,100.)
C-
      CALL HBOOK1(11,'ETL G-SPL.$',100,0.,100.)
      CALL HBOOK1(21,'ETL FUSION$',100,0.,100.)
      CALL HBOOK1(31,'ETL FL.EX.$',100,0.,100.)
C-
C-
+ADD,COJCONV,FILLH,15
C-
C-
C- SET UP FOR HADRONIC PRODUCTION OF BOTTOM QUARKS
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- ICH=1 - GLUON SPLITTING, ICH=2 - FUSION, ICH=3 FLAVOR EXCITATION
      ICH=1
      IF(LINV(1).GT.5) ICH=3
      IF((LINV(1).EQ.3.OR.LINV(1).EQ.5).AND.LFINV(1).GE.4) ICH=2
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(1,PTHARD,0.,1.)
      CALL HFILL(10*ICH+1,ETL,0.,1.)
1     CONTINUE
C-
C-
+ADD,COJCONV,OUTH,10
C-
C-
C- SET UP FOR HADRONIC PRODUCTION OF BOTTOM QUARKS
C-
      COMMON/CUTS/ETLC
C
      IF(NEVNTS.EQ.0) RETURN
      FACT=SIGTOT/FLOAT(NEVNTS)
      DO 20 ICH=1,3
      IST=10*ICH+1
      CALL HOPERA(IST,'+',IST,IST,FACT,0.)
20    CONTINUE
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-
