batgrbproduct analysis for trigger=143708

batgrbproduct Analysis for Trigger=143708

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): Wed Jan 22 21:49:13 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger143708
======================================================================

  Trigger: 00143708      Segment: 000

  BAT Trigger Time: 141910979.072   [s; MET]
      Trigger Stop: 141910980.096   [s; MET]
               UTC: 2005-07-01T11:42:59.472220    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 227.286665301381   Dec: -59.4186472467233  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 227.286665301381   Dec: -59.4186472467233  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 227.244707200737   Dec: -59.4002520490514  [deg; J2000]  
        { 15h 08m 58.7s , -59d 24' 00.9" }
         +/- 2.12745724446889 [arcmin] (estimated 90% radius based on SNR)
         +/- 1.05145821860101 [arcmin] (formal 1-sigma fit error)
    SNR: 10.345363389834
 
  Angular difference between Flight and Ground position is 1.691011 arcmin
 
  Partial Coding Fraction: 0.6796875    [ including projection effects ]

  Duration
     T90: 21.6960000097752 +/- 9.94378159478707
        Measured from: -3.64800000190735
                   to: 18.0480000078678     [s; relative to TRIGTIME]
     T50: 7.74400001764297 +/- 0.641760070223944
        Measured from: -0.0119999945163727
                   to: 7.7320000231266     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: 7.36399999260902     
                   to: 8.36399999260902     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -4.73199999332428     
                   to: 28.4079999923706     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.701489   0.925556   0.564948   0.165227
           0.031706   0.032800   0.028640   0.024377 [error]
 Peak      0.090507   0.124604   0.097765   0.039742
           0.006585   0.007245   0.006473   0.005329 [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 2.23e-07 3.66e-07 4.56e-07 1.12e-06 1.36e-06 2.17e-06   erg/cm2
              Error 1.73e-08 1.68e-08 2.89e-08 1.55e-07 6.74e-08 1.77e-07   erg/cm2
Peak flux Power-law 8.04e-01 8.90e-01 7.04e-01 9.19e-01 2.74e+00 3.32e+00   ph/cm2/s
              Error 9.58e-02 6.24e-02 5.51e-02 1.56e-01 1.74e-01 2.16e-01   ph/cm2/s
 
Fluence   Cutoff-PL 2.19e-07 3.70e-07 4.58e-07 9.35e-07 1.35e-06 1.98e-06   erg/cm2
              Error 1.55e-08 1.97e-08 2.98e-08 3.84e-07 7.78e-08 3.13e-07   erg/cm2
Peak flux Cutoff-PL 8.13e-01 8.94e-01 7.00e-01 8.97e-01 2.74e+00 3.30e+00   ph/cm2/s
              Error 9.65e-02 6.26e-02 5.55e-02 1.86e-01 1.31e-01 2.36e-01   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.684    0.084    -      -           4.43e-03 2.07e-04   61.470 57
CPL   1.573    0.260     186.1    0.0      5.06e-03 0.00e+00   61.110 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.339    0.116    -      -           2.28e-02 1.43e-03   56.930 57
CPL   1.335    0.161    9999.3    0.0      2.29e-02 3.23e-03   56.980 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: 8.064000 sec.
T50 in the 50-300 keV band: 4.608000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100   33.140   -4.732   28.408
T90    21.696   -3.648   18.048
T50     7.744   -0.012    7.732
Peak    1.000    7.364    8.364

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
  227.2447   -59.4003   0.0186    19.4777    82.1230         1.5764245     10.345 TRIG_00143708                 
 
Foreground time interval of the image:
 -4.732  28.408 (delta_t =  33.140 [sec])
Background time interval of the image:
-10.072  -4.732 (delta_t =   5.340 [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

8 ms binning

From T100_start-3 s to T100_start+3 s

2 ms binning

From T100_start-1 s to T100_start+1 s

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   -4.732 sec. to   28.408 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.68442     ( -0.0840815   0.0848213 )
Norm@50keV  : 4.43188E-03    ( -0.000207115   0.000206145 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          61.47 using 59 PHA bins.
# Reduced chi-squared =          1.079 for     57 degrees of freedom 
# Null hypothesis probability =   3.189229e-01
 
Photon flux (15-150 keV)    in 33.14 sec: 0.585436 ( -0.026887  0.0269378 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.36161e-06 ( -6.7189e-08  6.75642e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from   -4.732 sec. to   28.408 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.57294     ( -0.333402   0.185948 )
Epeak [keV] : 186.079      (   )
Norm@50keV  : 5.05546E-03    (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          61.11 using 59 PHA bins.
# Reduced chi-squared =          1.091 for     56 degrees of freedom 
# Null hypothesis probability =   2.974469e-01
 
Photon flux (15-150 keV)    in 33.14 sec: 0.581268 ( -0.029147  0.029081 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.34781e-06 ( -7.86313e-08  7.69643e-08 ) ergs/cm2

Band function

Time interval is from   -4.732 sec. to   28.408 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.368773   ( 0.95908   0.95908 )
beta        : -1.74606    ( 1.73635   0.107431 )
Epeak [keV] : 35.7848      ( -35.8359 221.307 )
Norm@50keV   : 3.55197E-02    ( -0.0365645   366.273 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          59.46 using 59 PHA bins.
# Reduced chi-squared =          1.081 for     55 degrees of freedom 
# Null hypothesis probability =   3.164805e-01
 
Photon flux (15-150 keV)    in 33.14 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]    : 13.2401  (  )
R^2/D10^2   : 9.56434E-02  (   )
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =         185.00 using 59 PHA bins.
# Reduced chi-squared =         3.2456 for     57 degrees of freedom 
# Null hypothesis probability =   1.943303e-15
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.178e-07   9.478e-09
25- 50  4.137e-07   2.071e-08
50-150  4.634e-07   4.557e-08
15-150  9.949e-07   5.286e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 147.618       (-31.6668      38.6801)
Norm        : 2.38237       (-0.0648377     0.157814)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          63.25 using 59 PHA bins.
# Reduced chi-squared =          1.110 for     57 degrees of freedom 
# Null hypothesis probability =   2.652248e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.030e-07   1.415e-08
25- 50  3.755e-07   3.369e-08
50-150  7.362e-07   2.089e-07
15-150  1.315e-06   2.522e-07

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 6.80600        (-0.923402      1.0127)
R1^2/D10^2  : 0.752713        (-0.282837      0.514771)
kT2 [keV]   : 25.9953        (-4.03541      6.33158)
R2^2/D10^2  : 6.53899E-03    (-0.00348317   0.00569175)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          62.69 using 59 PHA bins.
# Reduced chi-squared =          1.140 for     55 degrees of freedom 
# Null hypothesis probability =   2.223187e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.050e-07   2.701e-08
25- 50  3.775e-07   3.889e-08
50-150  7.521e-07   1.572e-07
15-150  1.335e-06   1.748e-07

Peak spectrum fit

Power-law model

Time interval is from    7.364 sec. to    8.364 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.33944     ( -0.116605   0.116288 )
Norm@50keV  : 2.27875E-02    ( -0.00143074   0.0014267 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.93 using 59 PHA bins.
# Reduced chi-squared =         0.9988 for     57 degrees of freedom 
# Null hypothesis probability =   4.776745e-01
 
Photon flux (15-150 keV)    in 1 sec: 2.73934 ( -0.174039  0.174364 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.23093e-07 ( -1.58984e-08  1.60087e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from    7.364 sec. to    8.364 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.33485     ( -0.202994   0.118172 )
Epeak [keV] : 9999.30      ( -9999.35 -9999.35 )
Norm@50keV  : 2.28723E-02    ( -0.0014368   0.00501687 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.98 using 59 PHA bins.
# Reduced chi-squared =          1.017 for     56 degrees of freedom 
# Null hypothesis probability =   4.385606e-01
 
Photon flux (15-150 keV)    in 1 sec: 2.74244 ( -0.09094  0.17116 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.21768e-07 ( -1.4692e-08  1.7018e-08 ) ergs/cm2

Band function

Time interval is from    7.364 sec. to    8.364 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 4.95857   ( -4.95857   -4.95857 )
beta        : -1.36600    ( -0.124874   0.123293 )
Epeak [keV] : 22.1389      ( -22.1389 5.12092 )
Norm@50keV   : 4046.02    ( -250.723   9.19341e+06 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.59 using 59 PHA bins.
# Reduced chi-squared =          1.029 for     55 degrees of freedom 
# Null hypothesis probability =   4.155650e-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]    : 16.6413  (  )
R^2/D10^2   : 0.214243  (   )
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =         149.43 using 59 PHA bins.
# Reduced chi-squared =         2.6216 for     57 degrees of freedom 
# Null hypothesis probability =   3.205446e-10
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.215e-08   1.595e-09
25- 50  5.285e-08   4.820e-09
50-150  1.038e-07   1.280e-08
15-150  1.688e-07   1.385e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 200.000       (      )
Norm        : 11.2574       (-0.704725     0.70467)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          70.49 using 59 PHA bins.
# Reduced chi-squared =          1.237 for     57 degrees of freedom 
# Null hypothesis probability =   1.081203e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.804e-08   1.915e-09
25- 50  5.393e-08   5.480e-09
50-150  1.183e-07   2.434e-08
15-150  2.003e-07   3.140e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 7.37109        (-1.07979      1.22034)
R1^2/D10^2  : 2.53785        (-1.03799      1.91754)
kT2 [keV]   : 39.3704        (-8.26166      15.4077)
R2^2/D10^2  : 1.16936E-02    (-0.00679329   0.0118173)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.25 using 59 PHA bins.
# Reduced chi-squared =         0.9136 for     55 degrees of freedom 
# Null hypothesis probability =   6.564427e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.566e-08   4.865e-09
25- 50  5.125e-08   9.075e-09
50-150  1.575e-07   5.176e-08
15-150  2.344e-07   6.205e-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.68442     ( -0.0840815   0.0848213 )
Norm@50keV  : 4.43188E-03    ( -0.000207115   0.000206145 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          61.47 using 59 PHA bins.
# Reduced chi-squared =          1.079 for     57 degrees of freedom 
# Null hypothesis probability =   3.189229e-01
 
Photon flux (15-150 keV)    in 33.14 sec: 0.585436 ( -0.026887  0.0269378 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.36161e-06 ( -6.7189e-08  6.75642e-08 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.57294     ( -0.333402   0.185948 )
Epeak [keV] : 186.079      (   )
Norm@50keV  : 5.05546E-03    (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          61.11 using 59 PHA bins.
# Reduced chi-squared =          1.091 for     56 degrees of freedom 
# Null hypothesis probability =   2.974469e-01
 
Photon flux (15-150 keV)    in 33.14 sec: 0.581268 ( -0.029147  0.029081 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.34781e-06 ( -7.86313e-08  7.69643e-08 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.368773   ( 0.95908   0.95908 )
beta        : -1.74606    ( 1.73635   0.107431 )
Epeak [keV] : 35.7848      ( -35.8359 221.307 )
Norm@50keV   : 3.55197E-02    ( -0.0365645   366.274 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          59.46 using 59 PHA bins.
# Reduced chi-squared =          1.081 for     55 degrees of freedom 
# Null hypothesis probability =   3.164805e-01
 
Photon flux (15-150 keV)    in 33.14 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)  = 3.66e-07
S(50-100 keV) = 4.56e-07

T90 vs. Hardness ratio plot

T90                                   = 21.6960000097752 sec.
Hardness ratio (energy fluence ratio) = 1.2459

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.31942
Count Ratio (50-100 keV) / (15-25 keV) = 0.805355

Mask shadow pattern

IMX = 4.847090538993031E-02, IMY = -3.503441877668790E-01

TIME vs. PHA plot around the trigger time

TIME vs. DetID plot around the trigger time

100 us light curve (15-350 keV)

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:
-10.072000 -4.732000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
135.6411 -40.6780   11.6  37.2  -93.1  -0.0500  -1.0    9.0	Vela X-1                      
170.2044 -60.7621   82.1   8.1  -96.2   0.0026   0.1    9.0	Cen X-3                       
186.2514 -62.8738    0.0   0.3  -61.6   0.0293   1.7   12.7	GX 301-2                      
227.2846 -59.4519    1.8  19.5   82.3   0.1483   6.3 ------	UNKNOWN                       
208.1890 -11.5160    2.5  53.8   10.8   0.3861   4.6 ------	UNKNOWN

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

Time interval of the image:
-4.732000 28.408000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
135.7644 -40.6168    0.0  37.2  -92.9   0.2718   0.8   11.4	Vela X-1                      
170.7263 -60.5729   22.5   8.0  -94.2   0.0617   0.5   12.5	Cen X-3                       
186.2514 -62.8738    0.0   0.3  -61.6  -0.2432  -2.2   12.7	GX 301-2                      
227.2447 -59.4002    1.1  19.5   82.1   1.5747  10.3 ------	UNKNOWN                       
 82.1541 -68.6268    2.4  37.9 -160.5   0.8296   4.9 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
-10.072000 32.998000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
135.5428 -40.6387    6.7  37.3  -93.2  -0.3250  -1.7    5.1	Vela X-1                      
170.6544 -60.5710    5.0   8.0  -94.3   0.1479   2.3   10.5	Cen X-3                       
186.4304 -62.5853    0.0   0.5  -32.0   0.0222   0.3   12.7	GX 301-2                      
356.2718 -66.4253    2.1  50.3  159.1   1.1453   5.4 ------	UNKNOWN                       
227.2623 -59.4255    0.3  19.5   82.2   2.7438  35.3 ------	UNKNOWN                       
 50.9465 -52.8983    3.4  59.0 -166.1   1.3023   3.4 ------	UNKNOWN                       
106.8834 -57.6921    2.9  37.2 -135.1   0.3700   4.0 ------	UNKNOWN                       
129.4186 -39.2781    3.8  41.1  -98.0   0.8555   3.0 ------	UNKNOWN

Plot creation:
Sun May 31 16:39:50 EDT 2015