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

From my_wiki
Jump to: navigation, search
(Created page with "=Analysis of LST-1 NGC 1068 data from observation cycle I= ==General information== * '''Name of the source''': NGC 1068 * '''Brief description of the source''': ** Object t...")
 
(→‎Dataset after good quality selection)
 
(15 intermediate revisions by the same user not shown)
Line 5: Line 5:
 
* '''Name of the source''': NGC 1068
 
* '''Name of the source''': NGC 1068
 
* '''Brief description of the source''':  
 
* '''Brief description of the source''':  
−
** Object type : Unknown, AGN, SNR, Pulsar, PWN, GRB, etc
+
** Object type : AGN, Seyfert II galaxy
−
** Redshift (z) :
+
** Redshift (z) : 0.0037
−
** Other relevant information  
+
** Other relevant information: neutrino excess found by IceCube (https://doi.org/10.1126/science.abg3395)
−
** RA: XX XX XX.XX (hh mm ss), Dec: +/-YY YY YY.YY (dd mm ss)
+
** RA, Dec in deg (ICRS): 40.66962153, -0.0132943
−
** RA, Dec in deg (ICRS): XX.XX, YY.YY
 
  
 
* Analysis by Name (Institute - email)
 
* Analysis by Name (Institute - email)
−
** Daniel Morcuende
+
** Daniel Morcuende (IAA-CSIC, dmorcuende@iaa.es)
−
** José Manuel
+
** José Manuel Gavilán (IAA-CSIC)
  
 
==Data-taking information==
 
==Data-taking information==
 +
 +
There are a total of 124 runs from 2023-08-19 to 2024-02-10
 +
 
* '''Dates of data-taking''':  
 
* '''Dates of data-taking''':  
 
** Date (time of observation h)
 
** Date (time of observation h)
−
***ZD range (deg): Zenith Range
+
*** ZD range (deg): 25-65
−
***AZ range (deg): Azimuth Range
 
−
***Runs:
 
  
 +
=== Dataset after good quality selection ===
  
−
==Monte Carlo information==
+
[[ File:Ngc1068_LST1_quality_selection_cycle_I_dmorcuende.png ]]
−
* Link to MC files used:
 
−
** Particle types:
 
−
** ZD (deg): Zenith Range
 
−
** AZ (deg): Azimuth Range
 
  
−
==DL1 data==
+
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
  
  
−
Reprocessed using DL1a files produced by LSTOSA (lstchain v0.7.3) and dl1ab script (v0.7.5)
+
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],
  
−
* original DL1a files
+
[[File:Ngc1068_LST1_rel_light_yield_cycle_I_dmorcuende.png]]
  
−
/fefs/aswg/data/real/DL1/20210808/v0.7.3/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
+
==Analysis information==
−
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/XXX/zenith_20deg/south_pointing/20210416_v0.7.3_prod5_trans_80_local_taicut_8_4/
+
* lstchain v0.10
 +
* source-independent analysis
 +
* Gammapy v1.1
  
−
* lstchain v0.7.5
+
==Monte Carlo information==
−
* real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning
 
  
−
  /fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/lstchain_standard_config.json
+
Closest declination line towards La Palma latitude (28.8 deg)
  
−
* MC: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF+NSB tuning
+
20240131_allsky_v0.10.5_all_dec_base/dec_931
−
  "image_modifier": {
 
−
    "increase_nsb": true,
 
−
    "extra_noise_in_dim_pixels": 0.74,
 
−
    "extra_bias_in_dim_pixels": 0.38,
 
−
    "transition_charge": 8,
 
−
    "extra_noise_in_bright_pixels": 0.052,
 
−
    "increase_psf": true,
 
−
    "smeared_light_fraction": 0.125
 
−
  },
 
  
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/lstchain_dl1ab_tune_MC_to_BLLac_config.json
+
==DL1 data==
  
−
* Produced DL1b files
+
Standard DL1b files produced by LSTOSA (lstchain v0.10, tailcut 8-4 with cleaning based on pedestal RMS, dynamic cleaning)
−
 
 
−
/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/data/
 
−
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/Run0XXXX
 
  
 +
/fefs/aswg/data/real/DL1/YYYYMMDD/v0.10/tailcut84/dl1_LST-1.Run?????.h5
  
 
==Random forest==
 
==Random forest==
 +
Models:
  
−
Please include your specific .json files used in producing your RFs and for which source analysis (source dependent or independent), and variables used in the RF.
+
/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931
  
 +
== DL2 ==
 +
Files:
 +
/fefs/aswg/workspace/daniel.morcuende/data/real/DL2/NGC1068
  
−
==DL2 data==
+
==IRF and DL3 data selection==
  
 +
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
  
−
* lstchain-v0.10
+
* IRF linear interpolation
 +
* intensity > 50 p.e.
 +
* gammaness energy-dependent cuts with 70% efficiency
 +
* theta energy-dependent cuts 70% efficiency
  
−
* source-indep
+
==High-level analysis==
−
/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcindep/DL2/data/
 
−
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/srcindep/DL2/data/dl2_LST-1.Run0XXXX_merged.h5
 
  
−
==DL3 data selection==
+
* Science Tool: Gammapy 1.1
 +
* point-like IRF, 1D analysis
  
−
Information about your DL3 data selection.
+
==Analysis Results==
  
−
'''Example'''
+
=== Theta2 plot ===
−
*intensity > 50
 
−
*r: [0,  1 ]
 
−
*wl: [0.1, 1 ]
 
−
*leakage_intensity_width_2: [0, 0.2 ]
 
  
−
* source-indep
+
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
  
−
==High-level analysis==
 
−
 
−
* lstchain to generate source-dep IRF, DL3
 
−
* Science Tool: Gammapy 1.1
 
−
* point-like IRF, 1D analysis
 
  
−
==Analysis Results==
+
[[ File:Ngc1068_LST1_theta2_HE_cycle_I_dmorcuende.png ]]
  
−
Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here.
+
Energy range: 0.2 - 1000000.0 TeV
−
 
+
Excess: 14.0; Off: 3798
−
=== Theta2 plot ===
+
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 108: Line 135:
  
 
=== Spectral results ===
 
=== Spectral results ===
 +
 +
* 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]]
 +
 +
== 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

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