This file summarises changes to the sort subroutines. 

Changes from 1.0-alpha to 1.0-beta
-----------------------------------

Bugfix: New return flags RD_ABORT and RD_SKIP in simple-c.srt
        (The previous ABORT and SKIP did not handle return values correctly.)

Scanga: Hitpat now defined for particles 3 to n (previously fixed to max of 12)

Scanga: Input deck (see sort-subs/simulations/eg_input.sca) now
        incorporates detector position resolution entered as a fraction 
        for x and y (in-detector coords). e.g. For a 50mm x 50mm PSSD 
        with 200 um resolution in x and 3 mm in y the resolution values 
        are entered as 0.004 (x) 0.06 (y).
	New options for simulation including Punchthrough calculation.

Scanga: evhead.pmult[det] now incremented properly for each detector hit
missing_p: Put reconstructed missing particle in next slot in event array.
           Previously it was written to the next slot after singles hits.

Error return flags applied consistently.

punchthru(): Recalculates kinematics for a punchthru event.






Changes from birth to 1.0-alpha
--------------------------------

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.
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.
nrp.c  same operation but modified to handle N>6 
nrhit.c recalculates all dependencies for new mass assignment

simulations   scanga/rps/mcarlo/efficiency simulations

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

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

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

reac_init now loads only m0,m1,m2
coords tidied into sort-subs and calib-routines
calib-routines tidied with an example calibration procedure
bug fix in event_dump.c
updated simulation code to write hits into event array in random sequence
wrtev.c generalised
missing_p.c  subroutine to create a particle from the missing momentum
ridgeproj  routines to project out ridges in basel coordinates
bin : links to all executables now in bin directory
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. 
updated simulation code for multiple l' values in Basel coordinates
Makefile updated
updated simulation code for different k,k' contributions for each l'
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. 
Scanga updated to allow for widths of excited states.
missing_p updated to calculate missing mass NB previously assumed hitmass
Fixed error message in scanga_input
Modified scanga reaction input to handle all breakup permutations 
Tidied files in simulations and calib-routines
File version.h contains the sort-subs version number
corr now called with (ev,#particle number).
Scanga detection simulation now has option : Energy (T)hreshold only
reac_init() now called with reac_init(REAC_SIM) or reac_init(REAC_SORT)
Tidied scanga output
Scanga option added...
Y       #Put Pseudodata in random order (Y)es,(N)o?

vcm added to reac structure
plot_spec.c  routines to plot spectra
detprocs.c contains template detector_init/event procedure calls
rewrite of example calib-inputs and sort files
FAQ/NEWS To be finished!
Standalone scanga updated