C++ Interface to Tauola
eli/demo.f
1 PROGRAM spintau
2C *****************
3 IMPLICIT DOUBLE PRECISION (a-h,o-z)
4 common/pysubs/msel,mselpd,msub(500),kfin(2,-40:40),ckin(200)
5 common/pypars/mstp(200),parp(200),msti(200),pari(200)
6 common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
7 common/pydat3/mdcy(500,3),mdme(8000,2),brat(8000),kfdp(8000,5)
8 COMMON / inout / inut,nout
9 EXTERNAL pydata !, LUDATA
10
11C to prevent tau decays in jetset
12 mdcy(15,1)=0
13c------------------------------------>>>>>>
14C-----output ident for TAUOLA.
15 nout2 = 6
16 nout = 6
17c.....JETSET output will be writen on NOUT2
18 mstu(11) = nout2
19c..... PHOTOS will write on 6
20CC-------------------------------------->>>>
21c OPEN(0 ,file='tauola.out')
22 OPEN(nout ,file='demo.out')
23
24c------------------------------------>>>>>>
25
26c------------------------------------>>>>>>
27C------number of requested events
28 nevt=9
29 keypol=1
30 CALL tauola(-1,keypol)
31
32
33 msel=0
34 msub(30)=1 !fg --> fZ
35 msub(15)=1 !ff --> gZ
36 DO i=174,189
37 mdme(i,1) = 0
38 ENDDO
39 mdme(186,1) = 1 ! only Z --> tau+ tau- decay
40 mstp(128)=1
41 call pyinit('cms','p','p',14000.d0)
42 call pylist(12)
43
44C>>>>>>>>>>>>>>
45 DO ii= 1,nevt
46 IF(mod(ii,50).EQ.1) WRITE(6,*)'event no=',ii
47 mdcy(15,1)=0
48 CALL pyevnt
49! if(ii.lt.5)CALL PYLIST(1)
50 CALL pyhepc(1)
51 n11=ihepdim(dum) ! no of entries in hepevt before tauola
52 CALL tauola(0,keypol)
53 n12=ihepdim(dum) ! no of entries in hepevt after tauola
54
55
56 if (n12.gt.n11) then
57 IF(ii.LT.101) CALL pylist(1)
58 CALL pyhepc(2)
59 IF(ii.LT.101) CALL pylist(1)
60 endif
61 ENDDO
62CC>>>>>>>>>>>>>>
63 CALL tauola(1,keypol)
64
65 CLOSE(nout2)
66 CLOSE(nout)
67 END
68
69