+PATCH,$CORR.
+DECK,CORR1.                     Updates version:  1.34/00  to   1.35/00
+REP,*TITLE*,TITLE,0-1.
&TITLE.
GTRAK     1.35/00  24/08/95  16.54.48  GLOBAL TRACKING SOFTWARE
+REP,$VERSION,V1_34,1-4.
*CMZU:  1.34/03 10/08/95  22.03.48  by  Christian Voelcker
*CMZ :  1.34/02 21/07/95  16.25.59  by  Christian Voelcker
*CMZ :  1.34/01 18/04/95  17.15.39  by  Christian Voelcker
*CMZ :  1.34/00 07/04/95  19.42.30  by  Christian Voelcker
*-- Author :    Christian Voelcker   07/04/95
*
*::>          VERSION  1.34/03   10/08/95 22.03.48
*         *>  Mark's version of install script allows for
*             using #include sequences
*
*::>          VERSION  1.34/02   21/07/95 16.25.59
*          *> install script available for gtrack:
*             UNIX: cmz -install gtrack <files_directory> [ <option> ... ]
*             VMS:  CMZ/INSTALL=GTRACK,<files_directory>,[<option>,...]
*          *> make code F77 consistent (to use gnu g77 compilers)
*
*::>          VERSION  1.34/01   18/04/95 17.15.39
*              typing error in TAXSTO, abs(logical), caused compiler problems
*              on some platforms
*
*::>          VERSION  1.34/00   07/04/95 19.42.30
+ADD,$VERSION,V1_35,*.
&DECK,V1_35.
*CMZU:  1.35/00 24/08/95  16.54.48  by  Christian Voelcker
*-- Author :    Christian Voelcker   24/08/95
*
*::>          VERSION  1.35/00   24/08/95 16.54.48
*          *> select $OS for SGI's and RISC's
*
+REP,GTPROGRM,GTBANK,1-22.
*CMZU:  1.34/02 06/07/95  18.49.34  by  Curtis A. Meyer
*CMZU:  1.34/00 06/04/95  20.07.07  by  Christian Voelcker
*CMZU:  1.33/04 29/03/95  15.31.37  by  Curtis A. Meyer
*CMZU:  1.33/02 13/02/95  10.39.48  by  Christian Voelcker
*CMZU:  1.33/00 18/01/95  20.34.30  by  Christian Voelcker
*CMZU:  1.30/02 02/07/94  20.27.41  by  Curtis A. Meyer
*CMZ :  1.30/01 01/07/94  17.18.41  by  Curtis A. Meyer
*CMZ :  1.21/11 19/05/94  11.42.54  by  Michael Doser
*CMZU:  1.21/07 12/01/92  14.44.47  by  Curtis A. Meyer
*CMZ :  1.21/05 04/12/91  17.17.10  by  Gunter Folger
*CMZ :  1.21/02 11/10/91  11.49.35  by  Mark Burchell
*CMZ :  1.21/01 20/08/91  15.21.11  by  Gunter Folger
*CMZ :  1.21/00 02/08/91  14.35.28  by  Mark Burchell
*CMZ :  1.20/09 08/05/91  15.07.22  by  Gunter Folger
*CMZU:  1.20/08 05/04/91  10.59.21  by  Mark Burchell
*CMZU:  1.20/05 15/02/91  17.43.05  by  Gunter Folger
*CMZ :  1.20/04 14/02/91  15.47.33  by  MARK J BURCHELL
*CMZU:  1.20/03 05/12/90  06.55.33  by  MARK BURCHELL
*-- Author : Mark Burchell
      SUBROUTINE GTBANK
*
*
*****************************************************************
+REP,GTPROGRM,GTBANK,122.
*---->   STRUCTURE OF BANK
+REP,GTPROGRM,GTBANK,131-142.
*-->    WORD  1 = PARTICLE NUMBER
*-->    WORD  2 = 0 TO FLAG BC AND 1 TO FLAG JDC DATA
*-->    WORD  3 = GEANT PARTICLE ID IF ANY FAVOURED GUESS EXISTS
*-->    WORD  4 = CHARGE
*-->    WORD  5 = # OF POINTS FOR DE/DX
*-->    WORD  6 = ID OF TRACK ENTRY IN TCTR  BANK
*-->    WORD  7 = FOR CHARGED TRACKS THIS IS
*                 ID OF MATCHED PED ENTRY IN TBTK BANK (IF MATCH EXISTS)
*-->    WORD  8 = ID OF PARTICLES ORIGINATING VERTEX
*                 For PEDs this is always vertex #1.
*-->    WORD  9 = ID OF PARTICLES DECAY VERTEX
*-->    WORD 10 = CHARGED TRACK RECONSTRUCTION QUALITY WORD :
+REP,GTPROGRM,GTBANK,148-195.
*-->    WORD 11 = PED RECONSTRUCION QUALITY WORD
*                 THIS IS (#PEDS IN THE CLUSTER)
*                 + 100* CENTRAL CRYSTAL COMPOUND INDEX
*-->    WORD 12 = more data about PEDS
*                 THIS IS (id of cluster) + 100*(number of crystals
*                 in the cluster)  and if in same cluster as a matched
*                 PED this is multiplied by -1
*   REMAINING WORDS ARE ALL REAL
*--> (CV not useful)   WORD 13 = Probability from LOCATER of charged track
*-->    WORD 13 = Chi2/ number of degrees of freedom. (good track < 1.2)
*-->    WORD 14 = MOMENTUM (P)
*-->    WORD 15 = ENERGY (FOR CHARGED TRACKS THIS ASSUMES THEY ARE PIONS)
*-->    WORD 16 = Px
*-->    WORD 17 = Py
*-->    WORD 18 = Pz
*-->    WORD 19 = ERROR ON MOMENTUM
*-->    WORD 20 = ERROR ON ENERGY
*-->    WORD 21 = ERROR ON PX
*-->    WORD 22 = ERROR ON Py
*-->    WORD 23 = ERROR ON Pz
*-->    WORD 24 = dE/dx
*-->    WORD 25 = ERROR ON dE/dx
*-->    WORD 26 = CHI**2(electron)
*-->    WORD 27 = CHI**2(muon)
*-->    WORD 28 = CHI**2(pion)
*-->    WORD 29 = CHI**2(kaon)
*-->    WORD 30 = CHI**2(proton)
*-->    WORD 31 = COSINE(ANGLE) between matched charged track and PED
*-->    WORD 32 = COSINE(ANGLE) between cloest unmatched track & PED
*-->    WORD 33 = E1/E9 for PEDs
*-->    WORD 34 = Second moment for PEDS
*-->    WORD 35 = Showermass (of parent cluster)
*-->    WORD 36 = Energy of parent cluster
*-->    WORD 37 = u directional cosine of parent cluster
*-->    WORD 37 = v directional cosine of parent cluster
*-->    WORD 37 = w directional cosine of parent cluster
*   THE LAST WORDS ARE THOSE NEEDED FOR KINEMATIC FITTING
*-->    WORD 40 = PSI ANGLE FOR CHARGED TRACKS
*-->    WORD 41 = 1/PT FOR CHARGED TRACKS
*-->    WORD 42 = TAN(LAMBDA) FOR PARTICLES
*-->    WORD 43 = PHI ANGLE FOR PEDS
*-->    WORD 44 = THETA ANGLE FOR PEDS
*-->    WORD 45 = SQUARE ROOT E FOR PEDS
*-->    WORD 46 = COVARIANT ERROR MATRIX IN LOWER TRIANGULAR FORM
*-->            = FOR CHARGED TRACKS IT IS FOR WORDS 40,41,42
*-->    WORD 51 = FOR PEDS           IT IS FOR WORDS 43,44,45
*
*************************************************************
+REP,GTPROGRM,GTBANK,213.
*  declaration of a multitude of previously undefined local variables
+REP,GTPROGRM,GTBANK,228-244.
*--->    NNTKS IS THE TOTAL NUMBER OF LIFTED SUB-BANKS
*--->    IVERS IS THE VERSION NUMBER OF THE BANK STRUCTURE
*
      INTEGER  I,IR,NNTKS,IVERSQ,ICD,IERH, NENT, NPBS, IERR
      SAVE     NTTKS,IVERSQ,ICD,NENT
      REAL     CRG,XVEC(5),PVEC(3),COVR(5,5),SIGP(3),AMP,TMP(4)
      SAVE     AMP
      INTEGER  MHOL,MTVTV(10),MHOP,MTVVP(10),MHOH,MTVVH(10)
      SAVE     MHOL,MHOP,MHOH
*
      INTEGER  IERTRK,IERTBK,IETCTR
      SAVE     IERTRK,IERTBK,IETCTR
*
      DATA     NTTKS/0/, IVERSQ/0/, ICD/1/, NENT/0/
      DATA     AMP/139.568/
      DATA     MHOL/4HTVTV/
      DATA     MHOP/4HTVVP/
      DATA     MHOH/4HTVVH/
+REP,GTPROGRM,GTBANK,256.
***>     WRITE(4,*) 'EVENT ',NENT,' LIFT NEW LHTRK BANK'
+REP,GTPROGRM,GTBANK,266-272.
***>     WRITE(4,*) 'EVENT ',NENT,' LIFT NEW LTCVX BANK'
         IF (LHTJD.LE.0) THEN
***>        WRITE(4,*) 'FIRST IT IS NECESSARY TO LIFT LHTJD'
            CALL MZLIFT(IXHTJD,LHTJD,LEVHD,-2,MHTJD,0)
         ENDIF
         MTCVX(5) = 2
***>     WRITE(4,*) (MTCVX(I),I=1,5)
+REP,GTPROGRM,GTBANK,322.
***>     WRITE(4,*) 'EVENT ',NENT,' LIFT NEW LTVTX BANK'
+REP,GTPROGRM,GTBANK,351-354.
*  1/4/95 CV
* do the same for the TCVH bank created recently. This bank has a
* homogenous structure, therefore we may just copy its info from
* MTCVH, including MTCVH(5)
+REP,GTPROGRM,GTBANK,390.
***>     WRITE(4,*) 'LIFT THE ITVVP BANK'
+REP,GTPROGRM,GTBANK,401.
***>     WRITE(4,*) ' TVVP BANK LIFTED OK'
+REP,GTPROGRM,GTBANK,408.
*  bank length is structure length * number of tracks at this vertex
+REP,GTPROGRM,GTBANK,418-419.
* cross the fingers that it works!! TVVH is the next bank in the linear
* structure (page 8, zebra manual)
+REP,GTPROGRM,GTBANK,428-434.
*  copy the TCVP bank
               ITCVP = LQ(ITCVT-1)
               ILP = (IK-1)*LENVP
               CALL UCOPY(IQ(ITCVP+1+ILP),IQ(ITVVP+1+ILP),2)
               CALL UCOPY(Q(ITCVP+3+ILP),Q(ITVVP+3+ILP),ILN)
400        CONTINUE
*  now the TCVH bank - only reals, copy at once!
+REP,GTPROGRM,GTBANK,441-446.
 500  CONTINUE
*
*
*        Always  create header ttks bank.
*
***>  WRITE(4,*) 'EVENT ',NENT,' LIFT NEW LTTKS BANK'
+REP,GTPROGRM,GTBANK,486.
***>           WRITE(4,*) 'TRK ',J,'LIFTS BANK FOR PARTICLE ',NNTKS
+REP,GTPROGRM,GTBANK,551-552.
*-->  IN WORD 5 ITCTK IS NOT A TYPING ERROR !
*-->  GET TRACK POINTER INTO TCTK BANK
+REP,GTPROGRM,GTBANK,577-585.
CV               IF (Q(ITCTR+17).GT.0.000001) THEN
CV                  ND = 3*IQ(ITCTR+1) - 5
CV                  Q(JTTKS+13) = PROB(Q(ITCTR+17),ND)
CV               ELSE
CV                  Q(JTTKS+13) = 0.0
CV               ENDIF
*
*  track probability was not useful. (and they did not substract PWC hits!)
*  now (18/01/95) fill the chi2/degree of freedom
+REP,GTPROGRM,GTBANK,592.
* a track should never have no hits ....
+REP,GTPROGRM,GTBANK,599.
*  GET PX FROM PT TIMES COS(ANGLE IN XY PLANE), PY SIMILARLY
+REP,GTPROGRM,GTBANK,605-607.
*  guard against floating point exception. Be careful reordering
*  the computation. The code has been written to try and
*  minimize overflow problems, aprticularly on VAXEN.
+REP,GTPROGRM,GTBANK,649-654.
*--> DE/Dx values.
*
               Q(JTTKS+24) = Q(ITCTR+9)
               Q(JTTKS+25) = Q(ITCTR+10)
*
*---> Next dE/dx. First set all probabilties to zero, then call gtdedx
+REP,GTPROGRM,GTBANK,665.
*-->   ANGLES BETWEEN TRACKS-PEDS
+REP,GTPROGRM,GTBANK,754.
 1000       CONTINUE
+REP,GTPROGRM,GTBANK,787.
***>           WRITE(4,*) 'PED ',NPEDS,'LIFTS BANK FOR PARTICLE ',NNTKS
+REP,GTPROGRM,GTBANK,918.
 1900        CONTINUE
+REP,GTPROGRM,GTBANK,928-932.
*******************************
* content of VERSQQ ...
*      VERSQQ = ' 1.31/00'
*      IVERSQ =  13100
*******************************
+REP,GTPROGRM,GTDEDX,1-5.
*CMZU:  1.34/02 06/07/95  18.50.12  by  Curtis A. Meyer
*CMZU:  1.21/11 14/09/93  14.19.53  by  Michael Doser
*CMZ :  1.20/08 08/05/91  15.02.25  by  Gunter Folger
*CMZ :  1.20/06 07/03/91  16.10.31  by  Mark Burchell
*CMZ :  1.20/03 05/12/90  06.55.33  by  MARK BURCHELL
*-- Author :
+REP,GTPROGRM,GTDEDX,26.
******************************************************************************
+REP,GTPROGRM,GTDEDX,35-37.
      EXTERNAL  PROB
      DATA    AMS/0.511,105.,139.57,493.67,938.5/
      DATA    AP1/1.5/
*
*     *************************************************
+REP,GTPROGRM,GTDEDX,97.
101   CONTINUE
+REP,GTPROGRM,GTDEDX,103.
102   CONTINUE
+REP,GTPROGRM,GTINIT,1-7.
*CMZU:  1.34/02 06/07/95  18.51.00  by  Curtis A. Meyer
*CMZU:  1.33/01 09/02/95  14.22.43  by  Christian Voelcker
*CMZU:  1.33/00 12/01/95  15.22.56  by  Christian Voelcker
*CMZU:  1.30/01 01/07/94  16.19.06  by  Curtis A. Meyer
*CMZ :  1.21/10 18/03/93  16.16.40  by  Michael Doser
*CMZ :  1.21/09 04/12/92  11.27.13  by  Michael Doser
*CMZ :  1.21/05 04/12/91  17.16.00  by  Gunter Folger
*-- Author :
+REP,GTPROGRM,GTINIT,18.
*********************************************************
+REP,GTPROGRM,GTINIT,52-65.
***************************************************************
      INTEGER IER
      INTEGER IDATQQ,ITIMQQ,I
      INTEGER NTAKS(5),NTNKS(5),NTRKS(5),NTVKS(5)
* dolby c initialization variables. (see below)
      REAL E, ASYM, CPSI, EPI, ENONG
*
      DATA NTAKS/4HHTRK,   3,  2,  3,  0/
      DATA NTNKS/4HTTKS,   0,  0,  3,  0/
      DATA NTVKS/4HTVTX,   1,  1,  1,  0/
      DATA NTRKS/4HTTKS,   0,  0, 58,  0/
*
*     __________________________________________________________
*
*--> SET ERROR CODE
+REP,GTPROGRM,GTINIT,81-83.
 6000 FORMAT(' WELCOME TO GLOBAL TRACKING ')
C
C--->   INITIALIZE THE COUNTERS FOR THE TOTAL NUMBER OF TRACKS AND
+REP,GTPROGRM,GTINIT,94-117.
C---> counters for splitoff packages
C
* DolBy-c
      NDBGGT = 0      ! 'Good gammas'
      NDBTGT = 0      ! 'subTreshold peds'
      NDBSGT = 0      ! 'Split-offs'
* SMart
      NSMGGT = 0      ! 'Good gammas'
      NSMAGT = 0      ! 'Ambiguous peds'
      NSMSGT = 0      ! 'Split-offs'
* BRain

* TAxi
      NTAGGT = 0      ! 'Good gammas'
      NTASGT = 0      ! 'Split-offs (hadronic)'
      NTACGT = 0      ! 'Charged tracks with some splitoffs'

* SIngle Xtal Peds
      NSIXGT = 0      ! 'peds(cluster) with only one Xtal'
* all together
      NSPLGT = 0
*
C
C--->  initialize variables to allow Splitoff recognition info being
+REP,GTPROGRM,GTINIT,126.
C--->  DOLBY C needs to be initialized. But careful: USINIT is called later,
+REP,GTPROGRM,GTINIT,144.
C--->  CALL MZFORM FOR THE GLOBAL TRACK BANKS
+REP,GTPROGRM,GTINIT,155.
 100  CONTINUE
+REP,GTPROGRM,GTMTCH,1-10.
*CMZU:  1.34/02 06/07/95  18.51.34  by  Curtis A. Meyer
*CMZU:  1.21/08 16/01/92  07.45.18  by  Curtis A. Meyer
*CMZ :  1.21/06 04/12/91  17.41.01  by  Gunter Folger
*CMZ :          04/12/91  17.37.49  by  Michael Doser
*                  declare a number of variables which were undeclared
*                  (fix necessitated by introduction of IMPLICIT NONE)
*CMZ :          30/11/91  11.35.37  by  Michael Doser
*                  fix to enforce the 50 PED limit consistently
*                  in *all* loops
*CMZ :          05/12/90  06.55.33  by  MARK BURCHELL
*-- Author :
+REP,GTPROGRM,GTMTCH,18.
*****************************************************************
+REP,GTPROGRM,GTMTCH,32-34.
*************************************************************
C
C--->   ARRAYS FOR CRYSTAL PED'S, ALLOW UP TO 50 ENTRIES
+REP,GTPROGRM,GTMTCH,40-47.
C-->    SET UP ARRAYS (MAX SIZE 20) FOR TRACK INTERSECTION AT BC
C        XTRK,YTRK,ZTRK,IFLG  (IFLG = 0 IF GOOD)
      REAL XTRK(20),YTRK(20),ZTRK(20)
      INTEGER IFLG(20)
C-->    VOU and VOV are some temporary storage arrays.
        REAL VOV(3),VOU(3)
C
C-->    ARRAYS FOR MATCHED TRACK/CLUSTER INDICES, FOR MATCH I
+REP,GTPROGRM,GTMTCH,54-62.
C-->  THIS IS THE NUMBER OF TRACKS THAT MATCH PEDS, NOTE IT IS A MANY
C     TO ONE MAPPING SO IT IS NOT THE NUMBER OF PEDS THAT ARE MATCHED
      INTEGER IMDCT
C-->  OCCSIONALLY THE TCTX BANK DOES NOT EXIST FOR A TRACK,
C     IBTK FLAGS SUCH TRACKS
      INTEGER IBTK(20)
C  M.Doser, Dec4, 1991
C  the following variables had to be defined for the routine to compile
*  after introducing the IMPLICIT NONE statement.
+REP,GTPROGRM,GTMTCH,69-78.
C-->   USE ALINIT FOR INITIALIZING REAL ARRAYS
      REAL ALINIT
C
C--->   SET CUT FOR MATCHING
      REAL CMAT
*
      DATA ALINIT/-1000./
      DATA CMAT/0.98/
C*********************************************************************
C
C-->   INITIALIZE THE # OF PEDS (NPEDS), # OF SHOWERS (NSHW),
+REP,GTPROGRM,GTMTCH,94-149.
C--->   START BY LOOKING AT CRYSTALS: PARTICLE ENERGY DEPOSITIONS
C              (SEE CB NOTE 92)
C--->   USE THE TBTK BANK (STILL UNDER DEVELOPMENT)
C--->   FIRST CHECK BANK POINTER IS THERE, AND THAT BANK EXISTS
      IF (LTBTK.GT.0 ) THEN
        IF (IQ(LTBTK+1).GT.0) THEN
          DO 90 NPEDS = 1,IQ(LTBTK+1)
C--->  WE ARE NOW IN  A LOOP OVER PEDS
C--->  SET 50 AS MAX NUMBER OF PEDS WE WILL CONSIDER
          IF (NPEDS.GT.50) GOTO 90
C--->   GET LINK AND CHECK
          LINK = LQ(LTBTK-IQ(LTBTK-2)-NPEDS)
          IF (LINK.LE.0) GOTO 90
C--->  STORE DIRECTIONAL COSINES
          CUC(NPEDS) = Q(LINK+5)
          CVC(NPEDS) = Q(LINK+6)
          CWC(NPEDS) = Q(LINK+7)
C--->   STORE ID OF CLUSTER CONTAINING THE PED
          NCST(NPEDS) = IQ(LINK+3)
C--->   FIND NUMBER OF CRYSTALS IN HIT CLUSTER
          IF (LTBCL.GT.0) THEN
            LOFFS = LQ(LTBCL-NCST(NPEDS))
            IF (LOFFS.GT.0) NXTCGT(NPEDS) = IQ(LOFFS)
          ENDIF
 90       CONTINUE
        ENDIF
      ENDIF
C
C--->   MAKE SURE NPEDS IS SET CORRECTLY
      NPEDS = 0
      IF (LTBTK.GT.0) NPEDS = IQ(LTBTK+1)
C
C--->     WE NOW WANT TO LOOK AT THE TRACKS
C--->  ZERO/SET APPROPRIATE COUNTERS
          CALL VFILL(DANG,20,ALINIT)
          CALL VZERO(D2NG,20)
C
C
C--->     TRY THE TCTR BANK AND SEE WHAT HAPPENS
C
C     LOOP OVER ALL TRACKS
C
C--->   GET THE NUMBER OF TRACKS
      IF (LTCTR.GT.0) THEN
        NTKS=IQ(LTCTR+1)
        IF (NTKS.GT.0) THEN
C--->    CHECK FOR X,Y,Z COORDS IF INTERSECT WITH BARREL
          CALL TCBARL(NTRKS,IFLG,XTRK,YTRK,ZTRK)
C
C--->    LOOP OVER TRACKS AND FLAG FOR USE AS MATCHING THOSE WITH
C        TCTR ERROR CODE .LT. 30, #HITS .GE.3 AND ONLY USE FIRST 20
          DO 40 J=1,NTKS
            IF (J.GT.20) GOTO 40
C--->  DON'T MATCH IF HAS ERROR FLAG SET BY TCBARL
            IF (IFLG(J).GT.0) GOTO 40
C--->    1 EVENT IN 2000 HAS A TRACK WITH LINK = 0, CHECK FOR THIS
+REP,GTPROGRM,GTMTCH,157.
C--->   DON'T MATCH IF THE #HITS IS LT 4
+REP,GTPROGRM,GTMTCH,168-172.
 40       CONTINUE
        ENDIF
      ENDIF
C
C---> HERE WE HAVE SOME ASSOCIATION OF TRACKS AND PEDS
+REP,GTPROGRM,GTMTCH,183-188.
C-->  FOR THE ANGLE DELIBERATELY KEEP DIS1 AS A COSINE
               VOU(1) = CUC(I)
               VOU(2) = CVC(I)
               VOU(3) = CWC(I)
               DIS1 = VDOTN(VOV,VOU,3)
C--> SELECT PAIR WITH BEST DIS1
+REP,GTPROGRM,GTMTCH,195-196.
 70         CONTINUE
C-->  TEST FOR MATCHING !!
+REP,GTPROGRM,GTMTCH,212-216.
 60      CONTINUE
      ENDIF
C
C
C-->     NOW FLAG IF THE CRYSTAL IS IN THE SAME CLUSTER AS A MATCHED PED
+REP,GTPROGRM,GTMTCH,224-227.
 162     CONTINUE
 161  CONTINUE
C
C-->     NOW FOR EACH PED FILL THE COS(ANGLE) TO THE CLOSEST TRACK
+REP,GTPROGRM,GTMTCH,246-252.
  140      CONTINUE
         ENDIF
  130 CONTINUE
C
C--->  FIND THE TOTAL NUMBER OF PARTICLES IN THIS EVENT
      NGTKGT = NTKS + NPEDS - IMSCGT
C--->  UPDATE CUMULATIVE TOTALS
+REP,GTPROGRM,GTRACK,1-7.
*CMZU:  1.34/02 06/07/95  18.51.56  by  Curtis A. Meyer
*CMZU:  1.33/04 29/03/95  15.10.31  by  Curtis A. Meyer
*CMZU:  1.33/00 12/01/95  17.27.19  by  Christian Voelcker
*CMZ :  1.21/08 17/03/92  17.16.43  by  Michael Doser
*CMZ :  1.21/05 04/12/91  17.09.27  by  Gunter Folger
*-- Author :
      SUBROUTINE GTRACK(IER)
*****************************************************************
+REP,GTPROGRM,GTRACK,35.
*************************************************************
+REP,GTPROGRM,GTRACK,43-62.
      SAVE     IER001,IER002,IER003,IER004,IER005,IER006
      DATA     IER001,IER002,IER003,IER004,IER005,IER006 /6*0/
C
C_____________________________________________________________
C
C-->  if retracking is required, wipe all global tracking banks
C     (MZWIPE is slower than MZDROP, but also cleaner...)

      IF (RTRKGT.AND.LHTRK.GE.0) CALL MZWIPE(IXHTRK)
C
C-->  IF MATCHING IS REQUESTED CALL GTMTCH
C
      IF (MTCHGT.AND.(LHTRK.EQ.0))  CALL GTMTCH
C
C
C-->  CREATE BANKS FOR THE EVENT
C
      IF (MTCHGT.AND.(LHTRK.EQ.0)) CALL GTBANK
*
*  do all splitoff recognition at once!
+REP,GTPROGRM,GTRACK,68-71.
C-->  SET ERROR CODES
      IER = 0
C
C-->    Verify that the split-off codes are good:
+REP,GTPROGRM,GTRACK,78.
C--> Check Dolby C codes:
+REP,GTPROGRM,GTRACK,85.
C--> Check TAXI codes:
+REP,GTPROGRM,GTRACK,92.
C--> Check SMART codes:
+REP,GTPROGRM,GTRACK,103.
C--> Check BRAIN codes:
+REP,GTPROGRM,GTRACK,114.
 1000 CONTINUE
+REP,GT_TAXI,TAROND,1-50.
*CMZU:  1.34/02 06/07/95  18.53.18  by  Curtis A. Meyer
*CMZU:  1.33/01 07/02/95  10.48.19  by  Christian Voelcker
*CMZU:  1.33/00 23/01/95  11.37.15  by  Christian Voelcker
*-- Author :    Benayoun 13/12/94
*234567890
      SUBROUTINE TAROND(IVEC,IERR)
      IMPLICIT NONE
**************************************************************
*                                                            *
*  Modifications                                             *
*  23/12/94 CV  put common blocks into COMMON area,          *
*               change name of subroutines                   *
*               define all variables the proper way          *
*               fill the TAXI information into TTKS bank.    *
*                                                            *
*  18/01/95 CV+BNY                                           *
*               correct for the external xtal rings ...      *
*  24/01/95 CV                                               *
*               call bitopt to calculate theta/phi safely!   *
*                                                            *
**************************************************************


C-----------------------------------------------------------------
C     STUDY OF THE CRYSTALS
C-------------------------------------------------------------
C   ROUTINE TO FIND IF UNMATCHED PEDs FIT MATCHED CLUSTERS
C
CC-------------------------------------------------------------

ccccc on exit IVEC contains the NXTAL compound indices
ccccc of the cluster Xtals
C-----------------------------------------------------------
c  error codes : IERR=0 normal termination
C                IERR=1 normal termination but
c                  1 track has no matched cluster
C                IERR=2 normal termination but
C                  2 tracks or more have no matched cluster
C  abnormal codes :   IERR=3 no TBEN bank
C                     IERR=4 no TTKS
C                     IERR=5 number of charged tracks>20
C                     IERR=6 no TBTK bank
C                     IERR=7 wrong xtal coordinates in theta/phi
C-----------------------------------------------------------
C result stored in IVEC
C IVEC(i),i=1,ncharged : nbr of unmatched peds
C                        matching track i
C IVEC(i), i=ncharged+1,npart : info on unmatched peds
C          IVEC(i)=0 true unmatched ped
C          ivec(i)=k unmatched ped matching track nbr=k
C-----------------------------------------------------------
+REP,GT_TAXI,TAROND,56.
C--- CONTAINS COMMON/BCENER/ENERBC(60,26)
+REP,GT_TAXI,TAROND,78-84.
      INTEGER ICALL
      DATA LPREM /.TRUE./
      DATA ICALL /0/
*
* do only at first call...
      IF(LPREM) THEN
* call taxido to zero the full CONFTA array
+REP,GT_TAXI,TAROND,91.
* set all crystal flags to zero (has to be done for each event!)
+REP,GT_TAXI,TAROND,97-111.
* set the return vector to zero
      DO I=1,50
          IVEC(I)=0.
      ENDDO
*
C--------------------------------
      IERR=3
      IF(LTBEN.EQ.0) RETURN             !<<<<<
C--------------------------------
C  load crystal informations when above ECUTTA into CONFTA
      NTOT=IQ(LTBEN+1)
      DO I=1,NTOT
         IPOS=IQ(LTBEN+2*I)
CV       IY=IPOS/60+1
CV       IX=IPOS-(IY-1)*60
+REP,GT_TAXI,TAROND,123.
C  complete all six external rings
+REP,GT_TAXI,TAROND,133.
C--------------------------------
+REP,GT_TAXI,TAROND,139-151.
C--------------------------------
      IERR=5
      IF(NPION.GT.20) RETURN   ! too much is too much
C--------------------------------
C--------------------------------
      IERR=6
      IF(LTBTK.EQ.0) RETURN
C--------------------------------
CV      IF(LECR) WRITE(LLOG,*) 'ICALL',ICALL,
CV     &                       'NPION,NPHOTON',NPION,NPHOT
************************************
* now loop over all charged tracks
************************************
+REP,GT_TAXI,TAROND,159-160.
CV            J=IW/60+1
CV            I=IW-60*(J-1)
+REP,GT_TAXI,TAROND,173.
* now load the cluster information into the IRINDX array
+REP,GT_TAXI,TAROND,179-192.
CV            IF(LECR) THEN
CV               WRITE(LLOG,*) 'ipion=',IP, 'nbre of Xtals',MAXC
CV               WRITE(LLOG,*) 'central Xtal',IW
CV               WRITE(LLOG,1000)(IRINDX(K,IP),K=1,MAXC)
CV 1000          FORMAT(5(2X,I10))
CV            ENDIF
   80       CONTINUE
         ENDIF
      ENDDO
C--------------------------------------------------------
C-- after the loop on IP, all hadronic clusters have been
C   identified ; now check the unmatched peds within
C   the hadronic clusters.
C--------------------------------------------------------
+REP,GT_TAXI,TAROND,200-204.
CV         IF(LECR) THEN
CV            IPQ=IL-NPION
CV            WRITE(LLOG,*)' photon nbr=',IPQ
CV            WRITE(LLOG,*)' central Xtal',IWR
CV         ENDIF
+REP,GT_TAXI,TAROND,216-218.
 90      CONTINUE
      ENDDO
c--   end clustering.
+REP,GT_TAXI,TAROND,227-229.
CV      IF(LECR) THEN
CV         WRITE(LLOG,1000) (IVEC(I),I=1,NPART)
CV      ENDIF
+REP,GT_TAXI,TAXSTO,1-35.
*CMZU:  1.34/02 06/07/95  18.53.55  by  Curtis A. Meyer
*CMZ :  1.34/01 18/04/95  09.59.10  by  Christian Voelcker
*CMZU:  1.33/04 01/04/95  20.22.12  by  Christian Voelcker
*CMZ :          15/03/95  20.46.27  by  Christian Voelcker
*CMZU:  1.33/03 07/03/95  17.04.24  by  F.-H.Heinsius
*CMZU:  1.33/00 19/01/95  18.05.50  by  Christian Voelcker
*-- Author :    Christian Voelcker   19/01/95
      SUBROUTINE TAXSTO(IVEC)
      IMPLICIT NONE
****************************************************
*
* Creation date: January, 1995
*
* References:    none
*
* Description:   write the taxi results into the TTKS bank as follows
*
* TTKS word 3  = negative, if the ped is found to be a splitoff.
*                set bit 4 if it is a hadronic splitoff identfied by TAROND.
* TTKS word 53 = charged tracks: number of associated Ped's in that cluster
*                neutral tracks: pointer to parent
*
* IEHDCB(12) should contain the number of *good* gammas, but we need to
*            know the total number of splitoffs before (DOLBYC, etc.).
*
* TTKS word 12 we should leave as it is. Otherwise we may cause confusion on
*              the definition of cluster and cluster numbering!
* TTKS words 35-39 give info on the parent cluster. we do not want to
*                  recalculate all those values, as they are only useful
*                  for merged pions
*
*  Modifications:
*  07/03/95 FHH: TTKS word 53 was written for wrong particle (shifted by -1).
*  14/03/95 CV:  set bit 32  in TTKS WORD 3 to mark as splitoff
*  29/03/95 CM : Add error logging on crazy codes -- just a check.
*  18/04/95 CV : Typo in last modification (bracket)
*
****************************************************
+REP,GT_TAXI,TAXSTO,44-48.
      SAVE    IEBDCD,IETOMC
      DATA    IEBDCD,IETOMC / 0,0 /
*
* only count charged tracks and unmatched PEDs
* but nothing else (e.g. TTKS entries from PI0FND)
+REP,GT_TAXI,TAXSTO,59-67.
         IF(ABS(IVEC(I)).GT.50) THEN
            CALL ERRLOG(IEBDCD,
     &' <TAXSTO> Crazy TAXI error code seen, (abs(c)>50) ')
         ENDIF
CFHH     JTTKS=LQ(LTTKS-I)
         JTTKS=LQ(LTTKS-IQ(LTTKS-2)-I)
         IQ(JTTKS+53) = IVEC(I)
*
* is this a neutral PED and a e-m splitoff?
+ADD,$KUMACS,00_PATCH,*.
&PATCH,$KUMACS.
*CMZU:  1.34/02 19/07/95  20.20.04  by  Curtis A. Meyer
+ADD,$KUMACS,HELP,*.
&DECK,HELP.
*CMZU:  1.34/02 19/07/95  20.27.46  by  Curtis A. Meyer
*-- Author :    Curtis A. Meyer  02.Nov.1994
   Macro HELP

   mess 'Command template:'
   mess 'UNIX: cmz -install gtrack <files_directory> [ <option> ... ]'
   mess 'VMS:  CMZ/INSTALL=GTRACK,<files_directory>,[<option>,...]'
   mess ' '
   mess ' The files_directory must exist before calling cmz'
   mess ' '
   mess ' Compile and selection options '
   mess ' DEBUG    - debug option'
   mess ' LIB      - create object library file only (Default) '
   mess ' SOURCE   - create compiled source code only '
   mess ' GNU      - use the g77 and gcc compilers. '
   mess ' '
   mess 'example: cmz -install cbdrop src    (for Unix machines)'
   mess '         CMZ /INSTALL=cbdrop,[.SRC] (for VMS machines)'
   Return
+ADD,$KUMACS,INSTALL,*.
&DECK,INSTALL.
*CMZU:  1.35/00 24/08/95  16.54.36  by  Christian Voelcker
*CMZ :  1.34/03 10/08/95  22.00.10  by  Mark Lakata
*CMZU:  1.34/02 19/07/95  20.27.45  by  Curtis A. Meyer
*-- Author :    Curtis A. Meyer 12 July 1995
*
*   This macro is similar to the install macro developed by
*   Mark Lakata for the Crystal Barrel Event Display
*   Ported to gtrack on 12 July, 1995
*
   Macro INSTALL SRC=X
*
* [SRC] = SRC directory
* [2-4] = options
*
* options: ALL       = create everything
*          DEBUG     = Create Debug version of code.
*          LIB       = create only library.
*          GNU       = use the g77 and gcc compilers.
*          SOURCE    = create only compiled source code.
*
   if ( [SRC]='X' .OR. [SRC]='?' ) Then
      goto help
   Endif
*
*  CODEDIR    Source code directory
*
   if ( .not. ($FEXIST([SRC]) =1)) then
     mess Directory [SRC] does not exist.
     exitm
   endif

*   CODEDIR  = $Lower([SRC])
   message Putting source code files in [SRC]
*
*  set defaults
*
   BUILDLEVEL = LIB
   DEBUGFLAG  = FALSE
   OBJECTSTAT = DELETE
   GNUFLAG    = FALSE
   INCFLAG    = FALSE
   Select  .
*
*  check options
*
   FAIL = FALSE
   Do I = 2,[#]
      Case $Upper([%I]) In
         (DEBUG)       DEBUGFLAG = TRUE
                       message ' Debug compile options selected. '
         (ALL)         BUILDLEVEL = LIB
                       OBJECTSTAT = KEEP
         (LIB)         BUILDLEVEL = LIB
         (CAR)         car gtrack.car
         (SOURCE)      BUILDLEVEL = SOURCE
         (OBJECT)      BUILDLEVEL = OBJECT
                       OBJECTSTAT = KEEP
* Note: use OBJECT LIB to keep object code AND generate LIB
*       (order is important!)
         (GNU)         GNUFLAG   = TRUE
                       message '* Selecting g77 and gcc compilers. '
         (INC)         INCFLAG   = TRUE
                       message '* Creating #include files'
         (*)     message '****************************************'
                 message Error: Option $Upper([%I]) not understood.
                 message '****************************************'
                 FAIL = TRUE
      Endcase
   Enddo

   Message * Build Level = [BUILDLEVEL]

   If [FAIL] = TRUE then
     goto help
   Endif

   Select $MACHINE
   Case $Machine In
     (ALPHA) Case $OS In
        (VMS)  Select VAX
        (UNIX) Select UNIX
               Select DECS
     (IBMRT) Select $OS
     (SGI)   Select $OS
     Endcase
   Endcase
*
   Mess ' '
   Mess Installation of GTRACK on $MACHINE $OS starting at $DATE $TIME .
   Select
   Mess ' '
*
*  Set up the FORTRAN compiling extensions
*
*      (IBM)    ext = fortran
   if $OS = VMS Then
      ext = for
   else
      ext = f
   endif
*
   Set FORTRAN -LAN
*
* Set up filenames for the extraction of the fortran code.
*
   if [OBJECTSTAT] = KEEP Then
      XACL = -XACL
   else
      XACL = -XADCL
   endif

   Case $OS In
     (UNIX) Set [SRC]/mkgtrack_f.sh  [XACL]
            Set libgtrack.a              -L
     (VMS)  Set [SRC]mkgtrack_f.com  [XACL]
            Set GTRACK.OLB               -L
   Endcase

   IF [INCFLAG] = TRUE Then
     Case $OS In
       (UNIX) set '#include "$SEQ.inc"'   -INC F77
              FileCase Keep
              ext = F
       (VMS)  set '      INCLUDE 'CB$INC:$SEQ.INC''   -INC F77
     Endcase
   Else
     set . -inc
   EndIf

   Case $OS In
       (UNIX) Set [SRC]/*.[ext]         -F   F77
              set [SRC]/*.c             -F   C
       (VMS)  Set [SRC]*.[ext]          -F   F77
              set [SRC]*.c              -F   C
   Endcase
*
*  Set up the debug options
*
   If [DEBUGFLAG] = TRUE Then
      Case $MACHINE In
         (APOLLO)                  DEBUGOPT=' -dba'
         (HPUX,SGI,SUN,DECS,IBMRT) DEBUGOPT=' -g'
         (VAX)                     DEBUGOPT='/NOOP/DEB'
         (ALPHA) Case $OS In
                   (VMS)  DEBUGOPT='/NOOP/DEB'
                   (UNIX) DEBUGOPT=' -g'
                 Endcase
      Endcase
   Else
      DEBUGOPT=' '
   Endif
*
* Generate sequences
*
   seq .
*
* read some sequences from CBOFF
*
   Case $OS In
       (UNIX) cboffcmz = $CBCMZ/cboff.cmz
              mess 'Reading SEQUENCE definitions from $CBCMZ/cboff.cmz'
              mess 'Where $CBCMZ is ' $CBCMZ
       (VMS)  cboffcmz = CB$CMZ:cboff.cmz
              mess Reading SEQUENCE definitions from CB$CMZ:cboff.cmz
   Endcase
   if ($FEXIST([cboffcmz]) = 1) then
     File [cboffcmz] -R
   else
     mess [cboffcmz] ' does not exist.'
     exitm
   endif

   Seq //cboff/commcb
   Release cboff

*
* choose fortran compiler options and other stuff
*
   If [GNUFLAG] = TRUE THEN
     Set '`which g77` -c'//[DEBUGOPT]//' -O -w $compfile' -C f77
     Set '`which gcc` -c'//[DEBUGOPT]//' -O -w $compfile' -C C
   else
     Case $MACHINE In
     (HPUX) Set 'f77 -c'//[DEBUGOPT]//' +ppu $compfile' -C f77
           Set 'cc -c -O -w '//[DEBUGOPT]//'        $compfile' -C C
     (SUN) Set 'f77 -c'//[DEBUGOPT]//' -fnonstd -O -Nl100 -w $compfile' -C f77
          Set '`which gcc` -c'//[DEBUGOPT]//' -O -w $compfile' -C C
     (IBMRT) Set 'xlf -c'//[DEBUGOPT]//' -qextname $compfile' -C f77
            Set 'cc -c -O -w '//[DEBUGOPT]//'        $compfile' -C C
     (SGI)   Set 'f77 -c'//[DEBUGOPT]//' $compfile G 3' -C f77
            Set 'cc -c -O -w '//[DEBUGOPT]//'        $compfile' -C C
     (DECS)  Set 'f77 -c -O -w -Nl99 '//[DEBUGOPT]//' $compfile' -C f77
            Set 'cc -c -O -w '//[DEBUGOPT]//'        $compfile' -C C
     (VAX)   Set 'FOR/OBJ=$compfile.OBJ'//[DEBUGOPT]//' $compfile' -C f77
            Set 'CC/OBJ=$compfile.obj'//[DEBUGOPT]//'  $compfile' -C C
*
     (ALPHA) Case $OS In
      (VMS) Set 'FOR/SEPAR/OBJ=$compfile.OBJ'//[DEBUGOPT]//' $compfile' -C f77
           Set 'CC/OBJ=$compfile.obj'//[DEBUGOPT]//' $compfile' -C C
      (UNIX) Set 'f77 -c -O -w -Nl99 '//[DEBUGOPT]//' $compfile' -C f77
            Set 'cc -c -O -w '//[DEBUGOPT]//' $compfile' -C C
      Endcase
*
     Endcase
   Endif
*
*
*
*  make the library...
*
   Cd //gtrack
   Case [BUILDLEVEL] In
   (NONE)
        message 'no code extracted'
   (SOURCE)
        set F77 -LAN
        Cxtract *
        set C -LAN
        Cxtract *
   (OBJECT)
        set F77 -LAN
        CCompile *
        set C -LAN
        CCompile *
   (LIB)
        set F77 -LAN
        CLib *
        set C -LAN
        CLib *
   EndCase
*
         Mess ************************** READ THIS ******************
         Mess *                                                     *
         Mess * The following files were created:                   *
         Mess *                                                     *
Case [BUILDLEVEL] In
(LIB)  CASE ( $OS ) IN
  (UNIX) Mess * libgtrack.a -  the object library                   *
  (VMS)  Mess * gtrack.olb  -  the object library                   *
  ENDCASE
         mess * - The source and object files                       *
(OBJECT) mess * - The source and object files                       *
(SOURCE) mess * - The source code files                             *
(NONE)   mess * - nothing                                           *
         mess * --> Perhaps you want to use the LIB option?         *
         mess * --> Use "cmz -install gtrack" for help              *
EndCase
If ([OBJECTSTAT] = DELETE .and. [BUILDLEVEL] = LIB )  Then
         Mess * Object code was not kept.                           *
Endif
If [INCFLAG] = TRUE Then
         Mess *                                                     *
         mess * You may want to delete *.inc from SRC dir, since    *
         mess * they are also in CBINC                              *
Endif
         Mess *                                                     *
         Mess ************************** READ THIS ******************
*
ExitM

Help:

   mess 'Command template:'
   mess 'UNIX: cmz -install gtrack <files_directory> [ <option> ... ]'
   mess 'VMS:  CMZ/INSTALL=GTRACK,<files_directory>,[<option>,...]'
   mess ' '
   mess ' The files_directory must exist before calling cmz'
   mess ' '
   mess ' Compile and selection options '
   mess ' DEBUG    - debug option'
   mess ' LIB      - create object library (default)'
   mess ' SOURCE   - create source code only'
   mess ' OBJECT   - create source and object code only'
   mess ' ALL      - create all the files (see LIB, SOURCE and OBJECT)'
   mess ' GNU      - Use the g77 and gcc compilers '
   mess ' INC      - Use cpp style #include commands for +SEQ'
   mess ' '
   mess 'example: cmz -install gtrack srcgtrack    (for Unix machines)'
   mess '         CMZ /INSTALL=gtrack,[.SRCGTRACK] (for VMS machines)'

   Return