aklog: Couldn't determine realm of user:aklog: No credentials cache found  while getting realm
/nfs/farm/g/glast/g/LCRepository/bin/miniconda2/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)


Waiting 223 seconds before starting...
Proceeding.

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 make4FGLxml version 01r04.
The default diffuse model files and names are for pass 8 and 4FGL and assume you have v11r5p3 of the Fermi Science Tools or higher.

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

/nfs/farm/g/glast/g/LCRepository/bin/miniconda2/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 = 74146
number of fields = 31
RA_SUN column already exists. File will not be modified.
Done.


Performing analysis on: 4FGL_J0348.8+4610_monthly_bin168
RA = 57.2185, Dec = 46.1695
tmin = 672451201, tmax = 675043201
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='614655', job_status='RUN', time_executed='2022-05-24 02:05:51.658523', time_epoch_executed='1653383151.66' WHERE source_name='4FGL J0348.8+4610' and cadence='monthly' and bin_id='168'

Database successfully updated.
Done.

Creating custom output directory:
/scratch/kocevski/614655

Creating custom pfile directory:
/scratch/kocevski/614655/pfiles/

Retrieving the FT1 file...

Issuing command:
/afs/slac.stanford.edu/g/glast/ground/bin/datacat find --group FT1 --filter 'nMetStop>=672451201.0 && nMetStart<675043201.0' /Data/Flight/Reprocess/P305

Issuing command:
/afs/slac.stanford.edu/g/glast/ground/bin/datacat find --group FT1 --filter 'nMetStart>564945418 && nMetStop>=672451201.0 && nMetStart<675043201.0' /Data/Flight/Level1/LPA

Retrieving 455 files...
Done.


Merging files...

Checksum keywords updated successfully.

Merged file saved to:
/scratch/kocevski/614655/ft1_4FGL_J0348.8+4610_monthly_bin168.fits


Retrieving the FT1 file...

Issuing command:
/afs/slac.stanford.edu/g/glast/ground/bin/datacat find --group FT2 --filter 'nMetStop>=672451201.0 && nMetStart<675043201.0' /Data/Flight/Level1/LPA

Retrieving 455 files...
Done.


Merging files...

Merged file saved to:
/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits


Selecting the photons:
time -p gtselect infile=/scratch/kocevski/614655/ft1_4FGL_J0348.8+4610_monthly_bin168.fits outfile=/scratch/kocevski/614655/ft1_filteredMKT_4FGL_J0348.8+4610_monthly_bin168.fits ra=57.2185 dec=46.1695 rad=12.0 tmin=672451201.0 tmax=675043201.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 28.77
user 27.93
sys 0.75

Adding sun position to ft2 file:
/nfs/farm/g/glast/g/LCRepository/Repository/bin/sunpos /scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits

Selecting the good time intervals:
time -p gtmktime scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits sctable="SC_DATA" filter="IN_SAA!=T && LIVETIME>0 && gtifilter('/nfs/farm/g/glast/g/LCRepository/Likelihood/GTIs/nogrb.gti',(START+STOP)/2) && (angsep(RA_SUN,DEC_SUN,RA_ZENITH,DEC_ZENITH)>115 || gtifilter('/nfs/farm/g/glast/g/LCRepository/Likelihood/GTIs/nosolarflares.gti',(START+STOP)/2))" roicut=no evfile=/scratch/kocevski/614655/ft1_filteredMKT_4FGL_J0348.8+4610_monthly_bin168.fits evtable="EVENTS" outfile="/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.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('/nfs/farm/g/glast/g/LCRepository/Likelihood/GTIs/nogrb.gti',(START+STOP)/2) && (angsep(RA_SUN,DEC_SUN,RA_ZENITH,DEC_ZENITH)>115 || gtifilter('/nfs/farm/g/glast/g/LCRepository/Likelihood/GTIs/nosolarflares.gti',(START+STOP)/2))
real 0.88
user 0.66
sys 0.20

Number of photons that the gtmktime cut: 18583

Generating the livetime cube:
time -p gtltcube evfile="/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits" evtable="EVENTS" scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits sctable="SC_DATA" outfile=/scratch/kocevski/614655/ltcube_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.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 /scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits
.....................!
real 136.91
user 136.24
sys 0.39

Generating the exposure map:
time -p gtexpmap evfile=/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits evtable="EVENTS" scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits sctable="SC_DATA" expcube=/scratch/kocevski/614655/ltcube_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits outfile=/scratch/kocevski/614655/expmap_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.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 /scratch/kocevski/614655/ltcube_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits
.....................!
real 192.26
user 186.22
sys 5.60

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 G150.3+4.5 in ROI with RadialDisk spatial model.
Extended source HB 9 in ROI, make sure /nfs/farm/g/glast/g/LCRepository/Likelihood/ExtendedSources/Extended_archive_v18/Templates/HB9.fits is the correct path to the extended template.
Extended source HB 3 in ROI with RadialDisk spatial model.
Extended source W 3 in ROI, make sure /nfs/farm/g/glast/g/LCRepository/Likelihood/ExtendedSources/Extended_archive_v18/Templates/W3.fits is the correct path to the extended template.
Added 169 point sources and 4 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 J0348.8+4610 @ 57.2185 46.1695

Creating the counts map:
time -p gtbin evfile=/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits outfile=/scratch/kocevski/614655/cmap_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.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=57.2185 yref=46.1695 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.82
user 0.63
sys 0.17

Creating the high resolution counts map:
time -p gtbin evfile=/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits outfile=/scratch/kocevski/614655/cmap_HighResolution_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.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=57.2185 yref=46.1695 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.82
user 0.64
sys 0.17

Computing the diffuse response:
time -p gtdiffrsp evfile=/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits evtable="EVENTS" scfile=/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits sctable="SC_DATA" srcmdl=/scratch/kocevski/614655/Model_4FGL_J0348.8+4610_monthly_bin168.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 G150.3+4.5
adding source GAL_V02
adding source HB 3
adding source HB 9
adding source W 3
Working on...
/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits.....................!
real 3415.36
user 3406.28
sys 1.95

Creating the likelihood object...
obs = UnbinnedObs('/scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits', '/scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits', expMap='/scratch/kocevski/614655/expmap_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits', expCube='/scratch/kocevski/614655/ltcube_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits', irfs='P8R3_SOURCE_V2')
like = UnbinnedAnalysis(obs, '/scratch/kocevski/614655/Model_4FGL_J0348.8+4610_monthly_bin168.xml', optimizer='MINUIT')

Using source of interest name: '4FGL J0348.8+4610'

4FGL J0348.8+4610
   Spectrum: LogParabola
299         norm:  8.952e+00  0.000e+00  1.000e-04  1.000e+04 ( 1.000e-13)
300        alpha:  2.791e+00  0.000e+00  0.000e+00  5.000e+00 ( 1.000e+00) fixed
301         beta:  7.194e-01  0.000e+00  0.000e+00  1.000e+01 ( 1.000e+00) fixed
302           Eb:  9.019e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed


Likelihood Fit Parameters:
Event file(s): /scratch/kocevski/614655/ft1_filtered_4FGL_J0348.8+4610_monthly_bin168.fits
Spacecraft file(s): /scratch/kocevski/614655/ft2_4FGL_J0348.8+4610_monthly_bin168.fits
Exposure map: /scratch/kocevski/614655/expmap_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits
Exposure cube: /scratch/kocevski/614655/ltcube_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.fits
IRFs: P8R3_SOURCE_V2
Source model file: /scratch/kocevski/614655/Model_4FGL_J0348.8+4610_monthly_bin168.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) = 178450.615643

Fit tolerance = 0.0001

Performing likelihood fit...
Done.

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

Fit tolerance = 1e-08

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

 PARAMETER DEFINITIONS:
    NO.   NAME         VALUE      STEP SIZE      LIMITS
     1 'Prefactor '    .10015E-03   1.0000         .10000E-03   10000.    
     2 'norm      '    1.8556       1.0000         .10000E-03   10000.    
     3 'Prefactor '    .66870E-03   1.0000         .10000E-03   10000.    
     4 'Prefactor '    .10000E-03   1.0000         .10000E-03   10000.    
     5 'norm      '    2.1200       1.0000         .10000E-03   10000.    
     6 'Prefactor '    40.965       1.0000         .10000E-03   10000.    
     7 'Prefactor '    22.504       1.0000         .10000E-03   10000.    
     8 'norm      '    .10001E-03   1.0000         .10000E-03   10000.    
     9 'Prefactor '    .50438E-02   1.0000         .10000E-03   10000.    
    10 'Prefactor '    2.2421       1.0000         .10000E-03   10000.    
    11 'norm      '    16.749       1.0000         .10000E-03   10000.    
    12 'Prefactor '    2.7982       1.0000         .10000E-03   10000.    
    13 'norm      '    .28284       1.0000         .10000E-03   10000.    
    14 'Normalizat'    1.4410       1.0000         .10000E-01   10.000    
    15 'Prefactor '    .95634       1.0000          .0000       10.000    
 **********
 **   15 **SET ERR    .5000    
 **********
 **********
 **   16 **SET GRAD    1.000    
 **********
 **********
 **   17 **MINIMIZE    .5940E+05   .2000E-04
 **********

 MIGRAD MINIMIZATION HAS CONVERGED.

 MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.

 MIGRAD MINIMIZATION HAS CONVERGED.

 FCN=   178450.5     FROM MIGRAD    STATUS=CONVERGED    323 CALLS     1282 TOTAL
                     EDM=   .11E-10  STRATEGY=1  ERROR MATRIX UNCERTAINTY=   .6%

  EXT PARAMETER                                   STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  Prefactor     .10000E-03    .13602       -.11968E-07** at limit **
   2  norm          1.8572        .98840       -.11762E-08   -.13583E-02
   3  Prefactor     .10002E-03    .76399        .74762E-07** at limit **
   4  Prefactor     .10000E-03    5.9948        .87072E-07** at limit **
   5  norm          2.1200        .27453       -.11561E-07   -.47485E-03
   6  Prefactor     40.955        24.783       -.86292E-05    .66757E-06
   7  Prefactor     22.504        11.366       -.15595E-06   -.62359E-05
   8  norm          .10000E-03    8.3241        .16089E-06** at limit **
   9  Prefactor     .10000E-03    1.5408       -.64020E-07** at limit **
  10  Prefactor     2.2420        .76629        .53766E-09   -.21150E-02
  11  norm          16.749        3.5971       -.24571E-07   -.26361E-04
  12  Prefactor     2.7980        2.6241        .69958E-08   -.44029E-03
  13  norm          .28752        .54068       -.22586E-07    .87893E-04
  14  Normalizat    1.4410        .11599        .75692E-06   -.38141E-04
  15  Prefactor     .95634        .12565E-01   -.27748E-07    .25131E-03
                               ERR DEF=  .500    
Final values: 
  Prefactor  = 0.0001
  norm       = 1.85723
  Prefactor  = 0.000100016
  Prefactor  = 0.0001
  norm       = 2.11999
  Prefactor  = 40.9548
  Prefactor  = 22.5042
  norm       = 0.0001
  Prefactor  = 0.0001
  Prefactor  = 2.24201
  norm       = 16.7493
  Prefactor  = 2.79797
  norm       = 0.287515
  Normalizat = 1.44102
  Prefactor  = 0.956342
 **********
 **   18 **HESSE 
 **********

 FCN=   178450.5     FROM HESSE     STATUS=OK           208 CALLS     1490 TOTAL
                     EDM=   .24E-08    STRATEGY= 1      ERROR MATRIX ACCURATE 

  EXT PARAMETER                                INTERNAL      INTERNAL  
  NO.   NAME        VALUE          ERROR       STEP SIZE       VALUE   
   1  Prefactor     .10000E-03    .43467        .30079E-02   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   2  norm          1.8572        1.1676        .15299E-02   -1.5435    
   3  Prefactor     .10002E-03    .94163        .10968       -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   4  Prefactor     .10000E-03    6.0039        .18142E-02   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   5  norm          2.1200        .27511        .73336E-04   -1.5417    
   6  Prefactor     40.955        24.787        .62629E-04   -1.4427    
   7  Prefactor     22.504        11.437        .84437E-04   -1.4759    
   8  norm          .10000E-03    8.3524        .18142E-02   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   9  Prefactor     .10000E-03    1.5216        .17525       -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
  10  Prefactor     2.2420        1.2058        .13724E-02   -1.5408    
  11  norm          16.749        3.5993        .97866E-04   -1.4889    
  12  Prefactor     2.7980        2.7742        .12160E-02   -1.5373    
  13  norm          .28752        2.0273        .39915E-02   -1.5601    
  14  Normalizat    1.4410        .11882        .86439E-03   -.79449    
  15  Prefactor     .95634        .13023E-01    .11833E-03   -.94199    
                               ERR DEF=  .500    
Minuit fit quality: 3   estimated distance: 2.40645e-09
Minuit parameter uncertainties:
  1  4.14993e-06
  2  1.16759
  3  0.00024928
  4  7.39671e-07
  5  0.275107
  6  24.7932
  7  11.438
  8  4.92645e-07
  9  1.72566e-05
  10  1.20581
  11  3.59936
  12  2.77435
  13  1.52674
  14  0.118844
  15  0.0130232
Done.

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

Extracting fit parameters...

4FGL J0348.8+4610
   Spectrum: LogParabola
299         norm:  1.000e-04  4.926e-07  1.000e-04  1.000e+04 ( 1.000e-13)
300        alpha:  2.791e+00  0.000e+00  0.000e+00  5.000e+00 ( 1.000e+00) fixed
301         beta:  7.194e-01  0.000e+00  0.000e+00  1.000e+01 ( 1.000e+00) fixed
302           Eb:  9.019e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed

GAL_V02
   Spectrum: PowerLaw
579    Prefactor:  9.563e-01  1.302e-02  0.000e+00  1.000e+01 ( 1.000e+00)
580        Index:  0.000e+00  0.000e+00 -1.000e+00  1.000e+00 ( 1.000e+00) fixed
581        Scale:  1.000e+02  0.000e+00  5.000e+01  2.000e+02 ( 1.000e+00) fixed

EG_v02
   Spectrum: FileFunction
574   Normalization:  1.441e+00  1.188e-01  1.000e-02  1.000e+01 ( 1.000e+00)


TS (Fixed Index): -0.00
Ra Dec: 57.218 46.169
Photon Flux Upper Limit: 1.56e-08 ph cm-2 s-1
Energy Flux Upper Limit: 6.40e-06 MeV cm-2 s-1

Saving results to:
/scratch/kocevski/614655/likelihoodResults_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.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): 102
logL (Free Index) = 178450.526107
delta logL = -0.000570752541535

Fit tolerance = 0.0001

Performing likelihood fit...
Done.

MINUIT Return Code (Free Index): 102
logL (Free Index) = 178450.525536
delta logL = 8.59785359353e-07

Fit tolerance = 1e-08

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

 PARAMETER DEFINITIONS:
    NO.   NAME         VALUE      STEP SIZE      LIMITS
     1 'norm      '    .10401E-03   1.0000         .10000E-03   10000.    
     2 'alpha     '    1.9132       1.0000          .0000       5.0000    
     3 'Normalizat'    1.4408       1.0000         .10000E-01   10.000    
     4 'Prefactor '    .95635       1.0000          .0000       10.000    
 **********
 **   15 **SET ERR    .5000    
 **********
 **********
 **   16 **SET GRAD    1.000    
 **********
 **********
 **   17 **MINIMIZE    .5940E+05   .2000E-04
 **********

 MIGRAD MINIMIZATION HAS CONVERGED.

 MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.

 MIGRAD MINIMIZATION HAS CONVERGED.

 MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
 MIGRAD TERMINATED WITHOUT CONVERGENCE.

 FCN=   178450.5     FROM MIGRAD    STATUS=FAILED       103 CALLS      283 TOTAL
                     EDM=   .24E-07    STRATEGY= 1      ERR MATRIX NOT POS-DEF

  EXT PARAMETER                APPROXIMATE        STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  norm          .10000E-03    6419.3        .18142E-02** at limit **
   2  alpha         1.6702        4.0202        .17074        .76707E-07
   3  Normalizat    1.4410        .62150E-01    .86438E-03   -.27218E-02
   4  Prefactor     .95634        .72119E-02    .11833E-03   -.20604E-01
                               ERR DEF=  .500    
 START SIMPLEX MINIMIZATION.    CONVERGENCE WHEN EDM .LT.   .20E-04
 SIMPLEX MINIMIZATION HAS CONVERGED.

 FCN=   178450.5     FROM SIMPLEX   STATUS=UNCHANGED    199 CALLS      482 TOTAL
                     EDM=   .17E-04    STRATEGY= 1      ERR MATRIX NOT POS-DEF

  EXT PARAMETER                APPROXIMATE        STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  norm          .10000E-03    6419.3        .76013E-03** at limit **
   2  alpha         1.6702        4.0202        .68238E-07    .76707E-07
   3  Normalizat    1.4410        .62150E-01    .44376E-08   -.27218E-02
   4  Prefactor     .95634        .72119E-02    .25023E-09   -.20604E-01
                               ERR DEF=  .500    

 MIGRAD MINIMIZATION HAS CONVERGED.

 MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.

 FCN=   178450.5     FROM MIGRAD    STATUS=CONVERGED     31 CALLS      513 TOTAL
                     EDM=   .40E-08    STRATEGY= 1      ERROR MATRIX ACCURATE 

  EXT PARAMETER                                   STEP         FIRST   
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE 
   1  norm          .10000E-03    16.984        .36285E-03** at limit **
   2  alpha         1.6673        4.1658        .68295E-02    .42615E-07
   3  Normalizat    1.4410        .10864        .34575E-04   -.27167E-02
   4  Prefactor     .95634        .12608E-01    .23667E-04   -.20549E-01
                               ERR DEF=  .500    
Final values: 
  norm       = 0.000100002
  alpha      = 1.66727
  Normalizat = 1.44102
  Prefactor  = 0.956342
 **********
 **   18 **HESSE 
 **********

 FCN=   178450.5     FROM HESSE     STATUS=OK            25 CALLS      538 TOTAL
                     EDM=   .40E-08    STRATEGY= 1      ERROR MATRIX ACCURATE 

  EXT PARAMETER                                INTERNAL      INTERNAL  
  NO.   NAME        VALUE          ERROR       STEP SIZE       VALUE   
   1  norm          .10000E-03    16.980        .72569E-04   -1.5708    
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   2  alpha         1.6673        3.2235        .13659E-01   -.33958    
   3  Normalizat    1.4410        .10862        .69151E-05   -.79450    
   4  Prefactor     .95634        .12607E-01    .47334E-05   -.94199    
                               ERR DEF=  .500    
Minuit fit quality: 3   estimated distance: 4.0412e-09
Minuit parameter uncertainties:
  1  0.000369891
  2  477.849
  3  0.108634
  4  0.0126066
shell-init: error retrieving current directory: getcwd: cannot access parent directories: No such file or directory
rmdir: failed to remove `/scratch/kocevski/614655': No such file or directory
shell-init: error retrieving current directory: getcwd: cannot access parent directories: No such file or directory
Done.

MINUIT Return Code (Free Index): 0
logL (Free Index) = 178450.525533
delta logL = 3.03838169202e-06

Extracting fit parameters...

4FGL J0348.8+4610
   Spectrum: LogParabola
299         norm:  1.000e-04  3.699e-04  1.000e-04  1.000e+04 ( 1.000e-13)
300        alpha:  1.667e+00  4.778e+02  0.000e+00  5.000e+00 ( 1.000e+00)
301         beta:  7.194e-01  0.000e+00  0.000e+00  1.000e+01 ( 1.000e+00) fixed
302           Eb:  9.019e+02  0.000e+00  3.000e+01  5.000e+05 ( 1.000e+00) fixed

GAL_V02
   Spectrum: PowerLaw
579    Prefactor:  9.563e-01  1.261e-02  0.000e+00  1.000e+01 ( 1.000e+00)
580        Index:  0.000e+00  0.000e+00 -1.000e+00  1.000e+00 ( 1.000e+00) fixed
581        Scale:  1.000e+02  0.000e+00  5.000e+01  2.000e+02 ( 1.000e+00) fixed

EG_v02
   Spectrum: FileFunction
574   Normalization:  1.441e+00  1.086e-01  1.000e-02  1.000e+01 ( 1.000e+00)


TS (Free Index): -0.00
Ra Dec: 57.218 46.169

Creating a high resolution counts map...
INFO: Auto-setting vmin to -3.000e-01 [aplpy.core]
INFO: Auto-setting vmax to  3.330e+00 [aplpy.core]
Saving /scratch/kocevski/614655/cmap_HighResolution_672451201_675043201_4FGL_J0348.8+4610_monthly_bin168.png

Cleaning up...

** Nuclear option initiated! **
rmdir /scratch/kocevski/614655

Updating database...

UPDATE lightcurve_data_v2 SET energy_flux_upper_limit='6.40e-06', job_status='DONE', photon_flux_upper_limit='1.56e-08', time_completed='2022-05-24 04:49:43.347724', job_id='614655', GAL2='0.96', fit_tolerance='1e-08', EG2='1.44', EG='1.44', time_epoch_completed='1653392983.35', return_code2='0', ts='-0.00', comment='SUCCESS', return_code='0', GAL='0.96', dlogl='3.04e-06', ts2='-0.00', photon_index_error='0.00', photon_index='2.79' WHERE bin_id='168' AND source_name='4FGL J0348.8+4610' AND cadence='monthly';

Database successfully updated.


Analysis Complete.


Tue May 24 04:49:43 PDT 2022


------------------------------------------------------------
Sender: LSF System <lsf@bullet0082>
Subject: Job 614655: <4FGL_J0348.8+4610_monthly_bin168> in cluster <slac> Done

Job <4FGL_J0348.8+4610_monthly_bin168> was submitted from host <rhel6-64q> by user <kocevski> in cluster <slac> at Tue May 24 02:00:36 2022
Job was executed on host(s) <bullet0082>, in queue <medium>, as user <kocevski> in cluster <slac> at Tue May 24 02:02:00 2022
</u/ki/kocevski> was used as the home directory.
</u/ki/kocevski> was used as the working directory.
Started at Tue May 24 02:02:00 2022
Terminated at Tue May 24 04:49:44 2022
Results reported at Tue May 24 04:49:44 2022

Your job looked like:

------------------------------------------------------------
# LSBATCH: User input
/nfs/farm/g/glast/g/LCRepository/Likelihood/Scripts/LikelihoodAnalysis.py 4FGL_J0348.8+4610_monthly_bin168 57.2185 46.1695 672451201 675043201 irfs=P8R3_SOURCE_V2 zmax=90 ROI=12 association=4FGL_J0348.8+4610 fixedModel=True performFreeIndexFit=True freeVariableSources=True nuke=True repositoryJob=True verbose=False
------------------------------------------------------------

Successfully completed.

Resource usage summary:

    CPU time :                                   9801.00 sec.
    Max Memory :                                 1101 MB
    Average Memory :                             567.14 MB
    Total Requested Memory :                     -
    Delta Memory :                               -
    Max Swap :                                   -
    Max Processes :                              5
    Max Threads :                                21
    Run time :                                   10064 sec.
    Turnaround time :                            10148 sec.

The output (if any) is above this job summary.

