This routine decides if a pair-production process occurs while the photon with momentum PP goes through the detector, and if so generates the process and updates the contents of LUJETS accordingly.
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 photon (in GeV/c) PMINPA : Minimum P of photon to consider pair-production. KKK : line-number of the incomming photon in LUJETS
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 2 : A interaction took place, the gamma marked as decayed, new electron and positron inserted into LUJETS.
From /ZXEVT/ ANC, ancestor of initial particle;
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.
(To get Track INFo at a layer)
(To get number of X0:s of a plane)
(To get number of X0:s of a cylinder)
(To get Plane Atom number)
(To get Cylinder Atom number)
(To get Lambda for PAiR-production)
(To Make a e+e- PAiR)
(To DELete LaYers)