ngc1068_d_aTDET

Sun Aug  8 17:18:14 1999 : obs_anal.pro started

Inputs


 Filename        = /nfs/spectra/d0/mta/ngc1068_d.fits
 Output_prefix   = ngc1068_d
 Output dir.     = /spectra/d0/Obs_out
 Grating         = HETG
   cd_width      =  5.00000 mm
   cd_offset     =  0.00000 mm
 Detector        = ACIS-S
   Spatial coord = TDET
   (FLIPY is not set)
 ASPECT keyword set
 LINE_ANALYSIS keyword set


Getting Data ...


 Source of data: /nfs/spectra/d0/mta/ngc1068_d.fits

 Reading FITS file: /nfs/spectra/d0/mta/ngc1068_d.fits

* STREAK_MARX keyword: called marx_timing_sim

 Tag names are:
TIME CCD_ID EXPNO CHIPX CHIPY TDETX TDETY DETX DETY X Y PHAS PHA ENERGY PI
 FLTGRADE GRADE SHELL ORDER ZCOS YCOS XCOS ZPOS YPOS XPOS MARX_ENERGY

* Looks like ENERGY is in eV, converting to keV
 Creating ftime from TIME

 Total Events =  18881
 Start time =  37037173. seconds
 Interval duration =  49993.119 seconds
 Ave Event Rate =  0.37767197 events/second


Coordinates and Pixel Size...


 (Not CC mode data)

 Spatial coordinates are : TDETXb, TDETYb
 pixel size  = 24.0000 microns
 Processing method = aTDET


Showing and Selecting Events...


 Keeping TIME-in-range events,  100.000 %
  (from  0.0000000 to  49993.119)
 Start time =  37037173. seconds
 Interval duration =  49993.119 seconds
 No STATUS selection
 No QUALCODE selection
 GRADE distribution:
    0 :    18881,   100.00 %
 Keeping GRADE = 0,2,3,4,6 events,  100.000 %
 Keeping ENERGY selected events,  99.7299 %
  (from  0.200000 to  12.0000 keV)

Linear  binning from  0.200000 to  12.0000 keV
 dE of one bin width =  0.0146000 keV



Timing Analysis...


 Selected Events =  18830
 Start time =  37037173. seconds
 Interval duration =  49993.119 seconds
 Ave Event Rate =  0.37665183 events/second

* Found  74 big ( >  15.0000 s) gaps, totalling  1346.0200 seconds

 Live Interval =  48647.099 seconds
 Revised Event Rate =  0.38707344 events/second



Determine Zero-order...


 Average TDETXb,TDETYb gives x,y =  4107.79,  526.682
 Average TDETXb,TDETYb gives x,y =  4150.18,  528.589
 Average TDETXb,TDETYb gives x,y =  4176.52,  530.200

==> MTA Monitor Task 5.2 <==

 Zero-order Location TDETXb,TDETYb ave-median =  4177.77,  531.033

 Events in Zero-order Region =  9304
 Ave. Event Rate in Zero-order Region =  0.18610561 events/seconds

 TDETXb,TDETYb RMS Radius =  22.2167 pixels
 Number of TDETXb LRF pixels GT 10% of peak =  44
 Number of TDETXb LRF pixels GT 50% of peak =  20



Creating Aspect Solution...


 Aspect solution from TDETXb,TDETYb data...
 Calculating TDETXb,TDETYb averages every  60.0000 seconds.
 Apodizing during the first and last  8332.1865 seconds of data.
 Requires  833 bins ...
 Min / Max number of events per time bin =  2 /  21


 Making aspect corrected coordinates, AX,AY ...

 Before aspect correction:

 TDETXb,TDETYb RMS Radius =  22.2167 pixels
 Number of TDETXb LRF width at 10% of peak =  44

 Aspect correction yields (including apodized events):

 AX,AY RMS Radius =  22.019822 pixels
 Width of AX LRF at 10% of peak =   37.8 pixels



Detailed Examine Zero-order Spatial...


 Total Events =  12582
 Start time =  37045623. seconds
 Interval duration =  33319.840 seconds
 Ave Event Rate =  0.37761286 events/second

 The spatial coordinates to be used for grating spectra are:
   AX, AY  =  aTDETXb, aTDETYb

 AX,AY RMS Radius =  21.975907 pixels
 Width of AX LRF at 10% of peak =   950.4 microns
 Width of AX LRF at 50% of peak =   292.8 microns
 Gaussian fit FWHM =   502.4 microns
        FWHM error =     7.2 microns

 Sky X,Y (Level 1) Coordinates are present... showing them too.

 X,Y RMS Radius =  21.833883 pixels
 Width of LRF at 10% of peak =   105.6 microns
 Width of LRF at 50% of peak =    48.0 microns
 Gaussian fit FWHM =   811.1 microns
        FWHM error =   206.4 microns

 Width of Streak LRF at 10% of peak =   220.8 microns
 Width of Streak LRF at 50% of peak =   220.8 microns
 Gaussian fit FWHM =   414.6 microns
        FWHM error =   427.9 microns
 Creating ISIS format: ngc1068_d_aTDET_Zero_isis.dat ...

Measure Grating Angle(s)...

==> MTA Monitor Task 5.5 <==

 Selecting events in the grating parts...
   cd_width =  5.00000 mm
   cd_offset =  0.00000 mm

 HEG_minus number of events =  800
 HEG_minus slope is  -0.090741588, -->  -5.1849103 degrees on detector
 MEG_plus number of events =  2138
 MEG_plus slope is  0.081787183, -->  4.6756536 degrees on detector
 MEG_minus number of events =  2471
 MEG_minus slope is  0.081734093, -->  4.6726320 degrees on detector
 HEG_plus number of events =  789
 HEG_plus slope is  -0.091657742, -->  -5.2369690 degrees on detector
 HEG_all number of events =  1589
 HEG_all slope is  -0.090484757, -->  -5.1703148 degrees on detector
 MEG_all number of events =  4609
 MEG_all slope is  0.083151024, -->  4.7532680 degrees on detector



Crude Dispersed Spectra...



 Rowland spacing used is  8634.00 mm

 Use ACIS ENERGY for Order selection...
 MEG minus order selection keeps  90.1255 %,   2227 events
 MEG plus  order selection keeps  92.4228 %,   1976 events
 HEG minus order selection keeps  90.0000 %,   720 events
 HEG plus  order selection keeps  89.8606 %,   709 events

 Constant d-Lambda binning from  0.400000 to  10.0000 keV

 E/dE of one bin width at 1 keV =  999.596



Compare Zero-order and Dispersed Spectra...

==> MTA Monitor Task 5.7 <==

 Smoothing the spectra by E/dE ~  40.1563



Outputing Histograms of Spectra...


==> MTA Monitor Task 5.4 <==

 Creating ngc1068_d_aTDET_MEGm1_hist.rdb ...
 Creating ISIS format: ngc1068_d_aTDET_MEGm1_isis.dat ...
Spectrum Line Fits E/dE vs E Line list .txt or .rdb

 Creating ngc1068_d_aTDET_MEGp1_hist.rdb ...
 Creating ISIS format: ngc1068_d_aTDET_MEGp1_isis.dat ...
Spectrum Line Fits E/dE vs E Line list .txt or .rdb

 Creating ngc1068_d_aTDET_HEGm1_hist.rdb ...
 Creating ISIS format: ngc1068_d_aTDET_HEGm1_isis.dat ...
Spectrum Line Fits E/dE vs E Line list .txt or .rdb

 Creating ngc1068_d_aTDET_HEGp1_hist.rdb ...
 Creating ISIS format: ngc1068_d_aTDET_HEGp1_isis.dat ...
Spectrum Line Fits E/dE vs E Line list .txt or .rdb


Check Level 1.5 Columns...


 No Level 1.5 columns found in the file.


Output Histogram files from L1.5...


 No Level 1.5 columns - cannot create L1.5 histograms.


Listing of _summary.rdb file ...
Sun Aug  8 17:19:01 1999 : obs_anal.pro finished