PKS1510-089 analysis notebook¶

Reference paper: https://www.aanda.org/articles/aa/pdf/2011/05/aa16128-10.pdf

In [1]:
from agilepy.api.AGAnalysis import AGAnalysis
In [2]:
# Interactive plots
# %matplotlib widget
! root -q
   ------------------------------------------------------------------
  | Welcome to ROOT 6.26/06                        https://root.cern |
  | (c) 1995-2021, The ROOT Team; conception: R. Brun, F. Rademakers |
  | Built for linuxx8664gcc on Jul 28 2022, 18:08:51                 |
  | From tags/v6-26-06@v6-26-06                                      |
  | With c++ (GCC) 8.5.0 20210514 (Red Hat 8.5.0-13)                 |
  | Try '.help', '.demo', '.license', '.credits', '.quit'/'.q'       |
   ------------------------------------------------------------------


Creating a configuration file¶

In [3]:
confFilePath = "$HOME/agilepy_conf.yaml"

Set the configuration and download the data from SSDC data server

In [4]:
AGAnalysis.getConfiguration(
    confFilePath = confFilePath,
    evtfile=None,
    logfile=None,
    userName = "username",
    sourceName = "PKS1510-089",
    tmin = 54891,
    tmax = 54921,
    timetype = "MJD",
    glon = 351.29,
    glat = 40.13,
    outputDir = "$HOME/agilepy_analysis",
    verboselvl = 0,
    userestapi=True,
    datapath="$HOME/agile_dataset"
)

Obtaining the AGAnalysis object¶

In [5]:
ag = AGAnalysis(confFilePath)
In [6]:
ag.setOptions(binsize=0.5)
In [7]:
#print all options of the configuration file
ag.printOptions()
{ 'ap': {'radius': 3, 'timeslot': 3600},
  'input': { 'datapath': '/root/agile_dataset',
             'evtfile': PosixPath('/root/agile_dataset/EVT.index'),
             'logfile': PosixPath('/root/agile_dataset/LOG.index'),
             'userestapi': True},
  'maps': { 'binsize': 0.5,
            'energybins': [[100, 10000]],
            'expstep': 4.0,
            'fovbinnumber': 1,
            'mapsize': 40,
            'offaxisangle': 30,
            'proj': 'ARC',
            'projtype': 'WCS',
            'spectralindex': 2.1,
            'timestep': 160,
            'useEDPmatrixforEXP': False},
  'mle': { 'contourpoints': 40,
           'edpcorrection': 1,
           'expratio_maxthr': 15,
           'expratio_minthr': 0,
           'expratio_size': 10,
           'expratioevaluation': True,
           'fluxcorrection': 0,
           'integratortype': 1,
           'loccl': 5.99147,
           'mindefaulttolerance': 0.01,
           'minimizeralg': 'Migrad',
           'minimizerdefstrategy': 2,
           'minimizertype': 'Minuit',
           'ranal': 10,
           'ulcl': 2},
  'model': { 'emax_sources': 10000,
             'emin_sources': 100,
             'galcoeff': [-1],
             'galmode': 1,
             'galmode2': 0,
             'galmode2fit': 0,
             'isocoeff': [-1],
             'isomode': 1,
             'isomode2': 0,
             'isomode2fit': 0,
             'modelfile': None},
  'output': { 'filenameprefix': 'analysis_product',
              'logfilenameprefix': 'analysis_log',
              'outdir': PosixPath('/root/agilepy_analysis/username_PKS1510-089_20221022-184850'),
              'sourcename': 'PKS1510-089',
              'username': 'username',
              'verboselvl': 0},
  'plotting': {'twocolumns': False},
  'selection': { 'albedorad': 80,
                 'bpointing': None,
                 'dq': 0,
                 'emax': 10000,
                 'emin': 100,
                 'filtercode': 5,
                 'fovradmax': 60,
                 'fovradmin': 0,
                 'glat': 40.13,
                 'glon': 351.290001,
                 'lonpole': 180,
                 'lpointing': None,
                 'maplistgen': 'None',
                 'phasecode': 18,
                 'proj': 'ARC',
                 'timelist': 'None',
                 'timetype': 'TT',
                 'tmax': 165542400.0,
                 'tmin': 162950400.0}}

Sources hypothesis¶

In [8]:
#sources = ag.loadSourcesFromCatalog("2AGL", rangeDist = (0, 5), show=True)

#aggiungeto una sorgente in questo modo non funziona più il selectSources
newSourceDict = {
    "glon" : 351.29,
    "glat": 40.13,
    "spectrumType" : "PowerLaw",
    "flux": 2.0e-07,
    "index": 2.1
}
newSource = ag.addSource("PKS1510-089", newSourceDict)

Selecting sources¶

In [9]:
sources = ag.selectSources("flux > 0", show = True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource)
  * Spectrum type: PowerLaw
  * Free parameters: none
  * Initial source parameters:
	- flux (ph/cm2s): 2.0000e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0

Free a source's parameter¶

Freeable params = [flux, index, index1, index2, cutoffEnergy, pivotEnergy, curvature, pos]

In [10]:
affectedSources = ag.freeSources('name == "PKS1510-089"', "flux", True, show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource)
  * Spectrum type: PowerLaw
  * Free parameters: flux
  * Initial source parameters:
	- flux (ph/cm2s): 2.0000e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0

In this example, only one source is affected.

Analysis¶

In [11]:
#Generate maps
maplistfile = ag.generateMaps()
286it [00:41,  6.94it/s]
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 33/33 [00:00<00:00, 21911.04it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
In [12]:
#Display counts map
ag.displayCtsSkyMaps(maplistFile=maplistfile, smooth=3, catalogRegions="2AGL", catalogRegionsColor="green")
Generating plot..please wait.
2022-10-22 18:50:52,143 [WARNING ] [AGAnalysis] singleMode has been turned off because only one map is going to be displayed.
WARNING: FITSFixedWarning: 'datfix' made the change 'Set DATEREF to '2004-01-01T00:01:06.184' from MJDREF.
Set MJD-OBS to 54891.000764 from DATE-OBS.
Set MJD-END to 54921.000764 from DATE-END'. [astropy.wcs.wcs]
Out[12]:
[None]
In [13]:
#Fix the diffuse emission coefficient to a default value. 
ag.setOptions(galcoeff=[0.7])
In [14]:
ag.printOptions("model")
{ 'emax_sources': 10000,
  'emin_sources': 100,
  'galcoeff': [0.7],
  'galmode': 1,
  'galmode2': 0,
  'galmode2fit': 0,
  'isocoeff': [-1],
  'isomode': 1,
  'isomode2': 0,
  'isomode2fit': 0,
  'modelfile': None}
In [15]:
#Perform a maximum likelihood estimator
ag.mle()
Out[15]:
['/root/agilepy_analysis/username_PKS1510-089_20221022-184850/mle/0/analysis_product_PKS1510-089.source']
In [16]:
#Display the results
ag.selectSources("sqrtTS > 0", show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 19.8115
  * Spectrum type: PowerLaw
  * Free parameters: flux
  * Initial source parameters:
	- flux (ph/cm2s): 2.0000e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:50:54.935770):
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- counts: 619.339 +/- 44.0936
	- upper limit(ph/cm2s): 1.98248e-06
	- ergLog(erg/cm2s): 2.79476e-10 +/- 1.98972e-11
	- galCoeff: [0.7]
	- isoCoeff: [11.0349] +/- [0.367294]
	- exposure(cm2s): 358153000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------
Out[16]:
[<agilepy.core.source.Source.PointSource at 0x7fc73a13ef70>]

Spectral index evaluation¶

In [17]:
#energybin=2 to calcolate spectral index
ag.setOptionEnergybin(2)
ag.setOptions(galcoeff=[0.7, 0.7, 0.7, 0.7])
ag.setOptions(isocoeff=[-1, -1, -1, -1])
In [18]:
maplistfile = ag.generateMaps()
286it [00:41,  6.82it/s]
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 33/33 [00:00<00:00, 43430.20it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
In [19]:
ag.displayCtsSkyMaps(maplistFile=maplistfile, smooth=3, normType="power", catalogRegions="2AGL", catalogRegionsColor="green")
Generating plot..please wait.
WARNING: FITSFixedWarning: 'datfix' made the change 'Set DATEREF to '2004-01-01T00:01:06.184' from MJDREF.
Set MJD-OBS to 54891.000764 from DATE-OBS.
Set MJD-END to 54921.000764 from DATE-END'. [astropy.wcs.wcs]
Out[19]:
[None]
In [20]:
affectedSources = ag.freeSources('name == "PKS1510-089"', "index", True, show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 19.8115
  * Spectrum type: PowerLaw
  * Free parameters: flux index
  * Initial source parameters:
	- flux (ph/cm2s): 2.0000e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:50:54.935770):
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- counts: 619.339 +/- 44.0936
	- upper limit(ph/cm2s): 1.98248e-06
	- ergLog(erg/cm2s): 2.79476e-10 +/- 1.98972e-11
	- galCoeff: [0.7]
	- isoCoeff: [11.0349] +/- [0.367294]
	- exposure(cm2s): 358153000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------
In [21]:
ag.mle()
Out[21]:
['/root/agilepy_analysis/username_PKS1510-089_20221022-184850/mle/1/analysis_product_PKS1510-089.source']
In [22]:
selectedSources = ag.selectSources('flux > 0', show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 20.0843
  * Spectrum type: PowerLaw
  * Free parameters: flux index
  * Initial source parameters:
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:56:59.510740):
	- flux (ph/cm2s): 1.5715e-06 +/- 1.4153e-07
	- index : 2.07581 +/- 0.0710626
	- counts: 563.679 +/- 50.7665
	- upper limit(ph/cm2s): 1.80527e-06
	- ergLog(erg/cm2s): 2.5398e-10 +/- 2.28741e-11
	- galCoeff: [0.7, 0.7, 0.7, 0.7]
	- isoCoeff: [6.05284, 3.10519, 0.877347, 0.103208] +/- [0.308163, 0.157711, 0.0575216, 0.023612]
	- exposure(cm2s): 358688000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------
In [23]:
ag.setOptionEnergybin(0)
ag.setOptions(galcoeff=[0.7])
affectedSources = ag.freeSources('name == "PKS1510-089"', "index", False, show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 20.0843
  * Spectrum type: PowerLaw
  * Free parameters: flux
  * Initial source parameters:
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:56:59.510740):
	- flux (ph/cm2s): 1.5715e-06 +/- 1.4153e-07
	- index : 2.07581 +/- 0.0710626
	- counts: 563.679 +/- 50.7665
	- upper limit(ph/cm2s): 1.80527e-06
	- ergLog(erg/cm2s): 2.5398e-10 +/- 2.28741e-11
	- galCoeff: [0.7, 0.7, 0.7, 0.7]
	- isoCoeff: [6.05284, 3.10519, 0.877347, 0.103208] +/- [0.308163, 0.157711, 0.0575216, 0.023612]
	- exposure(cm2s): 358688000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------

Light curve with default tmin and tmax¶

In [24]:
lightCurveData1 = ag.lightCurveMLE("PKS1510-089", binsize=86400)
Computing light curve bins..please wait.
286it [00:42,  6.73it/s]
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 33/33 [00:00<00:00, 21876.41it/s]
Temporal bin loop:   0%|          | 0/30 [00:00<?, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
0it [00:00, ?it/s]
1it [00:00,  5.60it/s]
2it [00:00,  6.06it/s]
3it [00:00,  6.47it/s]
4it [00:00,  6.55it/s]
5it [00:00,  6.66it/s]
6it [00:00,  6.78it/s]
7it [00:01,  6.79it/s]
8it [00:01,  6.87it/s]
9it [00:01,  7.02it/s]
10it [00:01,  7.12it/s]
11it [00:01,  7.23it/s]
12it [00:01,  7.43it/s]
13it [00:01,  7.28it/s]
14it [00:02,  7.08it/s]
15it [00:02,  6.53it/s]
16it [00:02,  6.54it/s]
17it [00:02,  6.24it/s]
19it [00:02,  6.90it/s]

100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 3/3 [00:00<00:00, 1609.69it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
In [25]:
cat $lightCurveData1
time_start_mjd time_end_mjd sqrt(ts) flux flux_err flux_ul gal gal_error iso iso_error l_peak b_peak dist_peak l b r ell_dist a b phi exposure ExpRatio counts counts_err Index Index_Err Par2 Par2_Err Par3 Par3_Err Erglog Erglog_Err Erglog_UL time_start_utc time_end_utc time_start_tt time_end_tt Fix index ULConfidenceLevel SrcLocConfLevel start_l start_b start_flux typefun par2 par3 galmode2 galmode2fit isomode2 isomode2fit edpcor fluxcor integratortype expratioEval expratio_minthr expratio_maxthr expratio_size Emin emax fovmin fovmax albedo binsize expstep phasecode fit_cts fit_fitstatus0 fit_fcn0 fit_edm0 fit_nvpar0 fit_nparx0 fit_iter0 fit_fitstatus1 fit_fcn1 fit_edm1 fit_nvpar1 fit_nparx1 fit_iter1 fit_Likelihood1
54891.0 54892.0 0 0 28.7323e-08 64.2527e-08 0.7 0 10.579 0.000710945 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 15327000.0 0 0 4.40381 2.07581 0 0 0 0 0 0 4.6436e-11 1.03843e-10 2009-03-01T00:00:00.000 2009-03-02T00:00:00.000 162950400.0 163036800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 171 0 358.338 1.4748e-15 1 8 29 0 358.338 3.68435e-10 2 8 29 514.089
54901.0 54902.0 4.05475 199.692e-08 71.3328e-08 363.143e-08 0.7 0 9.55789 0.38845 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 11588200.0 0 23.1407 8.26618 2.07581 0 0 0 0 0 3.22735e-10 1.15285e-10 5.86897e-10 2009-03-11T00:00:00.000 2009-03-12T00:00:00.000 163814400.0 163900800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 138 0 321.762 2.56683e-17 1 8 33 0 313.542 4.09373e-09 2 8 33 434.383
54902.0 54903.0 8.05892 418.403e-08 82.3936e-08 598.618e-08 0.7 0 10.5869 0.382312 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14133400.0 0 59.1346 11.645 2.07581 0 0 0 0 0 6.76207e-10 1.33161e-10 9.67464e-10 2009-03-12T00:00:00.000 2009-03-13T00:00:00.000 163900800.0 163987200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 212 0 409.112 5.83254e-16 1 8 31 0 376.638 6.19148e-09 2 8 31 553.274
54903.0 54904.0 7.79947 392.171e-08 81.0723e-08 570.322e-08 0.7 0 12.0593 0.358571 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 13751600.0 0 53.9297 11.1487 2.07581 0 0 0 0 0 6.33812e-10 1.31026e-10 9.21733e-10 2009-03-13T00:00:00.000 2009-03-14T00:00:00.000 163987200.0 164073600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 223 0 426.687 3.13691e-15 1 8 31 0 396.272 1.90977e-09 2 8 31 577.499
54904.0 54905.0 3.53833 142.073e-08 56.3876e-08 272.802e-08 0.7 0 9.8424 0.341388 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 13502200.0 0 19.1829 7.61356 2.07581 0 0 0 0 0 2.29613e-10 9.11315e-11 4.40893e-10 2009-03-14T00:00:00.000 2009-03-15T00:00:00.000 164073600.0 164160000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 159 0 333.967 2.42382e-15 1 8 31 0 327.707 1.2984e-06 2 8 31 481.161
54905.0 54906.0 3.56163 171.527e-08 65.8813e-08 322.337e-08 0.7 0 13.0983 0.354841 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 12904700.0 0 22.135 8.50177 2.07581 0 0 0 0 0 2.77216e-10 1.06475e-10 5.20949e-10 2009-03-15T00:00:00.000 2009-03-16T00:00:00.000 164160000.0 164246400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 196 0 383.827 1.56898e-14 1 8 43 0 377.484 5.15399e-10 2 8 43 555.802
54906.0 54907.0 2.62596 105.855e-08 55.3878e-08 238.8e-08 0.7 0 13.6628 0.312759 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 12575300.0 0 13.3116 6.96517 2.07581 0 0 0 0 0 1.71079e-10 8.95156e-11 3.85939e-10 2009-03-16T00:00:00.000 2009-03-17T00:00:00.000 164246400.0 164332800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 190 0 388.281 6.14462e-16 1 8 31 0 384.833 1.47379e-07 2 8 31 546.584
54907.0 54908.0 3.53062 141.11e-08 58.832e-08 280.348e-08 0.7 0 9.74173 0.335196 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 12076400.0 0 17.041 7.10479 2.07581 0 0 0 0 0 2.28057e-10 9.50819e-11 4.53087e-10 2009-03-17T00:00:00.000 2009-03-18T00:00:00.000 164332800.0 164419200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 141 0 319.023 6.71121e-16 1 8 31 0 312.79 1.19586e-06 2 8 31 444.086
54908.0 54909.0 5.21521 272.509e-08 76.7815e-08 444.481e-08 0.7 0 10.8178 0.371848 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 12029000.0 0 32.7799 9.23601 2.07581 0 0 0 0 0 4.40418e-10 1.24091e-10 7.18353e-10 2009-03-18T00:00:00.000 2009-03-19T00:00:00.000 164419200.0 164505600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 165 0 362.353 6.56599e-17 1 8 31 0 348.754 1.71679e-09 2 8 31 489.664
54909.0 54910.0 6.78882 420.475e-08 96.1237e-08 632.567e-08 0.7 0 11.7296 0.421565 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 11751700.0 0 49.413 11.2962 2.07581 0 0 0 0 0 6.79555e-10 1.55352e-10 1.02233e-09 2009-03-19T00:00:00.000 2009-03-20T00:00:00.000 164505600.0 164592000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 187 0 389.913 1.23565e-15 1 8 31 0 366.869 5.42187e-08 2 8 31 524.233
54910.0 54911.0 5.45734 277.593e-08 78.1827e-08 453.512e-08 0.7 0 7.7002 0.389604 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 10991400.0 0 30.5114 8.5934 2.07581 0 0 0 0 0 4.48635e-10 1.26356e-10 7.32949e-10 2009-03-20T00:00:00.000 2009-03-21T00:00:00.000 164592000.0 164678400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 119 0 294.622 3.78945e-16 1 8 31 0 279.731 8.51812e-09 2 8 31 384.345
54892.0 54893.0 0.118113 3.84498e-08 70.6507e-08 87.39e-08 0.7 0 10.1439 0.32329 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14312700.0 0 0.550319 10.112 2.07581 0 0 0 0 0 6.21411e-12 1.14183e-10 1.41236e-10 2009-03-02T00:00:00.000 2009-03-03T00:00:00.000 163036800.0 163123200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 154 0 338.839 1.75516e-17 1 8 31 0 338.832 9.00311e-08 2 8 31 478.446
54911.0 54912.0 1.73495 98.8951e-08 68.2961e-08 260.975e-08 0.7 0 13.5204 0.399269 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 11005500.0 0 10.8839 7.51635 2.07581 0 0 0 0 0 1.59831e-10 1.10378e-10 4.21778e-10 2009-03-21T00:00:00.000 2009-03-22T00:00:00.000 164678400.0 164764800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 160 0 343.754 3.76766e-16 1 8 31 0 342.249 7.00372e-07 2 8 31 491.159
54912.0 54913.0 3.95151 208.503e-08 77.093e-08 386.305e-08 0.7 0 11.8088 0.36921 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 10417300.0 0 21.7204 8.03099 2.07581 0 0 0 0 0 3.36975e-10 1.24595e-10 6.24331e-10 2009-03-22T00:00:00.000 2009-03-23T00:00:00.000 164764800.0 164851200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 147 0 329.361 9.64027e-16 1 8 33 0 321.554 1.91375e-09 2 8 33 456.077
54913.0 54914.0 2.71851 153.634e-08 74.8933e-08 329.501e-08 0.7 0 11.8345 0.382471 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 9797000.0 0 15.0516 7.3373 2.07581 0 0 0 0 0 2.48298e-10 1.2104e-10 5.32527e-10 2009-03-23T00:00:00.000 2009-03-24T00:00:00.000 164851200.0 164937600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 134 0 306.537 3.64352e-16 1 8 31 0 302.842 1.37689e-06 2 8 31 429.911
54914.0 54915.0 4.25687 291.499e-08 98.7378e-08 516.107e-08 0.7 0 12.911 0.379494 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 8249200.0 0.416419 24.0463 8.14508 2.07581 0 0 0 0 0 4.71109e-10 1.59576e-10 8.34113e-10 2009-03-24T00:00:00.000 2009-03-25T00:00:00.000 164937600.0 165024000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 132 0 307.296 6.21232e-16 1 8 31 0 298.236 1.56305e-09 2 8 31 420.532
54915.0 54916.0 7.08957 542.765e-08 119.967e-08 807.787e-08 0.7 0 14.2141 0.377602 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 8602120.0 1.36823 46.6893 10.3197 2.07581 0 0 0 0 0 8.77195e-10 1.93886e-10 1.30551e-09 2009-03-25T00:00:00.000 2009-03-26T00:00:00.000 165024000.0 165110400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 170 0 365.063 1.6272e-15 1 8 33 0 339.932 1.15734e-09 2 8 33 480.296
54916.0 54917.0 3.29673 250.002e-08 102.731e-08 486.709e-08 0.7 0 14.2527 0.392233 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 7576590.0 3.15289 18.9416 7.78352 2.07581 0 0 0 0 0 4.04043e-10 1.6603e-10 7.866e-10 2009-03-26T00:00:00.000 2009-03-27T00:00:00.000 165110400.0 165196800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 129 0 306.064 1.81549e-15 1 8 31 0 300.63 2.03001e-09 2 8 31 418.88
54917.0 54918.0 5.56527 353.984e-08 104.723e-08 593.889e-08 0.7 0 11.7234 0.357934 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 7598640.0 5.71089 26.898 7.95752 2.07581 0 0 0 0 0 5.72095e-10 1.69249e-10 9.5982e-10 2009-03-27T00:00:00.000 2009-03-28T00:00:00.000 165196800.0 165283200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 118 0 281.796 9.37738e-18 1 8 31 0 266.31 5.19027e-09 2 8 31 372.173
54918.0 54919.0 2.48055 178.839e-08 95.1463e-08 403.591e-08 0.7 0 12.6568 0.39534 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 7073150.0 9.04224 12.6496 6.72984 2.07581 0 0 0 0 0 2.89033e-10 1.53772e-10 6.52269e-10 2009-03-28T00:00:00.000 2009-03-29T00:00:00.000 165283200.0 165369600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 105 0 255.621 9.33425e-17 1 8 31 0 252.544 1.81754e-06 2 8 31 351.999
54919.0 54920.0 2.84601 207.07e-08 103.984e-08 456.135e-08 0.7 0 13.5054 0.396636 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 6560830.0 13.4444 13.5855 6.82224 2.07581 0 0 0 0 0 3.34659e-10 1.68056e-10 7.37188e-10 2009-03-29T00:00:00.000 2009-03-30T00:00:00.000 165369600.0 165456000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 103 0 257.922 6.0254e-13 1 8 41 0 253.872 3.85517e-09 2 8 41 348.554
54920.0 54921.0 3.29772 296.591e-08 119.671e-08 569.535e-08 0.7 0 10.9698 0.420335 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 5944460.0 18.4414 17.6308 7.11382 2.07581 0 0 0 0 0 4.79339e-10 1.93408e-10 9.20461e-10 2009-03-30T00:00:00.000 2009-03-31T00:00:00.000 165456000.0 165542400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 84 0 218.914 1.0323e-15 1 8 31 0 213.476 1.5114e-09 2 8 31 290.545
54893.0 54894.0 0 2164930.0e-08 12.1116e-08 41.7098e-08 0.7 0 9.49691 0.000335465 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 15551300.0 0 3.36676e-07 1.88351 2.07581 0 0 0 0 0 3.49888e-18 1.95742e-11 6.74098e-11 2009-03-03T00:00:00.000 2009-03-04T00:00:00.000 163123200.0 163209600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 157 0 341.689 1.8718e-15 1 8 29 0 341.689 2.1625e-07 2 8 29 483.44
54894.0 54895.0 0 66613400.0e-08 11.6348e-08 41.1691e-08 0.7 0 10.3525 0.000274007 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 15420500.0 0 1.02721e-07 1.79414 2.07581 0 0 0 0 0 1.07658e-18 1.88037e-11 6.65359e-11 2009-03-04T00:00:00.000 2009-03-05T00:00:00.000 163209600.0 163296000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 169 0 356.894 1.81457e-16 1 8 29 0 356.894 7.79199e-08 2 8 29 507.349
54895.0 54896.0 0 0 9.50689e-08 36.1909e-08 0.7 0 11.4341 0.000168804 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 15443200.0 0 0 1.46817 2.07581 0 0 0 0 0 0 1.53647e-11 5.84904e-11 2009-03-05T00:00:00.000 2009-03-06T00:00:00.000 163296000.0 163382400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 183 0 366.795 1.65408e-13 1 8 34 0 366.795 6.42006e-10 2 8 34 535.932
54896.0 54897.0 3.33681 121.106e-08 50.3652e-08 238.029e-08 0.7 0 9.34004 0.344954 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 15304200.0 0 18.5344 7.70802 2.07581 0 0 0 0 0 1.95727e-10 8.13984e-11 3.84693e-10 2009-03-06T00:00:00.000 2009-03-07T00:00:00.000 163382400.0 163468800.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 168 0 355.777 7.42986e-18 1 8 31 0 350.21 7.99511e-07 2 8 31 501.914
54897.0 54898.0 2.55625 103.173e-08 51.5861e-08 222.84e-08 0.7 0 9.29838 0.374695 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14884700.0 0 15.357 7.67845 2.07581 0 0 0 0 0 1.66744e-10 8.33715e-11 3.60145e-10 2009-03-07T00:00:00.000 2009-03-08T00:00:00.000 163468800.0 163555200.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 161 0 347.877 2.29795e-16 1 8 31 0 344.61 4.66805e-07 2 8 31 488.975
54898.0 54899.0 3.25476 133.796e-08 56.0206e-08 263.023e-08 0.7 0 8.99849 0.389475 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14806500.0 0 19.8105 8.29469 2.07581 0 0 0 0 0 2.16236e-10 9.05384e-11 4.25087e-10 2009-03-08T00:00:00.000 2009-03-09T00:00:00.000 163555200.0 163641600.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 160 0 338.051 1.25094e-16 1 8 31 0 332.754 1.8652e-06 2 8 31 484.437
54899.0 54900.0 4.08314 197.254e-08 67.691e-08 349.57e-08 0.7 0 10.5119 0.420623 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14619600.0 0 28.8378 9.89618 2.07581 0 0 0 0 0 3.18794e-10 1.094e-10 5.64962e-10 2009-03-09T00:00:00.000 2009-03-10T00:00:00.000 163641600.0 163728000.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 187 0 378.461 4.03457e-15 1 8 33 0 370.125 1.20358e-08 2 8 33 534.944
54900.0 54901.0 4.00797 169.194e-08 60.0694e-08 306.391e-08 0.7 0 10.2938 0.368614 351.29 40.13 0 -1 -1 -1 -1 -1 -1 -1 14455400.0 0 24.4576 8.68328 2.07581 0 0 0 0 0 2.73444e-10 9.70818e-11 4.95177e-10 2009-03-10T00:00:00.000 2009-03-11T00:00:00.000 163728000.0 163814400.0 1 2.07581 2 5.99147 351.29 40.13 1.5715e-06 0 0 0 0 0 0 0 1 0 1 1 0 15 10 100 10000 0 60 80 0.5 0 18 178 0 375.542 8.79846e-19 1 8 40 0 367.51 1.42583e-09 2 8 40 520.557
In [26]:
print(lightCurveData1)
/root/agilepy_analysis/username_PKS1510-089_20221022-184850/lc/0/light_curve_162950400_165542400.txt
In [27]:
ag.displayLightCurve("mle")
Generating plot..please wait.
In [28]:
affectedSources = ag.freeSources('name == "PKS1510-089"', "pos", True, show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 20.0843
  * Spectrum type: PowerLaw
  * Free parameters: flux pos
  * Initial source parameters:
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:56:59.510740):
	- flux (ph/cm2s): 1.5715e-06 +/- 1.4153e-07
	- index : 2.07581 +/- 0.0710626
	- counts: 563.679 +/- 50.7665
	- upper limit(ph/cm2s): 1.80527e-06
	- ergLog(erg/cm2s): 2.5398e-10 +/- 2.28741e-11
	- galCoeff: [0.7, 0.7, 0.7, 0.7]
	- isoCoeff: [6.05284, 3.10519, 0.877347, 0.103208] +/- [0.308163, 0.157711, 0.0575216, 0.023612]
	- exposure(cm2s): 358688000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------
In [29]:
lightCurveData1 = ag.lightCurveMLE("PKS1510-089", binsize=86400)
Computing light curve bins..please wait.
286it [00:41,  6.97it/s]
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 33/33 [00:00<00:00, 19141.48it/s]
Temporal bin loop:   0%|          | 0/30 [00:00<?, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
0it [00:00, ?it/s]
1it [00:00,  4.19it/s]
2it [00:00,  5.40it/s]
3it [00:00,  5.90it/s]
4it [00:00,  6.26it/s]
5it [00:00,  6.22it/s]
6it [00:00,  6.36it/s]
7it [00:01,  6.62it/s]
8it [00:01,  6.88it/s]
9it [00:01,  6.54it/s]
10it [00:01,  6.31it/s]
11it [00:01,  6.14it/s]
12it [00:01,  6.09it/s]
13it [00:02,  6.12it/s]
14it [00:02,  6.20it/s]
15it [00:02,  6.37it/s]
16it [00:02,  6.44it/s]
17it [00:02,  6.60it/s]
19it [00:02,  6.56it/s]

100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 3/3 [00:00<00:00, 2003.33it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
In [30]:
ag.displayLightCurve("mle")
Generating plot..please wait.
In [31]:
lightCurveData4 = ag.lightCurveMLE("PKS1510-089", binsize=86400*4)
Computing light curve bins..please wait.
286it [00:42,  6.77it/s]
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 33/33 [00:00<00:00, 22780.12it/s]
Temporal bin loop:   0%|          | 0/8 [00:00<?, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
2022-10-22 19:17:57,087 [WARNING ] [AGAnalysis] [LC] The last bin [165369600, 165715200] of the light curve analysis falls outside the range of the available data [.. , 165585534.0]. The binsize is reduced to 215934.0 seconds, the new bin is [165369600, 165585534.0]
0it [00:00, ?it/s]
1it [00:00,  5.93it/s]
2it [00:00,  6.38it/s]
3it [00:00,  6.85it/s]
4it [00:00,  6.98it/s]
5it [00:00,  7.13it/s]
6it [00:00,  6.94it/s]
7it [00:01,  6.71it/s]
8it [00:01,  6.49it/s]
9it [00:01,  6.52it/s]
10it [00:01,  6.50it/s]
11it [00:01,  6.62it/s]
12it [00:01,  6.73it/s]
13it [00:01,  6.80it/s]
14it [00:02,  6.85it/s]
15it [00:02,  6.30it/s]
16it [00:02,  5.97it/s]
17it [00:02,  5.89it/s]
18it [00:02,  5.93it/s]
19it [00:02,  6.05it/s]
20it [00:03,  6.19it/s]
21it [00:03,  6.35it/s]
22it [00:03,  6.52it/s]
23it [00:03,  6.82it/s]
24it [00:03,  6.98it/s]
25it [00:03,  7.08it/s]
26it [00:03,  6.81it/s]
27it [00:04,  6.61it/s]

100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 4/4 [00:00<00:00, 1951.75it/s]
Generating maps..please wait.
Fov bins loop:   0%|          | 0/1 [00:00<?, ?it/s]
Energy bins loop: 0it [00:00, ?it/s]
In [32]:
ag.displayLightCurve("mle")
Generating plot..please wait.
In [33]:
ag.selectSources('name == "PKS1510-089"', show=True)
-----------------------------------------------------------
 Source name: PKS1510-089 (PointSource) => sqrt(ts): 20.0843
  * Spectrum type: PowerLaw
  * Free parameters: flux pos
  * Initial source parameters:
	- flux (ph/cm2s): 1.7293e-06 +/- 1.2311e-07
	- index : 2.1
	- Source position (l,b): (351.29, 40.13)
	- Distance from map center (deg): 0.0
  * Last MLE analysis (2022-10-22 18:56:59.510740):
	- flux (ph/cm2s): 1.5715e-06 +/- 1.4153e-07
	- index : 2.07581 +/- 0.0710626
	- counts: 563.679 +/- 50.7665
	- upper limit(ph/cm2s): 1.80527e-06
	- ergLog(erg/cm2s): 2.5398e-10 +/- 2.28741e-11
	- galCoeff: [0.7, 0.7, 0.7, 0.7]
	- isoCoeff: [6.05284, 3.10519, 0.877347, 0.103208] +/- [0.308163, 0.157711, 0.0575216, 0.023612]
	- exposure(cm2s): 358688000.0
	- exp-ratio: 0
	- L_peak: 351.29
	- B_peak: 40.13
	- Distance from start pos: 0
	- ellipse position:
	    - L: -1
	    - B: -1
	    - Distance from start pos: -1
	    - radius of circle: -1
	    - ellipse:
		  - a: -1
		  - b: -1
		  - phi: -1
-----------------------------------------------------------
Out[33]:
[<agilepy.core.source.Source.PointSource at 0x7fc73a13ef70>]

time_start_mjd time_end_mjd sqrt(ts) flux flux_err flux_ul gal gal_error iso iso_error l_peak b_peak dist_peak l b r ell_dist a b phi exposure ExpRatio counts counts_err Index Index_Err Par2 Par2_Err Par3 Par3_Err Erglog Erglog_Err Erglog_UL time_start_utc time_end_utc time_start_tt time_end_tt Fix index ULConfidenceLevel SrcLocConfLevel start_l start_b start_flux typefun par2 par3 galmode2 galmode2fit isomode2 isomode2fit edpcor fluxcor integratortype expratioEval expratio_minthr expratio_maxthr expratio_size Emin emax fovmin fovmax albedo binsize expstep phasecode

In [34]:
ag.displayGenericColumns(lightCurveData4, ["ExpRatio"], um=[""])

Cleaning up¶

In [35]:
ag.deleteAnalysisDir()
Out[35]:
True

¶

¶