New functions and constructors added and some other fixes.
[u/mrichter/AliRoot.git] / PYTHIA / jetset / lutaud.F
CommitLineData
fe4da5cc 1
2C*********************************************************************
3
4 SUBROUTINE LUTAUD(ITAU,IORIG,KFORIG,NDECAY)
5
6C...Dummy routine, to be replaced by user, to handle the decay of a
7C...polarized tau lepton.
8C...Input:
9C...ITAU is the position where the decaying tau is stored in /LUJETS/.
10C...IORIG is the position where the mother of the tau is stored;
11C... is 0 when the mother is not stored.
12C...KFORIG is the flavour of the mother of the tau;
13C... is 0 when the mother is not known.
14C...Note that IORIG=0 does not necessarily imply KFORIG=0;
15C... e.g. in B hadron semileptonic decays the W propagator
16C... is not explicitly stored but the W code is still unambiguous.
17C...Output:
18C...NDECAY is the number of decay products in the current tau decay.
19C...These decay products should be added to the /LUJETS/ common block,
20C...in positions N+1 through N+NDECAY. For each product I you must
21C...give the flavour codes K(I,2) and the five-momenta P(I,1), P(I,2),
22C...P(I,3), P(I,4) and P(I,5). The rest will be stored automatically.
23
24 COMMON/LUJETS/N,K(4000,5),P(4000,5),V(4000,5)
25 COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
26 SAVE /LUJETS/,/LUDAT1/
27
28C...Stop program if this routine is ever called.
29C...You should not copy these lines to your own routine.
30 NDECAY=ITAU+IORIG+KFORIG
31 WRITE(MSTU(11),5000)
32 IF(RLU(0).LT.10.) STOP
33
34C...Format for error printout.
35 5000 FORMAT(1X,'Error: you did not link your LUTAUD routine ',
36 &'correctly.'/1X,'Dummy routine in JETSET file called instead.'/
37 &1X,'Execution stopped!')
38
39
40 RETURN
41 END