pythia_main.f

Go to the documentation of this file.

*
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
*
      common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
      common/pydat3/mdcy(500,3),mdme(8000,2),brat(8000),kfdp(8000,5)
      common/pysubs/msel,mselpd,msub(500),kfin(2,-40:40),ckin(200)
      common/pypars/mstp(200),parp(200),msti(200),pari(200)
*
      EXTERNAL pydata
      CHARACTER frame*4,beam*10,target*10
      REAL win
      REAL pdt(500,5)
*
      nevnt  = 100
*
      msel = 0
      msub(102) = 1
      msub(123) = 1
      msub(124) = 1
*
      pmas(6,1) = 172.9
      pmas(25,1) = 125.
*
*     turn off all Higgs decays except the ZZ
*
      ih = pycomp(25)
      DO idc=mdcy(ih,2),mdcy(ih,2)+mdcy(ih,3)-1
        IF (kfdp(idc,1).NE.23.AND.mdme(idc,1).EQ.1) mdme(idc,1)=0
      ENDDO
*
*     turn off all Z decays except mumu or ee
*
      iz = pycomp(23)
      DO idc=mdcy(iz,2),mdcy(iz,2)+mdcy(iz,3)-1
        IF (mdme(idc,1).EQ.1) THEN
          IF ((iabs(kfdp(idc,1)).NE.13)
     >   .AND.(iabs(kfdp(idc,1)).NE.11)) THEN
            mdme(idc,1)=0
          ENDIF
        ENDIF
      ENDDO
***********************************************************
      CALL pyinit('CMS','p','p',14000d0)
*
      DO ievt = 1, nevnt
*
       CALL pyevnt
       CALL pyedit(11)
       CALL pyedit(12)
       CALL pyedit(15)
*
       CALL pyhepc(1)
*
       CALL hep2g4
*
      ENDDO
***********************************************************
      stop
      END
***********************************************************
      SUBROUTINE hep2g4
*
* Output /HEPEVT/ event structure to G4HEPEvtInterface
*
* M.Asai (asai@kekvax.kek.jp)  --  24/09/96
*
***********************************************************
      parameter(nmxhep=4000)
      common/hepevt/nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
     >jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),vhep(4,nmxhep)
      DOUBLE PRECISION phep,vhep
*
      WRITE(6,*) nhep
      DO ihep=1,nhep
       WRITE(6,10)
     >  isthep(ihep),idhep(ihep),jdahep(1,ihep),jdahep(2,ihep),
     >  phep(1,ihep),phep(2,ihep),phep(3,ihep),phep(5,ihep)
10    FORMAT(i4,i10,i5,i5,4(1x,d15.8))
      ENDDO
*
      RETURN
      END