/sdf/group/fermi/u/kocevski/software/conda/envs/fermitools_v1.0.5/lib/python2.7/site-packages/pyfits/__init__.py:22: PyFITSDeprecationWarning: PyFITS is deprecated, please use astropy.io.fits
  PyFITSDeprecationWarning)  # noqa

Likelihood Analysis Tool
Support Contact: Daniel Kocevski (daniel.kocevski@nasa.gov)

Importing modules...
This is make3FGLxml version 01r0.
The default diffuse model files and names are for pass 8 and assume you have v10r00p05 of the Fermi Science Tools or higher.
This is make4FGLDR3xml version 01r09.
The default diffuse model files and names are for P8R3, 4FGL DRs, and the default fermitools directory locations.

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

/sdf/group/fermi/u/kocevski/software/conda/envs/fermitools_v1.0.5/lib/python2.7/site-packages/matplotlib/cbook/deprecation.py:106: MatplotlibDeprecationWarning: The mpl_toolkits.axes_grid module was deprecated in version 2.1. Use mpl_toolkits.axes_grid1 and mpl_toolkits.axisartist provies the same functionality instead.
  warnings.warn(message, mplDeprecation, stacklevel=1)

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

number of points = 7550
number of fields = 31
RA_SUN column already exists. File will not be modified.
Done.


Performing analysis on: 4FGL_J1954.6-1122_daily_bin1843
RA = 298.669, Dec = -11.3815
tmin = 717033601, tmax = 717292801
emin = 100.0, emax = 100000.0
irfs = P8R3_SOURCE_V2
ROI = 12
zmax = 90

Fitting strategies:
Fix sources to their catalog values = True
Leave variable source normalizations free = True
Perform pre-fit where ALL sources are frozen to catalog values = False
Perform secondary fit with a free spectral index = True
Remove all catalog sources = False
Remove weak (TS < 2) catalog sources = False
Freeze Galactic component = False

Likelihood fit parameters:
Statistic: UNBINNED
Optimizer: MINUIT

Catalog: 4FGL

Repository job = True

LAT data source:
Data catalog

Updating database...

UPDATE lightcurve_data_v2 SET job_id='28935907', job_status='RUN', time_executed='2023-09-25 00:47:06.497773', time_epoch_executed='1695628026.5' WHERE source_name='4FGL J1954.6-1122' and cadence='daily' and bin_id='1843'

Database successfully updated.
Done.

Creating custom output directory:
/lscratch/kocevski/28935907

Creating custom pfile directory:
/lscratch/kocevski/28935907/pfiles/

Retrieving the FT1 file...

Issuing command:
datacat find --group FT1 --filter 'nMetStop>=717033601.0 && nMetStart<717292801.0' /Data/Flight/Reprocess/P305

Issuing command:
datacat find --group FT1 --filter 'nMetStart>564945418 && nMetStop>=717033601.0 && nMetStart<717292801.0' /Data/Flight/Level1/LPA

Retrieving 47 files...
Done.


Merging files...

Checksum keywords updated successfully.

Merged file saved to:
/lscratch/kocevski/28935907/ft1_4FGL_J1954.6-1122_daily_bin1843.fits


Retrieving the FT1 file...

Issuing command:
datacat find --group FT2 --filter 'nMetStop>=717033601.0 && nMetStart<717292801.0' /Data/Flight/Level1/LPA

Retrieving 47 files...
Done.


Merging files...

Merged file saved to:
/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits


Selecting the photons:
time -p gtselect infile=/lscratch/kocevski/28935907/ft1_4FGL_J1954.6-1122_daily_bin1843.fits outfile=/lscratch/kocevski/28935907/ft1_filteredMKT_4FGL_J1954.6-1122_daily_bin1843.fits ra=298.669 dec=-11.3815 rad=12.0 tmin=717033601.0 tmax=717292801.0 emin=100.0 emax=100000.0 zmin=0.0 zmax=90.0 evclass=128 evtype="INDEF" convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

Done.
real 2.71
user 2.45
sys 0.20

Adding sun position to ft2 file:
/sdf/group/fermi/g/LCRepository/Repository/bin/sunpos /lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits

Selecting the good time intervals:
time -p gtmktime scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits sctable="SC_DATA" filter="IN_SAA!=T && LIVETIME>0 && gtifilter('/sdf/group/fermi/g/LCRepository/Likelihood/GTIs/nogrb.gti',(START+STOP)/2) && (angsep(RA_SUN,DEC_SUN,RA_ZENITH,DEC_ZENITH)>115 || gtifilter('/sdf/group/fermi/g/LCRepository/Likelihood/GTIs/nosolarflares.gti',(START+STOP)/2))" roicut=no evfile=/lscratch/kocevski/28935907/ft1_filteredMKT_4FGL_J1954.6-1122_daily_bin1843.fits evtable="EVENTS" outfile="/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits" apply_filter=yes overwrite=yes header_obstimes=yes tstart=0.0 tstop=0.0 gtifile="default" chatter=3 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

This is gtmktime version HEAD
Applying GTI filter:
IN_SAA!=T && LIVETIME>0 && gtifilter('/sdf/group/fermi/g/LCRepository/Likelihood/GTIs/nogrb.gti',(START+STOP)/2) && (angsep(RA_SUN,DEC_SUN,RA_ZENITH,DEC_ZENITH)>115 || gtifilter('/sdf/group/fermi/g/LCRepository/Likelihood/GTIs/nosolarflares.gti',(START+STOP)/2))
real 0.41
user 0.21
sys 0.15

Number of photons that the gtmktime cut: 1471

Generating the livetime cube:
time -p gtltcube evfile="/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits" evtable="EVENTS" scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits sctable="SC_DATA" outfile=/lscratch/kocevski/28935907/ltcube_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits dcostheta=0.025 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=90.0 chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

Working on file /lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits
.....................!
real 8.36
user 7.98
sys 0.25

Generating the exposure map:
time -p gtexpmap evfile=/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits evtable="EVENTS" scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits sctable="SC_DATA" expcube=/lscratch/kocevski/28935907/ltcube_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits outfile=/lscratch/kocevski/28935907/expmap_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits irfs="P8R3_SOURCE_V2" evtype="INDEF" srcrad=30.0 nlong=88 nlat=88 nenergies=10 submap=no nlongmin=0 nlongmax=0 nlatmin=0 nlatmax=0 chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

The exposure maps generated by this tool are meant
to be used for *unbinned* likelihood analysis only.
Do not use them for binned analyses.
Computing the ExposureMap using /lscratch/kocevski/28935907/ltcube_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits
.....................!
real 67.07
user 62.34
sys 3.71

Generating the xml model:
radLim = 0.001
maxRad = 12
normsOnly = True
freeVariableSources = True
Fixing all sources to their catalog values, except variable sources
Creating file and adding sources from 4FGL
Extended source G24.7+0.6 in ROI with RadialDisk spatial model.
Extended source W 41 in ROI with RadialGaussian spatial model.
Extended source FGES J1836.5-0651 in ROI with RadialDisk spatial model.
Extended source FGES J1838.9-0704 in ROI with RadialDisk spatial model.
Extended source Kes 73 in ROI with RadialDisk spatial model.
Extended source HESS J1841-055 in ROI, make sure /sdf/group/fermi/g/LCRepository/Likelihood/ExtendedSources/Extended_8years/Templates/HESSJ1841-055.fits is the correct path to the extended template.
Extended source Kes 79 in ROI with RadialDisk spatial model.
Extended source W 44 in ROI, make sure /sdf/group/fermi/g/LCRepository/Likelihood/ExtendedSources/Extended_8years/Templates/W44.fits is the correct path to the extended template.
Extended source HESS J1857+026 in ROI with RadialDisk spatial model.
Added 280 point sources and 9 extended sources
If using unbinned likelihood you will need to rerun gtdiffrsp for the extended sources or rerun the makeModel function with optional argument psForce=True

Setting source of interest to:
4FGL J1954.6-1122 @ 298.669 -11.3815

Creating the counts map:
time -p gtbin evfile=/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits outfile=/lscratch/kocevski/28935907/cmap_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits algorithm="CMAP" ebinalg="LOG" emin=30.0 emax=200000.0 ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=160 nypix=160 binsz=0.25 coordsys="CEL" xref=298.669 yref=-11.3815 axisrot=0.0 rafield="RA" decfield="DEC" proj="AIT" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

This is gtbin version HEAD
real 0.39
user 0.23
sys 0.14

Creating the high resolution counts map:
time -p gtbin evfile=/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits outfile=/lscratch/kocevski/28935907/cmap_HighResolution_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits algorithm="CMAP" ebinalg="LOG" emin=30.0 emax=200000.0 ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=500 nypix=500 binsz=0.1 coordsys="CEL" xref=298.669 yref=-11.3815 axisrot=0.0 rafield="RA" decfield="DEC" proj="AIT" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

This is gtbin version HEAD
real 0.38
user 0.21
sys 0.14

Computing the diffuse response:
time -p gtdiffrsp evfile=/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits evtable="EVENTS" scfile=/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits sctable="SC_DATA" srcmdl=/lscratch/kocevski/28935907/Model_4FGL_J1954.6-1122_daily_bin1843.xml irfs="P8R3_SOURCE_V2" evclsmin=0 evclass="INDEF" evtype="INDEF" convert=yes chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

adding source EG_v02
adding source FGES J1836.5-0651
adding source FGES J1838.9-0704
adding source G24.7+0.6
adding source GAL_V02
adding source HESS J1841-055
adding source HESS J1857+026
adding source Kes 73
adding source Kes 79
adding source W 41
adding source W 44
Working on...
/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits.....................!
real 102.23
user 100.34
sys 0.96

Creating the likelihood object...
obs = UnbinnedObs('/lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits', '/lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits', expMap='/lscratch/kocevski/28935907/expmap_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits', expCube='/lscratch/kocevski/28935907/ltcube_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits', irfs='P8R3_SOURCE_V2')
like = UnbinnedAnalysis(obs, '/lscratch/kocevski/28935907/Model_4FGL_J1954.6-1122_daily_bin1843.xml', optimizer='MINUIT')

Using source of interest name: '4FGL J1954.6-1122'

4FGL J1954.6-1122
   Spectrum: PowerLaw
683    Prefactor:  2.246e+00  0.000e+00  0.000e+00  1.000e+06 ( 1.000e-12)
684        Index:  2.443e+00  0.000e+00  0.000e+00  1.000e+01 (-1.000e+00) fixed
685        Scale:  8.706e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed


Likelihood Fit Parameters:
Event file(s): /lscratch/kocevski/28935907/ft1_filtered_4FGL_J1954.6-1122_daily_bin1843.fits
Spacecraft file(s): /lscratch/kocevski/28935907/ft2_4FGL_J1954.6-1122_daily_bin1843.fits
Exposure map: /lscratch/kocevski/28935907/expmap_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits
Exposure cube: /lscratch/kocevski/28935907/ltcube_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.fits
IRFs: P8R3_SOURCE_V2
Source model file: /lscratch/kocevski/28935907/Model_4FGL_J1954.6-1122_daily_bin1843.xml
Optimizer: MINUIT
Statistic: UNBINNED

Performing an iterative likelihood fit...

Fit tolerance = 1

Performing likelihood fit...
Done.

MINUIT Return Code (Fixed Index): 0
logL (Fixed Index) = 14847.3451954

Fit tolerance = 0.0001

Performing likelihood fit...
Done.

MINUIT Return Code (Fixed Index): 0
logL (Fixed Index) = 14847.3425802

Fit tolerance = 1e-08

Performing likelihood fit...
 **********
 **   14 **SET NOWARN 
 **********

 PARAMETER DEFINITIONS:
    NO.   NAME         VALUE      STEP SIZE      LIMITS
     1 'Prefactor '    .44797       1.0000          .0000       .10000E+07
     2 'norm      '    3.0956       1.0000          .0000       .10000E+07
     3 'Prefactor '    84.657       1.0000          .0000       .10000E+07
     4 'Prefactor '    .34570E-03   1.0000          .0000       .10000E+07
     5 'norm      '    3.8295       1.0000          .0000       .10000E+07
     6 'Prefactor '    .45495E-02   1.0000          .0000       .10000E+07
     7 'Prefactor '    30.141       1.0000          .0000       .10000E+07
     8 'Prefactor '    .22750E-02   1.0000          .0000       .10000E+07
     9 'Prefactor '    .31378E-05   1.0000          .0000       .10000E+07
    10 'Prefactor '    75.742       1.0000          .0000       .10000E+07
    11 'Prefactor '    38.335       1.0000          .0000       .10000E+07
    12 'Prefactor '    1.5934       1.0000          .0000       .10000E+07
    13 'Prefactor '    4.7854       1.0000          .0000       .10000E+07
    14 'norm      '    18.556       1.0000          .0000       .10000E+07
    15 'Prefactor '    .10716E-02   1.0000          .0000       .10000E+07
    16 'Prefactor '    .33649       1.0000          .0000       .10000E+07
    17 'norm      '    10.720       1.0000          .0000       .10000E+07
    18 'norm      '    31.648       1.0000          .0000       .10000E+07
    19 'norm      '    1.2748       1.0000          .0000       .10000E+07
    20 'Prefactor '    .50363E-03   1.0000          .0000       .10000E+07
    21 'norm      '    .88696       1.0000          .0000       .10000E+07
    22 'Normalizat'    .79399       1.0000         .10000E-01   10.000    
    23 'Prefactor '    1.0030       1.0000          .0000       10.000    
 **********
 **   15 **SET ERR    .5000    
 **********
 **********
 **   16 **SET GRAD    1.000    
 **********
 **********
 **   17 **MINIMIZE    .1025E+06   .2000E-04
 **********

 MIGRAD MINIMIZATION HAS CONVERGED.

 FCN=   14847.34     FROM MIGRAD    STATUS=CONVERGED    134 CALLS     3093 TOTAL
                     EDM=   .10E-09  STRATEGY=1  ERROR MATRIX UNCERTAINTY=   .6%

  EXT PARAMETER                                   STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  Prefactor     .45009        2.5137        .55994E-07   -.27767E-05
   2  norm          3.1159        1.2357       -.65428E-08    .27478E-02
   3  Prefactor     84.174        27.781        .33262E-07   -.29580E-01
   4  Prefactor     .12801E-06    5.6049        .28482E-07** at limit **
   5  norm          3.8336        3.9169       -.32092E-07    .15252E-03
   6  Prefactor     .47740E-08    8.0406        .38107E-07** at limit **
   7  Prefactor     30.023        18.963        .33252E-08    .30696E-01
   8  Prefactor      .0000        7.1352        .27886E-07** at limit **
   9  Prefactor      .0000        1.3149        .25472E-07** at limit **
  10  Prefactor     75.284        27.976        .32643E-07    .44371E-01
  11  Prefactor     38.129        12.154        .16283E-07   -.38689E-01
  12  Prefactor     1.6039        1.5947       -.38529E-07    .22714E-03
  13  Prefactor     4.7656        6.4862       -.72681E-09   -.48247E-02
  14  norm          18.289        9.4167        .21441E-07    .82386E-04
  15  Prefactor     .75495E-07    8.2063        .12626E-07** at limit **
  16  Prefactor     .13117E-06    3.7564        .31354E-07** at limit **
  17  norm          10.735        9.9386        .89428E-07    .17225E-02
  18  norm          31.613        9.2399        .84444E-08   -.18456E-01
  19  norm          1.2732        .90294       -.35596E-09    .25744E-05
  20  Prefactor     .55511E-10    5.9900        .11799E-06** at limit **
  21  norm          .87382        1.8091        .76792E-07   -.23883E-03
  22  Normalizat    .79592        .26633        .75695E-06    .35317E-05
  23  Prefactor     1.0030        .62542E-01   -.45995E-06    .32791E-04
                               ERR DEF=  .500    
Final values: 
  Prefactor  = 0.450089
  norm       = 3.11594
  Prefactor  = 84.1739
  Prefactor  = 1.28009e-07
  norm       = 3.83356
  Prefactor  = 4.77396e-09
  Prefactor  = 30.0233
  Prefactor  = 0
  Prefactor  = 0
  Prefactor  = 75.2843
  Prefactor  = 38.1293
  Prefactor  = 1.6039
  Prefactor  = 4.76561
  norm       = 18.2891
  Prefactor  = 7.54952e-08
  Prefactor  = 1.31173e-07
  norm       = 10.7353
  norm       = 31.6125
  norm       = 1.27321
  Prefactor  = 5.55112e-11
  norm       = 0.873816
  Normalizat = 0.795924
  Prefactor  = 1.00303
 **********
 **   18 **HESSE 
 **********

 FCN=   14847.34     FROM HESSE     STATUS=OK           430 CALLS     3523 TOTAL
                     EDM=   .42E-06    STRATEGY= 1      ERROR MATRIX ACCURATE 

  EXT PARAMETER                                INTERNAL      INTERNAL  
  NO.   NAME        VALUE          ERROR       STEP SIZE       VALUE   
   1  Prefactor     .45009        1.1905        .56212E-02   -1.5695    
   2  norm          3.1159        14.367        .57614E-03   -1.5673    
   3  Prefactor     84.174        38.642        .10869E-03   -1.5524    
   4  Prefactor     .12801E-06    13.625        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   5  norm          3.8336        3.9280        .51550E-03   -1.5669    
   6  Prefactor     .47740E-08    8.6930        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   7  Prefactor     30.023        23.109        .18218E-03   -1.5598    
   8  Prefactor      .0000        7.1418        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   9  Prefactor      .0000        1.3117        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
  10  Prefactor     75.284        40.740        .11491E-03   -1.5534    
  11  Prefactor     38.129        31.995        .16153E-03   -1.5584    
  12  Prefactor     1.6039        1.6567        .84008E-03   -1.5683    
  13  Prefactor     4.7656        62.245        .45968E-03   -1.5664    
  14  norm          18.289        9.4647        .23223E-03   -1.5622    
  15  Prefactor     .75495E-07    13.105        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
  16  Prefactor     .13117E-06    46.466        .97248E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
  17  norm          10.735        9.9511        .30575E-03   -1.5642    
  18  norm          31.613        26.085        .17778E-03   -1.5596    
  19  norm          1.2732        .90997        .98404E-03   -1.5685    
  20  Prefactor     .55511E-10    5.9924        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
  21  norm          .87382        1.8232        .46666E-03   -1.5689    
  22  Normalizat    .79592        .26552        .14142       -1.0022    
  23  Prefactor     1.0030        .61022E-01    .30317E-01   -.92629    
                               ERR DEF=  .500    
Minuit fit quality: 3   estimated distance: 4.17154e-07
Minuit parameter uncertainties:
  1  1.19047
  2  13.3819
  3  38.6426
  4  0.0026412
  5  3.92802
  6  0.000407912
  7  23.1093
  8  1.3767e-05
  9  2.24866e-07
  10  40.7397
  11  31.9952
  12  1.6567
  13  34.4468
  14  9.46474
  15  0.00198931
  16  0.00493807
  17  9.95113
  18  26.0852
  19  0.909973
  20  3.8276e-05
  21  1.82318
  22  0.265957
  23  0.0610262
Done.

MINUIT Return Code (Fixed Index): 0
logL (Fixed Index) = 14847.3378462

Extracting fit parameters...

4FGL J1954.6-1122
   Spectrum: PowerLaw
683    Prefactor:  0.000e+00  2.249e-07  0.000e+00  1.000e+06 ( 1.000e-12)
684        Index:  2.443e+00  0.000e+00  0.000e+00  1.000e+01 (-1.000e+00) fixed
685        Scale:  8.706e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed

GAL_V02
   Spectrum: PowerLaw
1001    Prefactor:  1.003e+00  6.103e-02  0.000e+00  1.000e+01 ( 1.000e+00)
1002        Index:  0.000e+00  0.000e+00 -1.000e+00  1.000e+00 ( 1.000e+00) fixed
1003        Scale:  1.000e+02  0.000e+00  5.000e+01  2.000e+02 ( 1.000e+00) fixed

EG_v02
   Spectrum: FileFunction
990   Normalization:  7.959e-01  2.660e-01  1.000e-02  1.000e+01 ( 1.000e+00)


TS (Fixed Index): 0.00
Ra Dec: 298.669 -11.382
Photon Flux Upper Limit: 7.49495387849e-08 ph cm-2 s-1
Energy Flux Upper Limit: 0.00580802202526 MeV cm-2 s-1

Saving results to:
/lscratch/kocevski/28935907/likelihoodResults_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.txt

Performing a free index iterative likelihood fit...

Freezing all point sources, except source of interest...

Done.

Freeing candidate source...
Done.

Fit tolerance = 1

Performing likelihood fit...
Done.

MINUIT Return Code (Free Index): 0
logL (Free Index) = 14847.3393971
delta logL = -0.00155083808204

Fit tolerance = 0.0001

Performing likelihood fit...
Done.

MINUIT Return Code (Free Index): 102
logL (Free Index) = 14847.3378556
delta logL = -9.30680653255e-06

Fit tolerance = 1e-08

Performing likelihood fit...
 **********
 **   14 **SET NOWARN 
 **********

 PARAMETER DEFINITIONS:
    NO.   NAME         VALUE      STEP SIZE      LIMITS
     1 'Prefactor '    .32257E-05   1.0000          .0000       .10000E+07
     2 'Index     '    1.5886       1.0000          .0000       10.000    
     3 'Normalizat'    .79496       1.0000         .10000E-01   10.000    
     4 'Prefactor '    1.0033       1.0000          .0000       10.000    
 **********
 **   15 **SET ERR    .5000    
 **********
 **********
 **   16 **SET GRAD    1.000    
 **********
 **********
 **   17 **MINIMIZE    .1025E+06   .2000E-04
 **********

 MIGRAD MINIMIZATION HAS CONVERGED.

 MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
  MNHESS FAILS AND WILL RETURN DIAGONAL MATRIX. 

 FCN=   14847.34     FROM MIGRAD    STATUS=CONVERGED     38 CALLS      244 TOTAL
                     EDM=   .27E-12  STRATEGY=1  ERROR MATRIX UNCERTAINTY=100.0%

  EXT PARAMETER                APPROXIMATE        STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  Prefactor      .0000        1.1086        .42118E-03** at limit **
   2  Index         1.6864        1.0661       -.53288E-05     .0000    
   3  Normalizat    .79592        1.1345       -.64837E-05    .12306E-04
   4  Prefactor     1.0030        1.2959        .64020E-06    .56050E-04
                               ERR DEF=  .500    
Final values: 
  Prefactor  = 0
  Index      = 1.68641
  Normalizat = 0.795924
  Prefactor  = 1.00303
 **********
 **   18 **HESSE 
 **********
  MNHESS FAILS AND WILL RETURN DIAGONAL MATRIX. 

 FCN=   14847.34     FROM HESSE     STATUS=FAILED         9 CALLS      253 TOTAL
                     EDM=   .27E-12  STRATEGY=1  ERROR MATRIX UNCERTAINTY=100.0%

  EXT PARAMETER                APPROXIMATE     INTERNAL      INTERNAL  
  NO.   NAME        VALUE          ERROR       STEP SIZE       VALUE   
   1  Prefactor      .0000        1.1086        .42118E-03   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   2  Index         .28018        1.5434       -.53288E-05   -.72444    
   3  Normalizat    .79592        1.1345       -.64837E-05   -1.0022    
   4  Prefactor     1.0030        1.2959        .64020E-06   -.92629    
                               ERR DEF=  .500    
Minuit fit quality: 1   estimated distance: 2.65513e-13
Minuit parameter uncertainties:
  1  6.75955e-08
  2  1.08105
  3  1.17115
  4  1.33987
shell-init: error retrieving current directory: getcwd: cannot access parent directories: No such file or directory
Done.

MINUIT Return Code (Free Index): 101
logL (Free Index) = 14847.3378462
delta logL = 0.0

Extracting fit parameters...

4FGL J1954.6-1122
   Spectrum: PowerLaw
683    Prefactor:  0.000e+00  6.760e-08  0.000e+00  1.000e+06 ( 1.000e-12)
684        Index:  1.686e+00  1.081e+00  0.000e+00  1.000e+01 (-1.000e+00)
685        Scale:  8.706e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed

GAL_V02
   Spectrum: PowerLaw
1001    Prefactor:  1.003e+00  1.340e+00  0.000e+00  1.000e+01 ( 1.000e+00)
1002        Index:  0.000e+00  0.000e+00 -1.000e+00  1.000e+00 ( 1.000e+00) fixed
1003        Scale:  1.000e+02  0.000e+00  5.000e+01  2.000e+02 ( 1.000e+00) fixed

EG_v02
   Spectrum: FileFunction
990   Normalization:  7.959e-01  1.171e+00  1.000e-02  1.000e+01 ( 1.000e+00)


TS (Free Index): 0.00
Ra Dec: 298.669 -11.382

Creating a high resolution counts map...
INFO: Auto-setting vmin to -1.000e-01 [aplpy.core]
INFO: Auto-setting vmax to  1.110e+00 [aplpy.core]
Saving /lscratch/kocevski/28935907/cmap_HighResolution_717033601_717292801_4FGL_J1954.6-1122_daily_bin1843.png

Cleaning up...

** Nuclear option initiated! **

Updating database...

UPDATE lightcurve_data_v2 SET energy_flux_upper_limit='5.81e-03', job_status='DONE', photon_flux_upper_limit='7.49e-08', time_completed='2023-09-25 01:01:53.820800', job_id='28935907', GAL2='1.00', fit_tolerance='1e-08', EG2='0.80', EG='0.80', time_epoch_completed='1695628913.82', return_code2='101', ts='0.00', comment='SUCCESS', return_code='0', GAL='1.00', dlogl='0.00e+00', ts2='0.00', photon_index_error='0.00', photon_index='2.44' WHERE bin_id='1843' AND source_name='4FGL J1954.6-1122' AND cadence='daily';

Database successfully updated.


Analysis Complete.


Mon Sep 25 01:01:53 PDT 2023

