eli/demo.f

00001       PROGRAM SPINTAU
00002 C     *****************
00003       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
00004       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
00005       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
00006       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
00007       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
00008       COMMON / INOUT / INUT,NOUT
00009       EXTERNAL PYDATA !, LUDATA
00010  
00011 C to prevent tau decays in jetset
00012        mdcy(15,1)=0
00013 c------------------------------------>>>>>>
00014 C-----output ident for TAUOLA.
00015          NOUT2 = 6
00016          NOUT  = 6
00017 c.....JETSET output will be writen on NOUT2
00018        MSTU(11) = NOUT2
00019 c..... PHOTOS will write on 6
00020 CC-------------------------------------->>>>
00021 c      OPEN(0      ,file='tauola.out')
00022       OPEN(NOUT   ,file='demo.out')
00023 
00024 c------------------------------------>>>>>>
00025     
00026 c------------------------------------>>>>>>
00027 C------number of requested events
00028        NEVT=9
00029          KEYPOL=1
00030          CALL TAUOLA(-1,KEYPOL)
00031 
00032 
00033          msel=0
00034          msub(30)=1  !fg --> fZ
00035          msub(15)=1  !ff --> gZ
00036          DO I=174,189
00037            MDME(I,1) = 0
00038          ENDDO
00039            MDME(186,1) = 1 ! only Z --> tau+ tau- decay
00040          mstp(128)=1
00041          call pyinit('cms','p','p',14000.D0)
00042          call PYlist(12)
00043 
00044 C>>>>>>>>>>>>>>
00045          DO II= 1,NEVT
00046            IF(MOD(II,50).EQ.1) WRITE(6,*)'event no=',II
00047            mdcy(15,1)=0
00048            CALL PYEVNT
00049 !           if(ii.lt.5)CALL PYLIST(1)
00050            CALL PYHEPC(1)
00051            n11=ihepdim(dum) ! no of entries in hepevt before tauola
00052            CALL TAUOLA(0,KEYPOL)
00053            n12=ihepdim(dum) ! no of entries in hepevt after  tauola
00054 
00055 
00056            if (n12.gt.n11) then
00057              IF(II.LT.101) CALL PYLIST(1)
00058                            CALL PYHEPC(2)
00059              IF(II.LT.101) CALL PYLIST(1)
00060            endif
00061          ENDDO
00062 CC>>>>>>>>>>>>>>
00063          CALL TAUOLA(1,KEYPOL)
00064 
00065       CLOSE(NOUT2)
00066       CLOSE(NOUT)
00067       END
00068 
00069