batgrbproduct Analysis for Trigger=318170

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Sun Jan 19 23:42:13 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger318170
======================================================================

  Trigger: 00318170      Segment: 000

  BAT Trigger Time: 238892858.368   [s; MET]
      Trigger Stop: 238892859.392   [s; MET]
               UTC: 2008-07-27T23:07:35.394100    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 32.6448768441636   Dec: 64.1515714708512  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 32.6448768441636   Dec: 64.1515714708512  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 32.639747754007   Dec: 64.1300985461938  [deg; J2000]  
        { 02h 10m 33.5s , +64d 07' 48.4" }
         +/- 0.769389733470983 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.292686512664784 [arcmin] (formal 1-sigma fit error)
    SNR: 44.2350373901105
 
  Angular difference between Flight and Ground position is 1.295349 arcmin
 
  Partial Coding Fraction: 0.7734375    [ including projection effects ]

  Duration
     T90: 77.363999992609 +/- 7.03936359051065
        Measured from: 1.68000000715256
                   to: 79.0439999997616     [s; relative to TRIGTIME]
     T50: 23.3200000226498 +/- 1.43458984090877
        Measured from: 14.9799999892712
                   to: 38.3000000119209     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: 15.2319999933243     
                   to: 16.2319999933243     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -3.30000001192093     
                   to: 110.247999995947     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     1.650692   2.958452   2.611177   0.606654
           0.057142   0.061597   0.056963   0.050429 [error]
 Peak      0.057938   0.103691   0.110149   0.022878
           0.006335   0.007377   0.007079   0.005600 [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 5.29e-07 1.15e-06 1.98e-06 7.85e-06 5.43e-06 1.15e-05   erg/cm2
              Error 2.81e-08 3.56e-08 5.83e-08 5.42e-07 1.48e-07 5.70e-07   erg/cm2
Peak flux Power-law 5.32e-01 7.01e-01 6.77e-01 1.17e+00 2.30e+00 3.08e+00   ph/cm2/s
              Error 7.43e-02 6.05e-02 5.50e-02 1.99e-01 1.69e-01 2.48e-01   ph/cm2/s
 
Fluence   Cutoff-PL 4.84e-07 1.18e-06 2.04e-06 4.94e-06 5.30e-06 8.64e-06   erg/cm2
              Error 3.58e-08 3.86e-08 6.75e-08 1.16e-06 1.63e-07 1.16e-06   erg/cm2
Peak flux Cutoff-PL 4.41e-01 7.39e-01 7.39e-01 5.60e-01 2.23e+00 2.48e+00   ph/cm2/s
              Error 9.16e-02 6.96e-02 6.89e-02 2.57e-01 1.75e-01 3.07e-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.221    0.047    -      -           4.74e-03 1.21e-04   68.620 57
CPL   0.864    0.193     178.9   86.7      6.85e-03 1.34e-03   57.750 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.050    0.121    -      -           1.99e-02 1.44e-03   77.440 57
CPL   0.251    0.550     121.3   70.9      4.64e-02 2.84e-02   69.910 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: 63.104000 sec.
T50 in the 50-300 keV band: 20.608000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100  113.548   -3.300  110.248
T90    77.364    1.680   79.044
T50    23.320   14.980   38.300
Peak    1.000   15.232   16.232
  • 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
       32.6397    64.1301   0.0044    20.2392   -64.6463         2.7711897     44.235 TRIG_00318170                 
     
    Foreground time interval of the image:
     -3.300  18.533 (delta_t =  21.833 [sec])
    Background time interval of the image:
    -239.368  -3.300 (delta_t = 236.068 [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   -3.300 sec. to  110.248 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.22108     ( -0.0473229   0.0471193 )
    Norm@50keV  : 4.74126E-03    ( -0.000120603   0.000120567 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          68.62 using 59 PHA bins.
    # Reduced chi-squared =          1.204 for     57 degrees of freedom 
    # Null hypothesis probability =   1.392519e-01
     
    Photon flux (15-150 keV)    in 113.5 sec: 0.558234 ( -0.0151255  0.0151522 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 5.42763e-06 ( -1.47903e-07  1.48171e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -3.300 sec. to  110.248 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.864335     ( -0.197145   0.188406 )
    Epeak [keV] : 178.906      ( -43.8743 129.495 )
    Norm@50keV  : 6.85075E-03    ( -0.00117953   0.00149286 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          57.75 using 59 PHA bins.
    # Reduced chi-squared =          1.031 for     56 degrees of freedom 
    # Null hypothesis probability =   4.104455e-01
     
    Photon flux (15-150 keV)    in 113.5 sec: 0.545664 ( -0.016358  0.016387 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 5.30465e-06 ( -1.6369e-07  1.62793e-07 ) ergs/cm2
     
    
    Band function
    Time interval is from   -3.300 sec. to  110.248 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 1.36234   ( -1.87319   -1.36309 )
    beta        : -1.31152    ( -0.0926621   0.0669797 )
    Epeak [keV] : 40.2881      ( -8.18006 49.2611 )
    Norm@50keV   : 0.233168    ( -0.234328   56.6652 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          52.82 using 59 PHA bins.
    # Reduced chi-squared =         0.9604 for     55 degrees of freedom 
    # Null hypothesis probability =   5.582491e-01
     
    Photon flux (15-150 keV)    in 113.5 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]    : 19.4562  (  )
    R^2/D10^2   : 2.73638E-02  (   )
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         287.10 using 59 PHA bins.
    # Reduced chi-squared =         5.0368 for     57 degrees of freedom 
    # Null hypothesis probability =   2.032451e-32
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.291e-07   1.266e-08
    25- 50  1.120e-06   4.495e-08
    50-150  3.130e-06   1.356e-07
    15-150  4.479e-06   1.300e-07
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 199.363       (      )
    Norm        : 2.36061       (     )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         205.50 using 59 PHA bins.
    # Reduced chi-squared =         3.6052 for     57 degrees of freedom 
    # Null hypothesis probability =   1.189941e-18
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  6.672e-07   1.856e-08
    25- 50  1.283e-06   3.689e-08
    50-150  2.815e-06   7.548e-08
    15-150  4.765e-06   1.379e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 8.51092        (-0.873063      0.959507)
    R1^2/D10^2  : 0.240004        (-0.0655194      0.0966244)
    kT2 [keV]   : 31.8162        (-3.06687      4.2261)
    R2^2/D10^2  : 4.49903E-03    (-0.00153031   0.00188145)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          50.92 using 59 PHA bins.
    # Reduced chi-squared =         0.9259 for     55 degrees of freedom 
    # Null hypothesis probability =   6.310303e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  4.501e-07   4.444e-08
    25- 50  1.215e-06   8.104e-08
    50-150  3.663e-06   3.047e-07
    15-150  5.329e-06   3.604e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from   15.232 sec. to   16.232 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.04968     ( -0.121956   0.120358 )
    Norm@50keV  : 1.98789E-02    ( -0.00144378   0.00144335 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          77.44 using 59 PHA bins.
    # Reduced chi-squared =          1.359 for     57 degrees of freedom 
    # Null hypothesis probability =   3.717234e-02
     
    Photon flux (15-150 keV)    in 1 sec: 2.29589 ( -0.169255  0.169425 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 2.11296e-07 ( -1.70347e-08  1.71941e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   15.232 sec. to   16.232 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.251374     ( -0.587278   0.513448 )
    Epeak [keV] : 121.347      ( -29.7206 112.134 )
    Norm@50keV  : 4.64149E-02    ( -0.0190896   0.0376419 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          69.91 using 59 PHA bins.
    # Reduced chi-squared =          1.248 for     56 degrees of freedom 
    # Null hypothesis probability =   1.001374e-01
     
    Photon flux (15-150 keV)    in 1 sec: 2.23139 ( -0.17492  0.17467 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 2.02601e-07 ( -1.8344e-08  1.8317e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   15.232 sec. to   16.232 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.251410   ( -0.457963   0.503122 )
    beta        : -10.0000    ( 10   20 )
    Epeak [keV] : 121.336      ( -22.0377 112.077 )
    Norm@50keV   : 4.64140E-02    ( -0.0190708   0.0376897 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          69.91 using 59 PHA bins.
    # Reduced chi-squared =          1.271 for     55 degrees of freedom 
    # Null hypothesis probability =   8.495413e-02
     
    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]    : 20.4455  (-1.79213  2.02688)
    R^2/D10^2   : 0.104890  (-0.0273815   0.035967)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          86.33 using 59 PHA bins.
    # Reduced chi-squared =          1.515 for     57 degrees of freedom 
    # Null hypothesis probability =   7.325413e-03
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  8.377e-09   1.255e-09
    25- 50  4.234e-08   4.360e-09
    50-150  1.316e-07   1.590e-08
    15-150  1.823e-07   1.765e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 199.363       (      )
    Norm        : 9.22161       (     )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         118.24 using 59 PHA bins.
    # Reduced chi-squared =         2.0744 for     57 degrees of freedom 
    # Null hypothesis probability =   3.546734e-06
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.296e-08   1.860e-09
    25- 50  4.417e-08   1.019e-08
    50-150  9.688e-08   3.358e-08
    15-150  1.640e-07   4.379e-08
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 7.40889        (-2.54818      4.16223)
    R1^2/D10^2  : 1.04939        (-0.726044      3.25699)
    kT2 [keV]   : 26.6522        (-4.28608      10.117)
    R2^2/D10^2  : 3.92365E-02    (-0.0270321   0.033807)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          69.03 using 59 PHA bins.
    # Reduced chi-squared =          1.255 for     55 degrees of freedom 
    # Null hypothesis probability =   9.672467e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.445e-08   6.661e-09
    25- 50  4.228e-08   1.651e-08
    50-150  1.455e-07   5.060e-08
    15-150  2.022e-07   6.736e-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.26889     ( -0.0597858   0.059601 )
    Norm@50keV  : 8.77444E-03    ( -0.000297209   0.000297208 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          70.15 using 59 PHA bins.
    # Reduced chi-squared =          1.231 for     57 degrees of freedom 
    # Null hypothesis probability =   1.133030e-01
     
    Photon flux (15-150 keV)    in 21.83 sec: 1.04105 ( -0.0367133  0.0367585 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.9076e-06 ( -7.10032e-08  7.11996e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 0.594978     ( -0.27164   0.255263 )
    Epeak [keV] : 111.904      ( -17.8179 35.4033 )
    Norm@50keV  : 1.81957E-02    ( -0.00427944   0.00600486 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.74 using 59 PHA bins.
    # Reduced chi-squared =         0.8524 for     56 degrees of freedom 
    # Null hypothesis probability =   7.760968e-01
     
    Photon flux (15-150 keV)    in 21.83 sec: 1.00449 ( -0.038957  0.03898 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.8218e-06 ( -7.89853e-08  7.88041e-08 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -0.594953   ( -0.179838   0.346606 )
    beta        : -10.0000    ( 10   5.32907e-15 )
    Epeak [keV] : 111.886      ( -10.6497 35.406 )
    Norm@50keV   : 1.81959E-02    ( -0.00428678   0.00600774 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          47.74 using 59 PHA bins.
    # Reduced chi-squared =         0.8679 for     55 degrees of freedom 
    # Null hypothesis probability =   7.457933e-01
     
    Photon flux (15-150 keV)    in 21.83 sec: 1.00826 (    ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.83343e-06 ( 0  0 ) ergs/cm2
     
    

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

    S(25-50 keV)  = 1.15e-06
    S(50-100 keV) = 1.98e-06
    

    T90 vs. Hardness ratio plot

    T90                                   = 77.363999992609 sec.
    Hardness ratio (energy fluence ratio) = 1.72174
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 1.79225
    Count Ratio (50-100 keV) / (15-25 keV) = 1.58187
    

    Mask shadow pattern

    IMX = 1.578814225197023E-01, IMY = 3.331922850841051E-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.368000 -3.300000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    299.5857  35.1879    0.7  40.7  135.9   8.7618  17.4    0.9	Cyg X-1                       
    308.1612  40.9372    4.2  32.0  133.4   0.8160   2.8    2.5	Cyg X-3                       
    326.0996  38.2689    3.9  26.0  106.9   0.8061   2.9    4.6	Cyg X-2                       
    285.6341  33.0528    2.3  49.6  146.9   4.3002   4.9 ------	UNKNOWN                       
    323.3445  31.6332    2.4  32.9  104.6   1.7250   4.9 ------	UNKNOWN                       
    106.4056  73.0572    4.3  41.1  -98.2   3.1772   2.7 ------	UNKNOWN                       
    

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

    Time interval of the image:
    -3.300000 18.532610
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    299.4654  35.0703    4.4  40.9  135.9   0.3782   2.6   10.0	Cyg X-1                       
    308.3157  40.8133   11.3  32.0  133.1  -0.0909  -1.0   12.5	Cyg X-3                       
    326.3772  38.4336    0.0  25.7  106.6  -0.1144  -1.5   11.8	Cyg X-2                       
     37.4393  41.9851    3.0  33.2  -22.3   0.2569   3.8 ------	UNKNOWN                       
     32.6397  64.1301    0.3  20.2  -64.6   2.7713  44.2 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    61.532600 362.774100
    420.632000 482.693600
    600.632000 662.724100
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    270.1649  59.0226    3.4  49.4  158.2   2.8548   3.4 ------	UNKNOWN                       
     32.7203  64.1684    1.9   0.0 -105.4   1.5542   6.0 ------	UNKNOWN                       
     39.2630  44.4543    1.9  20.0   -0.8   1.3837   6.0 ------	UNKNOWN                       
     74.7641  19.5421    3.1  52.7  -39.5   4.3553   3.7 ------	UNKNOWN                       
    
    Plot creation:
    Sun May 31 06:48:32 EDT 2015