batgrbproduct Analysis for Trigger=683906

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Wed Jun  1 15:37:22 EDT 2016
             Host Name: bat1
   Current Working Dir: /local/data/bat1/alien/Swift_3rdBATcatalog/event/batevent_reproc/trigger683906
======================================================================

  Trigger: 00683906      Segment: 000

  BAT Trigger Time: 483191363.52   [s; MET]
      Trigger Stop: 483191364.544   [s; MET]
               UTC: 2016-04-24T11:49:07.348300    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 319.480565124725   Dec: -60.4052836164155  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 319.480565124725   Dec: -60.4052836164155  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 319.48462934057   Dec: -60.4081088330894  [deg; J2000]  
        { 21h 17m 56.3s , -60d 24' 29.2" }
         +/- 1.13082706229561 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.419566901776941 [arcmin] (formal 1-sigma fit error)
    SNR: 25.5174988547371
 
  Angular difference between Flight and Ground position is 0.207934 arcmin
 
  Partial Coding Fraction: 0.3359375    [ including projection effects ]

  Duration
     T90: 6.4640000462532 +/- 0.951697413545641
        Measured from: -0.879999995231628
                   to: 5.58400005102158     [s; relative to TRIGTIME]
     T50: 1.98400002717972 +/- 0.354175118295394
        Measured from: 0.0960000157356262
                   to: 2.08000004291534     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: 0.0680000185966492     
                   to: 1.06800001859665     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -1.42399996519089     
                   to: 6.73600000143051     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.392838   0.687849   0.568343   0.183438
           0.034348   0.037512   0.034184   0.027809 [error]
 Peak      0.116218   0.200956   0.232535   0.085746
           0.013395   0.015500   0.015920   0.012445 [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.50e-07 2.93e-07 4.44e-07 1.46e-06 1.25e-06 2.35e-06   erg/cm2
              Error 1.68e-08 1.98e-08 3.23e-08 2.33e-07 8.06e-08 2.54e-07   erg/cm2
Peak flux Power-law 1.28e+00 1.59e+00 1.44e+00 2.27e+00 5.09e+00 6.58e+00   ph/cm2/s
              Error 1.93e-01 1.47e-01 1.20e-01 4.01e-01 4.01e-01 5.21e-01   ph/cm2/s
 
Fluence   Cutoff-PL 1.41e-07 3.00e-07 4.57e-07 9.79e-07 1.22e-06 1.88e-06   erg/cm2
              Error 2.04e-08 2.29e-08 3.73e-08 4.76e-07 8.50e-08 4.80e-07   erg/cm2
Peak flux Cutoff-PL 6.17e-01 4.81e-01 2.63e-01 2.51e-01 1.46e+00 1.68e+00   ph/cm2/s
              Error 2.86e-01 1.95e-01 7.38e-02 2.61e-02 5.63e-01 5.18e-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.400    0.104    -      -           1.58e-02 9.53e-04   50.790 57
CPL   1.105    0.391     163.5  -47.9      2.18e-02 9.80e-03   49.220 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.143    0.130    -      -           4.37e-02 3.31e-03   62.750 57
CPL   0.879    0.439     248.2    0.0      5.72e-02 2.85e-02   61.940 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: 3.072000 sec.
T50 in the 50-300 keV band: 1.024000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100    8.160   -1.424    6.736
T90     6.464   -0.880    5.584
T50     1.984    0.096    2.080
Peak    1.000    0.068    1.068
  • Click here for the BAT GCN circular template
  • Special notes of this burst
  • Estimated redshift from machine learning (Ukwatta et al. 2016) - may take a few minute to calcaulte

    Image

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

    #   RAcent    DECcent   POSerr    Theta     Phi              Peak Cts       SNR   Name
      319.4846   -60.4081   0.0076    44.3598   172.6921         1.8508401     25.517 TRIG_00683906                 
     
    Foreground time interval of the image:
     -1.424   6.736 (delta_t =   8.160 [sec])
    Background time interval of the image:
    -239.520  -1.424 (delta_t = 238.096 [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)

    Spectral Evolution

    Spectra

    Notes:

    Time averaged spectrum fit using the average DRM

    Power-law model
    Time interval is from   -1.424 sec. to    6.736 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.40043     ( -0.104163   0.103983 )
    Norm@50keV  : 1.57901E-02    ( -0.000952744   0.00095257 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          50.79 using 59 PHA bins.
    # Reduced chi-squared =         0.8911 for     57 degrees of freedom 
    # Null hypothesis probability =   7.057564e-01
     
    Photon flux (15-150 keV)    in 8.16 sec: 1.92316 ( -0.12274  0.12301 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.24506e-06 ( -8.04519e-08  8.06845e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -1.424 sec. to    6.736 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.10461     ( -0.439562   0.341897 )
    Epeak [keV] : 163.538      ( -67.6566 -163.538 )
    Norm@50keV  : 2.18320E-02    ( -0.0065832   0.0130122 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.22 using 59 PHA bins.
    # Reduced chi-squared =         0.8789 for     56 degrees of freedom 
    # Null hypothesis probability =   7.273753e-01
     
    Photon flux (15-150 keV)    in 8.16 sec: 1.89858 ( -0.12706  0.12726 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.22693e-06 ( -8.53128e-08  8.47571e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   -1.424 sec. to    6.736 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.193739   ( -1.0817   0.466058 )
    beta        : -1.53048    ( 1.52973   0.170736 )
    Epeak [keV] : 54.7027      ( -28.6934 226.496 )
    Norm@50keV   : 8.35509E-02    ( -0.087373   531.067 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.84 using 59 PHA bins.
    # Reduced chi-squared =         0.8698 for     55 degrees of freedom 
    # Null hypothesis probability =   7.422521e-01
     
    Photon flux (15-150 keV)    in 8.16 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.9525  (-1.29618  1.43032)
    R^2/D10^2   : 0.146476  (-0.036019   0.0469037)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         105.94 using 59 PHA bins.
    # Reduced chi-squared =         1.8586 for     57 degrees of freedom 
    # Null hypothesis probability =   8.907396e-05
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  7.001e-08   8.654e-09
    25- 50  3.091e-07   2.538e-08
    50-150  6.336e-07   6.654e-08
    15-150  1.013e-06   7.385e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 200.000       (-19.3956      -200)
    Norm        : 8.02955       (-0.485924     0.485924)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          56.42 using 59 PHA bins.
    # Reduced chi-squared =         0.9898 for     57 degrees of freedom 
    # Null hypothesis probability =   4.968706e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.632e-07   1.183e-08
    25- 50  3.138e-07   2.301e-08
    50-150  6.883e-07   1.565e-07
    15-150  1.165e-06   1.506e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 8.72445        (-1.77385      1.99535)
    R1^2/D10^2  : 0.935258        (-0.421544      0.955291)
    kT2 [keV]   : 33.2451        (-7.63739      18.0674)
    R2^2/D10^2  : 1.12273E-02    (-0.00818063   0.0170195)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.22 using 59 PHA bins.
    # Reduced chi-squared =         0.8948 for     55 degrees of freedom 
    # Null hypothesis probability =   6.942592e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.260e-07   3.481e-08
    25- 50  3.202e-07   6.230e-08
    50-150  7.873e-07   3.461e-07
    15-150  1.234e-06   3.753e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from    0.068 sec. to    1.068 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.14321     ( -0.130175   0.129078 )
    Norm@50keV  : 4.36651E-02    ( -0.00331107   0.0033054 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          62.75 using 59 PHA bins.
    # Reduced chi-squared =          1.101 for     57 degrees of freedom 
    # Null hypothesis probability =   2.797050e-01
     
    Photon flux (15-150 keV)    in 1 sec: 5.08819 ( -0.40038  0.40112 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 4.50474e-07 ( -3.65792e-08  3.67291e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from    0.068 sec. to    1.068 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.878850     ( -0.551403   0.327243 )
    Epeak [keV] : 248.160      ( -248.16 -248.16 )
    Norm@50keV  : 5.71722E-02    ( -0.0160036   0.0409074 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          61.94 using 59 PHA bins.
    # Reduced chi-squared =          1.106 for     56 degrees of freedom 
    # Null hypothesis probability =   2.724556e-01
     
    Photon flux (15-150 keV)    in 1 sec: 1.46068 ( -0.700342  0.42559 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 4.47208e-07 ( -3.75891e-08  3.74551e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from    0.068 sec. to    1.068 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.890464   ( 0.733515   0.733515 )
    beta        : -4.09615    ( 8.79332   18.7531 )
    Epeak [keV] : 259.635      ( -124.293 -259.581 )
    Norm@50keV   : 5.64053E-02    ( -0.00658276   0.0140523 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          61.95 using 59 PHA bins.
    # Reduced chi-squared =          1.126 for     55 degrees of freedom 
    # Null hypothesis probability =   2.421543e-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]    : 21.1635  (-1.99673  2.2779)
    R^2/D10^2   : 0.201905  (-0.057095   0.0768741)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          92.64 using 59 PHA bins.
    # Reduced chi-squared =          1.625 for     57 degrees of freedom 
    # Null hypothesis probability =   1.984153e-03
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.703e-08   2.600e-09
    25- 50  8.795e-08   9.365e-09
    50-150  2.938e-07   3.800e-08
    15-150  3.987e-07   3.610e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 200.000       (-14.161      -200)
    Norm        : 21.1753       (-1.64372     1.64372)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          90.84 using 59 PHA bins.
    # Reduced chi-squared =          1.594 for     57 degrees of freedom 
    # Null hypothesis probability =   2.913115e-03
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  5.274e-08   4.975e-09
    25- 50  1.015e-07   2.094e-08
    50-150  2.225e-07   7.138e-08
    15-150  3.767e-07   9.220e-08
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 6.78074        (-1.89489      2.70301)
    R1^2/D10^2  : 4.22597        (-2.79853      10.1947)
    kT2 [keV]   : 29.4370        (-4.76117      8.6989)
    R2^2/D10^2  : 6.09540E-02    (-0.0352503   0.0504421)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          61.54 using 59 PHA bins.
    # Reduced chi-squared =          1.119 for     55 degrees of freedom 
    # Null hypothesis probability =   2.535340e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  3.852e-08   1.795e-08
    25- 50  8.997e-08   2.435e-08
    50-150  3.211e-07   8.060e-08
    15-150  4.496e-07   1.098e-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.40043     ( -0.104163   0.103983 )
    Norm@50keV  : 1.57901E-02    ( -0.000952744   0.00095257 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          50.79 using 59 PHA bins.
    # Reduced chi-squared =         0.8911 for     57 degrees of freedom 
    # Null hypothesis probability =   7.057564e-01
     
    Photon flux (15-150 keV)    in 8.16 sec: 1.92312 ( -0.122726  0.123001 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.2449e-06 ( -8.04928e-08  8.07029e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.10461     ( -0.439562   0.341897 )
    Epeak [keV] : 163.538      ( -67.6566 -163.538 )
    Norm@50keV  : 2.18320E-02    ( -0.0065832   0.0130122 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.22 using 59 PHA bins.
    # Reduced chi-squared =         0.8789 for     56 degrees of freedom 
    # Null hypothesis probability =   7.273753e-01
     
    Photon flux (15-150 keV)    in 8.16 sec: 1.89839 ( -0.12703  0.12721 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.22741e-06 ( -8.54434e-08  8.4815e-08 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.193739   ( -1.0817   0.466058 )
    beta        : -1.53048    ( 1.52973   0.170736 )
    Epeak [keV] : 54.7027      ( -28.6981 218.391 )
    Norm@50keV   : 8.35509E-02    ( -0.087373   536.621 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.84 using 59 PHA bins.
    # Reduced chi-squared =         0.8698 for     55 degrees of freedom 
    # Null hypothesis probability =   7.422521e-01
     
    Photon flux (15-150 keV)    in 8.16 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.93e-07
    S(50-100 keV) = 4.44e-07
    

    T90 vs. Hardness ratio plot

    T90                                   = 6.4640000462532 sec.
    Hardness ratio (energy fluence ratio) = 1.51536
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 1.75097
    Count Ratio (50-100 keV) / (15-25 keV) = 1.44676
    

    Mask shadow pattern

    IMX = -9.699563817779852E-01, IMY = -1.243905168047467E-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.520000 -1.424000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    135.6687 -40.6234    6.0  46.4   46.1   1.6008   1.9    7.6	Vela X-1                      
    104.6806 -52.9156    2.4  23.1   33.3   1.1585   4.8 ------	UNKNOWN                       
     99.8563 -41.7114    2.3  29.6   11.1   1.2329   5.1 ------	UNKNOWN                       
    320.3496 -56.2678    3.0  47.3  176.8   1.8980   3.8 ------	UNKNOWN                       
     93.7130 -10.5784    3.0  56.2  -14.9   3.0955   3.8 ------	UNKNOWN                       
      5.0720 -32.0317    3.3  48.0 -129.3   3.8927   3.5 ------	UNKNOWN                       
     54.2818 -19.5420    2.8  43.8  -64.2   4.5285   4.0 ------	UNKNOWN                       
    

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

    Time interval of the image:
    -1.424000 6.736000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    135.3633 -40.4235   13.2  46.4   45.7  -0.1302  -0.9   10.9	Vela X-1                      
    319.4846 -60.4081    0.5  44.4  172.7   1.8509  25.5 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    709.680610 876.780600
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     83.6379  22.0151    0.3   6.7 -114.7   6.9352  44.2    0.3	Crab                          
     78.3784  77.4458    2.6  53.5  -10.1   2.1659   4.5 ------	UNKNOWN                       
     50.3976  79.1184    2.3  57.7  -15.1   3.7129   5.0 ------	UNKNOWN                       
     30.2752  70.9636    3.4  57.6  -26.4   3.3169   3.4 ------	UNKNOWN                       
    
    Plot creation:
    Wed Jun  1 17:46:36 EDT 2016