DPMJET/flukapro/(DCDRBS)

   1 *$ CREATE DCDRBS.ADD
   2 *COPY DCDRBS
   3 *
   4 *=== Dcdrbs ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     DeCay DiRection BiaSing:                                         *
   9 *                                                                      *
  10 *     Created on   30 march 1998   by    Alfredo Ferrari & Paola Sala  *
  11 *                                                   Infn - Milan       *
  12 *                                                                      *
  13 *     Last change on  13-apr-98    by    Alfredo Ferrari               *
  14 *                                                                      *
  15 *          Aldcdr = Lamda for eps=1-cos(theta). The polar direction    *
  16 *                   around the direction corresponding to the one      *
  17 *                   wished by the user is sample from                  *
  18 *                      P(eps) = exp (-eps/Aldcdr)                      *
  19 *      U,V,Wdcdrb = Lab direction the neutrino should go along         *
  20 *          Sdcdrb = Sin(acos(Wdcdrb))                                  *
  21 *          Spdcdr = Vdcdrb / Sdcdrb                                    *
  22 *          Cpdcdr = Udcdrb / Sdcdrb                                    *
  23 *          Aldcdc = current Lamda for eps=1-cos(theta)                 *
  24 *          Sth0dc = Sin(Theta_0), where (Theta_0,Phi_0) is the direc-  *
  25 *                   tion corresponding to U,V,Wdcdrb in the decaying   *
  26 *                   particle CMS                                       *
  27 *          Cth0dc = Cos(Theta_0)                                       *
  28 *          Sph0dc = Sin(Phi_0)                                         *
  29 *          Cph0dc = Cos(Phi_0)                                         *
  30 *          Kpdcdr = number (1,2,3...) of the decay product whose       *
  31 *                   direction is biased                                *
  32 *       Ldcdrb(i) = flag for decay direction biasing for i_th particle *
  33 *                   type (Paprop numbering scheme)                     *
  34 *          Ldcdbc = current flag for decay direction biasing           *
  35 *                                                                      *
  36 *----------------------------------------------------------------------*
  37 *
  38       LOGICAL LDCDRB, LDCDBC
  39 *
  40       COMMON / DCDRBS / ALDCDR, UDCDRB, VDCDRB, WDCDRB, SDCDRB, SPDCDR,
  41      &                  CPDCDR, ALDCDC, STH0DC, CTH0DC, SPH0DC, CPH0DC,
  42      &                  KPDCDR, LDCDRB (-6:NALLWP), LDCDBC
  43