Reference paper: https://www.aanda.org/articles/aa/pdf/2011/05/aa16128-10.pdf
from agilepy.api.AGAnalysis import AGAnalysis
# 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' | ------------------------------------------------------------------
confFilePath = "$HOME/agilepy_conf.yaml"
Set the configuration and download the data from SSDC data server
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"
)
ag = AGAnalysis(confFilePath)
ag.setOptions(binsize=0.5)
#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 = 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)
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
Freeable params = [flux, index, index1, index2, cutoffEnergy, pivotEnergy, curvature, pos]
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.
#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]
#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]
[None]
#Fix the diffuse emission coefficient to a default value.
ag.setOptions(galcoeff=[0.7])
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}
#Perform a maximum likelihood estimator
ag.mle()
['/root/agilepy_analysis/username_PKS1510-089_20221022-184850/mle/0/analysis_product_PKS1510-089.source']
#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 -----------------------------------------------------------
[<agilepy.core.source.Source.PointSource at 0x7fc73a13ef70>]
#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])
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]
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]
[None]
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 -----------------------------------------------------------
ag.mle()
['/root/agilepy_analysis/username_PKS1510-089_20221022-184850/mle/1/analysis_product_PKS1510-089.source']
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 -----------------------------------------------------------
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 -----------------------------------------------------------
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]
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
print(lightCurveData1)
/root/agilepy_analysis/username_PKS1510-089_20221022-184850/lc/0/light_curve_162950400_165542400.txt
ag.displayLightCurve("mle")
Generating plot..please wait.
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 -----------------------------------------------------------
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]
ag.displayLightCurve("mle")
Generating plot..please wait.
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]
ag.displayLightCurve("mle")
Generating plot..please wait.
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 -----------------------------------------------------------
[<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
ag.displayGenericColumns(lightCurveData4, ["ExpRatio"], um=[""])
ag.deleteAnalysisDir()
True