C++ Interface to Tauola
eli/demo.f
1  PROGRAM spintau
2 C *****************
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 
11 C to prevent tau decays in jetset
12  mdcy(15,1)=0
13 c------------------------------------>>>>>>
14 C-----output ident for TAUOLA.
15  nout2 = 6
16  nout = 6
17 c.....JETSET output will be writen on NOUT2
18  mstu(11) = nout2
19 c..... PHOTOS will write on 6
20 CC-------------------------------------->>>>
21 c OPEN(0 ,file='tauola.out')
22  OPEN(nout ,file='demo.out')
23 
24 c------------------------------------>>>>>>
25 
26 c------------------------------------>>>>>>
27 C------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 
44 C>>>>>>>>>>>>>>
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
62 CC>>>>>>>>>>>>>>
63  CALL tauola(1,keypol)
64 
65  CLOSE(nout2)
66  CLOSE(nout)
67  END
68 
69