pythia_main.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pythia_main.f

*
      COMMON /ludat2/ kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
      COMMON /ludat3/ mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
      COMMON /pysubs/ msel,msub(200),kfin(2,-40:40),ckin(200)
      COMMON /pypars/ mstp(200),parp(200),msti(200),pari(200)
***********      COMMON /LUJETS/ N,K(4000,5),P(4000,5),V(4000,5)
      EXTERNAL ludata,pydata
      CHARACTER frame*4,beam*10,target*10
      REAL win
      INTEGER lucomp
      REAL pdt(500,5)
*
      frame  = 'CMS'
      beam   = 'P'
      TARGET = 'P'
      win    = 14000.
      nevnt  = 3
*
      msel = 0
      msub(102) = 1
      msub(123) = 1
      msub(124) = 1
*
      pmas(6,1) = 176.
      pmas(25,1) = 500.
      ckin(1) = 470.
      ckin(2) = 530.
*
*     turn off all Higgs decays except the ZZ
*
      ih = lucomp(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 = lucomp(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(frame,beam,TARGET,win)
*
      DO ievt = 1, nevnt
*
       CALL pyevnt
       CALL luedit(11)
       CALL luedit(12)
       CALL luedit(15)
*
       CALL luhepc(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=2000)
      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(4i5,4(1x,d15.8))
      ENDDO
*
      RETURN
      END