]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | |
2 | C********************************************************************* | |
3 | ||
4 | SUBROUTINE LUTAUD(ITAU,IORIG,KFORIG,NDECAY) | |
5 | ||
6 | C...Dummy routine, to be replaced by user, to handle the decay of a | |
7 | C...polarized tau lepton. | |
8 | C...Input: | |
9 | C...ITAU is the position where the decaying tau is stored in /LUJETS/. | |
10 | C...IORIG is the position where the mother of the tau is stored; | |
11 | C... is 0 when the mother is not stored. | |
12 | C...KFORIG is the flavour of the mother of the tau; | |
13 | C... is 0 when the mother is not known. | |
14 | C...Note that IORIG=0 does not necessarily imply KFORIG=0; | |
15 | C... e.g. in B hadron semileptonic decays the W propagator | |
16 | C... is not explicitly stored but the W code is still unambiguous. | |
17 | C...Output: | |
18 | C...NDECAY is the number of decay products in the current tau decay. | |
19 | C...These decay products should be added to the /LUJETS/ common block, | |
20 | C...in positions N+1 through N+NDECAY. For each product I you must | |
21 | C...give the flavour codes K(I,2) and the five-momenta P(I,1), P(I,2), | |
22 | C...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 | ||
28 | C...Stop program if this routine is ever called. | |
29 | C...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 | ||
34 | C...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 |