batgrbproduct analysis for trigger=526642

batgrbproduct Analysis for Trigger=526642

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 03:16:28 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger526642
======================================================================

  Trigger: 00526642      Segment: 000

  BAT Trigger Time: 363993536   [s; MET]
      Trigger Stop: 363993600   [s; MET]
               UTC: 2012-07-14T21:18:46.573920    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 355.432702280664   Dec: -46.1858043791775  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 355.432702280664   Dec: -46.1858043791775  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 355.423600350055   Dec: -46.1957729965436  [deg; J2000]  
        { 23h 41m 41.7s , -46d 11' 44.8" }
         +/- 1.73190423869308 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.632058150241676 [arcmin] (formal 1-sigma fit error)
    SNR: 13.8793774623185
 
  Angular difference between Flight and Ground position is 0.707579 arcmin
 
  Partial Coding Fraction: 0.9296875    [ including projection effects ]

  Duration
     T90: 157.311999976635 +/- 23.9863005360993
        Measured from: -14.8080000281334
                   to: 142.503999948502     [s; relative to TRIGTIME]
     T50: 69.5239999890327 +/- 36.6918177621148
        Measured from: -2.38800001144409
                   to: 67.1359999775887     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -21.7080000042915     
                   to: -20.7080000042915     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -24.160000026226     
                   to: 151.248000025749     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.556728   0.890469   0.583334   0.032027
           0.069249   0.071346   0.062345   0.053216 [error]
 Peak      0.019426   0.014266   0.011642   0.001858
           0.005322   0.005469   0.004745   0.004165 [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.74e-07 3.13e-07 4.34e-07 1.25e-06 1.25e-06 2.18e-06   erg/cm2
              Error 2.64e-08 3.07e-08 5.84e-08 3.48e-07 1.41e-07 4.00e-07   erg/cm2
Peak flux Power-law 1.52e-01 1.18e-01 6.16e-02 4.63e-02 3.53e-01 3.78e-01   ph/cm2/s
              Error 7.11e-02 4.13e-02 3.81e-02 6.35e-02 1.21e-01 1.36e-01   ph/cm2/s
 
Fluence   Cutoff-PL 1.57e-07 3.46e-07 4.39e-07 3.62e-07 1.15e-06 1.30e-06   erg/cm2
              Error 3.00e-08 4.17e-08 6.69e-08 -------- 1.60e-07 4.33e-07   erg/cm2
Peak flux Cutoff-PL 1.54e-01 1.18e-01 6.05e-02 4.40e-02 3.53e-01 3.77e-01   ph/cm2/s
              Error 7.15e-02 3.40e-02 3.79e-02 6.21e-02 1.20e-01 1.36e-01   ph/cm2/s
 
Fluence        Band 1.51e-07 3.61e-07 4.11e-07 8.83e-07 1.19e-06 1.81e-06   erg/cm2
              Error 3.35e-08 4.62e-08 6.38e-08 -------- 1.47e-07 --------   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.529    0.163    -      -           7.54e-04 7.72e-05   51.990 57
CPL   0.688    0.754      71.2   59.4      2.16e-03 2.54e-03   47.400 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.941    0.730    -      -           2.42e-03 1.14e-03   54.930 57
CPL   1.928    1.880     181.2    0.0      2.45e-03 0.00e+00   54.930 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: 32.000000 sec.
T50 in the 50-300 keV band: 19.456000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100  175.408  -24.160  151.248
T90   157.312  -14.808  142.504
T50    69.524   -2.388   67.136
Peak    1.000  -21.708  -20.708

Click here for the BAT GCN circular template

Image

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

#   RAcent    DECcent   POSerr    Theta     Phi              Peak Cts       SNR   Name
  355.4236   -46.1958   0.0139    10.7754   106.7815         1.5977091     13.879 TRIG_00526642                 
 
Foreground time interval of the image:
-24.160  72.501 (delta_t =  96.661 [sec])
Background time interval of the image:
-394.799 -379.858 (delta_t =  14.941 [sec])
-239.000 -24.160 (delta_t = 214.840 [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 64.000000

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 64.000000

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 64.000000

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  -24.160 sec. to  151.248 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.52949     ( -0.162173   0.164338 )
Norm@50keV  : 7.54213E-04    ( -7.75211e-05   7.68799e-05 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          51.99 using 59 PHA bins.
# Reduced chi-squared =         0.9120 for     57 degrees of freedom 
# Null hypothesis probability =   6.631963e-01
 
Photon flux (15-150 keV)    in 175.4 sec: 9.49219E-02 ( -0.00919938  0.00920334 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.24914e-06 ( -1.40072e-07  1.42091e-07 ) ergs/cm2

Cutoff power-law model

Time interval is from  -24.160 sec. to  151.248 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.687663     ( -0.823942   0.68349 )
Epeak [keV] : 71.2384      ( -17.1039 101.72 )
Norm@50keV  : 2.15941E-03    ( -0.00122209   0.00385852 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          47.40 using 59 PHA bins.
# Reduced chi-squared =         0.8465 for     56 degrees of freedom 
# Null hypothesis probability =   7.864817e-01
 
Photon flux (15-150 keV)    in 175.4 sec: 9.18824E-02 ( -0.0095198  0.0095196 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.14956e-06 ( -1.57774e-07  1.62494e-07 ) ergs/cm2

Band function

Time interval is from  -24.160 sec. to  151.248 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 4.50821   ( -5.32965   -4.50821 )
beta        : -1.82077    ( -0.261045   0.231418 )
Epeak [keV] : 30.0867      ( -4.20645 4.2721 )
Norm@50keV   : 13.2309    ( -1.25874   88.536 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          44.81 using 59 PHA bins.
# Reduced chi-squared =         0.8147 for     55 degrees of freedom 
# Null hypothesis probability =   8.350908e-01
 
Photon flux (15-150 keV)    in 175.4 sec: 9.25059E-02 ( 6.90749  -0.0094009 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.1907e-06 ( -1.18947e-06  -1.45999e-07 ) ergs/cm2

Single BB

 
Spectral model blackbody:
------------------------------------------------------------
Parameters  : value    Lower 90%  Upper 90%
kT [keV]    : 13.1008  (-1.48115  1.72167)
R^2/D10^2   : 1.75254E-02  (-0.00625663   0.00940548)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          61.76 using 59 PHA bins.
# Reduced chi-squared =          1.084 for     57 degrees of freedom 
# Null hypothesis probability =   3.099405e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.119e-07   2.016e-08
25- 50  3.885e-07   4.639e-08
50-150  4.235e-07   9.182e-08
15-150  9.239e-07   1.164e-07

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 186.066       (      )
Norm        : 0.389889       (-0.0381887     0.0419013)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          49.19 using 59 PHA bins.
# Reduced chi-squared =         0.8629 for     57 degrees of freedom 
# Null hypothesis probability =   7.596233e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.707e-07   9.358e-08
25- 50  3.276e-07   1.750e-07
50-150  7.142e-07   3.814e-07
15-150  1.213e-06   6.510e-07

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 8.02669        (-2.17934      2.40725)
R1^2/D10^2  : 6.54180E-02        (-0.0346517      0.100029)
kT2 [keV]   : 25.9164        (-7.99578      28.1781)
R2^2/D10^2  : 1.04076E-03    (-0.000945762   0.00339848)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.13 using 59 PHA bins.
# Reduced chi-squared =         0.8206 for     55 degrees of freedom 
# Null hypothesis probability =   8.261713e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.549e-07   8.216e-08
25- 50  3.566e-07   1.685e-07
50-150  6.600e-07   3.550e-07
15-150  1.171e-06   5.160e-07

Peak spectrum fit

Power-law model

Time interval is from  -21.708 sec. to  -20.708 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.94096     ( -0.65744   0.803314 )
Norm@50keV  : 2.41851E-03    ( -0.00119053   0.00109948 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.93 using 59 PHA bins.
# Reduced chi-squared =         0.9637 for     57 degrees of freedom 
# Null hypothesis probability =   5.531860e-01
 
Photon flux (15-150 keV)    in 1 sec: 0.353275 ( -0.120499  0.120578 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.2251e-08 ( -9.63095e-09  1.06932e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from  -21.708 sec. to  -20.708 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.92831     ( -2.98646   0.774244 )
Epeak [keV] : 181.174      (   )
Norm@50keV  : 2.44850E-03    (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.93 using 59 PHA bins.
# Reduced chi-squared =         0.9810 for     56 degrees of freedom 
# Null hypothesis probability =   5.152650e-01
 
Photon flux (15-150 keV)    in 1 sec: 0.352962 ( -0.120282  0.120021 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.2088e-08 ( -1.04695e-08  1.0678e-08 ) ergs/cm2

Band function

Time interval is from  -21.708 sec. to  -20.708 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 4.97997   ( -4.97997   -4.97997 )
beta        : -2.15302    ( -1.84495   0.792592 )
Epeak [keV] : 21.6910      ( -21.691 9.15404 )
Norm@50keV   : 973.837    ( -327.466   5.55995e+06 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.45 using 59 PHA bins.
# Reduced chi-squared =         0.9899 for     55 degrees of freedom 
# Null hypothesis probability =   4.957786e-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]    : 7.67890  (-3.25109  7.16597)
R^2/D10^2   : 0.433008  (-0.415168   3.52835)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          58.20 using 59 PHA bins.
# Reduced chi-squared =          1.021 for     57 degrees of freedom 
# Null hypothesis probability =   4.307940e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.471e-09   2.917e-09
25- 50  7.279e-09   4.787e-09
50-150  1.667e-09   1.492e-09
15-150  1.342e-08   8.940e-09

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 70.8949       (-70.9474      -70.9474)
Norm        : 1.72376       (-0.725034     1.84986)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          55.32 using 59 PHA bins.
# Reduced chi-squared =         0.9705 for     57 degrees of freedom 
# Null hypothesis probability =   5.385027e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.428e-09   2.854e-09
25- 50  7.070e-09   4.488e-09
50-150  9.166e-09   8.075e-09
15-150  2.066e-08   1.494e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 4.99797        (-2.3597      4.33162)
R1^2/D10^2  : 1.98059        (-1.79671      30.6975)
kT2 [keV]   : 25.6691        (-14.7466      -25.6695)
R2^2/D10^2  : 3.90412E-03    (-0.00388985   0.0145949)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.17 using 59 PHA bins.
# Reduced chi-squared =         0.9849 for     55 degrees of freedom 
# Null hypothesis probability =   5.062751e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  5.250e-09   3.121e-09
25- 50  5.671e-09   3.342e-09
50-150  1.245e-08   8.220e-09
15-150  2.337e-08   1.254e-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: 1.48181     ( -0.158811   0.160078 )
Norm@50keV  : 1.07115E-03    ( -0.000105385   0.000104694 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          57.45 using 59 PHA bins.
# Reduced chi-squared =          1.008 for     57 degrees of freedom 
# Null hypothesis probability =   4.584637e-01
 
Photon flux (15-150 keV)    in 96.66 sec: 0.133078 ( -0.0125721  0.012579 ) ph/cm2/s
Energy fluence (15-150 keV)             : 9.85197e-07 ( -1.06787e-07  1.08193e-07 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.741558     ( -0.797822   0.665041 )
Epeak [keV] : 78.7982      ( -21.0461 229.043 )
Norm@50keV  : 2.67681E-03    ( -0.00147849   0.00444451 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          53.78 using 59 PHA bins.
# Reduced chi-squared =         0.9604 for     56 degrees of freedom 
# Null hypothesis probability =   5.592065e-01
 
Photon flux (15-150 keV)    in 96.66 sec: 0.129199 ( -0.013077  0.013059 ) ph/cm2/s
Energy fluence (15-150 keV)             : 9.15287e-07 ( -1.22598e-07  1.25045e-07 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 5.00000   ( -7.81706   -5 )
beta        : -1.72145    ( -0.272243   0.211833 )
Epeak [keV] : 28.7203      ( -4.67074 4.01196 )
Norm@50keV   : 50.7716    ( -4.70549   98.1373 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.19 using 59 PHA bins.
# Reduced chi-squared =         0.9125 for     55 degrees of freedom 
# Null hypothesis probability =   6.586163e-01
 
Photon flux (15-150 keV)    in 96.66 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 9.47609e-07 ( 0  0 ) ergs/cm2

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

S(25-50 keV)  = 3.13e-07
S(50-100 keV) = 4.34e-07

T90 vs. Hardness ratio plot

T90                                   = 157.311999976635 sec.
Hardness ratio (energy fluence ratio) = 1.38658

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.59947
Count Ratio (50-100 keV) / (15-25 keV) = 1.04779

Mask shadow pattern

IMX = -5.494803448911557E-02, IMY = -1.822095744744837E-01

TIME vs. PHA plot around the trigger time

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

TIME vs. DetID plot around the trigger time

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

100 us light curve (15-350 keV)

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

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:
-394.799390 -379.858100
-239.000000 -24.160000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
316.3440 -38.7086    2.4  38.9   81.6   3.6680   4.8 ------	UNKNOWN                       
  7.9791  -3.5970    2.1  47.7  171.7   2.1680   5.5  524.9	UNKNOWN

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

Time interval of the image:
-24.160000 72.500600
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
290.9333 -53.8374    2.6  45.5   48.5   2.3690   4.4 ------	UNKNOWN                       
355.4236 -46.1958    0.8  10.8  106.8   1.5977  13.9 ------	UNKNOWN                       
 64.0563 -45.0636    2.2  35.8  -83.9   1.2758   5.2 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
108.700600 166.300600
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
275.5987 -69.3635    3.1  44.1   24.9   0.5831   3.7 ------	UNKNOWN                       
326.0847 -61.1066    2.2  22.6   33.1   0.4418   5.3 ------	UNKNOWN                       
111.3841 -72.1484    2.6  53.6  -25.0   1.0720   4.4 ------	UNKNOWN                       
 66.6312 -52.2193    2.6  45.0  -60.1   1.5893   4.5 ------	UNKNOWN

Plot creation:
Sat May 30 19:10:10 EDT 2015