#!/bin/tcsh
# subsequent generation of "phase space" events with
# same vertex coordinates. These events can then be
# used for mixing with the original event file
#
# Maxim Potekhin, potekhin@bnl.gov

#*****************************************************
# Input for gkine
set geometry=year2003
set pid=5
set tracks_per_event=1
set ptlow=10.0
set pthigh=10.0
set ylow=-1.0
set yhigh=1.0
set philow=0.0
set phihigh=6.28
#*****************************************************
echo '******************************************************************************************************'
echo Commencing the generation of phase space tracks: geometry=$geometry, particle ID=$pid
echo Kinematic parameters: Pt_Low=$ptlow, Pt_High=$pthigh, Y_low=$ylow, Y_high=$yhigh, Phi_low=$philow, Phi_high=$phihigh
echo                 To change these parameters please edit this script, $0
echo '******************************************************************************************************'
# Form filenames:
set zfile=`basename $1`                      # grab the filename
set zfile=`echo $zfile | sed -e 's/\..*//g'` # lose extenstion
set vertex_file=$zfile.vertex                # file with the z coordinates
set phase_file=$zfile.phase                  # output file (phase space events)

echo Starting event production
# Generate events as defined by the user variables
staf -w 0 -b phase_gen $vertex_file $geometry $tracks_per_event $pid $ptlow $pthigh $ylow $yhigh $philow $phihigh $phase_file