tauola-F/jetset-F/demo.f

00001       PROGRAM SPINTAU
00002 C     *****************
00003       IMPLICIT REAL(A-H,O-Z)
00004       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) 
00005       COMMON/LUDAT3/MDCY(500,3),MDME(2000,2),BRAT(2000),KFDP(2000,5) 
00006       COMMON / INOUT / INUT,NOUT
00007  
00008 C to prevent tau decays in jetset
00009        mdcy(15,1)=0
00010 c------------------------------------>>>>>>
00011 C-----output ident for TAUOLA. 
00012          NOUT2 = 6
00013          NOUT  = 6
00014 c.....JETSET output will be writen on NOUT2
00015        MSTU(11) = NOUT2
00016 c..... PHOTOS will write on 6
00017 CC-------------------------------------->>>>
00018 c      OPEN(0      ,file='tauola.out')
00019       OPEN(NOUT   ,file='demo.out')
00020 
00021 c------------------------------------>>>>>>
00022     
00023 c------------------------------------>>>>>>
00024 C------number of requested events
00025        NEVT=99
00026 C polarization switch
00027        KEYPOL=1
00028 
00029        CALL TAUOLA(-1,KEYPOL) 
00030 
00031 C>>>>>>>>>>>>>>
00032          DO II= 1,NEVT
00033            IF(MOD(II,50).EQ.1) WRITE(6,*)'event no=',II
00034            mdcy(15,1)=0 
00035            CALL LUEEVT(5,200.0)
00036            CALL LUHEPC(1)
00037            n11=ihepdim(dum) ! no of entries in hepevt before tauola
00038            CALL TAUOLA(0,KEYPOL) 
00039            n12=ihepdim(dum) ! no of entries in hepevt after  tauola
00040 
00041            if (n12.gt.n11) then
00042              IF(II.LT.101) CALL LULIST(1)
00043                            CALL LUHEPC(2)
00044              IF(II.LT.101) CALL LULIST(1)
00045            endif
00046          ENDDO
00047 CC>>>>>>>>>>>>>>
00048          CALL TAUOLA(1,KEYPOL) 
00049 
00050       CLOSE(NOUT2)
00051       CLOSE(NOUT)
00052 
00053       END 
00054