C++ Interface to Tauola
tauola-fortran
tauola-BBB
jetset-F
tauola-BBB/jetset-F/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
Generated by
1.11.0