Fast Offline Analysis OP313 May 2024 P064

From my_wiki
Revision as of 08:35, 6 May 2024 by Lea (talk | contribs) (→‎Light curve)
Jump to: navigation, search

General information

  • Name of the source: OP 313
  • Brief description of the source:
    • Object type : FSRQ
    • Redshift (z) : 0.997
    • Other relevant information
    • RA: 13 10 28.6638 (hh mm ss), Dec: +32 20 43.783 (dd mm ss)
    • RA, Dec in deg (ICRS): 197.619432, 32.345495
  • Analysis by:
    • Irene Jiménez Martínez (MPP - irenejm@mpp.mpg.de or irene.jimenez@cta-consortium.org)

Data-taking information

  • Dates of data-taking:
    • Date (time of observation h)
      • ZD range (deg): 0 - 10 deg
      • AZ range (deg): -64 - 45 deg
      • Transmission9km:
      • Observations:
'2024-05-01': [17220, 17221, 17222, 17223] 0.48 h (after cuts)

Diagnostic op313 20240501 irenejm.png 

'2024-05-02': [17234, 17235, 17236, 17237] 0.82 h (no cuts)
  • Data quality selection applied on above data sample: Data quality selection could not be performed due to problem with run summary. During DL3 analysis an error ocurred for run 17237: Effective area is all zero at [197d37m09.957813s 32d20m43.782834s]. No safe energy band can be defined for the dataset '17237': setting `mask_safe` to all False.
  • Joint observations with MAGIC? : No
  • Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/05/run_summary_2024_05_02.pdf
'2024-05-04': [17281, 17282, 17283] h (all data cut but still analyzed to have an idea of the significance)

Diagnostic op313 20240504 leaheck.png

  • Data quality selection applied on above data sample: all 3 runs are excluded because of dR/dI rate and intensity threshold. Effective area is all zero at [197d37m09.957813s 32d20m43.782834s]. No safe energy band can be defined for the dataset '17282': setting `mask_safe` to all False.

General analysis information

  • source-independent analysis
  • lstchain version: v0.10.7 (software school)
  • DL1b standard production (v0.10/tailcut84), run-wise merged files in:
/fefs/aswg/data/real/DL1/20240501/v0.10/tailcut84/dl1_LST-1.Run17220.h5
/fefs/aswg/data/real/DL1/20240501/v0.10/tailcut84/dl1_LST-1.Run17221.h5
/fefs/aswg/data/real/DL1/20240501/v0.10/tailcut84/dl1_LST-1.Run17222.h5
/fefs/aswg/data/real/DL1/20240501/v0.10/tailcut84/dl1_LST-1.Run17223.h5

/fefs/aswg/data/real/DL1/20240502/v0.10/tailcut84/dl1_LST-1.Run17234.h5
/fefs/aswg/data/real/DL1/20240502/v0.10/tailcut84/dl1_LST-1.Run17235.h5
/fefs/aswg/data/real/DL1/20240502/v0.10/tailcut84/dl1_LST-1.Run17236.h5
/fefs/aswg/data/real/DL1/20240502/v0.10/tailcut84/dl1_LST-1.Run17237.h5

/fefs/aswg/data/real/DL1/20240504/v0.10/tailcut84/dl1_LST-1.Run17281.h5
/fefs/aswg/data/real/DL1/20240504/v0.10/tailcut84/dl1_LST-1.Run17282.h5
/fefs/aswg/data/real/DL1/20240504/v0.10/tailcut84/dl1_LST-1.Run17283.h5

Monte Carlo information

Write below the MC production id among the existing MC

  • All-sky MC production used: 20240131_allsky_v0.10.5_all_dec_base
    • Declination line (training): dec_3476


Random forest

Path to the models used: 

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

DL2 data

Config file: 

/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_3476/lstchain_config_2024-01-31.json

Files in: 

/fefs/aswg/workspace/irene.jimenez/FOA/OP313_20240501_data/DL2

/fefs/aswg/workspace/irene.jimenez/FOA/OP313_20240502_data/DL2

/fefs/aswg/workspace/lea.heckmann/FOA/OP313_20240504/DL2

DL3 data & IRF

 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.01, 1],
     "leakage_intensity_width_2": [0, 1],
     "event_type": [32, 32]
   }
 },
 "DL3Cuts": {
   "min_event_p_en_bin": 100,
   "min_gh_cut": 0.1,
   "max_gh_cut": 0.98,
   "min_theta_cut": 0.1,
   "max_theta_cut": 0.32,
   "fill_theta_cut": 0.32,
   "allowed_tels": [1]
 },
  • 70 % efficiency for both gammaness and theta gamma-ray selection cuts
  • Point-like IRF

High-level analysis

  • Gammapy 1.1
  • 1D analysis (point-like IRF)

Analysis Results

Theta2 plot

Settings:

  • Theta2 cut for signal extraction: 0.04 deg2
  • Normalization range: 0.1 - 0.4 deg2
  • Number of off positions: 3
  • Energy range: 100 GeV - 10 TeV

Daily results

  • 2024-05-01

Theta2 op313 20240501 irenejm.png.png

  • 2024-05-02

Theta2 op313 20240502 irenejm.png


  • 2024-05-04

Theta2 op313 20240504 leaheck.png


Stacked observations

  • 2024-05-01 to 2024-05-04


Theta2 op313 2024050104 leaheck.png


Significance and Excess sky maps

Spectral results

  • 2024-05-01

Spectral parameters: PL+EBL(dominguez) 

  type      name         value         unit        error       min     max frozen is_norm link
--------  ----------  ----------  --------------  --------- ---------  --- ------ ------- ----
spectral       index  1.9394e+00                  7.889e-01       nan  nan  False   False 
spectral   amplitude  9.2076e-09  TeV-1 s-1 cm-2  2.605e-09       nan  nan  False    True  
spectral   reference  1.0000e-01             TeV  0.000e+00       nan  nan   True   False 
spectral  alpha_norm  1.0000e+00                  0.000e+00       nan  nan   True   False 
spectral    redshift  9.9700e-01                  0.000e+00       nan  nan   True   False

Integral flux E > 100 GeV: 1.3565e-10 s-1 cm-2

SED op313 20240501 irenejm.png

  • 2024-05-02

Spectral parameters: PL+EBL(dominguez) 

  type      name         value         unit        error       min     max frozen is_norm link
--------  ----------  ----------  --------------  --------- ---------  --- ------ ------- ----
spectral       index  3.9883e+00                  1.507e+00       nan  nan  False   False 
spectral   amplitude  2.8111e-09  TeV-1 s-1 cm-2  3.443e-09       nan  nan  False    True  
spectral   reference  1.0000e-01             TeV  0.000e+00       nan  nan   True   False 
spectral  alpha_norm  1.0000e+00                  0.000e+00       nan  nan   True   False 
spectral    redshift  9.9700e-01                  0.000e+00       nan  nan   True   False

Integral flux E > 100 GeV: 2.7110e-11 s-1 cm-2

SED op313 20240502 irenejm.png


  • 2024-05-04

Spectral parameters: PL+EBL(dominguez) 

  type      name         value         unit        error       min     max frozen is_norm link
--------  ----------  ----------  --------------  --------- ---------  --- ------ ------- ----
spectral       index  3.6192                       1.4352       nan  nan  False   False 
spectral   amplitude  2.18903-09  TeV-1 s-1 cm-2  3.94539e-09       nan  nan  False    True  
spectral   reference  1.0000e-01             TeV  0.000e+00       nan  nan   True   False 
spectral  alpha_norm  1.0000e+00                  0.000e+00       nan  nan   True   False 
spectral    redshift  9.9700e-01                  0.000e+00       nan  nan   True   False

Integral flux E > 100 GeV: XXe-11 s-1 cm-2

SED op313 20240504 leaheck.png

  • Stacked 2024-05-01 to 2024-05-04


Spectral parameters: PL+EBL(dominguez) 

  type      name         value         unit        error       min     max frozen is_norm link
--------  ----------  ----------  --------------  --------- ---------  --- ------ ------- ----

spectral index 2.3038e+00 4.329e-01 nan nan False False spectral amplitude 5.7638e-09 TeV-1 s-1 cm-2 1.572e-09 nan nan False True spectral reference 1.0000e-01 TeV 0.000e+00 nan nan True False spectral alpha_norm 1.0000e+00 0.000e+00 nan nan True False spectral redshift 9.9700e-01 0.000e+00 nan nan True False

Integral flux E > 100 GeV: 7.7934464e-11 s-1 cm-2

SED op313 2024050104 leaheck.png

Light curve

  • 2024-05-01

LC op313 20240501 irenejm.png

  • 2024-05-02

LC op313 20240502 irenejm.png


  • 2024-05-04

LC op313 20240504 leaheck.png


  • Stacked 2024-05-01 to 2024-05-04

LC op313 2024050104 leaheck.png

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