GEANT321/erdecks/eufilp.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/03/06 15:37:35  mclareni
   6 * Add geane321 source directories
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.49  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE EUFILP (N, EIN, PLI, PLF)
  13 *
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *    User routine to fill the input values of the commons :      *
  17 C.    *               /EROPTS/, /EROPTC/ and /ERTRIO/ for CHOPT = 'P'  *
  18 C.    *         N     Number of predictions where to store results     *
  19 C.    *         EIN   Input error matrix (in the 'Plane' system )      *
  20 C.    *         PLI   Defines the start plane                          *
  21 C.    *                      PLI(3,1) - and                            *
  22 C.    *                      PLI(3,2) - 2 unit vectors in the plane    *
  23 C.    *         PLF   Defines the end plane                            *
  24 C.    *                      PLF(3,1,I) - and                          *
  25 C.    *                      PLF(3,2,I) - 2 unit vectors in the plane  *
  26 C.    *                      PLF(3,3,I) - point on the plane           *
  27 C.    *                                   at intermediate point I      *
  28 C.    *                                                                *
  29 C.    *    ==>Called by : USER (before calling ERTRAK)                 *
  30 C.    *       Author    M.Maire, E.Nagy  *********                     *
  31 C.    *                                                                *
  32 C.    ******************************************************************
  33 C.
  34 #include "geant321/ertrio.inc"
  35       DIMENSION EIN(15), PLF(3,3,N), PLI(3,2),PWORK(3)
  36 *
  37       NEPRED = MIN0 (N, MXPRED)
  38 *
  39       CALL UCOPY (EIN(1), ERRIN(1), 15)
  40       CALL UCOPY (PLI(1,1), ERPLI(1,1), 6)
  41       DO 5 I = 1, NEPRED
  42          CALL UCOPY (PLF(1,1,I), ERPLO(1,1,I), 9)
  43          CALL CROSS (PLF(1,1,I), PLF(1,2,I), PWORK(1))
  44          PNORM = VMOD(PWORK,3)
  45          CALL VSCALE(PWORK,PNORM,ERPLO(1,4,I),3)
  46     5 CONTINUE
  47 *
  48 *                                                           END EUFILP
  49   99  END