batgrbproduct analysis for trigger=312958

batgrbproduct Analysis for Trigger=312958

Analysis Summary

Notes:

1) No attempt has been made to eliminate the non-significant decimal places.
2) This batgrbproduct analysis is done using the HEASARC data.
3) This batgrbproduct analysis is using the flight position for the mask-weighting calculation.

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

  Trigger: 00312958      Segment: 000

  BAT Trigger Time: 234063030.912   [s; MET]
      Trigger Stop: 234063032.96   [s; MET]
               UTC: 2008-06-02T01:30:28.077080    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 19.1726774847547   Dec: -9.25077306570324  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 19.1726774847547   Dec: -9.25077306570324  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 19.1738896892608   Dec: -9.23794765140966  [deg; J2000]  
        { 01h 16m 41.7s , -09d 14' 16.6" }
         +/- 1.67606817425502 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.826878743819133 [arcmin] (formal 1-sigma fit error)
    SNR: 14.544597081511
 
  Angular difference between Flight and Ground position is 0.772866 arcmin
 
  Partial Coding Fraction: 0.1972656    [ including projection effects ]

  Duration
     T90: 74.2880000174046 +/- 7.14318726944463
        Measured from: -9.75600001215935
                   to: 64.5320000052452     [s; relative to TRIGTIME]
     T50: 58.9760000109673 +/- 1.19843230514097
        Measured from: 1.20800000429153
                   to: 60.1840000152588     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: 59.5959999859333     
                   to: 60.5959999859333     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -11.4880000054836     
                   to: 69.2439999878407     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.923140   1.512023   1.557406   0.445529
           0.103868   0.110691   0.105159   0.101678 [error]
 Peak      0.032527   0.110120   0.153191   0.061038
           0.007375   0.008660   0.009172   0.008094 [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 3.40e-07 7.16e-07 1.18e-06 4.41e-06 3.26e-06 6.65e-06   erg/cm2
              Error 4.77e-08 6.00e-08 9.99e-08 8.39e-07 2.53e-07 8.95e-07   erg/cm2
Peak flux Power-law 4.55e-01 7.39e-01 9.10e-01 2.21e+00 2.73e+00 4.32e+00   ph/cm2/s
              Error 7.80e-02 7.78e-02 6.83e-02 3.64e-01 2.09e-01 3.82e-01   ph/cm2/s
 
Fluence   Cutoff-PL 3.07e-07 7.38e-07 1.23e-06 2.57e-06 3.19e-06 4.85e-06   erg/cm2
              Error 5.96e-08 6.72e-08 1.18e-07 -------- 2.68e-07 1.62e-06   erg/cm2
Peak flux Cutoff-PL 3.84e-01 7.52e-01 9.76e-01 1.33e+00 2.68e+00 3.44e+00   ph/cm2/s
              Error 9.77e-02 8.43e-02 9.18e-02 6.44e-01 2.13e-01 6.66e-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.277    0.126    -      -           4.06e-03 2.92e-04   58.360 57
CPL   0.844    0.516     149.5  -47.4      6.39e-03 3.63e-03   56.120 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    0.700    0.132    -      -           2.36e-02 1.89e-03   72.470 57
CPL   0.135    0.567     201.3  695.2      4.07e-02 2.30e-02   69.100 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: 74.048000 sec.
T50 in the 50-300 keV band: 59.456000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100   80.732  -11.488   69.244
T90    74.288   -9.756   64.532
T50    58.976    1.208   60.184
Peak    1.000   59.596   60.596

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
   19.1739    -9.2379   0.0133    51.2871    -0.8423         2.7714487     14.545 TRIG_00312958                 
 
Foreground time interval of the image:
-11.488  29.589 (delta_t =  41.077 [sec])
Background time interval of the image:
-280.411 -11.488 (delta_t = 268.923 [sec])

Lightcurves

Notes:

1) The mask-weighted light curves are using the flight position.
2) Multiple plots of different time binning/intervals are shown to cover all scenarios of short/long GRBs, rate/image triggers, and real/false positives.
3) For all the mask-weighted lightcurves, the y-axis units are counts/sec/det where a det is 0.4 x 0.4 = 0.16 sq cm.
4) The verticle lines correspond to following: green dotted lines are T50 interval, black dotted lines are T90 interval, blue solid line(s) is a spacecraft slew start time, and orange solid line(s) is a spacecraft slew end time.
5) Time of each bin is in the middle of the bin.

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 2.048000

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 2.048000

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 2.048000

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:

1) The fitting includes the systematic errors.
2) When the burst includes telescope slew time periods, the fitting uses an average response file made from multiple 5-s response files through out the slew time plus single time preiod for the non-slew times, and weighted by the total counts in the corresponding time period using addrmf. An average response file is needed becuase a spectral fit using the pre-slew DRM will introduce some errors in both a spectral shape and a normalization if the PHA file contains a lot of the slew/post-slew time interval.
3) For fits to more complicated models (e.g. a power-law over a cutoff power-law), the BAT team has decided to require a chi-square improvement of more than 6 for each extra dof.

Time averaged spectrum fit using the average DRM

Power-law model

Time interval is from  -11.488 sec. to   69.244 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.27732     ( -0.126987   0.126004 )
Norm@50keV  : 4.05925E-03    ( -0.000292525   0.000292455 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          58.36 using 59 PHA bins.
# Reduced chi-squared =          1.024 for     57 degrees of freedom 
# Null hypothesis probability =   4.252667e-01
 
Photon flux (15-150 keV)    in 80.73 sec: 0.482316 ( -0.0366921  0.0367995 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.2566e-06 ( -2.52184e-07  2.53089e-07 ) ergs/cm2

Cutoff power-law model

Time interval is from  -11.488 sec. to   69.244 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.844447     ( -0.556998   0.475937 )
Epeak [keV] : 149.464      ( -54.7415 -149.464 )
Norm@50keV  : 6.39394E-03    ( -0.00243873   0.0048207 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.12 using 59 PHA bins.
# Reduced chi-squared =          1.002 for     56 degrees of freedom 
# Null hypothesis probability =   4.705032e-01
 
Photon flux (15-150 keV)    in 80.73 sec: 0.471987 ( -0.038518  0.038559 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.18925e-06 ( -2.68589e-07  2.6699e-07 ) ergs/cm2

Band function

Time interval is from  -11.488 sec. to   69.244 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.844408   ( -0.181372   0.741835 )
beta        : -10.0000    ( 10   0 )
Epeak [keV] : 149.435      ( -74.117 468.917 )
Norm@50keV   : 6.39379E-03    ( -0.00244574   0.00482225 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.12 using 59 PHA bins.
# Reduced chi-squared =          1.020 for     55 degrees of freedom 
# Null hypothesis probability =   4.328073e-01
 
Photon flux (15-150 keV)    in 80.73 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.9420  (-1.68868  1.90769)
R^2/D10^2   : 2.60631E-02  (-0.00720414   0.00969036)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          86.38 using 59 PHA bins.
# Reduced chi-squared =          1.515 for     57 degrees of freedom 
# Null hypothesis probability =   7.258410e-03
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.486e-07   2.309e-08
25- 50  7.130e-07   7.944e-08
50-150  1.879e-06   2.365e-07
15-150  2.740e-06   2.414e-07

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 200.000       (      )
Norm        : 2.02849       (-0.147889     0.147863)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          70.99 using 59 PHA bins.
# Reduced chi-squared =          1.245 for     57 degrees of freedom 
# Null hypothesis probability =   1.007329e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.079e-07   3.605e-08
25- 50  7.847e-07   1.022e-07
50-150  1.721e-06   4.622e-07
15-150  2.913e-06   7.149e-07

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 7.53529        (-2.63053      3.18719)
R1^2/D10^2  : 0.298461        (-0.190907      0.606799)
kT2 [keV]   : 27.9037        (-5.55417      12.3666)
R2^2/D10^2  : 5.99195E-03    (-0.0042951   0.00752225)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          57.02 using 59 PHA bins.
# Reduced chi-squared =          1.037 for     55 degrees of freedom 
# Null hypothesis probability =   3.997539e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.973e-07   1.268e-07
25- 50  7.389e-07   2.520e-07
50-150  2.154e-06   8.708e-07
15-150  3.190e-06   1.018e-06

Peak spectrum fit

Power-law model

Time interval is from   59.596 sec. to   60.596 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.700165     ( -0.134205   0.12927 )
Norm@50keV  : 2.36289E-02    ( -0.00189881   0.00187873 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          72.47 using 59 PHA bins.
# Reduced chi-squared =          1.271 for     57 degrees of freedom 
# Null hypothesis probability =   8.126769e-02
 
Photon flux (15-150 keV)    in 1 sec: 2.73125 ( -0.208438  0.208855 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.88395e-07 ( -2.25847e-08  2.26604e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from   59.596 sec. to   60.596 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.135262     ( -0.608663   0.525537 )
Epeak [keV] : 201.336      ( -73.5645 1316.9 )
Norm@50keV  : 4.06863E-02    ( -0.0157157   0.0302213 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          69.10 using 59 PHA bins.
# Reduced chi-squared =          1.234 for     56 degrees of freedom 
# Null hypothesis probability =   1.122168e-01
 
Photon flux (15-150 keV)    in 1 sec: 2.68391 ( -0.21292  0.21325 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.82946e-07 ( -2.4025e-08  2.3744e-08 ) ergs/cm2

Band function

Time interval is from   59.596 sec. to   60.596 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 5.00000   ( -5   -5 )
beta        : -0.702475    ( -0.129268   0.134491 )
Epeak [keV] : 18.5969      ( -17.7129 -17.7129 )
Norm@50keV   : 6419.68    ( -469.886   5.23971e+06 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          72.40 using 59 PHA bins.
# Reduced chi-squared =          1.316 for     55 degrees of freedom 
# Null hypothesis probability =   5.789469e-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]    : 26.7240  (-2.73097  3.26539)
R^2/D10^2   : 5.88075E-02  (-0.0167428   0.0221471)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          82.37 using 59 PHA bins.
# Reduced chi-squared =          1.445 for     57 degrees of freedom 
# Null hypothesis probability =   1.561824e-02
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  7.035e-09   1.064e-09
25- 50  4.124e-08   4.730e-09
50-150  2.157e-07   2.550e-08
15-150  2.639e-07   2.400e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 199.363       (      )
Norm        : 10.3997       (     )
------------------------------------------------------------
#Fit statistic : Chi-Squared =         183.39 using 59 PHA bins.
# Reduced chi-squared =         3.2173 for     57 degrees of freedom 
# Null hypothesis probability =   3.435843e-15
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.590e-08   3.065e-09
25- 50  4.981e-08   9.515e-09
50-150  1.093e-07   5.064e-08
15-150  1.850e-07   8.363e-08

Double BB

------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 13.6563        (-4.30066      5.60508)
R1^2/D10^2  : 0.236891        (-0.129909      0.361616)
kT2 [keV]   : 56.0957        (-42.7543      -55.2335)
R2^2/D10^2  : 5.84154E-03    (   )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          66.38 using 59 PHA bins.
#Fit statistic : Chi-Squared =          66.36 using 59 PHA bins.
# Reduced chi-squared =          1.207 for     55 degrees of freedom 
# Reduced chi-squared =          1.206 for     55 degrees of freedom 
# Null hypothesis probability =   1.399104e-01
# Null hypothesis probability =   1.403624e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.129e-08   5.823e-09
25- 50  4.737e-08   2.238e-08
50-150  2.258e-07   1.137e-07
15-150  2.845e-07   1.420e-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.35635     ( -0.17661   0.17626 )
Norm@50keV  : 4.70761E-03    ( -0.000479037   0.000478993 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          59.55 using 59 PHA bins.
# Reduced chi-squared =          1.045 for     57 degrees of freedom 
# Null hypothesis probability =   3.830944e-01
 
Photon flux (15-150 keV)    in 41.08 sec: 0.567880 ( -0.0616567  0.0619854 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.8861e-06 ( -2.03517e-07  2.04299e-07 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.655265     ( -0.911496   0.746678 )
Epeak [keV] : 110.058      ( -33.2416 -110.056 )
Norm@50keV  : 9.63616E-03    ( -0.0038871   0.0137868 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          57.04 using 59 PHA bins.
# Reduced chi-squared =          1.019 for     56 degrees of freedom 
# Null hypothesis probability =   4.362605e-01
 
Photon flux (15-150 keV)    in 41.08 sec: 0.546719 ( -0.06578  0.065933 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.83645e-06 ( -2.14865e-07  2.13743e-07 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.656492   ( -0.632316   0.991501 )
beta        : -10.0000    ( 10   1.77636e-15 )
Epeak [keV] : 110.110      ( -33.3183 170.559 )
Norm@50keV   : 9.62568E-03    ( -0.00503731   0.0138243 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          57.04 using 59 PHA bins.
# Reduced chi-squared =          1.037 for     55 degrees of freedom 
# Null hypothesis probability =   3.991944e-01
 
Photon flux (15-150 keV)    in 41.08 sec: 0.556750 (    ) 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)  = 7.16e-07
S(50-100 keV) = 1.18e-06

T90 vs. Hardness ratio plot

T90                                   = 74.2880000174046 sec.
Hardness ratio (energy fluence ratio) = 1.64804

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.63791
Count Ratio (50-100 keV) / (15-25 keV) = 1.68707

Mask shadow pattern

IMX = 1.247494626952796E+00, IMY = 1.833962613966093E-02

TIME vs. PHA plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 2.048000

TIME vs. DetID plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 2.048000

100 us light curve (15-350 keV)

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 2.048000

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:
-280.411400 -11.488000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
326.3740  38.2981    5.8  34.3  108.7   0.8660   2.0    9.6	Cyg X-2                       
 12.4578   8.0656    2.2  33.5    7.1   1.2156   5.2 ------	UNKNOWN                       
 38.5410  -6.7742    2.4  54.3  -24.9   3.2234   4.8 ------	UNKNOWN                       
 40.2756  -0.0737    2.3  49.1  -30.7   2.3824   4.9 ------	UNKNOWN

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

Time interval of the image:
-11.488000 29.588600
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
326.1098  38.2747   11.2  34.5  108.8  -0.1548  -1.0    4.1	Cyg X-2                       
231.7057  78.5165    1.9  57.6 -179.8   2.5613   6.1 ------	UNKNOWN                       
 19.1739  -9.2379    0.8  51.3   -0.8   2.7717  14.5 ------	UNKNOWN                       
 85.7396  85.9484    2.4  47.5 -165.8   0.5946   4.8 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
100.588610 962.161900
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
  0.0561 -42.4540    2.9  37.2   56.7   3.8860   3.9 ------	UNKNOWN                       
 28.8173 -30.1185    2.9  22.7   11.1   1.6516   3.9 ------	UNKNOWN                       
347.0243  35.5001    2.9  54.1 -179.3   6.1157   4.0 ------	UNKNOWN

Plot creation:
Sun May 31 07:24:05 EDT 2015