batgrbproduct analysis for trigger=1038247

batgrbproduct Analysis for Trigger=1038247

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.48
    Process Start Time (local time): Fri Jul 29 23:44:07 EDT 2022
             Host Name: batproc1
   Current Working Dir: /local/data/bat1/alien/Swift_3rdBATcatalog/event/batevent_reproc/trigger1038247
======================================================================

  Trigger: 01038247      Segment: 000

  BAT Trigger Time: 638229765.184   [s; MET]
      Trigger Stop: 638229765.248   [s; MET]
               UTC: 2021-03-23T22:02:18.464260    [includes UTCF correction] 
       Where From?: (*SHORT BURST*) TDRSS scaled map message

  BAT 
    RA: 317.952404831532   Dec: 25.3779613350552  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 317.952404831532   Dec: 25.3779613350552  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 317.944940772893   Dec: 25.3523417898211  [deg; J2000]  
        { 21h 11m 46.8s , +25d 21' 08.4" }
         +/- 1.44634869399025 [arcmin] (estimated 90% radius based on SNR)
         +/- 0.603476226560251 [arcmin] (formal 1-sigma fit error)
    SNR: 17.9538230409482
 
  Angular difference between Flight and Ground position is 1.589546 arcmin
 
  Partial Coding Fraction: 0.765625    [ including projection effects ]

  Duration
     T90: 1.12000000476837 +/- 0.322893082552845
        Measured from: -0.0720000267028809
                   to: 1.04799997806549     [s; relative to TRIGTIME]
     T50: 0.247999906539917 +/- 0.109178796885355
        Measured from: 0.00800001621246338
                   to: 0.25599992275238     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.0959999561309814     
                   to: 0.904000043869019     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: -0.0920000076293945     
                   to: 1.25600004196167     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.105003   0.140407   0.110229   0.028012
           0.011194   0.011800   0.011049   0.009198 [error]
 Peak      0.096267   0.127363   0.107437   0.026438
           0.010015   0.010766   0.010214   0.008364 [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.13e-08 5.90e-08 8.61e-08 2.68e-07 2.44e-07 4.45e-07   erg/cm2
              Error 4.99e-09 5.67e-09 1.08e-08 7.28e-08 2.60e-08 8.11e-08   erg/cm2
Peak flux Power-law 8.99e-01 9.64e-01 7.34e-01 9.11e-01 2.94e+00 3.51e+00   ph/cm2/s
              Error 1.47e-01 9.39e-02 8.45e-02 2.23e-01 2.66e-01 3.27e-01   ph/cm2/s
 
Fluence   Cutoff-PL 3.14e-08 5.90e-08 8.61e-08 2.74e-07 2.44e-07 4.52e-07   erg/cm2
              Error 5.44e-09 6.23e-09 1.12e-08 1.04e-07 2.52e-08 1.10e-07   erg/cm2
Peak flux Cutoff-PL 7.12e-01 4.54e-01 1.96e-01 1.40e-01 1.42e+00 1.49e+00   ph/cm2/s
              Error 2.01e-01 1.39e-01 1.12e-01 6.20e-02 4.24e-01 4.80e-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.456    0.166    -      -           1.90e-02 1.79e-03   53.630 57
CPL   1.453    0.721    9998.4    0.0      1.90e-02 1.25e-02   53.630 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.393    0.164    -      -           2.42e-02 2.20e-03   55.520 57
CPL   1.390    0.325    10000.0    0.0      2.43e-02 2.19e-03   55.530 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: 0.192000 sec.
T50 in the 50-300 keV band: 0.128000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100    1.348   -0.092    1.256
T90     1.120   -0.072    1.048
T50     0.248    0.008    0.256
Peak    1.000   -0.096    0.904

Click here for the BAT GCN circular template

Estimated redshift from machine learning (Ukwatta et al. 2016) - may take a few minute to calcaulte

Image

Pre-slew 15.0-350.0 keV image (Event data)

#   RAcent    DECcent   POSerr    Theta     Phi              Peak Cts       SNR   Name
  317.9449    25.3523   0.0107    30.3162   -22.4353         0.4002490     17.954 TRIG_01038247                 
 
Foreground time interval of the image:
 -0.092   1.256 (delta_t =   1.348 [sec])
Background time interval of the image:
-19.083  -0.092 (delta_t =  18.991 [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 0.064000

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 0.064000

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 0.064000

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   -0.092 sec. to    1.256 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.45562     ( -0.165649   0.166856 )
Norm@50keV  : 1.89549E-02    ( -0.00179892   0.00178719 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   53.63     using 59 bins.
 Reduced chi-squared = 0.940877
# Null hypothesis probability of 6.02e-01 with 57 degrees of freedom
 
Photon flux (15-150 keV)    in 1.348 sec: 2.33926 ( -0.21723  0.2175 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.44287e-07 ( -2.5793e-08  2.61085e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from   -0.092 sec. to    1.256 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.45289     ( -0.574258   0.868206 )
Epeak [keV] : 9998.38      ( -9998.49 -9998.49 )
Norm@50keV  : 1.90168E-02    ( -0.00558584   0.0193319 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   53.63     using 59 bins.
 Reduced chi-squared = 0.957679
# Null hypothesis probability of 5.65e-01 with 56 degrees of freedom
 
Photon flux (15-150 keV)    in 1.348 sec: 1.51509 ( -0.37074  0.31911 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.44647e-07 ( -2.88976e-08  2.14498e-08 ) ergs/cm2

Band function

Time interval is from   -0.092 sec. to    1.256 sec.
 
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       :    (     )
beta        :     (     )
Epeak [keV] :       (   )
Norm@50keV   :     (     )
------------------------------------------------------------
 
Photon flux (15-150 keV)    in  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]    : 15.3465  (-1.89803  2.22318)
R^2/D10^2   : 0.242023  (-0.0898997   0.140834)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   86.67     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 6.85e-03 with 57 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.602e-08   2.627e-09
25- 50  6.500e-08   8.695e-09
50-150  1.052e-07   2.025e-08
15-150  1.863e-07   2.386e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 200.000       (-50.8925      -200)
Norm        : 9.53581       (-0.876541     0.876541)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   55.28     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 5.40e-01 with 57 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.213e-08   7.549e-09
25- 50  6.181e-08   1.794e-08
50-150  1.356e-07   7.174e-08
15-150  2.295e-07   8.680e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 6.14391        (-1.8282      2.5421)
R1^2/D10^2  : 3.77926        (-2.57894      6.82577)
kT2 [keV]   : 26.0667        (-6.04711      15.0953)
R2^2/D10^2  : 3.30498E-02    (-0.0256391   0.051664)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   55.07     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 4.72e-01 with 55 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.156e-08   1.555e-08
25- 50  5.751e-08   1.854e-08
50-150  1.518e-07   6.766e-08
15-150  2.408e-07   9.355e-08

Peak spectrum fit

Power-law model

Time interval is from   -0.096 sec. to    0.904 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.39327     ( -0.163945   0.164686 )
Norm@50keV  : 2.41883E-02    ( -0.00220385   0.00219379 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   55.52     using 59 bins.
 Reduced chi-squared = 0.974035
# Null hypothesis probability of 5.31e-01 with 57 degrees of freedom
 
Photon flux (15-150 keV)    in 1 sec: 2.94127 ( -0.26605  0.26642 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.34078e-07 ( -2.40535e-08  2.43324e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from   -0.096 sec. to    0.904 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.39025     ( -0.449063   0.200291 )
Epeak [keV] : 1.00000E+04      ( -10000 -10000 )
Norm@50keV  : 2.42752E-02    ( -0.00221766   0.00215589 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   55.53     using 59 bins.
 Reduced chi-squared = 0.991607
# Null hypothesis probability of 4.93e-01 with 56 degrees of freedom
 
Photon flux (15-150 keV)    in 1 sec: 1.42409 ( -0.40211  0.44632 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.36331e-07 ( -2.68764e-08  2.17702e-08 ) ergs/cm2

Band function

Time interval is from   -0.096 sec. to    0.904 sec.
 
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       :    (     )
beta        :     (     )
Epeak [keV] :       (   )
Norm@50keV   :     (     )
------------------------------------------------------------
 
Photon flux (15-150 keV)    in  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]    : 16.1426  (-2.04649  2.41004)
R^2/D10^2   : 0.256677  (-0.0960363   0.15125)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   92.77     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 1.93e-03 with 57 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  1.380e-08   2.365e-09
25- 50  5.852e-08   7.245e-09
50-150  1.069e-07   1.876e-08
15-150  1.793e-07   2.220e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 200.000       (-39.3555      -200)
Norm        : 12.0553       (-1.07686     1.07686)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   59.45     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 3.87e-01 with 57 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.002e-08   4.635e-09
25- 50  5.774e-08   1.096e-08
50-150  1.266e-07   6.657e-08
15-150  2.144e-07   7.544e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 6.36847        (-1.96137      2.71772)
R1^2/D10^2  : 4.12874        (-2.82564      12.8653)
kT2 [keV]   : 28.4161        (-6.94063      19.9183)
R2^2/D10^2  : 3.27748E-02    (-0.0262833   0.0533242)
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   57.72     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 3.75e-01 with 55 degrees of freedom
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  2.840e-08   1.229e-08
25- 50  5.310e-08   2.060e-08
50-150  1.531e-07   6.815e-08
15-150  2.346e-07   8.739e-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.45562     ( -0.165649   0.166856 )
Norm@50keV  : 1.89549E-02    ( -0.00179892   0.00178719 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   53.63     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 6.02e-01 with 57 degrees of freedom
 
Photon flux (15-150 keV)    in 1.348 sec: 2.33921 ( -0.217213  0.217491 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.44256e-07 ( -2.58373e-08  2.61144e-08 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.45289     ( -0.574867   0.867274 )
Epeak [keV] : 9998.43      ( -9998.54 -9998.54 )
Norm@50keV  : 1.90168E-02    ( -0.00558544   0.0193282 )
------------------------------------------------------------
#Fit statistic  : Chi-Squared                   53.63     using 59 bins.
 Reduced chi-squared = 
# Null hypothesis probability of 5.65e-01 with 56 degrees of freedom
 
Photon flux (15-150 keV)    in 1.348 sec: 2.33875 ( -0.21383  0.21664 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.442e-07 ( -2.74992e-08  2.60245e-08 ) ergs/cm2

Band function

 
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       :    (     )
beta        :     (     )
Epeak [keV] :       (   )
Norm@50keV   :     (     )
------------------------------------------------------------
 
Photon flux (15-150 keV)    in  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)  = 5.90e-08
S(50-100 keV) = 8.61e-08

T90 vs. Hardness ratio plot

T90                                   = 1.12000000476837 sec.
Hardness ratio (energy fluence ratio) = 1.45932

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.33717
Count Ratio (50-100 keV) / (15-25 keV) = 1.04977

Mask shadow pattern

IMX = 5.404735181931878E-01, IMY = 2.231571670947622E-01

TIME vs. PHA plot around the trigger time

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

TIME vs. DetID plot around the trigger time

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

100 us light curve (15-350 keV)

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

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:
-19.083400 -0.092000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
288.7594  10.9292   31.4  53.7    9.6   0.0798   0.4    2.5	GRS 1915+105                  
299.5880  35.3590    3.2  27.7   16.3   0.2406   3.6    9.4	Cyg X-1                       
307.8575  40.8811    0.0  19.4   11.7   0.2149   3.8   12.2	Cyg X-3                       
326.2987  38.4315    5.4  16.3  -35.6   0.1248   2.1    8.9	Cyg X-2                       
283.0477  13.5075    2.4  54.1   17.2   1.3120   4.8 ------	UNKNOWN                       
306.2217  -1.8534    2.2  59.4  -15.4   1.9746   5.2 ------	UNKNOWN                       
 62.9687  41.4907    2.6  59.6 -161.4   1.7531   4.5 ------	UNKNOWN

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

Time interval of the image:
-0.092000 1.256000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
288.7847  11.0568    8.5  53.6    9.7   0.1004   1.3    6.7	GRS 1915+105                  
299.5374  35.3408    4.6  27.7   16.4  -0.0558  -2.5    8.7	Cyg X-1                       
308.1924  40.7536    0.0  19.3   10.9   0.0081   0.4   12.5	Cyg X-3                       
325.9519  38.2698    4.6  16.5  -34.7   0.0458   2.5   10.8	Cyg X-2                       
317.9450  25.3524    0.6  30.3  -22.4   0.4002  18.0 ------	UNKNOWN                       
 57.2925  38.0522    1.8  59.2 -154.9   0.7243   6.5 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
69.316600 962.859600
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
288.9129  11.0574    6.6  31.0   46.5   2.4087   1.7    9.5	GRS 1915+105                  
299.6121  35.1663    2.4  18.6  104.8   5.0177   4.9    2.4	Cyg X-1                       
308.1315  40.9070    3.2  17.6  133.4   3.2739   3.6    2.9	Cyg X-3                       
326.1594  38.2187    5.9  14.6 -174.9   1.7566   2.0    6.2	Cyg X-2                       
280.7769  22.2548    3.0  34.0   71.4   7.3121   3.9 ------	UNKNOWN                       
310.6977  23.0244    2.7   7.0   50.7   3.2004   4.4 ------	UNKNOWN                       
302.4285  -3.6034    2.5  32.7    8.4   4.7651   4.7 ------	UNKNOWN                       
 19.1401  43.7174    3.4  52.3 -148.2   8.9097   3.4 ------	UNKNOWN

Plot creation:
Sat Jul 30 07:08:48 EDT 2022