+PATCH,$CORR.
+DECK,CORR3.                     Updates version:  5.03/00  to   5.03/02
+REP,*TITLE*,TITLE,1.
CBGEANT  5.03/02  04/08/95  20.51.01
+REP,$VERSION,V5_03,1-2.
*CMZ :  5.03/02 04/08/95  20.51.01  by  Curtis A. Meyer
*CMZ :  5.03/01 17/03/95  10.00.09  by  Curtis A. Meyer
*CMZ :  5.03/00 15/03/95  15.26.48  by  Curtis A. Meyer
*-- Author :    Curtis A. Meyer   15/03/95
*
*::>          VERSION  5.03/02   04/08/95 20.51.01
*          *> Problem in BIDCAY that caused beam momentum to
*             be treated as -beam momentum has been fixed.
*             This was apparently ok in an earlier version.
*
*::>          VERSION  5.03/01   17/03/95 10.00.08
*          *> Inclusion of Gas Target details from
*             Ch. Strassburger.
*          *> Clean up the automatic intsallation.
*          *> Other small changes.
+ADD,$KUMACS,INSTALL,0.
*CMZ :  5.03/01 16/03/95  10.47.07  by  Curtis A. Meyer
+REP,$KUMACS,INSTALL,52.
   FXGEANT   = FALSE
+ADD,CBKINE,BIDCAY,0.
*CMZ :  5.03/02 04/08/95  18.53.23  by  Curtis A. Meyer
+ADD,CBKINE,BIDCAY,50.
C.    *                  C.Meyer        Aug 1995                       *
C.    *                  After looking at the GLOREN code, I belive    *
C.    *                  that VELO(3) should have a negative sign in   *
C.    *                  it. The GLOREN routine seems to be backwards  *
C.    *                  of common sense!!                             *
+REP,CBKINE,BIDCAY,88-90.
C
C-- C.A.Meyer, insert negative sign into VELO(3) term.
C
      VELO(1) = 0.
      VELO(2) = 0.
      VELO(3) = -(PBEAM+PTARG)/E12
+ADD,CBKINE,KIVERT,0.
*CMZU:  5.03/01 13/03/95  16.33.42  by  Christoph Strassburger
+ADD,CBKINE,KIVERT,49.
C *
C *         Modified 13 March 95 C. Strassburger
C *            Put the right Crystal Barrel TGHP-stop-distribution in
C *
+REP,CBKINE,KIVERT,62-70.
*        SIGXHP,Y,ZCB Crystal Barrel Stop distribution for high presssure
*                                           gas target
*
      REAL SIGX,SIGY,SIGZ
      REAL SIGXCB,SIGYCB,SIGZCB
      REAL SIGXHPCB,SIGYHPCB,SIGZHPCB
*
      SAVE SIGX,SIGY,SIGZ
      SAVE SIGXCB,SIGYCB,SIGZCB
      SAVE SIGXHPCB,SIGYHPCB,SIGZHPCB
*
      DATA SIGX,SIGY,SIGZ       / 2.300,2.300,5.710 /
      DATA SIGXCB,SIGYCB,SIGZCB / 0.280,0.280,0.400 /
      DATA SIGXHPCB,SIGYHPCB,SIGZHPCB / 0.400,0.400,0.700 /
+ADD,CBKINE,KIVERT,126.
*
      ELSEIF (IP.EQ.6) THEN
*
*          p-bar stop distribution for Crystal Barrel high press. gas target.
*
        VERTEX(1)=SIGXHPCB*XRN1
        VERTEX(2)=SIGYHPCB*XRN2
        VERTEX(3)=SIGZHPCB*XRN3
*
+ADD,CBMATE,MATERI,0.
*CMZU:  5.03/01 15/03/95  10.17.28  by  Christoph Strassburger
+REP,CBMATE,MATERI,462-464.
* H2 gas to H2 gas at 1 atm density :
*
       DATA PFACT /12.0/
+REP,CBMATE,MATERI,525.
      H2DENS = PFACT * 0.0838
+ADD,PW,GETGHP,0.
*CMZU:  5.03/01 14/03/95  11.52.56  by  Christoph Strassburger
+ADD,PW,GETGHP,21.
*.    *       modified : Ulrike Thoma and Christoph Strassburger       *
*.    *                  implementation of new high pressure           *
*.    *                  gas target according to a sketch              *
*.    *                  of Wolfgang Duennweber                        *
*.    *                                                                *
+REP,PW,GETGHP,28-34.
      CALL ERRWRN(0,' TGHP GEOMETRY/MATERIALS from a sketch of WD$')
*
*  Target
*
      PAR(1) = 0.
      PAR(2) = 0.14
      PAR(3) = 5.9
+REP,PW,GETGHP,45-50.
*--  outer carbon-mylar vessel
*
      PAR(1) = 0.14
      PAR(2) = 0.15
      PAR(3) = 6.0
      CALL GSVOLU('THPO','TUBE',54,PAR,3,IVOL)
+REP,PW,GETGHP,57-62.
      PAR(2) = 0.15
      PAR(3) = 0.05
      CALL GSVOLU('THPE','TUBE',53,PAR,3,IVOL)
      IF(IVOL.LE.0) CALL ERRMES(IVOL,'GETGHP : GSVOLU THPE failed$')
      CALL  GSPOS('THPE',1,'TGHP',0.,0.,ZTARG-5.95,0,'ONLY')
      CALL  GSPOS('THPE',2,'TGHP',0.,0.,ZTARG+5.95,0,'ONLY')
+REP,OBSOLETE,00_PATCH,0-1.
&PATCH,OBSOLETE,IF=NEVER.
*CMZ :  5.03/01 16/03/95  08.59.09  by  Curtis A. Meyer
+ADD,DFILES,CARDFILE,0.
*CMZ :  5.03/01 16/03/95  09.20.27  by  Curtis A. Meyer
+REP,DFILES,CARDFILE,8-21.
RNDM     14  0

PRINT    'VERT' 'KINE'
DEBUG    1 10

KINE     3   1.  0.0  3.  7. 9. 8.
SETV     91  0. 0. 0.0
BWID     0.1 0.1 0.1
RAWD     20 'EVHD' 'RPWC' 'RJDC' 'RBCL' 'RBCF' 'RMCB'

JDCD     1
JDC      'jdc1_15.dat'
LUTF     'mcfera.lut'
LUT2     'mc2282.lut'

END

STOP
STOP  <===Explanation Below This Line ===

RUNG     1    ==  Run Number 1
TRIG   100    ==  Produce 100 Events
RNDM   14 0   ==  Use internal random sequence 14

PRINT 'VERT' 'KINE' == Debug info is VERT and KINE banks.
DEBUG  1 10   ==  Print out Debug info on events 1 to 10.

SETV   91  0.  0. 0. == Vertex type 91 at (0,0,0)
BWID   0.1 0.1 0.1   == Vertex (sig_x,sig_y,sig_z)