2001.05.16  D.Huenemoerder

 Test new QE ARD - w/ adjusted filter, as per HLM.

 ARD:	     acisD1997-04-17qeN0004.fits
 Made with:  acis_qbar_ard.pro

Test gARFs made w/ mkgarf and new QE, for obsid 1198:

  acisf01198_001N001LEG_-1_garf.fits
  acisf01198_001N001LEG_1_garf.fits
  acisf01198_002N001LEG_-1_garf.fits
  acisf01198_002N001LEG_1_garf.fits

Requires psets to ardlib.par for mkgarf:

#---------------------------------------------------------------------------
# AXAF ACIS Mean QE Files.  These MUST include the FILTER transmission QE.
#---------------------------------------------------------------------------
AXAF_ACIS0_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE1]",,,"Enter ACIS-0 Mean QE File"
AXAF_ACIS1_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE2]",,,"Enter ACIS-1 Mean QE File"
AXAF_ACIS2_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE3]",,,"Enter ACIS-2 Mean QE File"
AXAF_ACIS3_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE4]",,,"Enter ACIS-3 Mean QE File"
AXAF_ACIS4_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE5]",,,"Enter ACIS-4 Mean QE File"
AXAF_ACIS5_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE6]",,,"Enter ACIS-5 Mean QE File"
AXAF_ACIS6_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE7]",,,"Enter ACIS-6 Mean QE File"
AXAF_ACIS7_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE8]",,,"Enter ACIS-7 Mean QE File"
AXAF_ACIS8_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE9]",,,"Enter ACIS-8 Mean QE File"
AXAF_ACIS9_QE_FILE,f,a,"./Ard/acisD1997-04-17qeN0004.fits[AXAF_QE10]",,,"Enter ACIS-9 Mean QE File"



ISIS/Slang script to apply to spectrum is Cmp_arf_fits.sl


Ps files:         Obs_1198_p1_oldqe.ps      
		  Obs_1198_m1_oldqe.ps      
		  Obs_1198_p1_newqe.ps      
		  Obs_1198_m1_newqe.ps      

		  Garf_aciss_letg-qeN0004.ps
                  Garf_aciss_letg-qeN0004.png

  First 4 show fits in ISIS of wabs*powerlaw (using xspec module for
  wabs) +/-1 fit simultaneously (made w/ ISIS Slang script)

  Last 2 (same figure) are comparison of old and new ARFs. (made w/
  IDL Plt_newqe.pro) 

  (quality/accuracy requires cal review)


Fit results:

% New QE fit:

isis> fit_counts;
 Parameters[Variable] = 3[3]
            Data bins = 2060
  Best fit chi-square = 1744
   Reduced chi-square = 0.848

isis> list_par;
wabs(1)*powerlaw(1)
 idx  param          tie-to  freeze         value          min          max
  1  wabs[1].NH_22       0     0        0.01230491            0            1
  2  powerlaw[1].norm    0     0        0.02296741            0            1
  3  powerlaw[1].alpha   0     0         -1.718404          -10            0

isis> conf(1,1);  % 90% limits
0.0134436
0.011494

isis> conf(2,1)(2,1);
0.0231878
0.022773

isis> conf(3,1);
-1.70709
-1.73111





% Old QE fit:

isis> fit_counts;
 Parameters[Variable] = 3[3]
            Data bins = 2060
  Best fit chi-square = 1389
   Reduced chi-square = 0.6752

isis> list_par;

wabs(1)*powerlaw(1)
 idx  param          tie-to  freeze         value          min          max
  1  wabs[1].NH_22       0     0        0.01238499            0            1
  2  powerlaw[1].norm    0     0        0.02200437            0            1
  3  powerlaw[1].alpha   0     0         -1.662157          -10            0

isis> conf(1,1);  % 90% limits
0.0134581
0.0113228

isis> conf(2,1);
0.0222007
0.0218083

isis> conf(3,1);
-1.65073
-1.67333