batgrbproduct analysis for trigger=452685

batgrbproduct Analysis for Trigger=452685

Analysis Summary

Notes:

1) No attempt has been made to eliminate the non-significant decimal places.
2) This batgrbproduct analysis is done using the HEASARC data.
3) This batgrbproduct analysis is using the flight position for the mask-weighting calculation.

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Mon Jan 13 22:32:32 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger452685
======================================================================

  Trigger: 00452685      Segment: 000

  BAT Trigger Time: 326136952.448   [s; MET]
      Trigger Stop: 326136953.472   [s; MET]
               UTC: 2011-05-03T17:35:45.530200    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 132.79984224966   Dec: 52.2156973681742  [deg; J2000]

     Catalogged Source?: NO
     Point Source?:      YES
     GRB Indicated?:     YES    [ by BAT flight software ]
     Image S/N Ratio:    20.4
     Image Trigger?:     NO
     Rate S/N Ratio:     134.029847422132  [ if not an image trigger ]
     Image S/N Ratio:    20.4

  Flight Position:   [ source = BAT ]
    RA: 132.79984224966   Dec: 52.2156973681742  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 132.799490954221   Dec: 52.211078499856  [deg; J2000]  
        { 08h 51m 11.9s , +52d 12' 39.9" }
         +/- 0.979525171125444 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.348119802711741 [arcmin] (formal 1-sigma fit error)
    SNR: 31.329483576278
 
  Angular difference between Flight and Ground position is 0.277433 arcmin
 
  Partial Coding Fraction: 0.1191406    [ including projection effects ]

  Duration
     T90: 58.6999999284744 +/- 47.0071018572133
        Measured from: -1.64800000190735
                   to: 57.0519999265671     [s; relative to TRIGTIME]
     T50: 3.91199994087219 +/- 0.540592264284623
        Measured from: -0.00400000810623169
                   to: 3.90799993276596     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.120000004768372     
                   to: 0.879999995231628     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -6.57600003480911     
                   to: 203.587999999523     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     2.994171   5.089695   5.349838   1.786912
           0.162915   0.183378   0.185515   0.195359 [error]
 Peak      0.499534   1.103982   1.507427   0.511213
           0.034995   0.047989   0.054829   0.042377 [error]

     [ fluence units of on-axis counts / fully illuminated detector ] 

======================================================================
 
=== Fluence/Peak Flux Summary ===
 
              Model    Band1    Band2    Band3    Band4    Band5    Total      
                       15-25    25-50   50-100  100-350   15-150   15-350   keV
Fluence   Power-law 1.16e-06 2.47e-06 4.12e-06 1.56e-05 1.13e-05 2.34e-05   erg/cm2
              Error 9.03e-08 1.12e-07 1.87e-07 1.64e-06 4.72e-07 1.73e-06   erg/cm2
Peak flux Power-law 6.74e+00 9.11e+00 9.07e+00 1.63e+01 3.02e+01 4.12e+01   ph/cm2/s
              Error 5.71e-01 4.86e-01 4.18e-01 1.53e+00 1.35e+00 1.92e+00   ph/cm2/s
 
Fluence   Cutoff-PL 1.10e-06 2.51e-06 4.22e-06 1.14e-05 1.12e-05 1.92e-05   erg/cm2
              Error 1.16e-07 1.23e-07 2.19e-07 4.11e-06 5.06e-07 3.51e-06   erg/cm2
Peak flux Cutoff-PL 5.34e+00 9.50e+00 1.02e+01 8.49e+00 2.96e+01 3.35e+01   ph/cm2/s
              Error 7.21e-01 5.46e-01 5.87e-01 1.94e+00 1.36e+00 2.33e+00   ph/cm2/s
 
Fluence        Band -------- -------- -------- -------- -------- --------   erg/cm2
              Error -------- -------- -------- -------- -------- --------   erg/cm2
Peak flux      Band -------- -------- -------- -------- -------- --------   ph/cm2/s
              Error -------- -------- -------- -------- -------- --------   ph/cm2/s

 
=== Time-averaged spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.262    0.070    -      -           5.41e-03 2.08e-04   68.770 57
CPL   1.022    0.277     229.4  -66.8      6.94e-03 1.90e-03   66.580 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.007    0.069    -      -           2.62e-01 1.15e-02   77.430 57
CPL   0.163    0.332     130.2   32.0      6.08e-01 1.99e-01   54.750 56
 
Model      : spectral model: simple power-law (PL) and cutoff power-law (CPL)
Ph_index   : photon index: E^-{Ph_index} (minus sign is not included in the index)
Ph_90err   : 90% error of the photon index
Epeak      : Epeak in keV
Epeak_90err: 90% error of Epeak
Norm       : normalization at 50 keV in PL model and at 1 keV in CPL model
Nomr_90err : 90% error of the normalization
chi2       : Chi-Squared
dof        : degree of freedom

 
=== T90/T50 calculation in the 50-300 keV band (BATSE band) using 64 ms lc ===
 
T90 in the 50-300 keV band: 7.616000 sec.
T50 in the 50-300 keV band: 1.984000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100  210.164   -6.576  203.588
T90    58.700   -1.648   57.052
T50     3.912   -0.004    3.908
Peak    1.000   -0.120    0.880

Click here for the BAT GCN circular template

Special notes of this burst

Image

Pre-slew 15.0-350.0 keV image (Event data)

#   RAcent    DECcent   POSerr    Theta     Phi              Peak Cts       SNR   Name
  132.7995    52.2111   0.0062    48.6614    23.5343        14.6476680     31.329 TRIG_00452685                 
 
Foreground time interval of the image:
 -6.576  15.453 (delta_t =  22.029 [sec])
Background time interval of the image:
-179.448  -6.576 (delta_t = 172.872 [sec])

Lightcurves

Notes:

1) The mask-weighted light curves are using the flight position.
2) Multiple plots of different time binning/intervals are shown to cover all scenarios of short/long GRBs, rate/image triggers, and real/false positives.
3) For all the mask-weighted lightcurves, the y-axis units are counts/sec/det where a det is 0.4 x 0.4 = 0.16 sq cm.
4) The verticle lines correspond to following: green dotted lines are T50 interval, black dotted lines are T90 interval, blue solid line(s) is a spacecraft slew start time, and orange solid line(s) is a spacecraft slew end time.
5) Time of each bin is in the middle of the bin.

1 s binning

From T0-310 s to T0+310 s

From T100_start - 20 s to T100_end + 30 s

Full time interval

64 ms binning

Full time interval

From T100_start to T100_end

From T0-5 sec to T0+5 sec

16 ms binning

From T100_start-5 s to T100_start+5 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

8 ms binning

From T100_start-3 s to T100_start+3 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

2 ms binning

From T100_start-1 s to T100_start+1 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

SN=5 or 10 sec. binning (whichever binning is satisfied first)

SN=5 or 10 sec. binning (T < 200 sec)

Quad-rate summed light curves (from T0-300s to T0+1000s)

Spectral Evolution

Spectra

Notes:

1) The fitting includes the systematic errors.
2) When the burst includes telescope slew time periods, the fitting uses an average response file made from multiple 5-s response files through out the slew time plus single time preiod for the non-slew times, and weighted by the total counts in the corresponding time period using addrmf. An average response file is needed becuase a spectral fit using the pre-slew DRM will introduce some errors in both a spectral shape and a normalization if the PHA file contains a lot of the slew/post-slew time interval.
3) For fits to more complicated models (e.g. a power-law over a cutoff power-law), the BAT team has decided to require a chi-square improvement of more than 6 for each extra dof.

Time averaged spectrum fit using the average DRM

Power-law model

Time interval is from   -6.576 sec. to  203.588 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.26180     ( -0.070416   0.0702077 )
Norm@50keV  : 5.40543E-03    ( -0.00020801   0.000207966 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          68.77 using 59 PHA bins.
# Reduced chi-squared =          1.207 for     57 degrees of freedom 
# Null hypothesis probability =   1.365087e-01
 
Photon flux (15-150 keV)    in 210.2 sec: 0.640564 ( -0.0261521  0.0262081 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.13362e-05 ( -4.71123e-07  4.72221e-07 ) ergs/cm2

Cutoff power-law model

Time interval is from   -6.576 sec. to  203.588 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.02157     ( -0.29482   0.258557 )
Epeak [keV] : 229.438      ( -95.8073 -229.438 )
Norm@50keV  : 6.93861E-03    ( -0.00139312   0.00240489 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          66.58 using 59 PHA bins.
# Reduced chi-squared =          1.189 for     56 degrees of freedom 
# Null hypothesis probability =   1.573656e-01
 
Photon flux (15-150 keV)    in 210.2 sec: 0.631773 ( -0.027912  0.027954 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.11913e-05 ( -5.08011e-07  5.04816e-07 ) ergs/cm2

Band function

Time interval is from   -6.576 sec. to  203.588 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.890752   ( 0.925538   1.25484 )
beta        : -1.43170    ( 1.41923   0.186875 )
Epeak [keV] : 153.318      ( -164.435 5055.69 )
Norm@50keV   : 8.14712E-03    ( -0.00216495   178.353 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          66.32 using 59 PHA bins.
# Reduced chi-squared =          1.206 for     55 degrees of freedom 
# Null hypothesis probability =   1.410943e-01
 
Photon flux (15-150 keV)    in 210.2 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

Single BB

 
Spectral model blackbody:
------------------------------------------------------------
Parameters  : value    Lower 90%  Upper 90%
kT [keV]    : 19.0034  (  )
R^2/D10^2   : 3.39356E-02  (   )
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =         180.88 using 59 PHA bins.
# Reduced chi-squared =         3.1733 for     57 degrees of freedom 
# Null hypothesis probability =   8.288099e-15
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  5.063e-07   3.857e-08
25- 50  2.436e-06   1.419e-07
50-150  6.464e-06   4.351e-07
15-150  9.406e-06   4.246e-07

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 199.363       (      )
Norm        : 2.69356       (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =         119.40 using 59 PHA bins.
# Reduced chi-squared =         2.0948 for     57 degrees of freedom 
# Null hypothesis probability =   2.576195e-06
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.410e-06   6.180e-08
25- 50  2.712e-06   1.146e-07
50-150  5.948e-06   2.649e-07
15-150  1.007e-05   4.109e-07

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 7.74999        (-1.39232      1.71865)
R1^2/D10^2  : 0.377148        (-0.168512      0.35248)
kT2 [keV]   : 29.8864        (-4.06145      6.94555)
R2^2/D10^2  : 6.34075E-03    (-0.00317255   0.00422818)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          71.74 using 59 PHA bins.
# Reduced chi-squared =          1.304 for     55 degrees of freedom 
# Null hypothesis probability =   6.422633e-02
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.027e-06   1.743e-07
25- 50  2.540e-06   2.852e-07
50-150  7.593e-06   1.132e-06
15-150  1.116e-05   1.163e-06

Peak spectrum fit

Power-law model

Time interval is from   -0.120 sec. to    0.880 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.00651     ( -0.0694166   0.0689762 )
Norm@50keV  : 0.262199    ( -0.0114741   0.0114609 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          77.43 using 59 PHA bins.
# Reduced chi-squared =          1.358 for     57 degrees of freedom 
# Null hypothesis probability =   3.720030e-02
 
Photon flux (15-150 keV)    in 1 sec: 30.1981 ( -1.34535  1.34493 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.82841e-06 ( -1.30377e-07  1.31092e-07 ) ergs/cm2

Cutoff power-law model

Time interval is from   -0.120 sec. to    0.880 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.162546     ( -0.344629   0.319419 )
Epeak [keV] : 130.196      ( -21.4148 42.5596 )
Norm@50keV  : 0.608028    ( -0.161537   0.237442 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.75 using 59 PHA bins.
# Reduced chi-squared =         0.9777 for     56 degrees of freedom 
# Null hypothesis probability =   5.221646e-01
 
Photon flux (15-150 keV)    in 1 sec: 29.6269 ( -1.3595  1.3613 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.79426e-06 ( 0  0 ) ergs/cm2

Band function

Time interval is from   -0.120 sec. to    0.880 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.142342   ( -0.302402   0.702159 )
beta        : -2.05101    ( 2.05101   0.742159 )
Epeak [keV] : 127.657      ( -43.7502 41.0838 )
Norm@50keV   : 0.622809    ( -0.146503   0.810112 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.75 using 59 PHA bins.
# Reduced chi-squared =         0.9954 for     55 degrees of freedom 
# Null hypothesis probability =   4.843029e-01
 
Photon flux (15-150 keV)    in 1 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

Single BB

 
Spectral model blackbody:
------------------------------------------------------------
Parameters  : value    Lower 90%  Upper 90%
kT [keV]    : 22.1738  (-1.12266  1.20583)
R^2/D10^2   : 1.11507  (-0.17337   0.201004)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          96.44 using 59 PHA bins.
# Reduced chi-squared =          1.692 for     57 degrees of freedom 
# Null hypothesis probability =   8.564212e-04
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.011e-07   8.300e-09
25- 50  5.371e-07   3.345e-08
50-150  1.973e-06   1.265e-07
15-150  2.611e-06   1.335e-07

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 199.363       (      )
Norm        : 121.935       (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =         232.54 using 59 PHA bins.
# Reduced chi-squared =         4.0797 for     57 degrees of freedom 
# Null hypothesis probability =   4.583793e-23
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.036e-07   1.410e-08
25- 50  5.841e-07   2.965e-08
50-150  1.281e-06   6.900e-08
15-150  2.169e-06   1.070e-07

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 8.73238        (-3.15548      4.94968)
R1^2/D10^2  : 7.61937        (-4.63553      23.0504)
kT2 [keV]   : 29.1082        (-4.07003      12.8984)
R2^2/D10^2  : 0.401704    (-0.294266   0.296505)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          57.35 using 59 PHA bins.
# Reduced chi-squared =          1.043 for     55 degrees of freedom 
# Null hypothesis probability =   3.881235e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.655e-07   4.713e-08
25- 50  5.558e-07   1.055e-07
50-150  2.078e-06   3.035e-07
15-150  2.799e-06   4.285e-07

Pre-slew spectrum (Pre-slew PHA with pre-slew DRM)

Power-law model

 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.28588     ( -0.05521   0.0550844 )
Norm@50keV  : 4.92201E-02    ( -0.00152273   0.00152241 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          63.08 using 59 PHA bins.
# Reduced chi-squared =          1.107 for     57 degrees of freedom 
# Null hypothesis probability =   2.701389e-01
 
Photon flux (15-150 keV)    in 22.03 sec: 5.85709 ( -0.194328  0.194713 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.07526e-05 ( -3.52654e-07  3.52806e-07 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.876181     ( -0.243769   0.230941 )
Epeak [keV] : 156.316      ( -37.7916 121.86 )
Norm@50keV  : 7.49929E-02    ( -0.0154969   0.0207206 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          53.59 using 59 PHA bins.
# Reduced chi-squared =         0.9569 for     56 degrees of freedom 
# Null hypothesis probability =   5.667269e-01
 
Photon flux (15-150 keV)    in 22.03 sec: 5.72366 ( -0.20641  0.20694 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.05477e-05 ( -3.77506e-07  3.75766e-07 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.682440   ( -0.356974   0.737121 )
beta        : -1.53957    ( 1.5393   0.171068 )
Epeak [keV] : 112.557      ( -51.3063 108.806 )
Norm@50keV   : 9.54556E-02    ( -0.0322085   0.185882 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          52.26 using 59 PHA bins.
# Reduced chi-squared =         0.9501 for     55 degrees of freedom 
# Null hypothesis probability =   5.801463e-01
 
Photon flux (15-150 keV)    in 22.03 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

S(25-50 keV) vs. S(50-100 keV) plot

S(25-50 keV)  = 2.47e-06
S(50-100 keV) = 4.12e-06

T90 vs. Hardness ratio plot

T90                                   = 58.6999999284744 sec.
Hardness ratio (energy fluence ratio) = 1.66802

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.69987
Count Ratio (50-100 keV) / (15-25 keV) = 1.78675

Mask shadow pattern

IMX = 1.042179667613151E+00, IMY = -4.538943176827844E-01

TIME vs. PHA plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

TIME vs. DetID plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

100 us light curve (15-350 keV)

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 1.024000

Spacecraft aspect plot

Bright source in the pre-burst/pre-slew/post-slew images

Pre-burst image of 15-350 keV band

Time interval of the image:
-179.448000 -6.576000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
 83.6230  22.0246    0.5  11.4  -58.6   4.6671  21.7    0.8	Crab                          
 87.1016 -22.1631    3.1  37.7 -178.7   1.0736   3.7 ------	UNKNOWN                       
 65.1533 -27.7141    2.3  50.6 -155.3   2.6267   4.9 ------	UNKNOWN                       
 54.9539 -25.8545    4.6  55.2 -145.2   3.5998   2.5 ------	UNKNOWN                       
 48.3354  42.2708    4.2  47.0  -53.3   2.6966   2.7 ------	UNKNOWN

Pre-slew background subtracted image of 15-350 keV band

Time interval of the image:
-6.576000 15.452600
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
 83.5862  21.9888    2.4  11.4  -58.9   0.4281   4.9    3.0	Crab                          
132.7995  52.2111    0.4  48.7   23.5  14.6455  31.3 ------	UNKNOWN                       
 75.2414  -3.1360    3.7  25.4 -142.1   0.2844   3.1 ------	UNKNOWN                       
 54.9816 -25.9982    2.3  55.3 -145.4   2.2284   5.0 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
82.252600 962.601700
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
135.2676  29.6390    3.3  22.7  163.0   2.3257   3.5 ------	UNKNOWN                       
 98.8527  62.2860    2.9  20.6  -58.9   2.5806   3.9 ------	UNKNOWN                       
104.8658   0.7831    3.9  56.5 -157.2  10.0416   3.0 ------	UNKNOWN                       
 96.1684   5.8413    7.2  55.3 -145.2   7.7016   1.6 ------	UNKNOWN

Plot creation:
Sat May 30 22:20:42 EDT 2015