[u/mrichter/AliRoot.git] / GEANT321 / gtrak / gstord.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:43  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.23  by  S.Giani
*-- Author :
      SUBROUTINE GSTORD
C.
C.    ******************************************************************
C.    *                                                                *
C.    *   Orders the particles in the stack in ascending order by      *
C.    *   the user word.                                               *
C.    *                                                                *
C.    *   Called by : GLTRAC                                           *
C.    *   Authors   : F.Carminati                                      *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gcbank.inc"
#include "geant321/gcjloc.inc"
#include "geant321/gcstak.inc"
#include "geant321/gctrak.inc"
#include "geant321/gcnum.inc"
      PARAMETER (NWPOLA=3,ISTPOS=12)
*
* ----------------------------------------------------------------------
*
      JSTPOS = JSTAK+3+ISTPOS
      VALMAX = Q(JSTPOS)
      INDMAX = 1
      NPSTAK = IQ(JSTAK+1)
      DO 10 J=2,IQ(JSTAK+1)
         JSTPOS = JSTPOS+NWSTAK
         VALTST = Q(JSTPOS)
         IF(VALTST.GT.VALMAX) THEN
            VALMAX = VALTST
            INDMAX = J
         ENDIF
   10 CONTINUE
      IF(INDMAX.NE.NPSTAK) THEN
         JSMOVE = JSTAK+3+(INDMAX-1)*NWSTAK
         JSLAST = JSTAK+3+(NPSTAK-1)*NWSTAK
         DO 20 J=1,NWSTAK
            TEMP = Q(JSMOVE+J)
            Q(JSMOVE+J) = Q(JSLAST+J)
            Q(JSLAST+J) = TEMP
   20    CONTINUE
         IF(JPOLAR.NE.0) THEN
            JPMOVE = JPOLAR+(INDMAX-1)*NWPOLA
            JPLAST = JPOLAR+(NPSTAK-1)*NWPOLA
            DO 30 J=1,NWPOLA
               TEMP = Q(JPMOVE+J)
               Q(JPMOVE+J) = Q(JPLAST+J)
               Q(JPLAST+J) = TEMP
   30       CONTINUE
        ENDIF
      ENDIF
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:21:43 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.23 by S.Giani
	11	*-- Author :
	12	SUBROUTINE GSTORD
	13	C.
	14	C. ******************************************************************
	15	C. * *
	16	C. * Orders the particles in the stack in ascending order by *
	17	C. * the user word. *
	18	C. * *
	19	C. * Called by : GLTRAC *
	20	C. * Authors : F.Carminati *
	21	C. * *
	22	C. ******************************************************************
	23	C.
	24	#include "geant321/gcbank.inc"
	25	#include "geant321/gcjloc.inc"
	26	#include "geant321/gcstak.inc"
	27	#include "geant321/gctrak.inc"
	28	#include "geant321/gcnum.inc"
	29	PARAMETER (NWPOLA=3,ISTPOS=12)
	30	*
	31	* ----------------------------------------------------------------------
	32	*
	33	JSTPOS = JSTAK+3+ISTPOS
	34	VALMAX = Q(JSTPOS)
	35	INDMAX = 1
	36	NPSTAK = IQ(JSTAK+1)
	37	DO 10 J=2,IQ(JSTAK+1)
	38	JSTPOS = JSTPOS+NWSTAK
	39	VALTST = Q(JSTPOS)
	40	IF(VALTST.GT.VALMAX) THEN
	41	VALMAX = VALTST
	42	INDMAX = J
	43	ENDIF
	44	10 CONTINUE
	45	IF(INDMAX.NE.NPSTAK) THEN
	46	JSMOVE = JSTAK+3+(INDMAX-1)*NWSTAK
	47	JSLAST = JSTAK+3+(NPSTAK-1)*NWSTAK
	48	DO 20 J=1,NWSTAK
	49	TEMP = Q(JSMOVE+J)
	50	Q(JSMOVE+J) = Q(JSLAST+J)
	51	Q(JSLAST+J) = TEMP
	52	20 CONTINUE
	53	IF(JPOLAR.NE.0) THEN
	54	JPMOVE = JPOLAR+(INDMAX-1)*NWPOLA
	55	JPLAST = JPOLAR+(NPSTAK-1)*NWPOLA
	56	DO 30 J=1,NWPOLA
	57	TEMP = Q(JPMOVE+J)
	58	Q(JPMOVE+J) = Q(JPLAST+J)
	59	Q(JPLAST+J) = TEMP
	60	30 CONTINUE
	61	ENDIF
	62	ENDIF
	63	END