C++ Interface to Tauola
tauola-BBB/jetset-F/demo.f
1 PROGRAM spintau
2C *****************
3 IMPLICIT REAL(A-H,O-Z)
4 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
5 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
6 COMMON / inout / inut,nout
7
8C to prevent tau decays in jetset
9 mdcy(15,1)=0
10c------------------------------------>>>>>>
11C-----output ident for TAUOLA.
12 nout2 = 6
13 nout = 6
14c.....JETSET output will be writen on NOUT2
15 mstu(11) = nout2
16c..... PHOTOS will write on 6
17CC-------------------------------------->>>>
18c OPEN(0 ,file='tauola.out')
19 OPEN(nout ,file='demo.out')
20
21c------------------------------------>>>>>>
22
23c------------------------------------>>>>>>
24C------number of requested events
25 nevt=99
26C polarization switch
27 keypol=1
28
29 CALL tauola(-1,keypol)
30
31C>>>>>>>>>>>>>>
32 DO ii= 1,nevt
33 IF(mod(ii,50).EQ.1) WRITE(6,*)'event no=',ii
34 mdcy(15,1)=0
35 CALL lueevt(5,200.0)
36 CALL luhepc(1)
37 n11=ihepdim(dum) ! no of entries in hepevt before tauola
38 CALL tauola(0,keypol)
39 n12=ihepdim(dum) ! no of entries in hepevt after tauola
40
41 if (n12.gt.n11) then
42 IF(ii.LT.101) CALL lulist(1)
43 CALL luhepc(2)
44 IF(ii.LT.101) CALL lulist(1)
45 endif
46 ENDDO
47CC>>>>>>>>>>>>>>
48 CALL tauola(1,keypol)
49
50 CLOSE(nout2)
51 CLOSE(nout)
52
53 END
54