batgrbproduct Analysis for Trigger=586680

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Tue Mar  4 17:50:10 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger586680
======================================================================

  Trigger: 00586680      Segment: 000

  BAT Trigger Time: 414130041.92   [s; MET]
      Trigger Stop: 414130042.944   [s; MET]
               UTC: 2014-02-15T04:07:10.273340    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 104.135942727015   Dec: 41.7803604777544  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 104.135942727015   Dec: 41.7803604777544  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 104.131035415408   Dec: 41.7870063941203  [deg; J2000]  
        { 06h 56m 31.4s , +41d 47' 13.2" }
         +/- 1.11504409092991 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.474871573306095 [arcmin] (formal 1-sigma fit error)
    SNR: 26.0350443650305
 
  Angular difference between Flight and Ground position is 0.455202 arcmin
 
  Partial Coding Fraction: 0.1542969    [ including projection effects ]

  Duration
     T90: 25.6320000290871 +/- 5.64181532462767
        Measured from: -0.304000020027161
                   to: 25.3280000090599     [s; relative to TRIGTIME]
     T50: 8.78400003910065 +/- 0.840875777678445
        Measured from: 1.74399995803833
                   to: 10.527999997139     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.0120000243186951     
                   to: 0.987999975681305     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -1.47200000286102     
                   to: 35.8560000061989     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     1.647878   3.072008   3.280514   1.500788
           0.108582   0.122947   0.118974   0.099929 [error]
 Peak      0.140164   0.302490   0.404192   0.142714
           0.021060   0.027269   0.030847   0.023297 [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 6.08e-07 1.43e-06 2.67e-06 1.20e-05 7.24e-06 1.67e-05   erg/cm2
              Error 5.37e-08 7.46e-08 1.15e-07 1.22e-06 3.00e-07 1.26e-06   erg/cm2
Peak flux Power-law 1.64e+00 2.33e+00 2.44e+00 4.74e+00 7.90e+00 1.12e+01   ph/cm2/s
              Error 2.96e-01 2.61e-01 2.28e-01 9.26e-01 7.17e-01 1.09e+00   ph/cm2/s
 
Fluence   Cutoff-PL 6.10e-07 1.43e-06 2.67e-06 1.19e-05 7.24e-06 1.66e-05   erg/cm2
              Error 5.49e-08 7.42e-08 1.15e-07 1.49e-06 2.98e-07 1.63e-06   erg/cm2
Peak flux Cutoff-PL 1.50e+00 2.36e+00 2.56e+00 3.58e+00 7.85e+00 1.00e+01   ph/cm2/s
              Error 3.72e-01 2.78e-01 3.07e-01 1.53e+00 7.21e-01 1.76e+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.099    0.073    -      -           1.86e-02 7.64e-04   53.370 57
CPL   1.096    0.103    9999.4    0.0      1.86e-02 2.05e-03   53.430 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    0.928    0.146    -      -           6.88e-02 6.20e-03   56.570 57
CPL   0.612    0.497     268.5    0.0      9.39e-02 5.06e-02   55.600 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: 13.632000 sec.
T50 in the 50-300 keV band: 6.592000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100   37.328   -1.472   35.856
T90    25.632   -0.304   25.328
T50     8.784    1.744   10.528
Peak    1.000   -0.012    0.988
  • 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
      104.1310    41.7870   0.0074    46.6078   155.6688         8.1692408     26.035 TRIG_00586680                 
     
    Foreground time interval of the image:
     -1.472  14.981 (delta_t =  16.453 [sec])
    Background time interval of the image:
    -239.920  -1.472 (delta_t = 238.448 [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   -1.472 sec. to   35.856 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.09874     ( -0.0733605   0.0724535 )
    Norm@50keV  : 1.85515E-02    ( -0.000765325   0.000762752 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          53.37 using 59 PHA bins.
    # Reduced chi-squared =         0.9363 for     57 degrees of freedom 
    # Null hypothesis probability =   6.120264e-01
     
    Photon flux (15-150 keV)    in 37.33 sec: 2.15158 ( -0.0937347  0.0940828 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.24348e-06 ( -2.99514e-07  2.99971e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -1.472 sec. to   35.856 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.09636     ( -0.134624   0.070617 )
    Epeak [keV] : 9999.36      ( -9999.36 -9999.36 )
    Norm@50keV  : 1.86464E-02    ( -0.000764677   0.00334072 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          53.43 using 59 PHA bins.
    # Reduced chi-squared =         0.9541 for     56 degrees of freedom 
    # Null hypothesis probability =   5.727744e-01
     
    Photon flux (15-150 keV)    in 37.33 sec: 2.15079 (    ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.24198e-06 ( -2.94944e-07  3.00768e-07 ) ergs/cm2
     
    
    Band function
    Time interval is from   -1.472 sec. to   35.856 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 4.99996   ( -4.99996   -4.99996 )
    beta        : -1.11625    ( -0.0758599   0.0760599 )
    Epeak [keV] : 23.9489      ( -23.9489 4.38679 )
    Norm@50keV   : 1743.91    ( -69.0446   3.50235e+06 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          52.23 using 59 PHA bins.
    # Reduced chi-squared =         0.9496 for     55 degrees of freedom 
    # Null hypothesis probability =   5.812187e-01
     
    Photon flux (15-150 keV)    in 37.33 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.9566  (  )
    R^2/D10^2   : 6.43175E-02  (   )
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         202.25 using 59 PHA bins.
    # Reduced chi-squared =         3.5482 for     57 degrees of freedom 
    # Null hypothesis probability =   3.916469e-18
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.295e-07   2.087e-08
    25- 50  1.244e-06   8.101e-08
    50-150  4.897e-06   2.800e-07
    15-150  6.371e-06   2.930e-07
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 199.363       (      )
    Norm        : 9.24132       (     )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         177.65 using 59 PHA bins.
    # Reduced chi-squared =         3.1167 for     57 degrees of freedom 
    # Null hypothesis probability =   2.553258e-14
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  8.590e-07   3.864e-08
    25- 50  1.652e-06   7.149e-08
    50-150  3.624e-06   1.689e-07
    15-150  6.136e-06   2.613e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 8.86303        (-1.25433      1.41761)
    R1^2/D10^2  : 0.866533        (-0.319441      0.559923)
    kT2 [keV]   : 41.9955        (-6.41499      10.7303)
    R2^2/D10^2  : 8.57112E-03    (-0.0040229   0.00560561)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          52.37 using 59 PHA bins.
    # Reduced chi-squared =         0.9522 for     55 degrees of freedom 
    # Null hypothesis probability =   5.756338e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  5.610e-07   7.597e-08
    25- 50  1.504e-06   1.820e-07
    50-150  5.298e-06   6.757e-07
    15-150  7.363e-06   9.295e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from   -0.012 sec. to    0.988 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.927833     ( -0.147941   0.144951 )
    Norm@50keV  : 6.87777E-02    ( -0.00620998   0.00618281 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          56.57 using 59 PHA bins.
    # Reduced chi-squared =         0.9924 for     57 degrees of freedom 
    # Null hypothesis probability =   4.912277e-01
     
    Photon flux (15-150 keV)    in 1 sec: 7.89735 ( -0.716438  0.717196 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.63668e-07 ( -7.22602e-08  7.24054e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -0.012 sec. to    0.988 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.612403     ( -0.609691   0.384718 )
    Epeak [keV] : 268.456      ( -268.455 -268.455 )
    Norm@50keV  : 9.38684E-02    ( -0.0295047   0.0717171 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.60 using 59 PHA bins.
    # Reduced chi-squared =         0.9928 for     56 degrees of freedom 
    # Null hypothesis probability =   4.901357e-01
     
    Photon flux (15-150 keV)    in 1 sec: 7.85076 ( -0.72054  0.72164 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.62158e-07 ( -7.3434e-08  7.3278e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   -0.012 sec. to    0.988 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 5.00000   ( -5   -5 )
    beta        : -1.00330    ( -0.16968   0.165143 )
    Epeak [keV] : 29.7016      ( -29.7016 6.72912 )
    Norm@50keV   : 1652.02    ( -144.708   9.35018e+06 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          54.48 using 59 PHA bins.
    # Reduced chi-squared =         0.9905 for     55 degrees of freedom 
    # Null hypothesis probability =   4.944746e-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]    : 24.9472  (-2.79143  3.27864)
    R^2/D10^2   : 0.202327  (-0.0635896   0.0889605)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          75.53 using 59 PHA bins.
    # Reduced chi-squared =          1.325 for     57 degrees of freedom 
    # Null hypothesis probability =   5.079151e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.190e-08   3.795e-09
    25- 50  1.241e-07   1.815e-08
    50-150  5.729e-07   7.575e-08
    15-150  7.189e-07   8.445e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 200.000       (      )
    Norm        : 31.2215       (-2.94924     2.94917)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         103.75 using 59 PHA bins.
    # Reduced chi-squared =         1.8202 for     57 degrees of freedom 
    # Null hypothesis probability =   1.529928e-04
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  7.772e-08   2.210e-08
    25- 50  1.495e-07   7.747e-08
    50-150  3.279e-07   1.761e-07
    15-150  5.551e-07   2.036e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 7.30873        (-1.96528      2.8171)
    R1^2/D10^2  : 4.43137        (-2.87962      8.81203)
    kT2 [keV]   : 34.1354        (-5.93538      10.639)
    R2^2/D10^2  : 6.89921E-02    (-0.0386772   0.0724488)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          52.44 using 59 PHA bins.
    # Reduced chi-squared =         0.9535 for     55 degrees of freedom 
    # Null hypothesis probability =   5.729114e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  5.126e-08   2.433e-08
    25- 50  1.340e-07   5.100e-08
    50-150  5.920e-07   1.907e-07
    15-150  7.773e-07   2.622e-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: 0.978906     ( -0.0687717   0.067874 )
    Norm@50keV  : 3.50744E-02    ( -0.00140291   0.00139649 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          64.66 using 59 PHA bins.
    # Reduced chi-squared =          1.134 for     57 degrees of freedom 
    # Null hypothesis probability =   2.269292e-01
     
    Photon flux (15-150 keV)    in 16.45 sec: 4.03415 ( -0.165256  0.165615 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 6.28632e-06 ( -3.69383e-07  2.46225e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.972492     ( -0.265432   0.0754755 )
    Epeak [keV] : 9999.99      (   )
    Norm@50keV  : 3.52897E-02    ( -0.00140628   0.00132466 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          64.65 using 59 PHA bins.
    # Reduced chi-squared =          1.154 for     56 degrees of freedom 
    # Null hypothesis probability =   2.002191e-01
     
    Photon flux (15-150 keV)    in 16.45 sec: 4.03225 ( -0.1645  0.16593 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 6.29757e-06 ( -2.36633e-07  2.38821e-07 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 5.00000   ( -4.99996   -4.99996 )
    beta        : -0.987836    ( -0.0698851   0.0703277 )
    Epeak [keV] : 23.1876      ( -23.1876 5.82681 )
    Norm@50keV   : 3561.42    ( -133.672   5.25589e+06 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          64.30 using 59 PHA bins.
    # Reduced chi-squared =          1.169 for     55 degrees of freedom 
    # Null hypothesis probability =   1.831021e-01
     
    Photon flux (15-150 keV)    in 16.45 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.43e-06
    S(50-100 keV) = 2.67e-06
    

    T90 vs. Hardness ratio plot

    T90                                   = 25.6320000290871 sec.
    Hardness ratio (energy fluence ratio) = 1.86713
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 1.86422
    Count Ratio (50-100 keV) / (15-25 keV) = 1.99075
    

    Mask shadow pattern

    IMX = -9.638086568856981E-01, IMY = -4.358076949019931E-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.920000 -1.472000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    109.7571  30.5394    2.8  46.5  138.9   4.2201   4.1 ------	UNKNOWN                       
    161.1466  15.7843    3.4  22.8   53.1   0.8800   3.4 ------	UNKNOWN                       
    194.8764  -7.4623    2.3  52.9    3.9   3.3668   5.0 ------	UNKNOWN                       
    223.8811  84.9261    3.0  49.2 -133.7   4.2546   3.9 ------	UNKNOWN                       
    223.5093  71.1299    2.5  43.8 -116.2   5.1476   4.7 ------	UNKNOWN                       
    227.0921  79.6737    2.7  47.4 -127.1   4.5142   4.2 ------	UNKNOWN                       
    

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

    Time interval of the image:
    -1.472000 14.980600
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    182.0660 -13.4622    3.2  53.8   21.3   1.2487   3.6 ------	UNKNOWN                       
    104.1310  41.7870    0.4  46.6  155.7   8.1688  26.0 ------	UNKNOWN                       
    198.9244  31.1300    2.7  27.7  -44.6   0.3223   4.3 ------	UNKNOWN                       
    222.3708  52.2314    2.4  40.9  -88.3   1.0883   4.8 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    105.380610 962.150000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     83.6251  22.0158    0.2  26.2 -156.1  21.0297  46.5    0.5	Crab                          
    136.1658  12.7169    3.9  40.2  103.4   5.0623   3.0 ------	UNKNOWN                       
     95.5681  -2.9443    2.7  45.4  168.6   3.9109   4.2 ------	UNKNOWN                       
     43.1048  55.0778    2.6  41.1  -73.1   7.2845   4.4 ------	UNKNOWN                       
     51.8114  31.1664    3.4  42.7 -110.3   5.8578   3.4 ------	UNKNOWN                       
    
    Plot creation:
    Mon Jun  8 22:33:03 EDT 2015