batgrbproduct Analysis for Trigger=431128

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Fri Jan 17 16:47:06 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger431128
======================================================================

  Trigger: 00431128      Segment: 000

  BAT Trigger Time: 302865199.68   [s; MET]
      Trigger Stop: 302865200.704   [s; MET]
               UTC: 2010-08-07T09:13:13.622520    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 55.324189348915   Dec: 67.6787633437716  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 55.324189348915   Dec: 67.6787633437716  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 55.2995122773925   Dec: 67.6648710966458  [deg; J2000]  
        { 03h 41m 11.9s , +67d 39' 53.5" }
         +/- 1.80938584818832 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.773966194921213 [arcmin] (formal 1-sigma fit error)
    SNR: 13.0381712529179
 
  Angular difference between Flight and Ground position is 1.005581 arcmin
 
  Partial Coding Fraction: 0.3710938    [ including projection effects ]

  Duration
     T90: 7.74400001764297 +/- 1.43822114437232
        Measured from: -5.18400001525879
                   to: 2.56000000238419     [s; relative to TRIGTIME]
     T50: 2.75199997425079 +/- 1.11771195772738
        Measured from: -2.17599999904633
                   to: 0.575999975204468     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.340000033378601     
                   to: 0.659999966621399     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -6.0799999833107     
                   to: 3.13599997758865     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.274524   0.270979   0.093802   0.032216
           0.025555   0.026140   0.022505   0.019499 [error]
 Peak      0.060233   0.087174   0.058129   0.014330
           0.009229   0.009676   0.008161   0.006682 [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 9.72e-08 1.09e-07 8.71e-08 1.16e-07 3.36e-07 4.09e-07   erg/cm2
              Error 1.41e-08 1.30e-08 2.03e-08 5.13e-08 4.65e-08 7.78e-08   erg/cm2
Peak flux Power-law 6.54e-01 6.11e-01 3.96e-01 3.96e-01 1.82e+00 2.06e+00   ph/cm2/s
              Error 1.32e-01 8.00e-02 6.54e-02 1.28e-01 2.34e-01 2.61e-01   ph/cm2/s
 
Fluence   Cutoff-PL 9.43e-08 1.17e-07 7.71e-08 2.94e-08 3.09e-07 3.18e-07   erg/cm2
              Error 1.52e-08 1.91e-08 2.97e-08 7.89e-08 5.86e-08 1.07e-07   erg/cm2
Peak flux Cutoff-PL 5.79e-01 6.62e-01 4.21e-01 1.83e-01 1.78e+00 1.85e+00   ph/cm2/s
              Error 1.55e-01 1.00e-01 7.56e-02 1.73e-01 2.39e-01 2.92e-01   ph/cm2/s
 
Fluence        Band -------- -------- -------- -------- -------- --------   erg/cm2
              Error -------- -------- -------- -------- -------- --------   erg/cm2
Peak flux      Band 5.86e-01 6.55e-01 4.20e-01 1.93e-01 1.79e+00 1.86e+00   ph/cm2/s
              Error 1.52e-01 -------- 7.27e-02 -------- 2.38e-01 2.86e-01   ph/cm2/s
 
=== Time-averaged spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    2.321    0.234    -      -           3.80e-03 6.43e-04   49.100 57
CPL   1.587    0.000      20.4    0.0      1.10e-02 6.67e-02   47.950 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.627    0.208    -      -           1.41e-02 1.86e-03   57.320 57
CPL   0.865    0.912      72.9  -26.4      3.50e-02 4.99e-02   54.820 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: 1.152000 sec.
T50 in the 50-300 keV band: 0.640000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100    9.216   -6.080    3.136
T90     7.744   -5.184    2.560
T50     2.752   -2.176    0.576
Peak    1.000   -0.340    0.660
  • 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
       55.2995    67.6649   0.0148    32.9709  -110.1038         0.6500968     13.038 TRIG_00431128                 
     
    Foreground time interval of the image:
     -6.080   3.136 (delta_t =   9.216 [sec])
    Background time interval of the image:
    -239.680  -6.080 (delta_t = 233.600 [sec])
    

    Lightcurves

    Notes:

    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)

    The mask-weighted light curves based on cleaned DPHs

    Spectral Evolution

    Spectra

    Notes:

    Time averaged spectrum fit using the average DRM

    Power-law model
    Time interval is from   -6.080 sec. to    3.136 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 2.32094     ( -0.224106   0.243463 )
    Norm@50keV  : 3.79747E-03    ( -0.000648869   0.000636276 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.10 using 59 PHA bins.
    # Reduced chi-squared =         0.8615 for     57 degrees of freedom 
    # Null hypothesis probability =   7.622584e-01
     
    Photon flux (15-150 keV)    in 9.216 sec: 0.671458 ( -0.0742514  0.0742526 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.3576e-07 ( -4.57455e-08  4.71879e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -6.080 sec. to    3.136 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.58739     (     )
    Epeak [keV] : 20.3693      (   )
    Norm@50keV  : 1.09848E-02    ( -0.0108887   0.122526 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.95 using 59 PHA bins.
    # Reduced chi-squared =         0.8563 for     56 degrees of freedom 
    # Null hypothesis probability =   7.692493e-01
     
    Photon flux (15-150 keV)    in 9.216 sec: 0.655124 ( -0.078692  0.078516 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.09257e-07 ( -5.48251e-08  6.22909e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   -6.080 sec. to    3.136 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.970309   ( 0.81885   0.81885 )
    beta        : -2.73464    ( 2.79009   0.550737 )
    Epeak [keV] : 24.6113      ( -25.4544 7.43161 )
    Norm@50keV   : 2.97927E-02    ( -0.0427724   2299.69 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          46.76 using 59 PHA bins.
    # Reduced chi-squared =         0.8502 for     55 degrees of freedom 
    # Null hypothesis probability =   7.776114e-01
     
    Photon flux (15-150 keV)    in 9.216 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]    : 7.75697  (-0.876953  0.986013)
    R^2/D10^2   : 0.832092  (-0.321724   0.528818)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.54 using 59 PHA bins.
    # Reduced chi-squared =         0.9744 for     57 degrees of freedom 
    # Null hypothesis probability =   5.299443e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  8.147e-08   1.325e-08
    25- 50  1.351e-07   2.212e-08
    50-150  3.201e-08   1.253e-08
    15-150  2.485e-07   3.677e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 38.3726       (-9.64603      14.5539)
    Norm        : 4.45606       (-0.921139     1.22208)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          48.02 using 59 PHA bins.
    # Reduced chi-squared =         0.8425 for     57 degrees of freedom 
    # Null hypothesis probability =   7.956813e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  9.295e-08   1.303e-08
    25- 50  1.201e-07   1.668e-08
    50-150  8.871e-08   3.183e-08
    15-150  3.017e-07   4.649e-08
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 6.73359        (-1.59193      1.16077)
    R1^2/D10^2  : 1.33947        (-0.609286      0.833629)
    kT2 [keV]   : 33.1377        (-20.0112      -33.1348)
    R2^2/D10^2  : 1.49428E-03    (-0.00144485   0.0221046)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          46.76 using 59 PHA bins.
    # Reduced chi-squared =         0.8502 for     55 degrees of freedom 
    # Null hypothesis probability =   7.775290e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  8.844e-08   2.414e-08
    25- 50  1.236e-07   3.617e-08
    50-150  1.213e-07   6.408e-08
    15-150  3.333e-07   1.052e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from   -0.340 sec. to    0.660 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.62671     ( -0.205778   0.211092 )
    Norm@50keV  : 1.40772E-02    ( -0.00187235   0.00185385 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          57.32 using 59 PHA bins.
    # Reduced chi-squared =          1.006 for     57 degrees of freedom 
    # Null hypothesis probability =   4.631439e-01
     
    Photon flux (15-150 keV)    in 1 sec: 1.82427 ( -0.233919  0.234233 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.31245e-07 ( -1.83973e-08  1.85784e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -0.340 sec. to    0.660 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.864920     ( -1.00529   0.818519 )
    Epeak [keV] : 72.8581      ( -20.1085 -72.8581 )
    Norm@50keV  : 3.49781E-02    ( -0.0214037   0.0783594 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          54.82 using 59 PHA bins.
    # Reduced chi-squared =         0.9790 for     56 degrees of freedom 
    # Null hypothesis probability =   5.194040e-01
     
    Photon flux (15-150 keV)    in 1 sec: 1.78171 ( -0.23931  0.23912 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.25639e-07 ( -1.9395e-08  1.9637e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   -0.340 sec. to    0.660 sec.
    XSPEC12>exit
     XSPEC: quit
     
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       :    (     )
    beta        :     (     )
    Epeak [keV] :       (   )
    Norm@50keV   :     (     )
    ------------------------------------------------------------
     
    Photon flux (15-150 keV)    in 1 sec: 1.78838 ( 5.21162  -0.24451 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.26044e-07 ( 0  0 ) ergs/cm2
     
    
    Single BB
     
    Spectral model blackbody:
    ------------------------------------------------------------
    Parameters  : value    Lower 90%  Upper 90%
    kT [keV]    : 13.8776  (-1.99531  2.35764)
    R^2/D10^2   : 0.274705  (-0.119143   0.20592)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          65.92 using 59 PHA bins.
    # Reduced chi-squared =          1.157 for     57 degrees of freedom 
    # Null hypothesis probability =   1.955639e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.118e-08   2.697e-09
    25- 50  4.123e-08   7.365e-09
    50-150  5.205e-08   1.246e-08
    15-150  1.045e-07   1.974e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 161.862       (-68.7288      -160.312)
    Norm        : 7.54193       (-1.01083     1.3014)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.47 using 59 PHA bins.
    # Reduced chi-squared =         0.9732 for     57 degrees of freedom 
    # Null hypothesis probability =   5.325744e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.924e-08   1.042e-08
    25- 50  3.604e-08   1.948e-08
    50-150  7.334e-08   4.231e-08
    15-150  1.286e-07   7.025e-08
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 6.86186        (-3.38249      3.6413)
    R1^2/D10^2  : 1.93015        (-1.37818      23.4431)
    kT2 [keV]   : 21.8868        (-6.29645      19.4069)
    R2^2/D10^2  : 4.00492E-02    (-0.0364625   0.122913)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          54.92 using 59 PHA bins.
    # Reduced chi-squared =         0.9985 for     55 degrees of freedom 
    # Null hypothesis probability =   4.777323e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.785e-08   9.895e-09
    25- 50  3.759e-08   1.947e-08
    50-150  7.004e-08   3.989e-08
    15-150  1.255e-07   6.645e-08
    

    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: 2.32094     ( -0.224106   0.243463 )
    Norm@50keV  : 3.79747E-03    ( -0.000648869   0.000636276 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.10 using 59 PHA bins.
    # Reduced chi-squared =         0.8615 for     57 degrees of freedom 
    # Null hypothesis probability =   7.622584e-01
     
    Photon flux (15-150 keV)    in 9.216 sec: 0.671458 ( -0.0742514  0.0742526 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.3576e-07 ( -4.57455e-08  4.71879e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.58739     (     )
    Epeak [keV] : 20.3693      (   )
    Norm@50keV  : 1.09848E-02    ( -0.0108887   0.122526 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.95 using 59 PHA bins.
    # Reduced chi-squared =         0.8563 for     56 degrees of freedom 
    # Null hypothesis probability =   7.692493e-01
     
    Photon flux (15-150 keV)    in 9.216 sec: 0.655124 ( -0.078692  0.078516 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.09257e-07 ( -5.48251e-08  6.22909e-08 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.970309   ( 0.81885   0.81885 )
    beta        : -2.73464    ( 2.79009   0.550737 )
    Epeak [keV] : 24.6113      ( -25.4544 7.43161 )
    Norm@50keV   : 2.97927E-02    ( -0.0427724   2299.69 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          46.76 using 59 PHA bins.
    # Reduced chi-squared =         0.8502 for     55 degrees of freedom 
    # Null hypothesis probability =   7.776114e-01
     
    Photon flux (15-150 keV)    in 9.216 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)  = 1.09e-07
    S(50-100 keV) = 8.71e-08
    

    T90 vs. Hardness ratio plot

    T90                                   = 7.74400001764297 sec.
    Hardness ratio (energy fluence ratio) = 0.799083
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 0.987087
    Count Ratio (50-100 keV) / (15-25 keV) = 0.34169
    

    Mask shadow pattern

    IMX = -2.229676086232162E-01, IMY = 6.091630567530077E-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:
    -239.680000 -6.080000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    299.5777  35.1874    2.4  47.4  134.9   3.7875   4.9    1.1	Cyg X-1                       
    308.1237  40.9250    2.6  38.7  136.7   1.8325   4.4    1.8	Cyg X-3                       
    326.0509  38.3141    3.8  27.9  120.0   0.7734   3.0    5.7	Cyg X-2                       
     10.2984  10.2226    2.4  38.1   23.6   1.2393   4.8 ------	UNKNOWN                       
    

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

    Time interval of the image:
    -6.080000 3.136000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    299.5885  35.2000   10.8  47.4  134.9   0.1445   1.1    0.1	Cyg X-1                       
    308.1020  40.9323   10.3  38.7  136.7  -0.0640  -1.1    1.3	Cyg X-3                       
    326.3474  38.2922   13.9  27.8  119.7   0.0341   0.8    8.5	Cyg X-2                       
    333.6237  46.7562    2.2  19.5  134.0   0.1638   5.3 ------	UNKNOWN                       
    342.4993  58.2971    3.0  15.9  175.1   0.1060   3.8 ------	UNKNOWN                       
      2.7460  47.8787    1.9   0.2  -94.5   0.1709   5.9 ------	UNKNOWN                       
    313.0377  68.5979    2.2  31.9 -175.5   0.1856   5.3 ------	UNKNOWN                       
     43.6996   4.5294    2.5  55.7  -16.5   0.5140   4.7 ------	UNKNOWN                       
     55.2995  67.6649    0.9  33.0 -110.1   0.6500  13.0 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    70.020600 962.353600
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     83.6370  22.0202    0.5  48.9  -36.6  33.7629  25.3    0.4	Crab                          
    337.1348  56.4012    2.6  35.5  110.5   4.5130   4.4 ------	UNKNOWN                       
     96.3086  27.3740    3.6  47.2  -53.5   7.5974   3.3 ------	UNKNOWN                       
    
    Plot creation:
    Sun May 31 00:13:30 EDT 2015