batgrbproduct Analysis for Trigger=509337

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Mon May 25 18:19:02 EDT 2015
             Host Name: bat1
   Current Working Dir: /local/data/bat1/alien/Swift_3rdBATcatalog/event/batevent_reproc/trigger509337
======================================================================

  Trigger: 00509337      Segment: 000

  BAT Trigger Time: 345107960   [s; MET]
      Trigger Stop: 345108024   [s; MET]
               UTC: 2011-12-09T07:19:12.342960    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 14.362573622164   Dec: -46.8016976759424  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 14.362573622164   Dec: -46.8016976759424  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 14.3499054195286   Dec: -46.7991139590406  [deg; J2000]  
        { 00h 57m 24.0s , -46d 47' 56.8" }
         +/- 0.479293136934408 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.116581972297158 [arcmin] (formal 1-sigma fit error)
    SNR: 86.9765877008079
 
  Angular difference between Flight and Ground position is 0.542918 arcmin
 
  Partial Coding Fraction: 1.015625    [ including projection effects ]

  Duration
     T90: 810.96799993515 +/- 52.1094840475323
        Measured from: 23.1440001130104
                   to: 834.112000048161     [s; relative to TRIGTIME]
     T50: 243.275999963284 +/- 16.8582879558111
        Measured from: 172.344000101089
                   to: 415.620000064373     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -64.3840000033379     
                   to: -63.3840000033379     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -66.9520000219345     
                   to: 914.212000012398     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total    17.292129  21.892395  17.291401   4.119940
           0.178111   0.186194   0.167361   0.145194 [error]
 Peak      0.020490   0.060550   0.031330   0.014640
           0.013633   0.013278   0.012589   0.011378 [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 4.85e-06 8.98e-06 1.28e-05 3.84e-05 3.65e-05 6.51e-05   erg/cm2
              Error 1.67e-07 1.76e-07 2.11e-07 1.55e-06 5.40e-07 1.59e-06   erg/cm2
Peak flux Power-law 3.99e-01 3.62e-01 2.26e-01 2.16e-01 1.08e+00 1.20e+00   ph/cm2/s
              Error 2.21e-01 1.28e-01 1.01e-01 1.97e-01 3.73e-01 4.04e-01   ph/cm2/s
 
Fluence   Cutoff-PL 4.76e-06 8.99e-06 1.29e-05 3.58e-05 3.64e-05 6.24e-05   erg/cm2
              Error 2.27e-07 1.80e-07 2.45e-07 2.31e-06 6.00e-07 3.92e-06   erg/cm2
Peak flux Cutoff-PL 3.98e-01 3.62e-01 2.26e-01 2.14e-01 1.08e+00 1.20e+00   ph/cm2/s
              Error 2.32e-01 1.25e-01 1.07e-01 2.32e-01 3.74e-01 4.40e-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.489    0.029    -      -           3.92e-03 5.92e-05   36.410 57
CPL   1.488    0.000    9720.2    0.0      3.93e-03 3.60e-04   36.360 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.678    0.617    -      -           8.18e-03 2.91e-03   74.580 57
CPL   1.658    2.347     814.6    0.0      8.39e-03 0.00e+00   74.580 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: 935.680000 sec.
T50 in the 50-300 keV band: 257.728000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100  981.164  -66.952  914.212
T90   810.968   23.144  834.112
T50   243.276  172.344  415.620
Peak    1.000  -64.384  -63.384
  • 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
       14.3499   -46.7991   0.0022     0.0366  -106.8900        57.8054139     86.977 TRIG_00509337                 
     
    Foreground time interval of the image:
    -16.299 164.560 (delta_t = 180.859 [sec])
    166.000 914.212 (delta_t = 748.212 [sec])
    Background time interval of the image:
    

    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 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:

    Time averaged spectrum fit using the average DRM

    Power-law model
    Multiple GTIs in the spectrum
    tstart tstop
    -66.952000 164.560000
    166.000000 914.212000
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.48925     ( -0.029508   0.0294588 )
    Norm@50keV  : 3.91935E-03    ( -5.92695e-05   5.91967e-05 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          36.41 using 59 PHA bins.
    # Reduced chi-squared =         0.6387 for     57 degrees of freedom 
    # Null hypothesis probability =   9.847144e-01
     
    Photon flux (15-150 keV)    in 979.7 sec: 0.487888 ( -0.00900603  0.00904641 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.6491e-05 ( -5.39986e-07  5.40067e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Multiple GTIs in the spectrum
    tstart tstop
    -66.952000 164.560000
    166.000000 914.212000
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.48780     (     )
    Epeak [keV] : 9720.16      ( -9646.41 -9646.41 )
    Norm@50keV  : 3.93203E-03    ( -5.98322e-05   0.00065929 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          36.36 using 59 PHA bins.
    # Reduced chi-squared =         0.6493 for     56 degrees of freedom 
    # Null hypothesis probability =   9.806091e-01
     
    Photon flux (15-150 keV)    in 979.7 sec: 0.484700 ( -0.009765  0.011484 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.63791e-05 ( -5.63425e-07  6.36903e-07 ) ergs/cm2
     
    
    Band function
    Multiple GTIs in the spectrum
    tstart tstop
    -66.952000 164.560000
    166.000000 914.212000
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -1.45042   ( 1.11455   5.94631 )
    beta        : -1.53905    ( 1.59954   0.137876 )
    Epeak [keV] : 855.055      ( -863.918 -863.918 )
    Norm@50keV   : 4.06610E-03    ( -6.03856e-05   0.000228922 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          36.06 using 59 PHA bins.
    # Reduced chi-squared =         0.6557 for     55 degrees of freedom 
    # Null hypothesis probability =   9.773661e-01
     
    Photon flux (15-150 keV)    in 979.7 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]    : 18.9910  (  )
    R^2/D10^2   : 2.30909E-02  (   )
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         943.90 using 59 PHA bins.
    # Reduced chi-squared =         16.560 for     57 degrees of freedom 
    # Null hypothesis probability =  6.843204e-161
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.604e-06   5.192e-08
    25- 50  7.712e-06   1.793e-07
    50-150  2.044e-05   4.507e-07
    15-150  2.976e-05   5.143e-07
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 200.000       (      )
    Norm        : 2.02430       (-0.0298736     0.0298742)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =         102.55 using 59 PHA bins.
    # Reduced chi-squared =         1.7991 for     57 degrees of freedom 
    # Null hypothesis probability =   2.048368e-04
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  4.940e-06   9.356e-08
    25- 50  9.503e-06   1.705e-07
    50-150  2.084e-05   3.527e-07
    15-150  3.528e-05   6.221e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 7.68779        (-0.492792      0.521975)
    R1^2/D10^2  : 0.345523        (-0.0716344      0.0952025)
    kT2 [keV]   : 29.7429        (-1.5877      1.88171)
    R2^2/D10^2  : 4.06472E-03    (-0.00079706   0.000893325)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          67.73 using 59 PHA bins.
    # Reduced chi-squared =          1.231 for     55 degrees of freedom 
    # Null hypothesis probability =   1.163409e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  4.061e-06   2.371e-07
    25- 50  9.126e-06   2.479e-07
    50-150  2.267e-05   6.417e-07
    15-150  3.585e-05   8.964e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from  -64.384 sec. to  -63.384 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.67788     ( -0.594956   0.638697 )
    Norm@50keV  : 8.18168E-03    ( -0.00297239   0.0028546 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          74.58 using 59 PHA bins.
    # Reduced chi-squared =          1.308 for     57 degrees of freedom 
    # Null hypothesis probability =   5.901308e-02
     
    Photon flux (15-150 keV)    in 1 sec: 1.07835 ( -0.371952  0.374742 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.58931e-08 ( -2.77618e-08  2.84288e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from  -64.384 sec. to  -63.384 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.65759     ( -3.75932   0.935636 )
    Epeak [keV] : 814.644      (   )
    Norm@50keV  : 8.38679E-03    (     )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          74.58 using 59 PHA bins.
    # Reduced chi-squared =          1.332 for     56 degrees of freedom 
    # Null hypothesis probability =   4.913258e-02
     
    Photon flux (15-150 keV)    in 1 sec: 1.07826 ( -0.375659  0.372 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 7.55397e-08 ( -4.45477e-08  2.87013e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from  -64.384 sec. to  -63.384 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 4.95003   ( -4.95003   -4.95003 )
    beta        : -2.03382    ( -2.20114   0.81367 )
    Epeak [keV] : 30.8726      ( -30.8726 55.7524 )
    Norm@50keV   : 263.753    ( -87.3087   1.11098e+07 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          74.02 using 59 PHA bins.
    # Reduced chi-squared =          1.346 for     55 degrees of freedom 
    # Null hypothesis probability =   4.450198e-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]    : 10.4495  (-3.49155  8.62948)
    R^2/D10^2   : 0.429613  (-0.426938   1.6008)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          77.14 using 59 PHA bins.
    # Reduced chi-squared =          1.353 for     57 degrees of freedom 
    # Null hypothesis probability =   3.905761e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  9.683e-09   6.035e-09
    25- 50  2.559e-08   1.588e-08
    50-150  1.515e-08   1.112e-08
    15-150  5.042e-08   3.275e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 136.521       (      )
    Norm        : 4.52823       (-1.63942     2.86443)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          74.68 using 59 PHA bins.
    # Reduced chi-squared =          1.310 for     57 degrees of freedom 
    # Null hypothesis probability =   5.807016e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.101e-08   6.780e-09
    25- 50  2.083e-08   1.252e-08
    50-150  4.358e-08   2.754e-08
    15-150  7.542e-08   4.504e-08
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 7.84028        (-3.33105      5.00179)
    R1^2/D10^2  : 1.05128        (-0.664271      2.10334)
    kT2 [keV]   : 43.1662        (-28.8713      -43.5794)
    R2^2/D10^2  : 2.27147E-03    (-0.00220121   0.0852505)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          73.30 using 59 PHA bins.
    # Reduced chi-squared =          1.333 for     55 degrees of freedom 
    # Null hypothesis probability =   5.008550e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.202e-08   7.420e-09
    25- 50  2.304e-08   1.358e-08
    50-150  4.244e-08   2.597e-08
    15-150  7.751e-08   4.026e-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.49585     ( -0.0296012   0.0295545 )
    Norm@50keV  : 3.89799E-03    ( -5.90774e-05   5.90092e-05 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          36.55 using 59 PHA bins.
    # Reduced chi-squared =         0.6412 for     57 degrees of freedom 
    # Null hypothesis probability =   9.839909e-01
     
    Photon flux (15-150 keV)    in 929.1 sec: 0.486081 ( -0.00899577  0.00903561 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.43799e-05 ( -5.10913e-07  5.10989e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.43193     ( -0.125035   0.0737666 )
    Epeak [keV] : 513.350      ( -321.871 -549.842 )
    Norm@50keV  : 4.14210E-03    ( -0.000293424   0.000494677 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          35.91 using 59 PHA bins.
    # Reduced chi-squared =         0.6413 for     56 degrees of freedom 
    # Null hypothesis probability =   9.831629e-01
     
    Photon flux (15-150 keV)    in 929.1 sec: 0.482915 ( -0.010647  0.010655 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 3.42699e-05 ( -5.648e-07  5.0924e-07 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : -1.43430   ( 1.10157   5.99891 )
    beta        : -10.0000    ( 10   1.77636e-15 )
    Epeak [keV] : 552.675      ( -177.816 -552.746 )
    Norm@50keV   : 4.13262E-03    ( -0.00022499   0.000496974 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          35.77 using 59 PHA bins.
    # Reduced chi-squared =         0.6504 for     55 degrees of freedom 
    # Null hypothesis probability =   9.792906e-01
     
    Photon flux (15-150 keV)    in 929.1 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)  = 8.98e-06
    S(50-100 keV) = 1.28e-05
    

    T90 vs. Hardness ratio plot

    T90                                   = 810.96799993515 sec.
    Hardness ratio (energy fluence ratio) = 1.42539
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 1.26603
    Count Ratio (50-100 keV) / (15-25 keV) = 0.999958
    

    Mask shadow pattern

    IMX = -1.855401688589516E-04, IMY = 6.110664473265062E-04
    

    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:
     
    

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

    Time interval of the image:
    -16.299410 164.560000
    166.000000 914.212000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     59.1157 -41.4584    2.4  32.1  133.3   5.3548   4.8 ------	UNKNOWN                       
    114.4612 -56.7569    2.5  57.1  176.2  19.4703   4.6 ------	UNKNOWN                       
    115.3314 -59.5268    2.9  55.8  179.2  14.3795   3.9 ------	UNKNOWN                       
     55.8419 -63.9388    2.3  28.4  178.4   3.7535   5.0 ------	UNKNOWN                       
     14.3499 -46.7991    0.1   0.0 -106.9  57.8057  87.0 ------	UNKNOWN                       
    125.2951 -65.1611    2.8  56.0 -172.0  11.0640   4.1 ------	UNKNOWN                       
    329.6172  -3.0686    2.7  58.4  -19.5  16.1338   4.2 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    -239.000000 164.560000
    166.000000 963.057900
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     59.1179 -41.4696    2.4  32.1  133.3   6.2729   4.9 ------	UNKNOWN                       
    349.0415   3.8965    3.6  55.4    4.9  11.8919   3.2 ------	UNKNOWN                       
    115.3235 -59.5208    3.0  55.8  179.2  16.6512   3.9 ------	UNKNOWN                       
     55.8509 -63.9284    2.3  28.4  178.4   4.5011   5.1 ------	UNKNOWN                       
     14.3506 -46.8004    0.2   0.0 -108.8  59.6262  75.8 ------	UNKNOWN                       
    125.3151 -65.1603    2.9  56.1 -172.0  12.5274   4.0 ------	UNKNOWN                       
    329.6091  -3.0697    2.5  58.4  -19.5  19.9374   4.6 ------	UNKNOWN                       
    324.3848 -56.5684    3.5  31.7  -90.4   4.6904   3.3 ------	UNKNOWN                       
    
    Plot creation:
    Sat May 30 20:57:25 EDT 2015