C++ Interface to Tauola
tauola/demo-jetset/demo.f
1  PROGRAM spintau
2 C *****************
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 
8 C to prevent tau decays in jetset
9  mdcy(15,1)=0
10 c------------------------------------>>>>>>
11 C-----output ident for TAUOLA.
12  nout2 = 6
13  nout = 6
14 c.....JETSET output will be writen on NOUT2
15  mstu(11) = nout2
16 c..... PHOTOS will write on 6
17 CC-------------------------------------->>>>
18 c OPEN(0 ,file='tauola.out')
19  OPEN(nout ,file='demo.out')
20 
21 c------------------------------------>>>>>>
22 
23 c------------------------------------>>>>>>
24 C------number of requested events
25  nevt=99
26 C polarization switch
27  keypol=1
28 
29  CALL tauola(-1,keypol)
30 
31 C>>>>>>>>>>>>>>
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
47 CC>>>>>>>>>>>>>>
48  CALL tauola(1,keypol)
49 
50  CLOSE(nout2)
51  CLOSE(nout)
52 
53  END
54