batgrbproduct Analysis for Trigger=533268

Analysis Summary

Notes:

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Mon Jan 13 01:01:29 EST 2014
             Host Name: batevent1
   Current Working Dir: /local/data/batevent/batevent_reproc/trigger533268
======================================================================

  Trigger: 00533268      Segment: 000

  BAT Trigger Time: 369040123.392   [s; MET]
      Trigger Stop: 369040124.416   [s; MET]
               UTC: 2012-09-11T07:08:33.754960    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 357.960763192049   Dec: 63.0987077305092  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 357.960763192049   Dec: 63.0987077305092  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 357.966505198214   Dec: 63.0901179230837  [deg; J2000]  
        { 23h 51m 52.0s , +63d 05' 24.4" }
         +/- 1.5861721238515 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.894022123582005 [arcmin] (formal 1-sigma fit error)
    SNR: 15.7363369772024
 
  Angular difference between Flight and Ground position is 0.538452 arcmin
 
  Partial Coding Fraction: 0.2363281    [ including projection effects ]

  Duration
     T90: 21.4679999351501 +/- 3.84755763774157
        Measured from: -1.12800002098083
                   to: 20.3399999141693     [s; relative to TRIGTIME]
     T50: 5.77999991178513 +/- 7.84400228456224
        Measured from: 0.495999991893768
                   to: 6.27599990367889     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.300000011920929     
                   to: 0.699999988079071     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -1.6240000128746     
                   to: 24.0999999642372     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.536009   0.793181   0.370379   0.205111
           0.056045   0.057683   0.050343   0.042397 [error]
 Peak      0.079126   0.160280   0.076247   0.022147
           0.012863   0.014351   0.011685   0.008417 [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.80e-07 2.90e-07 3.52e-07 8.42e-07 1.06e-06 1.66e-06   erg/cm2
              Error 2.58e-08 2.71e-08 4.77e-08 2.32e-07 1.13e-07 2.77e-07   erg/cm2
Peak flux Power-law 9.66e-01 9.01e-01 5.83e-01 5.82e-01 2.69e+00 3.03e+00   ph/cm2/s
              Error 1.77e-01 1.12e-01 8.74e-02 1.65e-01 3.27e-01 3.63e-01   ph/cm2/s
 
Fluence   Cutoff-PL 1.68e-07 3.14e-07 3.52e-07 3.00e-07 9.97e-07 1.13e-06   erg/cm2
              Error 2.92e-08 3.91e-08 5.48e-08 4.09e-07 1.36e-07 5.58e-07   erg/cm2
Peak flux Cutoff-PL 7.80e-01 1.12e+00 6.16e-01 9.88e-02 2.60e+00 2.61e+00   ph/cm2/s
              Error 2.13e-01 1.76e-01 1.19e-01 1.19e-01 3.37e-01 3.49e-01   ph/cm2/s
 
Fluence        Band 1.56e-07 3.33e-07 3.31e-07 5.69e-07 1.01e-06 1.39e-06   erg/cm2
              Error 3.22e-08 4.12e-08 5.33e-08 -------- 1.19e-07 --------   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.717    0.161    -      -           4.47e-03 4.54e-04   57.870 57
CPL   1.068    0.628      63.4  -23.5      1.03e-02 1.44e-02   55.550 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.628    0.182    -      -           2.07e-02 2.57e-03   75.460 57
CPL  -0.248    1.347      53.3   13.3      2.27e-01 8.35e-01   66.140 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: 4.160000 sec.
T50 in the 50-300 keV band: 1.984000 sec.
 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100   25.724   -1.624   24.100
T90    21.468   -1.128   20.340
T50     5.780    0.496    6.276
Peak    1.000   -0.300    0.700
  • 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
      357.9665    63.0901   0.0122    34.0713    99.8536         1.9677411     15.736 TRIG_00533268                 
     
    Foreground time interval of the image:
     -1.624  24.100 (delta_t =  25.724 [sec])
    Background time interval of the image:
    -85.691  -1.624 (delta_t =  84.067 [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.624 sec. to   24.100 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.71683     ( -0.159191   0.162762 )
    Norm@50keV  : 4.46748E-03    ( -0.000457013   0.000451684 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          57.87 using 59 PHA bins.
    # Reduced chi-squared =          1.015 for     57 degrees of freedom 
    # Null hypothesis probability =   4.428182e-01
     
    Photon flux (15-150 keV)    in 25.72 sec: 0.596863 ( -0.0562073  0.0562531 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.06253e-06 ( -1.128e-07  1.14074e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -1.624 sec. to   24.100 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.06815     ( -0.902637   0.353193 )
    Epeak [keV] : 63.4344      ( -16.5125 -63.4669 )
    Norm@50keV  : 1.02805E-02    ( -0.00596633   0.0228307 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.55 using 59 PHA bins.
    # Reduced chi-squared =         0.9919 for     56 degrees of freedom 
    # Null hypothesis probability =   4.919071e-01
     
    Photon flux (15-150 keV)    in 25.72 sec: 0.583013 ( -0.058742  0.058484 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 9.96853e-07 ( -1.36033e-07  1.36203e-07 ) ergs/cm2
     
    
    Band function
    Time interval is from   -1.624 sec. to   24.100 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 1.29086   ( -1.43214   -0.235649 )
    beta        : -2.05254    ( -0.628429   0.349726 )
    Epeak [keV] : 34.9427      ( -14.3484 24.6812 )
    Norm@50keV   : 0.420323    ( -0.438186   92.2594 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          51.42 using 59 PHA bins.
    # Reduced chi-squared =         0.9349 for     55 degrees of freedom 
    # Null hypothesis probability =   6.121561e-01
     
    Photon flux (15-150 keV)    in 25.72 sec:  (    ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.01101e-06 ( -1.01083e-06  -1.19482e-07 ) ergs/cm2
     
    
    Single BB
     
    Spectral model blackbody:
    ------------------------------------------------------------
    Parameters  : value    Lower 90%  Upper 90%
    kT [keV]    : 11.4099  (-1.16913  1.34174)
    R^2/D10^2   : 0.177208  (-0.0603283   0.0882814)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          70.55 using 59 PHA bins.
    # Reduced chi-squared =          1.238 for     57 degrees of freedom 
    # Null hypothesis probability =   1.072156e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.248e-07   2.175e-08
    25- 50  3.695e-07   4.900e-08
    50-150  2.792e-07   6.871e-08
    15-150  7.735e-07   9.619e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 117.627       (-37.0564      67.0077)
    Norm        : 2.55640       (-0.294873     0.371872)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.64 using 59 PHA bins.
    # Reduced chi-squared =         0.9762 for     57 degrees of freedom 
    # Null hypothesis probability =   5.261006e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.714e-07   2.294e-08
    25- 50  3.056e-07   4.249e-08
    50-150  5.399e-07   2.284e-07
    15-150  1.017e-06   2.692e-07
    
    Double BB
     
     
    Spectral model: bbodyrad<1> + bbodyrad<2>
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 8.93548        (-1.73778      1.54256)
    R1^2/D10^2  : 0.371893        (-0.126755      0.11057)
    kT2 [keV]   : 47.8033        (-50.2741      -77.3541)
    R2^2/D10^2  : 7.21057E-04    (-0.000198357   0.00425331)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          49.48 using 59 PHA bins.
    # Reduced chi-squared =         0.8997 for     55 degrees of freedom 
    # Null hypothesis probability =   6.845546e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  1.529e-07   3.922e-08
    25- 50  3.456e-07   8.575e-08
    50-150  5.164e-07   2.557e-07
    15-150  1.015e-06   3.404e-07
    

    Peak spectrum fit

    Power-law model
    Time interval is from   -0.300 sec. to    0.700 sec.
     
    Spectral model in power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.62838     ( -0.180427   0.182805 )
    Norm@50keV  : 2.07442E-02    ( -0.00257668   0.00256138 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          75.46 using 59 PHA bins.
    # Reduced chi-squared =          1.324 for     57 degrees of freedom 
    # Null hypothesis probability =   5.136634e-02
     
    Photon flux (15-150 keV)    in 1 sec: 2.68968 ( -0.32728  0.327667 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.9337e-07 ( -2.51313e-08  2.53066e-08 ) ergs/cm2
     
    
    Cutoff power-law model
    Time interval is from   -0.300 sec. to    0.700 sec.
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: -0.247503     ( -1.55931   1.1357 )
    Epeak [keV] : 53.3143      ( -8.21262 18.4031 )
    Norm@50keV  : 0.227234    ( -0.228083   1.44226 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          66.14 using 59 PHA bins.
    # Reduced chi-squared =          1.181 for     56 degrees of freedom 
    # Null hypothesis probability =   1.666039e-01
     
    Photon flux (15-150 keV)    in 1 sec: 2.60111 ( -0.33818  0.33585 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.71392e-07 ( -2.6845e-08  2.8415e-08 ) ergs/cm2
     
    
    Band function
    Time interval is from   -0.300 sec. to    0.700 sec.
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 1.22552   ( -1.53245   -1.24859 )
    beta        : -2.52508    ( -1.30982   0.56923 )
    Epeak [keV] : 45.4729      ( -10.2737 12.3067 )
    Norm@50keV   : 0.940679    ( -1.0528   166.661 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          62.28 using 59 PHA bins.
    # Reduced chi-squared =          1.132 for     55 degrees of freedom 
    # Null hypothesis probability =   2.330193e-01
     
    Photon flux (15-150 keV)    in 1 sec:  (    ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.81052e-07 ( -1.81045e-07  -2.6314e-08 ) ergs/cm2
     
    
    Single BB
     
    Spectral model blackbody:
    ------------------------------------------------------------
    Parameters  : value    Lower 90%  Upper 90%
    kT [keV]    : 12.5121  (-1.36361  1.56398)
    R^2/D10^2   : 0.628589  (-0.225166   0.334453)
    (R is the radius in km and D10 is the distance to the source in units of 10 kpc)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          68.37 using 59 PHA bins.
    # Reduced chi-squared =          1.199 for     57 degrees of freedom 
    # Null hypothesis probability =   1.439936e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.086e-08   4.480e-09
    25- 50  6.889e-08   1.027e-08
    50-150  6.663e-08   1.746e-08
    15-150  1.564e-07   2.290e-08
    
    Thermal bremsstrahlung (OTTB)
     
     
    Spectral model: thermal bremsstrahlung
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT [keV]    : 143.131       (-52.106      -139.793)
    Norm        : 11.3572       (-1.52916     1.90872)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          71.44 using 59 PHA bins.
    # Reduced chi-squared =          1.253 for     57 degrees of freedom 
    # Null hypothesis probability =   9.438272e-02
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.928e-08   4.975e-09
    25- 50  5.390e-08   1.726e-08
    50-150  1.043e-07   5.737e-08
    15-150  1.875e-07   7.526e-08
    
    Double BB
    ------------------------------------------------------------
    Parameters  : value         Lower 90%    Upper 90%
    kT1 [keV]   : 10.9845        (-1.36387      1.69715)
    R1^2/D10^2  : 0.925076        (-0.358424      0.56304)
    kT2 [keV]   : 200.000        (      )
    R2^2/D10^2  : 1.91527E-04    (-9.94884e-05   0.00122635)
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          59.19 using 59 PHA bins.
    #Fit statistic : Chi-Squared =          58.75 using 59 PHA bins.
    # Reduced chi-squared =          1.076 for     55 degrees of freedom 
    # Reduced chi-squared =          1.068 for     55 degrees of freedom 
    # Null hypothesis probability =   3.253602e-01
    # Null hypothesis probability =   3.397884e-01
     
    Energy  Fluence     90% Error
    [keV]   [erg/cm2]   [erg/cm2]
    15- 25  2.362e-08   6.550e-09
    25- 50  6.770e-08   1.849e-08
    50-150  9.657e-08   4.557e-08
    15-150  1.879e-07   7.190e-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.71683     ( -0.159191   0.162762 )
    Norm@50keV  : 4.46748E-03    ( -0.000457013   0.000451684 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          57.87 using 59 PHA bins.
    # Reduced chi-squared =          1.015 for     57 degrees of freedom 
    # Null hypothesis probability =   4.428182e-01
     
    Photon flux (15-150 keV)    in 25.72 sec: 0.596863 ( -0.0562073  0.0562531 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.06253e-06 ( -1.128e-07  1.14074e-07 ) ergs/cm2
     
    
    Cutoff power-law model
    Spectral model in the cutoff power-law:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    Photon index: 1.06815     ( -0.902637   0.353193 )
    Epeak [keV] : 63.4344      ( -16.5125 -63.4669 )
    Norm@50keV  : 1.02805E-02    ( -0.00596633   0.0228307 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          55.55 using 59 PHA bins.
    # Reduced chi-squared =         0.9919 for     56 degrees of freedom 
    # Null hypothesis probability =   4.919071e-01
     
    Photon flux (15-150 keV)    in 25.72 sec: 0.583013 ( -0.058742  0.058484 ) ph/cm2/s
    Energy fluence (15-150 keV)             : 9.96853e-07 ( -1.36033e-07  1.36203e-07 ) ergs/cm2
     
    
    Band function
     
    Spectral model in the Band function:
    ------------------------------------------------------------
    Parameters  : value       lower 90%   higher 90%
    alpha       : 1.29086   ( -1.43214   -0.235649 )
    beta        : -2.05254    ( -0.628429   0.349726 )
    Epeak [keV] : 34.9427      ( -14.3484 24.6812 )
    Norm@50keV   : 0.420323    ( -0.438186   92.2594 )
    ------------------------------------------------------------
    #Fit statistic : Chi-Squared =          51.42 using 59 PHA bins.
    # Reduced chi-squared =         0.9349 for     55 degrees of freedom 
    # Null hypothesis probability =   6.121561e-01
     
    Photon flux (15-150 keV)    in 25.72 sec:  (    ) ph/cm2/s
    Energy fluence (15-150 keV)             : 1.01101e-06 ( -1.01083e-06  -1.19482e-07 ) ergs/cm2
     
    

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

    S(25-50 keV)  = 2.90e-07
    S(50-100 keV) = 3.52e-07
    

    T90 vs. Hardness ratio plot

    T90                                   = 21.4679999351501 sec.
    Hardness ratio (energy fluence ratio) = 1.21379
    

    Color vs. Color plot

    Count Ratio  (25-50 keV) / (15-25 keV) = 1.47979
    Count Ratio (50-100 keV) / (15-25 keV) = 0.690994
    

    Mask shadow pattern

    IMX = -1.157392562803334E-01, IMY = -6.663447067642172E-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:
    -85.691400 -1.624000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
     31.1866  68.6732    2.2  24.1   77.9   0.6053   5.2 ------	UNKNOWN                       
    268.7684  44.5958    3.1  54.2  171.1   1.5834   3.7 ------	UNKNOWN                       
    260.0376  50.4404    2.4  47.6  177.5   0.8944   4.7 ------	UNKNOWN                       
    

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

    Time interval of the image:
    -1.624000 24.100000
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    357.9665  63.0901    0.7  34.1   99.9   1.9670  15.7 ------	UNKNOWN                       
    

    Post-slew image of 15-350 keV band

    Time interval of the image:
    556.108590 962.654600
     
    # RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
    #  [deg]    [deg]    ['] [deg]  [deg]                   [']
    135.6009 -40.6368    5.6  55.7 -169.5  -2.1923  -2.1    5.9	Vela X-1                      
    170.2169 -60.5401    6.4  38.1  158.2   1.0415   1.8    5.8	Cen X-3                       
    186.6509 -62.7958    1.8  34.9  145.7   3.4711   6.5    1.5	GX 301-2                      
    244.9831 -15.6433    0.4  39.9   30.0  20.4882  32.7    0.3	Sco X-1                       
    212.6614 -33.6105    2.9   8.0   71.7   1.0907   3.9 ------	UNKNOWN                       
    167.7587 -33.1577    2.9  31.2 -154.0   1.4303   4.0 ------	UNKNOWN                       
    134.7336 -31.6246    4.8  58.8 -159.3   6.2995   2.4 ------	UNKNOWN                       
    
    Plot creation:
    Sat May 30 18:31:49 EDT 2015