Difference between revisions of "Dragonfly Pulsar analysis CristianPozo"

From my_wiki
Jump to: navigation, search
 
(8 intermediate revisions by the same user not shown)
Line 109: Line 109:
 
* Performed with gammapy-v0.20.1
 
* 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
−
     0.01, 10, nbin=40, per_decade=False, unit="TeV", name="energy"
+
  real:
−
)
+
     emin: 0.003
−
energy_axis_true = MapAxis.from_energy_bounds(
+
    emax: 1
−
     0.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  
 
  spectral_points_binning  
Line 124: Line 129:
 
* Phase Regions:
 
* 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
 
* Spectral fitting of stacked LST dataset
Line 142: Line 147:
  
 
* Phaseogram with all data
 
* 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
+
* Phaseogram in different energy bins (without zd cut):
 +
[[File:energy_bin_1_all_zenit.jpg|500px]]
 +
[[File:energy_bin_2_all_zenit.jpg|500px]]
  
−
[[File:phaseogram_low50.png|1000px]]
+
* Phaseogram in different energy bins (zd<20):
−
 
+
[[File:energy_bin_1_zd_20.jpg|500px]]
−
* Phaseogram with zd>35
+
[[File:energy_bin_2_zd_20.jpg|500px]]
−
[[File:phaseogram_above50.png|1000px]]
 
−
 
 
−
* Phaseogram in different energy bins:
 
−
[[File:phaseogram_bin1_1.png|500px]]
 
−
[[File:phaseogram_bin2_1.png|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):
 
* 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 ===
 
=== 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]]
 

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

Phaseograma all zenit.jpg

  • Phaseogram with zd<20

Phaseograme 20 zd.jpg

  • Phaseogram in different energy bins (without zd cut):

Energy bin 1 all zenit.jpg Energy bin 2 all zenit.jpg

  • Phaseogram in different energy bins (zd<20):

Energy bin 1 zd 20.jpg Energy bin 2 zd 20.jpg

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

Time vs significance all zenit.jpg

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

Time vs significace zd 20.jpg


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

    • Using a PowerLaw model (without zd cut):

Fit fermi all zd.png

    • Using a PowerLaw model (zd<20):

Fit fermi zd 20.jpg

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Dragonfly_Pulsar_analysis_CristianPozo&oldid=10945"