This file summarises changes to the sort subroutines. 

Changes up to 4.07-beta
------------------------

CHANGES to sort-subs

24/6/97	
Corrected error in c.m. momentum calc in corr-sub.f and rp-sub.f
Started conversion to C
Moved fortran files to directory oldfortran (These are no longer supported)
Fortran callable copies of functions messy and mostly untested to date. 
Event-structure stuff reworked. PARTIAL file auto generation in a 
sub directory.

1/7/97
Can set and read the event structure from the fortran by using the
include file subs.i which points to ev.i.


New files...
rp.c    subroutine ``rp'' reconstructs resonant particle kinematics
nrp.c   subroutine ``nrp'' runs rp correlating particles in each n-hit event
rhit.c  subroutine ``rhit'' recalculates hit for new mass assignment
files.c handles input of data files
corr.c  subroutine ``corr'' reconstructs correlation phase space in various 
                            co-ordinate systems.

20/7/97
nrp.c  same operation but modified to handle N>6 
nrhit.c recalculates all dependencies for new mass assignment

17/9/97
simulations   scanga/rps/mcarlo/efficiency simulations

21/10/97
event_dump.c subroutine to format and dump sunsort event data
`z' charge of particle added to event structure

12/12/97
scanga simulation linked in with sunsort!
coords        coordinate definitions and routines 

1/2/98
``Pandoras box'' of next generation routines linked in
structure.c
properties.c
loadadcs.c
markhits.c
fortprop.c

17/3/98
reac_init now loads only m0,m1,m2

21/4/98
coords tidied into sort-subs and calib-routines
calib-routines tidied with an example calibration procedure

30/4/98
bug fix in event_dump.c

7/5/98
updated simulation code to write hits into event array in random sequence
wrtev.c generalised

19/5/98
missing_p.c  subroutine to create a particle from the missing momentum

26/5/98
ridgeproj  routines to project out ridges in basel coordinates
bin : links to all executables now in bin directory

9/6/98
nrp() now called as nrp(fold) to generate up to fold resonant particles only. 
The call "nrp(evhead.n[1])" will give the same operation as before.
chit(old,new): copies hit in ev array from old to new.
crhit(newmass,old,new): copies hit then recalculates the copy for newmass. 

22/6/98
updated simulation code for multiple l' values in Basel coordinates

24/6/98
Makefile updated
updated simulation code for different k,k' contributions for each l'

6/7/98: rp() q calc now independent of daughter grouping into hi and li
        implemented using new variable in event structure "ev[].etot" to
        keep running sum of daughter particle energies. 
