Difference between revisions of "NGC1068 LST1 analysis Cycle I DanielMorcuende"

From my_wiki
Jump to: navigation, search
(Spectral results)
(Dataset after good quality selection)
 
(5 intermediate revisions by the same user not shown)
Line 7: Line 7:
 
** Object type : AGN, Seyfert II galaxy
 
** Object type : AGN, Seyfert II galaxy
 
** Redshift (z) : 0.0037
 
** Redshift (z) : 0.0037
** Other relevant information: neutrino excess found by IceCube from this source location
+
** Other relevant information: neutrino excess found by IceCube (https://doi.org/10.1126/science.abg3395)
 
** RA, Dec in deg (ICRS): 40.66962153, -0.0132943
 
** RA, Dec in deg (ICRS): 40.66962153, -0.0132943
  
Line 23: Line 23:
  
 
=== Dataset after good quality selection ===
 
=== Dataset after good quality selection ===
 +
 +
[[ File:Ngc1068_LST1_quality_selection_cycle_I_dmorcuende.png ]]
  
 
  Number of runs (% is w.r.t. those in Sky region & zenith range):
 
  Number of runs (% is w.r.t. those in Sky region & zenith range):
Line 59: Line 61:
 
  '2024-01-31': [16445, 16446, 16447, 16448],
 
  '2024-01-31': [16445, 16446, 16447, 16448],
 
  '2024-02-10': [16660, 16661, 16662, 16663],
 
  '2024-02-10': [16660, 16661, 16662, 16663],
 +
 +
[[File:Ngc1068_LST1_rel_light_yield_cycle_I_dmorcuende.png]]
  
 
==Analysis information==
 
==Analysis information==
Line 86: Line 90:
 
  /fefs/aswg/workspace/daniel.morcuende/data/real/DL2/NGC1068
 
  /fefs/aswg/workspace/daniel.morcuende/data/real/DL2/NGC1068
  
==DL3 data selection==
+
==IRF and DL3 data selection==
Files:  
+
 
 +
DL3 files:  
 
  /fefs/aswg/workspace/daniel.morcuende/data/real/DL3/NGC1068/20240131_allsky_v0.10.5_all_dec_base/dec_931/gh_eff_0.7_th_cont_0.7
 
  /fefs/aswg/workspace/daniel.morcuende/data/real/DL3/NGC1068/20240131_allsky_v0.10.5_all_dec_base/dec_931/gh_eff_0.7_th_cont_0.7
  
 +
* IRF linear interpolation
 
* intensity > 50 p.e.
 
* intensity > 50 p.e.
 
* gammaness energy-dependent cuts with 70% efficiency
 
* gammaness energy-dependent cuts with 70% efficiency
 
* theta energy-dependent cuts 70% efficiency
 
* theta energy-dependent cuts 70% efficiency
 
  
 
==High-level analysis==
 
==High-level analysis==
Line 103: Line 108:
  
 
=== Theta2 plot ===
 
=== Theta2 plot ===
 +
 +
No detection
 +
 +
[[ File:Ngc1068_LST1_theta2_LE_cycle_I_dmorcuende.png ]]
 +
 +
Energy range: 0.0 - 0.2 TeV
 +
Excess: -5324.0; Off: 450085
 +
Gamma rate: -6.0 events / minute
 +
Off rate: 510.9 events / minute
 +
alpha (backg normalization): 1.00
 +
Li & Ma Significance: 0.0 standard deviations
 +
 +
 +
[[ File:Ngc1068_LST1_theta2_HE_cycle_I_dmorcuende.png ]]
 +
 +
Energy range: 0.2 - 1000000.0 TeV
 +
Excess: 14.0; Off: 3798
 +
Gamma rate: 0.0 events / minute
 +
Off rate: 4.3 events / minute
 +
alpha (backg normalization): 0.33
 +
Li & Ma Significance: 0.3 standard deviations
  
 
=== Significance map ===
 
=== Significance map ===
Line 110: Line 136:
 
=== Spectral results ===
 
=== Spectral results ===
  
For comparison 125-h long MAGIC results (V. A. Acciari et al 2019 ApJ 883 135)
+
* Fermi-LAT 3FHL and 3FGL SED shown.
 +
** Spectral model plotted corresponds to the Fermi-LAT data, extended to VHE.
 +
* For comparison 125-h long MAGIC results (V. A. Acciari et al 2019 ApJ 883 135)
  
 
[[File:NGC1068_LST1_SED_ULs_cycle_I_dmorcuende.png]]
 
[[File:NGC1068_LST1_SED_ULs_cycle_I_dmorcuende.png]]
 +
 +
== Light curve ==
 +
[[ File:Ngc1068_LST1_LC_cycle_I_dmorcuende.png ]]

Latest revision as of 14:50, 24 July 2024

Analysis of LST-1 NGC 1068 data from observation cycle I[edit]

General information[edit]

  • Name of the source: NGC 1068
  • Brief description of the source:
    • Object type : AGN, Seyfert II galaxy
    • Redshift (z) : 0.0037
    • Other relevant information: neutrino excess found by IceCube (https://doi.org/10.1126/science.abg3395)
    • RA, Dec in deg (ICRS): 40.66962153, -0.0132943
  • Analysis by Name (Institute - email)
    • Daniel Morcuende (IAA-CSIC, dmorcuende@iaa.es)
    • José Manuel Gavilán (IAA-CSIC)

Data-taking information[edit]

There are a total of 124 runs from 2023-08-19 to 2024-02-10

  • Dates of data-taking:
    • Date (time of observation h)
      • ZD range (deg): 25-65

Dataset after good quality selection[edit]

Ngc1068 LST1 quality selection cycle I dmorcuende.png

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

    In the requested Sky region and range of dates:	 128
  + zenith in requested range:				 127
  + NSB in requested range:				 127 (100.0%) 

  + FF and pedestal interleaved events are present:	 125 (98.4%)
  + Stable pointing:					 123 (96.9%)
  + dR/dI fit P-value ok:				 123 (96.9%)
  + dR/dI LS periodogram ok:				 118 (92.9%)
  + dR/dI index ok:					 94 (74.0%)
  + dR/dI rate ok:					 60 (47.2%)
  + intensity threshold ok:				 60 (47.2%)
Total observation time: 15.46 h


Total number of runs: 60 (in 16 nights)
Total observation time: 15.46 hours 

'2023-09-22': [14696, 14697, 14698, 14699],
'2023-10-06': [14927, 14929],
'2023-11-03': [15256, 15257, 15258, 15259],
'2023-11-04': [15263],
'2023-11-10': [15387, 15388, 15389],
'2023-11-12': [15474, 15475, 15476, 15477],
'2023-11-16': [15615, 15616, 15618],
'2023-11-17': [15656, 15657, 15658, 15659, 15672, 15673, 15674, 15675],
'2023-11-18': [15706, 15707, 15708, 15709, 15710, 15713],
'2023-12-08': [15883, 15884, 15886, 15887, 15888],
'2023-12-09': [15917, 15918, 15919, 15920],
'2023-12-14': [16085, 16086, 16087],
'2023-12-17': [16195, 16198],
'2024-01-29': [16427, 16428, 16429],
'2024-01-31': [16445, 16446, 16447, 16448],
'2024-02-10': [16660, 16661, 16662, 16663],

Ngc1068 LST1 rel light yield cycle I dmorcuende.png

Analysis information[edit]

  • lstchain v0.10
  • source-independent analysis
  • Gammapy v1.1

Monte Carlo information[edit]

Closest declination line towards La Palma latitude (28.8 deg)

20240131_allsky_v0.10.5_all_dec_base/dec_931

DL1 data[edit]

Standard DL1b files produced by LSTOSA (lstchain v0.10, tailcut 8-4 with cleaning based on pedestal RMS, dynamic cleaning)

/fefs/aswg/data/real/DL1/YYYYMMDD/v0.10/tailcut84/dl1_LST-1.Run?????.h5

Random forest[edit]

Models:

/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931

DL2[edit]

Files:

/fefs/aswg/workspace/daniel.morcuende/data/real/DL2/NGC1068

IRF and DL3 data selection[edit]

DL3 files:

/fefs/aswg/workspace/daniel.morcuende/data/real/DL3/NGC1068/20240131_allsky_v0.10.5_all_dec_base/dec_931/gh_eff_0.7_th_cont_0.7
  • IRF linear interpolation
  • intensity > 50 p.e.
  • gammaness energy-dependent cuts with 70% efficiency
  • theta energy-dependent cuts 70% efficiency

High-level analysis[edit]

  • Science Tool: Gammapy 1.1
  • point-like IRF, 1D analysis

Analysis Results[edit]

Theta2 plot[edit]

No detection

Ngc1068 LST1 theta2 LE cycle I dmorcuende.png

Energy range: 0.0 - 0.2 TeV
Excess: -5324.0; Off: 450085
Gamma rate: -6.0 events / minute
Off rate: 510.9 events / minute
alpha (backg normalization): 1.00
Li & Ma Significance: 0.0 standard deviations


Ngc1068 LST1 theta2 HE cycle I dmorcuende.png

Energy range: 0.2 - 1000000.0 TeV
Excess: 14.0; Off: 3798
Gamma rate: 0.0 events / minute
Off rate: 4.3 events / minute
alpha (backg normalization): 0.33
Li & Ma Significance: 0.3 standard deviations

Significance map[edit]

Excess map[edit]

Spectral results[edit]

  • Fermi-LAT 3FHL and 3FGL SED shown.
    • Spectral model plotted corresponds to the Fermi-LAT data, extended to VHE.
  • For comparison 125-h long MAGIC results (V. A. Acciari et al 2019 ApJ 883 135)

NGC1068 LST1 SED ULs cycle I dmorcuende.png

Light curve[edit]

Ngc1068 LST1 LC cycle I dmorcuende.png