Difference between revisions of "Fast Offline Analysis OP313 May 2024 P064"

From my_wiki
Jump to: navigation, search
(Spectral results)
(Data-taking information)
Line 61: Line 61:
  
 
[[File:Diagnostic_op313_20240516_leaheck.png | 600 px]]
 
[[File:Diagnostic_op313_20240516_leaheck.png | 600 px]]
*Data quality selection applied on above data sample: All runs excluded by standard cuts because of dR/dI rate and intensity. Analysis still presented below. But even with noff=1, it says: Found overlapping on/off regions, choose less off regions ReflectedRegionsBackgroundMaker failed. No OFF region found outside exclusion mask for dataset '17426'. ReflectedRegionsBackgroundMaker failed. Setting 17426 mask to False.
+
*Data quality selection applied on above data sample: All runs excluded by standard cuts because of several reasons (see plot below). Analysis still presented below. But even with noff=1, it says: Found overlapping on/off regions, choose less off regions ReflectedRegionsBackgroundMaker failed. No OFF region found outside exclusion mask for dataset '17426'. ReflectedRegionsBackgroundMaker failed. Setting 17426 mask to False.
 
*Joint observations with MAGIC? : No
 
*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_16.pdf
 
*Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/05/run_summary_2024_05_16.pdf

Revision as of 12:49, 17 May 2024

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)
    • Lea Heckmann (MPP - heckmann@mpp.mpg.de)

Data-taking information

  • Dates of data-taking:
    • Date (time of observation h)
      • ZD range (deg): 0 - 16 deg
      • AZ range (deg): -64 - 73 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.
'2024-05-07': [17352, 17353, 17354, 17355, 17356] 0.95 h (no cuts)

Diagnostic op313 20240507 irenejm.png


'2024-05-15': [17417] 0.2 h (no cuts)

Diagnostic op313 20240515 leaheck.png


'2024-05-16': [17425,17426,17427,17428] 0.9 h (before cuts, 0 after)

Diagnostic op313 20240516 leaheck.png

  • Data quality selection applied on above data sample: All runs excluded by standard cuts because of several reasons (see plot below). Analysis still presented below. But even with noff=1, it says: Found overlapping on/off regions, choose less off regions ReflectedRegionsBackgroundMaker failed. No OFF region found outside exclusion mask for dataset '17426'. ReflectedRegionsBackgroundMaker failed. Setting 17426 mask to 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_16.pdf

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
/fefs/aswg/data/real/DL1/20240507/v0.10/tailcut84/dl1_LST-1.Run17352.h5
/fefs/aswg/data/real/DL1/20240507/v0.10/tailcut84/dl1_LST-1.Run17353.h5
/fefs/aswg/data/real/DL1/20240507/v0.10/tailcut84/dl1_LST-1.Run17354.h5
/fefs/aswg/data/real/DL1/20240507/v0.10/tailcut84/dl1_LST-1.Run17355.h5
/fefs/aswg/data/real/DL1/20240507/v0.10/tailcut84/dl1_LST-1.Run17356.h5
/fefs/aswg/data/real/DL1/20240515/v0.10/tailcut84/dl1_LST-1.Run17417.h5
/fefs/aswg/data/real/DL1/20240516/v0.10/tailcut84/dl1_LST-1.Run17425.h5
/fefs/aswg/data/real/DL1/20240516/v0.10/tailcut84/dl1_LST-1.Run17426.h5
/fefs/aswg/data/real/DL1/20240516/v0.10/tailcut84/dl1_LST-1.Run17427.h5
/fefs/aswg/data/real/DL1/20240516/v0.10/tailcut84/dl1_LST-1.Run17428.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
/fefs/aswg/workspace/irene.jimenez/FOA/OP313_20240507_data/DL2
/fefs/aswg/workspace/lea.heckmann/FOA/OP313_20240515/DL2
/fefs/aswg/workspace/lea.heckmann/FOA/OP313_20240516/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


  • 2024-05-07

Theta2 op313 20240507 irenejm.png

  • 2024-05-15

Theta2 op313 20240515 leaheck.png


  • 2024-05-16

Theta2 op313 20240516 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


  • 2024-05-07

Not enough signal!

  • 2024-05-15

Not enough signal!

  • 2024-05-16

Not enough signal!

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


  • 2024-05-07

Not enough signal!


  • 2024-05-15

Not enough signal!


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

LC op313 2024050104 leaheck.png