Dear "beta-testers",

Finally....

Here is a trial version of the refinement of the rear segment response
matrices.  Right now it has some restrictions.  Now that the simulations
have all been run flexibly, most of these should be easy to change around.

RESTRICTIONS:
1) Single rear segments only -- no coincident modes
2) No front/rear anticoincidence (opposite of original DRMs)
3) As of now, no partial response elements -- i.e. the DRM you generate
     will contain Earth reflection, energy resolution effects, etc.  I
     haven't yet made it capable of shutting these off in the new mode

IMPROVEMENTS:
1) "input side" of the response matrix can be put up to 150MeV (soon to
    become 200 MeV and beyond if requested) reliably
2) Original normalization error in Compton component of direct
    photon interactions is fixed (this accounts for half of the
    "2.2 edge" problem, with coincidence mode accounting for the
    other half)
3) Default behavior is now *not* to have front/rear anticoincidence, which
    is suitable for large X flares with high deadtime in the fronts

STILL AVAILABLE:
1) You can still specify a date for radiation damage calculations OR specify
    a trapping parameter (lambda) by hand
2) You can still select what assortment of rear segments you want to use

NOT YET AVAILABLE (but available in a month or so if requested):
1) F/R anticoincident matrices (old mode but with other fixes included)
2) R/R sum pair modes (slight increase in sensitivity & photopeak)
3) Optimized matrices (removal of bits of response to correspond to
    removal of bits of spectrum in various detectors to get rid of
    artifacts) (with help from Albert once his dissertation is all done)

How to use:
1) Copy all the .pro files into ssw/hessi/idl/spectra
2) Copy all the .dat files intw ssw/hessi/dbase/spec_resp
3) Call hessi_build_srm with a /direct flag.  For example, this
produces a matrix using all the rear detectors except G2 and a big
energy binning scheme for nuclear spectroscopy, plus an extension to
very high energies (100 MeV) on the input side, with radiation
damage for 20 Jan 2005:

segs=[bytarr(9),1,0,bytarr(7)+1]
hsi_rd_ct_edges,20,es
ps=[es,11000.+1000.*findgen(90)]
hessi_build_srm,es,segs,srm,geo,/direct,ph_edges=ps,$
   time_wanted='2005-Jan-20 00:00:00'

Let me know if you have any questions.

Cheers,
David