PYTHIA/jetset/lutaud.F

   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