Difference between revisions of "Dragonfly Pulsar analysis CristianPozo"

Latest revision as of 10:40, 10 July 2024

General information[edit]

Name of the source: Dragonfly Pulsar (PSRJ2021+3651)
Brief description of the source:

- Object type : Pulsar

- RA, Dec in deg (ICRS): 305.27252, 36.85126

People involved in this analysis (LST-1 and new Fermi-LAT analysis)[edit]

Cristian Pozo González
Alvaro Mas-Aguilar
Ruben Lopez-Coto
Daniel Morcuende

Run selection[edit]

Run selection with the use of the notebook from Abelardo Moralejo [1]
Atmospheric transmission extracted from ELOG
In total 48 runs selected, obstime 13.51 hours

Summary of selection cuts:

- Zenith < 70 deg and Zenith < 20 deg

- elapsed_time > 2 min

- transmission_cut > 0.60

- pedestal charge std dev < 1.9 p.e.

- Cosmic rate > 3000 ev/s

- Cosmic rate( > 10 p.e.) > 25-15 ev/s (depends on zd)

- Cosmic rate( > 30 p.e.) > 3-4 ev/s (depends on zd)

Number of runs with zenith < 70 deg (% is w.r.t. those in Sky region & zenith range):

 + In the requested Sky region and range of dates:	 132
 + zenith in requested range:				 132
 + NSB in requested range:				 131 (99.2%)
 + FF and pedestal interleaved events are present:	 130 (98.5%)
 + Stable pointing:					 129 (97.7%)
 + dR/dI fit P-value ok:				 128 (97.0%)
 + dR/dI LS periodogram ok:				 127 (96.2%)
 + dR/dI index ok:					 121 (91.7%)
 + dR/dI rate ok:					 48 (36.4%)
 + intensity threshold ok:				 48 (36.4%)

Note: about 64% of all *dark-night* observations within ZD<80 deg fulfill all quality cuts. (in the stable, good-quality period 2022-11-18 - 2023-02-14, 92% of *dark-night* observations within ZD<80 deg do).

List of selected runs (Zenith < 70 deg):

1  :  2022-06-28  :  [8948 8949 8950 8951 8952 8953 8954 8955 8956 8957 8958 8959]
2  :  2022-06-29  :  [8980 8982 8983 8985 8986 8987]
3  :  2023-05-16  :  [13062 13063 13064 13066 13067]
4  :  2023-05-18  :  [13096 13097 13098 13099]
5  :  2023-06-22  :  [13519 13520 13521 13522]
6  :  2023-06-25  :  [13588 13589 13590 13591]
7  :  2023-10-05  :  [14910]
8  :  2023-10-06  :  [14918 14919 14920 14921 14922 14923 14924]
9  :  2023-11-13  :  [15489 15490 15491 15495 15497]

Number of runs with zenith < 20 deg (% is w.r.t. those in Sky region & zenith range):

 + In the requested Sky region and range of dates:	 132
 + zenith in requested range:				 77
 + NSB in requested range:				 76 (98.7%)
 + FF and pedestal interleaved events are present:	 75 (97.4%)
 + Stable pointing:					 74 (96.1%)
 + dR/dI fit P-value ok:				 73 (94.8%)
 + dR/dI LS periodogram ok:				 73 (94.8%)
 + dR/dI index ok:					 72 (93.5%)
 + dR/dI rate ok:					 31 (40.3%)
 + intensity threshold ok:				 31 (40.3%)

Note: about 64% of all *dark-night* observations within ZD<80 deg fulfill all quality cuts. (in the stable, good-quality period 2022-11-18 - 2023-02-14, 92% of *dark-night* observations within ZD<80 deg do).

List of selected runs (Zenith < 20 deg):

1  :  2022-06-28  :  [8951 8952 8953 8954 8955 8956 8957 8958 8959]
2  :  2022-06-29  :  [8980 8982 8983 8985 8986 8987]
3  :  2023-05-16  :  [13066 13067]
4  :  2023-05-18  :  [13098 13099]
5  :  2023-06-22  :  [13520 13521 13522]
6  :  2023-06-25  :  [13588 13589 13590 13591]
7  :  2023-10-05  :  [14910]
8  :  2023-10-06  :  [14918 14919 14920 14921]

DL1 data[edit]

- DL1a files produced by LSTOSA (lstchain v0.9.9)

- lstchain v0.9 tailcut8-4

- /fefs/aswg/data/real/DL1/{date}/v0.9/tailcut84/

DL2 data[edit]

- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL2

IRFs[edit]

- Point-like

- IRFs for each MC node and for each intensity cut

- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/IRF

DL3 data[edit]

- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL3

DL3 pulsar data[edit]

- Add phases using PINT-pulsar tools for LST-1 (https://github.com/alvmas/PulsarTimingAnalysis)

- Ephemeris file: /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021+3651/J2021+3651_58685_60244_cube.par

- DL3 pulsar data: /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL3_pulsar

High-level analysis[edit]

1D Spectral analysis[edit]

Performed with gammapy-v0.20.1

energy_geometry: #Geometry on true energy and reco energy for gammapy. True energy should be larger
 real:
   emin: 0.003
   emax: 1
   units: 'TeV'
   nbinning: 40
 reco:
   emin: 0.01
   emax: 0.5
   units: 'TeV'
   nbinning: 20

spectral_points_binning 
   e_min, e_max = 0.02, 0.7
   nbins = 8 
   logspace
   sqrt_ts>1.5

Phase Regions:

 Bkg: [0.73,1]
 P1: [0.105,0.13]
 P2: [0.57,0.615]
 P1+P2: [0.105,0.13,0.57,0.615]
 P3: None

Spectral fitting of stacked LST dataset
Performed for point-like assumption
Stacked analysis
LST-1 data alone: Power-law spectral model
Fermi + LST-1 data: Power-Law Spectral Model, ExponentialCutOff, LogParabola

Results[edit]

Phaseogram[edit]

Phaseogram with all data

Phaseogram with zd<20

Phaseogram in different energy bins (without zd cut):

Phaseogram in different energy bins (zd<20):

Evolution of the signal in time (without zd cut):

Evolution of the signal in time (zd<20):

P1 LST-1 + Fermi-LAT joint fit[edit]

- Using a PowerLaw model (without zd cut):

- Using a PowerLaw model (zd<20):

Difference between revisions of "Dragonfly Pulsar analysis CristianPozo"

Latest revision as of 10:40, 10 July 2024

Contents

General information[edit]

People involved in this analysis (LST-1 and new Fermi-LAT analysis)[edit]

Run selection[edit]

DL1 data[edit]

DL2 data[edit]

IRFs[edit]

DL3 data[edit]

DL3 pulsar data[edit]

High-level analysis[edit]

1D Spectral analysis[edit]

Results[edit]

Phaseogram[edit]

P1 LST-1 + Fermi-LAT joint fit[edit]

Navigation menu

Personal tools

Namespaces

Variants

Views

More

Search

Navigation

Tools

@@ Line 1: / Line 1: @@
 ==General information==
-* '''Name of the source''': Dragonfly Pulsar
+* '''Name of the source''': Dragonfly Pulsar (PSRJ2021+3651)
 * '''Brief description of the source''':
 :- Object type : Pulsar
-:- RA, Dec in deg (ICRS): 83.63, 22.01
+:- RA, Dec in deg (ICRS): 305.27252, 36.85126
 ==People involved in this analysis (LST-1 and new Fermi-LAT analysis)==
@@ Line 19: / Line 19: @@
 * '''Summary of selection cuts''':
-:- Zenith < 70 deg
+:- Zenith < 70 deg and Zenith < 20 deg
 :- elapsed_time > 2 min
 :- transmission_cut > 0.60
@@ Line 27: / Line 27: @@
 :- Cosmic rate( > 30 p.e.) > 3-4 ev/s (depends on zd)
+* '''Number of runs with zenith < 70 deg (% is w.r.t. those in Sky region & zenith range)''':
-* '''List of selected runs''':
+  + In the requested Sky region and range of dates:	 132
+  + zenith in requested range:				 132
+  + NSB in requested range:				 131 (99.2%)
+  + FF and pedestal interleaved events are present:	 130 (98.5%)
+  + Stable pointing:					 129 (97.7%)
+  + dR/dI fit P-value ok:				 128 (97.0%)
+  + dR/dI LS periodogram ok:				 127 (96.2%)
+  + dR/dI index ok:					 121 (91.7%)
+  + dR/dI rate ok:					 48 (36.4%)
+  + intensity threshold ok:				 48 (36.4%)
+Note: about 64% of all *dark-night* observations within ZD<80 deg fulfill all quality cuts.
+(in the stable, good-quality period 2022-11-18 - 2023-02-14, 92% of *dark-night* observations within ZD<80 deg do).
+* '''List of selected runs (Zenith < 70 deg)''':
   :  2022-06-28  :  [8948 8949 8950 8951 8952 8953 8954 8955 8956 8957 8958 8959]
@@ Line 39: / Line 54: @@
   :  2023-10-06  :  [14918 14919 14920 14921 14922 14923 14924]
   :  2023-11-13  :  [15489 15490 15491 15495 15497]
-== MC information ==
+* '''Number of runs with zenith < 20 deg (% is w.r.t. those in Sky region & zenith range)''':
-:- 20221027_v0.9.9_crab_tuned
+  + In the requested Sky region and range of dates:	 132
-:- AllSky MC using different test nodes for different runs depending on sky location.
+  + zenith in requested range:				 77
+  + NSB in requested range:				 76 (98.7%)
+  + FF and pedestal interleaved events are present:	 75 (97.4%)
+  + Stable pointing:					 74 (96.1%)
+  + dR/dI fit P-value ok:				 73 (94.8%)
+  + dR/dI LS periodogram ok:				 73 (94.8%)
+  + dR/dI index ok:					 72 (93.5%)
+  + dR/dI rate ok:					 31 (40.3%)
+  + intensity threshold ok:				 31 (40.3%)
-* '''PSF and NSB tuned''':
+Note: about 64% of all *dark-night* observations within ZD<80 deg fulfill all quality cuts.
-  "image_modifier": {
+(in the stable, good-quality period 2022-11-18 - 2023-02-14, 92% of *dark-night* observations within ZD<80 deg do).
-     "increase_nsb": true,
-     "extra_noise_in_dim_pixels": 1.27,
-     "extra_bias_in_dim_pixels": 0.665,
-     "transition_charge": 8,
-     "extra_noise_in_bright_pixels": 2.08,
-     "increase_psf": false,
-     "smeared_light_fraction": 0
-     },
-* '''Tuned DL1b MC''':
+* '''List of selected runs (Zenith < 20 deg)''':
-:- Training sample: /fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/TrainingDataset/
-:- Testing sample (AllSky MC): /fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/TestingDataset/
-* '''Random Forests - Source Dependent''':
+  :  2022-06-28  :  [8951 8952 8953 8954 8955 8956 8957 8958 8959]
-:- General path: /fefs/aswg/workspace/seiya.nozaki/Crab_performance_paper/20221027_v0.9.9_crab_tuned/std/RF/trained_models/
+  :  2022-06-29  :  [8980 8982 8983 8985 8986 8987]
-:- Configuration file: /fefs/aswg/workspace/seiya.nozaki/Crab_performance_paper/20221027_v0.9.9_crab_tuned/std/DL2/lstchain_src_dep_config.json
+  :  2023-05-16  :  [13066 13067]
+  :  2023-05-18  :  [13098 13099]
-* '''DL2 MC (testing only)''':
+  :  2023-06-22  :  [13520 13521 13522]
-:- Souce-dependent
+  :  2023-06-25  :  [13588 13589 13590 13591]
-:- AllSky MC depending on the sky location.
+  :  2023-10-05  :  [14910]
-:- General path: /fefs/aswg/workspace/seiya.nozaki/Crab_performance_paper/20221027_v0.9.9_crab_tuned/std/DL2_MC/
+  :  2023-10-06  :  [14918 14919 14920 14921]
 == DL1 data ==
@@ Line 76: / Line 88: @@
 == DL2 data ==
-:- Data before August 2021: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL2_preAug/
+:- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL2
-:- Data after August 2021: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL2_posAug/
 ===IRFs===
 :- Point-like
 :- IRFs for each MC node and for each intensity cut
-:- IRF path before Aug2021: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/IRF/IRF_cut80/eff_0.7/
+:- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/IRF
-:- IRF path after Aug2021: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/IRF/IRF_cut50/eff_0.7/
-* IRF for different MC efficiency for extra checks: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/IRF/IRF_cut{int_cut}/eff_{eff}/
 == DL3 data ==
-:- General path for standard MC efficiency: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL3/eff_0.7/
+:- /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL3
-:- General path for different MC efficiency: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL3/eff_{eff}/
-:- Distribution of runs in AllSKy MC test nodes:
-* mc_node = [10,102.199]
-  [2929,2930,2931,2950,2970,2971,2972,2973,2988,3095,3096,3273,3274,3275,3318,3319,3320,3338,9438,10001,
-,10880,10909,10981,11087,11120,11121,11374,11463,11464,11491,11492,11576,11577,11625,11626]
-* mc_node = [10,248.117]
-  [2914,2932,2933,2934,2974,2975,2976,2989,2990,2991,3004,3005,3006,3276,3277,3278,3279,3321,3328,3329,
-,3339,3340,3355,3356,3598,3632,3672,3673,3674,3675,6872,6873,6874,6892,6893,7133,7161,7195,7196,
-,7227,7228,7253,7274,10881,10882,10910,10911,10946,10947,10948,10982,10983,11014,11088,11089,11122,
-,11402,11465,11466,11467,11493,11578,11579,11627,11628,11667]
-* mc_node = [14.984,175.158]
-  [9436,9437,9907,9998,9999,10883,10912,10913,10949,10950,10984,10985,11015,11090,11091,11118,11119,11124,
-,11373,11432,11433,11468,11469,11489,11490,11575,11630,11631,11668,11669,11681,11682,11699]
-* mc_node = [23.63,100.758]
-  [2949,2968,2969,3093,3094,3271,3272,6304, 9434, 9435, 9596, 9905, 9996,9997, 11117,11430,11431,11487,11488]
-* mc_node = [23.63,259.265]
-  [2977,2992,3007,3008,3231,3243,3373,3599,3600,3615,3633,3634,3676,3677,3706,3707,4067,4086,6875,6894,7097,
-,7099,7199,7231,7254,7255,7256,7275,10914,10915,10986,10987,11016,11017,11092,11093,11470,11471,11560,
-,11562,11603,11604,11670,11671,11683,11700,11701,11738,11739,11763]
-* mc_node = [32.059,102.217]
-  [2925,2966,2967,3090,3268,3269,3270,6016,6045,6072,6073,6302,9903,9904,11428,11429,11486]
-* mc_node=[32.059,248.099]
-  [2959,3009,3010,3169,3232,3233,3244,3245,3601,3635,3636,3708,4068,4069,4087,4088,6895,6896,7100,7101,7136,7137,
-,7201,7232,7233,7234,7276,7277,7278,7301,7302,10916,10917,10918,10988,10989,11018,11019,11094,11095,11472,
-,11563,11605,11606,11673,11674,11684,11685,11702,11703,11740,11741,11742,11764,11765,11766,11923]
-* mc_node = [43.197,87.604]
-  [2923,2924,2965,3088,3089,3266,3267,6043,6044,6071,6279,6301,6070,9274,11427]
-* mc_node = [43.197,262.712]
-  [2919,2961,3011,3012,3143,3170,3171,3234,3382,3383,3384,3637,3638,3639,3711,3925,4007,4008,4031,4032,4070,4071,
-,4090,6897,7102,7138,7139,7168,7169,7235,7236,7237,7238,7279,7280,7303,7304,7465,10919,11510,11511,11675,
-,11743,11744,11767,11768,11828,11829,11890,11891,11926]
-* mc_node = [52.374,110.312]
-  [2759,3264,3265,2922,3087,6013,6014,6041,6042,6069]
-* mc_node = [52.374,301.217]
-  [3235,3236,3385,3712,3961,3980,4033,4034,4091,4125,4126,4073,4092,6898,7105,7106,7140,7170,7171,7281,7305,11512,11513,
-,11709,11745,11746,11769,11770,11830,11831,11892,11893,11927,11928]
-* mc_node = [60.528,99.126]
-  [6012]
-* mc_node = [60.528,251.19]
-  [3237,4074,7172]
-* Fixed cuts (before August 2021):
- "intensity": [80, Infinity],
- "width": [0, Infinity],
- "length": [0, Infinity],
- "r": [0, 1],
- "wl": [0.1, 1],
- "leakage_intensity_width_2": [0, 1.0],
- "event_type": [32, 32]
-* Fixed cuts (after August 2021):
- "intensity": [50, Infinity],
- "width": [0, Infinity],
- "length": [0, Infinity],
- "r": [0, 1],
- "wl": [0.1, 1],
- "leakage_intensity_width_2": [0, 1.0],
- "event_type": [32, 32]
-* Global gammaness and alpha cuts optimized on MC gamma efficiency of 70%, point-like source assumption
- gh-efficiency 0.7
- alpha-containment 0.7
 == DL3 pulsar data ==
 :- Add phases using PINT-pulsar tools for LST-1 (https://github.com/alvmas/PulsarTimingAnalysis)
-:- Ephemeris file: /home/alvaro.mas/software/Pulsar_analysis/data/all.gro
+:- Ephemeris file: /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021+3651/J2021+3651_58685_60244_cube.par
-:- General path for standard MC efficiency: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL3_pulsar/eff_0.7/
+:- DL3 pulsar data: /fefs/aswg/workspace/cristian.pozo/data/PSRJ2021_3651/DL3_pulsar
-:- General path for different MC efficiency: /fefs/aswg/workspace/alvaro.mas/data/real/Crab/v0.9.9_final/DL3_pulsar/eff_{eff}/
 ==High-level analysis==
@@ Line 178: / Line 109: @@
 * Performed with gammapy-v0.20.1
-  energy_axis = MapAxis.from_energy_bounds(
+  energy_geometry: #Geometry on true energy and reco energy for gammapy. True energy should be larger
-.01, 10, nbin=40, per_decade=False, unit="TeV", name="energy"
+  real:
- )
+     emin: 0.003
- energy_axis_true = MapAxis.from_energy_bounds(
+    emax: 1
-.05, 100, nbin=100, per_decade=False, unit="TeV", name="energy_true"
+    units: 'TeV'
- )
+    nbinning: 40
+  reco:
+    emin: 0.01
+     emax: 0.5
+    units: 'TeV'
+    nbinning: 20
   spectral_points_binning
@@ Line 193: / Line 129: @@
 * Phase Regions:
-   Bkg: [0.52,0.87]
+   Bkg: [0.73,1]
-   P1: [0,0.026,0.983,1]
+   P1: [0.105,0.13]
-   P2: [0.377,0.422]
+   P2: [0.57,0.615]
-   P1+P2: [0,0.026,0.377,0.422, 0.983,1]
+   P1+P2: [0.105,0.13,0.57,0.615]
-   P3: [0.08,0.24]
+   P3: None
 * Spectral fitting of stacked LST dataset
@@ Line 211: / Line 147: @@
 * Phaseogram with all data
+[[File:phaseograma_all_zenit.jpg|1000px]]
-[[File:phaseogram_all_0.7.png|1000px]]
+* Phaseogram with zd<20
+[[File:phaseograme_20_zd.jpg|1000px]]
-* Phaseogram with zd<35
-[[File:phaseogram_low50.png|1000px]]
-* Phaseogram with zd>35
+* Phaseogram in different energy bins (without zd cut):
-[[File:phaseogram_above50.png|1000px]]
+[[File:energy_bin_1_all_zenit.jpg|500px]]
+[[File:energy_bin_2_all_zenit.jpg|500px]]
-* Phaseogram in different energy bins:
+* Phaseogram in different energy bins (zd<20):
-[[File:phaseogram_bin1_1.png|500px]]
+[[File:energy_bin_1_zd_20.jpg|500px]]
-[[File:phaseogram_bin2_1.png|500px]]
+[[File:energy_bin_2_zd_20.jpg|500px]]
-[[File:phaseogram_bin3_1.png|500px]]
-[[File:phaseogram_bin4_1.png|500px]]
-[[File:phaseogram_bin5_1.png|500px]]
-[[File:phaseogram_bin6_1.png|500px]]
-[[File:phaseogram_bin7_1.png|500px]]
 * Evolution of the signal in time (without zd cut):
+[[File:time_vs_significance_all_zenit.jpg|1000px]]
-[[File:time_evolution_lst.png|1000px]]
+* Evolution of the signal in time (zd<20):
-* Fitting results
+[[File:time_vs_significace_zd_20.jpg|1000px]]
-* P1/P2 differential ratio (with Fermi-LAT). Until 200 GeV due to large uncertainties
-[[File:p1p2_differential_2.png|500px]]
-* P1/P2 integral ratio (with Fermi-LAT).
-[[File:p1p2_ratio_integral_2.png|500px]]
-=== LST-1 only fit and SED (zd<50) ===
-*P1
-[[File:SED_P1_1.png]]
-[[File:LST_fit_only_P1_fine_bin_fermi.png]]
-[[File:LST_fit_only_P1_fine_wide_fermi.png]]
-[[File:fitting_sed_p1.png|800px]]
-*P2
-[[File:SED_p2_1.png]]
-[[File:LST_fit_only_P2_fine_bin_fermi.png]]
-[[File:LST_fit_only_P2_wide_bin_fermi.png]]
-[[File:fitting_sed_p2.png|800px]]
 === P1 LST-1 + Fermi-LAT joint fit ===
-** Using a SubExponentialCutOffPL model:
+** Using a PowerLaw model (without zd cut):
+[[File:fit_fermi_all_zd.png|1000px]]
-[[File:P1_exponential_model_3.png|500px]]
-[[File:Exponential_result_fit_P1.png|500px]]
-'''Total stat: -2logL = 32.56'''
-** Using a SmoothBrokenPowerLaw model
-[[File:SmoothBrokenPL_P1_2.png|500px]]
-[[File:SmoothBPL_results_fit_P1.png|500px]]
-'''Total stat: -2logL = 24.44'''
-*P1 for E>10 GeV:
-** Using a LogParabola model
-[[File:P1_logParabolaModel_2.png|500px]]
-[[File:LogParabola_results_fit_P1.png|500px]]
-'''Total stat: -2logL = 7.34'''
-** Using a PowerLaw model
-[[File:P1_PLmodel_2.png|500px]]
-[[File:PL_results_fit_P1.png|500px]]
-'''Total stat: -2logL = 7.83'''
-=== P2 LST-1 + Fermi-LAT joint fit ===
-** Using a SubExponentialCutOffPL model:
-[[File:P2_Exponential_model.png|500px]]
-[[File:Exponential_results_fit_P2.png|500px]]
-'''Total stat: -2logL = 51.40'''
-** Using a SmoothBrokenPowerLaw model
-[[File:SmoothBrokenPL_P2.png|500px]]
-[[File:SmoothBPL_results_fit_P2.png|500px]]
-'''Total stat: -2logL = 34.74'''
-*P2 for E>10 GeV:
-** Using a LogParabola model
-[[File:P2_logParabolaModel.png|500px]]
-[[File:LogParabola_results_fit_P2.png|500px]]
-'''Total stat: -2logL = 17.85'''
-** Using a PowerLaw model
-[[File:P2_PLmodel.png|500px]]
-[[File:PL_results_fit_P2.png|500px]]
-'''Total stat: -2logL = 19.35'''
-=== LST-1+ Fermi-LAT energy calibration ===
-*P1 in left and P2 in right (Using SmoothBrokenPowerLaw model in all energy range)
-[[File:energy_factor_fermi_LST_P1_3.png|500px]]
-[[File:energy_factor_LST_fermi_P2_3.png|500px]]
-*P1 in left and P2 in right (Using PowerLaw model at E>10 GeV)
-[[File:energy_factor_fermi_LST_P1_PL_3.png|500px]]
-[[File:energy_factor_LST_fermi_P2_PL_3.png|500px]]
-'''In all the cases the energy factor that minimize the chi2 value is below 5%'''
-=== Lightcurve ===
-[[File:lightcurve.png|1000px]]
-[[File:lightcurve_together.png|1000px]]
-=== Systematics tests ===
-==== Phaseogram and SED for different samples ====
-*Phaseogram Before August 2021
-[[File:histo_preAug.png|600px]]
-'''sign(P1)/sqrt(Tobs) = 0.856 h^(-1/2)'''
-'''sign(P2)/sqrt(Tobs) = 1.059 h^(-1/2)'''
-*Phaseogram After August 2021
-[[File:histo_posAug.png|600px]]
-'''sign(P1)/sqrt(Tobs) = 1.063 h^(-1/2)'''
-'''sign(P2)/sqrt(Tobs) = 1.183 h^(-1/2)'''
-*SED for P1:
-[[File:spectra_P1_bothsamples.png|600px]]
-'''Difference in spectral index of ~10% and flux 60% (lacking some statistics for first sample)'''
-*SED for P2:
-[[File:spectra_P2_bothsamples.png|600px]]
-'''Difference in spectral index of ~2% and flux of 15%'''
-==== Phaseogram and SED with different efficiencies ====
-[[File:phaseogram_different_efficiencies.png|1000px]]
-* P1 SED:
-[[File:SED_different_efficiencies.png|1000px]]
-'''Relative maximum error of 3.07% in spectral index and 35.4% in flux.'''
-* P1 SED:
-[[File:SED_different_efficiencies_P2.png|1000px]]
-'''Relative maximum error of 7.36% in spectral index and 29.2% in flux.'''
-==== Effect on SED due to shift on MC true energy =====
-** P1:
-[[File:P1_energy_factor_SED_2.png|1000px]]
-'''Relative maximum error of 1.01% in spectral index and 24.9% in flux.'''
-**P2:
-[[File:energy_factor_LST_P2.png|1000px]]
-'''Relative maximum error of 0.45% in spectral index and 16.4% in flux.'''
-==== Final uncertainties ====
-**P1:
-[[File:relative_error_index_P1.png|300px]]
-[[File:relative_error_flux_P1.png|300px]]
-[[File:total_uncertainty_P1.png|500px]]
-**P2:
+** Using a PowerLaw model (zd<20):
-[[File:relative_error_index_P2.png|300px]]
+[[File:fit_fermi_zd_20.jpg|1000px]]
-[[File:relative_error_flux_P2.png|300px]]
-[[File:total_uncertainty_P2.png|500px]]