4.3.9 SUBROUTINE ZDBREM(LAYLIS, PP, VERT, PMINBR, PTLOSLIM, IPRF, IPZ, PHIIP, KKK, IGEO, DST, NNEW) generate BREMsstrahlung

This routine decides if a brems-strahlung process occurs while the electron with momentum PP goes through the detector, and if so generates the process and updates the contents of LUJETS accordingly.

ARGUMENTS:

INPUT :

     LAYLIS : List of detector layers hit by particle :
              LAYLIS(1,i) = number of the i:th hit layer.
                            LAYLIS(2,i) tells if it is a  plane
                            or a cylinder, if it measures or not,
                            and if it is the reference layer.
              LAYLIS(2,i) = 0  : mathematical layer : no material,
                                 no measurement. Ignored in all
                                 calculations
                            1  : seen by a cylinder
                            2  : seen by a plane
                            -1 : particle went through a cylinder, but
                                 was not seen by the detector.
                            -2 : particle went through a plane, but
                                 was not seen by the detector.
                            10 : the layer is a reference cylinder
                            20 : the layer is a reference plane
              LAYLIS(3,i) = No. of radiation lengths in layer.
              LAYLIS(1,i) = -1 -> end of list.
    PP : 3-momentum of the incomming electron (in GeV/c)
    VERT : Position where PP is given ( in cm).
    PMINBR : Minimum P of electron to consider brems.
    PTLOSLIM : Minimal pt(before)/pt(after) of the electron
               for the brems to be un-noticed by the track-
               reconstruction. If more pt than this is lost,
               a new electron will be created
    IPRP,IPZ : True values of Impact parameters in RPhi and z
    PHIIP : True phi-angle to perigee.
    KKK : line-number of the incomming electron in LUJETS
    IGEO : Index of the detector description to be used.

OUTPUT :

    LAYLIS : List of detector layers hit by particle, updated if there
             was an interaction
    DST : Distance travled from VERT to the point of interaction.
          0.0 if no interaction.
    NNEW : number of new particles generated :
             0 : no interaction
             1 : A interaction took place, but the loss of energy
                 of the electron was so small the the original
                 electron was kept (see PTLOSLIM), and only the
                 gamma inserted into LUJETS
             2 : A interaction took place, old electron marked as
                 decayed, new electron and the gamma inserted into
                 LUJETS.

IMPLICIT INPUTS:

    From /ZXEVT/ ANC, ancestor of initial particle;

IMPLICIT OUTPUTS:

    To /LUJETS/ P, momentum of created particles; V, production-point
       of created particles; K, id. and decay info of initial and
       created particles; NLUND, total number of particles.

CALLED BY :

ZDETSI, DETector SImulation (see 4.3.1).

CALLED ROUTINE(S):

ZTTINF (see 4.4.48)

(To get Track INFo at a layer)

ZTPX0 (see 4.4.45)

(To get number of X0:s of a plane)

ZTCX0 (see 4.4.44)

(To get number of X0:s of a cylinder)

ZTPANB (see 4.4.47)

(To get Plane Atom number)

ZTCANB (see 4.4.46)

(To get Cylinder Atom number)

ZDLBRM (see 4.3.11)

(To get Lambda for BRemsStrahlung)

ZTTISC (see 4.4.5)

(To get TRacK InterSections with Cylinder.)

ZTTISP (see 4.4.7)

(To get TRacK InterSections with Plane.)

ZDMBRM (see 4.3.13)

(To Make a BReMs gamma)

ZDELLY (see 4.3.16)

(To DELete LaYers)

VMOD

(CERNLIB routine)

SIDE EFFECTS:

none

AUTHORS:

M. Berggren, UFRJ

CREATION DATE:

05-may-1996