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

From my_wiki
Jump to: navigation, search
(Data-taking information)
(Light curve)
 
(35 intermediate revisions by 2 users not shown)
Line 10: Line 10:
 
* Analysis by:
 
* Analysis by:
 
**Irene Jiménez Martínez (MPP - irenejm@mpp.mpg.de or irene.jimenez@cta-consortium.org)
 
**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==
 
==Data-taking information==
 
* '''Dates of data-taking''':  
 
* '''Dates of data-taking''':  
 
** Date (time of observation h)
 
** Date (time of observation h)
***ZD range (deg): 0 - 10 deg
+
***ZD range (deg): 0 - 16 deg
***AZ range (deg): -64 - 45 deg
+
***AZ range (deg): -64 - 73 deg
 
***Transmission9km:  
 
***Transmission9km:  
 
***Observations:
 
***Observations:
Line 32: Line 33:
 
*Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/05/run_summary_2024_05_02.pdf
 
*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] 0 h (all data cut)
+
  '2024-05-04': [17281, 17282, 17283] h (all data cut but still analyzed to have an idea of the significance)
  
 
[[File:Diagnostic_op313_20240504_leaheck.png | 600 px]]
 
[[File:Diagnostic_op313_20240504_leaheck.png | 600 px]]
*Data quality selection applied on above data sample: all 3 runs are excluded because of  dR/dI rate and intensity threshold
+
*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.
 +
 
 
*Joint observations with MAGIC? : Yes
 
*Joint observations with MAGIC? : Yes
 
*Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/05/run_summary_2024_05_04.pdf
 
*Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/05/run_summary_2024_05_04.pdf
 +
 +
'2024-05-07': [17352, 17353, 17354, 17355, 17356] 0.95 h (no cuts)
 +
 +
[[File:Diagnostic_op313_20240507_irenejm.png | 600 px]]
 +
*Data quality selection applied on above data sample: No data survived the cuts, there were too many clouds. Fast analysis performed with no cuts. Run 17356 has less than 5 minutes of data.
 +
*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_07.pdf
 +
 +
 +
'2024-05-15': [17417] 0.2 h (no cuts)
 +
 +
[[File:Diagnostic_op313_20240515_leaheck.png | 600 px]]
 +
*Data quality selection applied on above data sample: All data survived the cuts
 +
*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_15.pdf
 +
 +
 +
 +
'2024-05-16': [17425,17426,17427,17428] 0.9 h (before cuts, 0 after)
 +
 +
[[File:Diagnostic_op313_20240516_leaheck.png | 600 px]]
 +
*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==
 
==General analysis information==
Line 44: Line 70:
 
* DL1b standard production (v0.10/tailcut84), run-wise merged files in:  
 
* 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.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.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.Run17234.h5
Line 52: Line 80:
 
  /fefs/aswg/data/real/DL1/20240502/v0.10/tailcut84/dl1_LST-1.Run17237.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==
 
==Monte Carlo information==
Line 76: Line 120:
  
 
  /fefs/aswg/workspace/lea.heckmann/FOA/OP313_20240504/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==
 
==DL3 data & IRF==
Line 120: Line 170:
  
 
=== Daily results ===
 
=== Daily results ===
 
 
=== Stacked observations ===
 
 
 
* 2024-05-01
 
* 2024-05-01
  
Line 131: Line 177:
  
 
[[File:Theta2_op313_20240502_irenejm.png|600 px]]
 
[[File:Theta2_op313_20240502_irenejm.png|600 px]]
 +
 +
 +
* 2024-05-04
 +
 +
[[File:Theta2_op313_20240504_leaheck.png|600 px]]
 +
 +
 +
* 2024-05-07
 +
 +
[[File:Theta2_op313_20240507_irenejm.png|600 px]]
 +
 +
* 2024-05-15
 +
 +
[[File:Theta2_op313_20240515_leaheck.png|600 px]]
 +
 +
 +
* 2024-05-16
 +
 +
[[File:Theta2_op313_20240516_leaheck.png|600 px]]
 +
 +
=== Stacked observations ===
 +
* 2024-05-01 to 2024-05-04
 +
 +
 +
[[File:Theta2_op313_2024050104_leaheck.png|600 px]]
  
  
Line 171: Line 242:
  
  
 +
* 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
 +
 +
[[File:SED_op313_20240504_leaheck.png|600 px]]
 +
 +
* 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
 +
 +
[[File:SED_op313_2024050104_leaheck.png|600 px]]
 +
 +
 +
* 2024-05-07
 +
 +
Not enough signal!
 +
 +
* 2024-05-15
 +
 +
Not enough signal!
 +
 +
* 2024-05-16
 +
 +
Not enough signal!
  
 
=== Light curve ===
 
=== Light curve ===
Line 181: Line 297:
  
 
[[File:LC_op313_20240502_irenejm.png|600 px]]
 
[[File:LC_op313_20240502_irenejm.png|600 px]]
 +
 +
 +
* 2024-05-04
 +
 +
[[File:LC_op313_20240504_leaheck.png|600 px]]
 +
 +
 +
* 2024-05-07
 +
 +
Not enough signal!
 +
 +
 +
* 2024-05-15
 +
 +
Not enough signal!
 +
 +
 +
* 2024-05-16
 +
 +
Not enough signal!
 +
 +
 +
 +
* Stacked 2024-05-01 to 2024-05-04
 +
 +
[[File:LC_op313_2024050104_leaheck.png|600 px]]

Latest revision as of 12:50, 17 May 2024

General information[edit]

  • 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[edit]

  • 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[edit]

  • 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[edit]

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[edit]

Path to the models used:

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

DL2 data[edit]

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[edit]

 "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[edit]

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

Analysis Results[edit]

Theta2 plot[edit]

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[edit]

  • 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[edit]

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


Theta2 op313 2024050104 leaheck.png


Significance and Excess sky maps[edit]

Spectral results[edit]

  • 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[edit]

  • 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!


  • 2024-05-16

Not enough signal!


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

LC op313 2024050104 leaheck.png