c   A main program for S(200GeV/c)+S

      include 'jam1.inc'
      include 'jam2.inc'
      character frame*8,proj*8,targ*8,cwin*15

c....Initialize JAM
      mstc(1) =48827      ! random seed.
      mevent=100          ! total simulation event
      bmin=0.0D0          ! minimum impact parameter
      bmax=-1.0D0         ! maximum impact parameter
      dt=100.0D0          ! collision time(fm/c)
      nstep=1             ! time step
      cwin='200gevc'      ! incident energy
      frame='nn'          ! comp. frame
      proj='32S'          ! projectile
      targ='32S'          ! target

c...Initialize JAM.
      call jaminit(mevent,bmin,bmax,dt,nstep,frame,proj,targ,cwin)

c...Initialize analysis.
      call anal1

c...Simulation start.
      do iev=1,mevent

c...Simulate one event.
        call jamevt(iev)

c...Data analysis.
        call anal2

c...List phase space data.
        if(ive.eq.1)  call jamlist(1)

      end do

c...Final output.
      call jamfin

c...Print analysis results.
      call anal3

      end

c***********************************************************************

      subroutine anal1

      include 'jam1.inc'
      include 'jam2.inc'
      save ylab,wy

c....Rapidity distribution.
      ymin=-3.0D0
      ymax=7.0D0
      nymx=30
      wy=(ymax-ymin)/nymx
      call vbook1(11,'dN/dy - proton',nymx,ymin,ymax)
      call vbook1(12,'dN/dy - pion- ',nymx,ymin,ymax)
      call vbook1(13,'dN/dy - pion+ ',nymx,ymin,ymax)
      call vbook1(14,'dN/dy - charged',nymx,ymin,ymax)
      call vbook1(15,'dN/dy - p-p-bar',nymx,ymin,ymax)
      call vbook1(16,'dN/dy - h-',nymx,ymin,ymax)
      ylab=pard(17)

      return

c***********************************************************************

      entry anal2

c...Loop over all particles.
      do i=1,nv

        rap=0.5D0*log( max(p(4,i)+p(3,i),1.D-8)/max(p(4,i)-p(3,i), 
     & 1.D-8) )
        if(mstc(4).eq.0) then
        else if(mstc(4).eq.3) then
        else
         rap=rap+ylab
        endif

        kch=jamchge(k(2,i))
c...Charged particles.
        if(kch.ne.0) then
          call vfill1(14,rap,1.D0/wy)
c...Negative charged particles.
          if(kch.lt.0) call vfill1(16,rap,1.D0/wy)

        endif

        kf=k(2,i)

c.......Protons.
        if(abs(kf).eq.2212) then
          call vfill1(15,rap,isign(1,kf)/wy) 
          if(kf.eq.2212) call vfill1(11,rap,1.D0/wy) 

c.......Pions.
        else if(kf.eq.-211) then
          call vfill1(12,rap,1.D0/wy) 
        else if(kf.eq.211) then
          call vfill1(13,rap,1.D0/wy) 
        else if(kf.eq.111) then
        endif

      end do

      return

c***********************************************************************

      entry anal3

c...Output of histograms.
c...Event weight
      fac=1.D0/real(mstc(2))
      mnorm=0
      mform=0
      do i=1,6
       call vscale(10+i,fac)
       call vprint(10+i,mnorm,mform)
      end do

      end