Difference between revisions of "Analysis OP313 Dec 2023 ToO"
Jorge.otero (talk | contribs) (→Analysis 1 (Jorge Otero)) |
Jorge.otero (talk | contribs) (→Analysis 2 (Jorge Otero)) |
||
(26 intermediate revisions by 2 users not shown) | |||
Line 26: | Line 26: | ||
== Meeting == | == Meeting == | ||
* Kickoff meeting (19 Jan 2024): https://indico.cta-observatory.org/event/5330/ | * Kickoff meeting (19 Jan 2024): https://indico.cta-observatory.org/event/5330/ | ||
+ | * meeting 02 February 2024: https://indico.cta-observatory.org/event/5379/ | ||
+ | * meeting 16 February 2024: https://indico.cta-observatory.org/event/5412/ | ||
− | ==Data-taking information== | + | ==Data-taking information and data quality selection== |
* observation summary | * observation summary | ||
{| class="wikitable" | {| class="wikitable" | ||
Line 69: | Line 71: | ||
|} | |} | ||
− | == | + | Final data quality selection (using the latest notebook, default settings + min_drdi_at_422pe = 1.2 instead of the default 1.5 to recover 1.7 h from January 2024 where we are limited by low statistic) |
+ | |||
+ | '2023-12-09': [15942, 15943, 15944, 15945], | ||
+ | '2023-12-10': [15975, 15976, 15977, 15978], | ||
+ | '2023-12-11': [16007, 16008, 16009, 16010], | ||
+ | '2023-12-12': [16042, 16043, 16044, 16045], | ||
+ | '2023-12-13': [16079, 16080, 16081, 16082, 16083], | ||
+ | '2023-12-14': [16114, 16116, 16117, 16118, 16119], | ||
+ | '2023-12-15': [16150, 16151, 16152, 16153, 16154, 16155, 16156], | ||
+ | '2023-12-16': [16188, 16189, 16191, 16192, 16193], | ||
+ | '2023-12-17': [16220, 16221, 16222, 16223, 16224, 16225], | ||
+ | '2023-12-18': [16240, 16241, 16242, 16243, 16244, 16246], | ||
+ | '2024-01-09': [16293, 16294, 16295], | ||
+ | '2024-01-13': [16350, 16351, 16352, 16353, 16354, 16356, 16357, 16358, 16359], | ||
+ | '2024-01-17': [16391, 16392], | ||
+ | '2024-01-19': [16412, 16413, 16414, 16415], | ||
+ | '2024-01-21': [16418, 16419, 16420] | ||
+ | |||
+ | == Analysis 1 (Daniel Morcuende) == | ||
[[File:OP313_20231213_Morcuende.png|thumb|centre|600px|Src-independent analysis, Gammaness cut efficiency of 40%, Theta containment 70%, 3 OFF positions, theta2 cut at 0.04 deg2]] | [[File:OP313_20231213_Morcuende.png|thumb|centre|600px|Src-independent analysis, Gammaness cut efficiency of 40%, Theta containment 70%, 3 OFF positions, theta2 cut at 0.04 deg2]] | ||
Line 77: | Line 97: | ||
Standard MC dark analysis | Standard MC dark analysis | ||
− | -Zd~ | + | Source independent analysis |
+ | |||
+ | -Zd~15-54 deg | ||
-MC prod: 20230901_v0.10.4_allsky_base_prod/…/dec_3476 | -MC prod: 20230901_v0.10.4_allsky_base_prod/…/dec_3476 | ||
+ | |||
+ | -3 OFFs | ||
+ | |||
+ | -Safe mask > 5% max. Aeff → Fit energy range: [70 GeV, 1 TeV] | ||
+ | |||
+ | A Crab Nebula analysis sanity check for the same nights as OP 313 can be checked here: https://www.lst1.iac.es/wiki/index.php/Crab_Sanity_check_JorgeOtero. | ||
* Theta2 plots | * Theta2 plots | ||
− | |||
− | |||
− | |||
− | + | Energy dependent cuts: 0.7gh efficiency, 0.7theta efficiency | |
− | Energy dependent cuts: 0. | + | (maximum theta cut restricted to 0.2 to avoid ON-OFF overlapping in post-DL3 analysis) |
{| align='center' | {| align='center' | ||
− | | [[File: | + | | [[File:OP313_20231209-20240121_endep_Otero.png|thumb|centre|400px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs (all Dec.+Jan. data)]] |
− | | [[File: | + | | [[File:OP313_20231209-20231220_endep_high_Otero.png|thumb|centre|400px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, high state (Dec. 2023)]] |
+ | | [[File:OP313_20240109-20240121_endep_low_Otero.png|thumb|centre|400px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, low state (Jan. 2024)]] | ||
|} | |} | ||
− | Different energy bins: | + | Different energy bins: |
− | [[File: | + | [[File:OP313_20231209-20240121_endep_binned_Otero.png|thumb|centre|800px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] |
+ | |||
+ | Nightwise theta2 plots: | ||
+ | {| align='center' | ||
+ | | [[File:OP313_20231209_endep_binned_Otero.png|thumb|centre|400px|2023-12-09, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231210_endep_binned_Otero.png|thumb|centre|400px|2023-12-10, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231211_endep_binned_Otero.png|thumb|centre|400px|2023-12-11, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | |} | ||
+ | {| align='center' | ||
+ | | [[File:OP313_20231212_endep_binned_Otero.png|thumb|centre|400px|2023-12-12, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231213_endep_binned_Otero.png|thumb|centre|400px|2023-12-13, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231214_endep_binned_Otero.png|thumb|centre|400px|2023-12-14, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | |} | ||
+ | {| align='center' | ||
+ | | [[File:OP313_20231215_endep_binned_Otero.png|thumb|centre|400px|2023-12-15, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231216_endep_binned_Otero.png|thumb|centre|400px|2023-12-16, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20231217_endep_binned_Otero.png|thumb|centre|400px|2023-12-17, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | |} | ||
+ | {| align='center' | ||
+ | | [[File:OP313_20231218_endep_binned_Otero.png|thumb|centre|400px|2023-12-18, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20240109_endep_binned_Otero.png|thumb|centre|400px|2024-01-09, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20240113_endep_binned_Otero.png|thumb|centre|400px|2024-01-13, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | |} | ||
+ | {| align='center' | ||
+ | | [[File:OP313_20240117_endep_binned_Otero.png|thumb|centre|400px|2024-01-17, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20240119_endep_binned_Otero.png|thumb|centre|400px|2024-01-19, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | | [[File:OP313_20240121_endep_binned_Otero.png|thumb|centre|400px|2024-01-21, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] | ||
+ | |} | ||
* SED | * SED | ||
− | |||
− | |||
− | -alpha= | + | Assuming a PWL+EBL spectral model (Saldana-Lopez 2021 EBL model): |
+ | {| align='center' | ||
+ | | [[File:OP313_20231209-20231218_SED_powerlaw_EBL_SL21_Otero.png|thumb|centre|400px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, high state (Dec. 2023)]] | ||
+ | | [[File:OP313_20240109-20240121_SED_powerlaw_Otero.png|thumb|centre|400px|energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, low state (Jan. 2024)]] | ||
+ | |} | ||
+ | |||
+ | Spectral parameters for the high state (December 2023): | ||
+ | F0=(1.15 土 0.13)・10e-11 TeV-1 cm-2 s-1 | ||
+ | alpha=2.07 土 0.51 | ||
+ | E0=100 GeV | ||
− | + | Spectral parameters for the low state (January 2024): | |
+ | F0=(2.6 土 2.5)・10e-12 TeV-1 cm-2 s-1 | ||
+ | alpha=2.4 土 1.4 | ||
+ | E0=450 GeV | ||
− | + | * Lightcurve | |
− | + | Complete lightcurve above 100 GeV assuming the SED obtained for the high state in December 2023: | |
− | [[File: | + | [[File:OP313_20231209-20240121_LC_DecJan_Otero.png|thumb|centre|600px|Daily binned lightcurve above 100GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] |
+ | Average flux: | ||
+ | F=(4.69 土 0.10)・10e-11 cm-2 s-1 = 10% Crab Nebula flux | ||
− | + | Separating high and low state: | |
− | |||
{| align='center' | {| align='center' | ||
− | | [[File: | + | | [[File:OP313_20231209-20231218_LC_Dec_Otero.png|thumb|centre|600px|Daily-binned lightcurve above 100GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, high state (Dec. 2023)]] |
− | | [[File: | + | | [[File:OP313_20240109-20240121_LC_Jan_Otero.png|thumb|centre|600px|Daily-binned lightcurve above 100GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs, low state (Jan. 2024)]] |
|} | |} | ||
− | == | + | Average flux for the high state (December 2023): |
+ | F=(1.35 土 0.19)・10e-10 cm-2 s-1 = 28.3% Crab Nebula flux | ||
+ | |||
+ | Average flux for the low state (January 2024) compatible with 0 flux. | ||
+ | |||
+ | == Analysis 3 (Seiya Nozaki) == | ||
+ | Standard MC dark analysis | ||
+ | |||
+ | Source dependent analysis | ||
+ | |||
+ | -Zd~15-54 deg | ||
+ | |||
+ | -MC prod: 20230901_v0.10.4_allsky_base_prod/…/dec_3476 | ||
+ | |||
* 2023-12-09 | * 2023-12-09 | ||
** run 15942, 15943, 15944, 15945 | ** run 15942, 15943, 15944, 15945 | ||
Line 204: | Line 281: | ||
| [[File:OP313_20240109_srcdep_g07_i50_stacked_Nozaki.png|thumb|centre|400px|intensity>50, gammaness>0.7 (stacked)]] | | [[File:OP313_20240109_srcdep_g07_i50_stacked_Nozaki.png|thumb|centre|400px|intensity>50, gammaness>0.7 (stacked)]] | ||
|} | |} | ||
+ | |||
+ | == Crosschecks of plots to be approved == | ||
+ | |||
+ | ADD HERE SEDs AND LCs OF DANIEL, SEIYA AND JORGE PLOTTED ALL TOGETHER | ||
== Fast Analysis (Joshua Baxter) == | == Fast Analysis (Joshua Baxter) == |
Latest revision as of 13:52, 19 February 2024
Contents
General information[edit]
- Name of the source: OP313
- Brief description of the source:
- - Object type : FSRQ
- - Redshift : z=0.997 (Schneider et al. 2010, AJ, 139, 2360)
Atel[edit]
The Large-Sized Telescope (LST-1) on La Palma has been monitoring the very distant Flat Spectrum Radio Quasar (FSRQ) OP 313 (z=0.997, Schneider et al. 2010, AJ, 139, 2360) since November 2023. Following the announcement of enhanced gamma-ray emission by Fermi-LAT (ATel #16356) and several optical facilities (ATel #16360) in early December, the Fermi- LAT emission of OP 313 has been closely monitored using the FlaapLUC pipeline (Astronomy and Computing, Volume 22, p. 9-15, 2018). This monitoring revealed the detection of renewed activity in the high-energy (HE, E>100 MeV) band and so, Target of Opportunity observations with LST-1 were triggered on December 10th 2023. OP 313 was detected by LST-1 with a preliminary offline analysis using data from 2023/12/11 to 2023/12/14. It was detected with a significance greater than 5 sigma and an integrated flux, above 100 GeV, at 15% flux of the Crab Nebula. LST-1 observations on OP 313 will continue during the next few nights and therefore multi-wavelength observations are highly encouraged. LST-1 is a prototype of the Large-Sized Telescope for the Cherenkov Telescope Array Observatory, and is located on the Canary island of La Palma, Spain. The telescope design is optimized for observation of gamma rays in the range from 20 GeV to 3 TeV. The preliminary offline analysis has been performed by Daniel Morcuende (dmorcuende@iaa.es), Jorge Otero-Santos (joteros@iaa.es) and Seiya Nozaki (nozaki@mpp.mpg.de). The LST-1 contact persons for these observations are Masahiro Teshima (mteshima@mpp.mpg.de) and Juan Cortina (juan.cortina@ciemat.es).
Meeting[edit]
- Kickoff meeting (19 Jan 2024): https://indico.cta-observatory.org/event/5330/
- meeting 02 February 2024: https://indico.cta-observatory.org/event/5379/
- meeting 16 February 2024: https://indico.cta-observatory.org/event/5412/
Data-taking information and data quality selection[edit]
- observation summary
Date | Runs | time (h) | Zd | Az | Transmission | Comments |
2023-12-09 | 15942, 15943, 15944, 15945 | 1.31 | 35-53 | - | - | 50-sec shifted timestamp -> solved! |
2023-12-10 | 15973, 15974, 15975, 15976, 15977, 15978 | 1.61 | 30-51 | - | - | |
2023-12-11 | 16007, 16008, 16009, 16010, 16011 | 1.39 | 30-54 | - | - | |
2023-12-12 | 16042, 16043, 16044, 16045, 16046 | 1.63 | 30-54 | - | - | |
2023-12-13 | 16079, 16080, 16081, 16082, 16083 | 1.69 | 30-53 | - | - | |
2023-12-14 | 16114, 16115, 16116, 16117, 16118, 16119 | 1.60 | 30-52 | - | - | |
2023-12-15 | 16150, 16151, 16152, 16153, 16154, 16155, 16156, 16157* | 1.75 | 30-52 | - | - | No tracking during run 16157 |
2023-12-16 | 16188, 16189, 16190*, 16191, 16192, 16193 | 1.67 | 25-52 | - | - | Drive issue during run 16190 (no tracking) |
2023-12-17 | 16220, 16221, 16222, 16223, 16224, 16225 | 1.72 | 25-50 | - | - | |
2023-12-18 | 16240, 16241, 16242, 16243, 16244, 16246 | 25-50 | - | - | online gain-selected data | |
2023-12-19 | 16258, 16259, 16260, 16261, 16262, 16263 | 25-50 | - | - | online gain-selected data | |
2023-12-20 | 16269, 16270 | 40-43 | - | - | online gain-selected data, humidity is rising during run 16270 | |
2024-01-09 | 16293, 16294, 16295 | 0.63 | 20-26 | - | - | |
2024-01-13 | 16350, 16351, 16352, 16353, 16354, 16355, 16356, 16357, 16358, 16359 | 2.08 | 15-43 | - | - | |
2024-01-17 | 16391, 16392 | 0.46 | 25-30 | - | - | Timestamp is shifted?? -> solved! |
2024-01-18 | 16407 | 0.10 | 4-6 | - | - | |
2024-01-19 | 16412, 16413, 16414, 16415 | 1.34 | 5-30 | - | - | |
2024-01-21 | 16418, 16419, 16420 | 0.70 | 4-8 | - | - |
Final data quality selection (using the latest notebook, default settings + min_drdi_at_422pe = 1.2 instead of the default 1.5 to recover 1.7 h from January 2024 where we are limited by low statistic)
'2023-12-09': [15942, 15943, 15944, 15945], '2023-12-10': [15975, 15976, 15977, 15978], '2023-12-11': [16007, 16008, 16009, 16010], '2023-12-12': [16042, 16043, 16044, 16045], '2023-12-13': [16079, 16080, 16081, 16082, 16083], '2023-12-14': [16114, 16116, 16117, 16118, 16119], '2023-12-15': [16150, 16151, 16152, 16153, 16154, 16155, 16156], '2023-12-16': [16188, 16189, 16191, 16192, 16193], '2023-12-17': [16220, 16221, 16222, 16223, 16224, 16225], '2023-12-18': [16240, 16241, 16242, 16243, 16244, 16246], '2024-01-09': [16293, 16294, 16295], '2024-01-13': [16350, 16351, 16352, 16353, 16354, 16356, 16357, 16358, 16359], '2024-01-17': [16391, 16392], '2024-01-19': [16412, 16413, 16414, 16415], '2024-01-21': [16418, 16419, 16420]
Analysis 1 (Daniel Morcuende)[edit]
Analysis 2 (Jorge Otero)[edit]
Standard MC dark analysis
Source independent analysis
-Zd~15-54 deg
-MC prod: 20230901_v0.10.4_allsky_base_prod/…/dec_3476
-3 OFFs
-Safe mask > 5% max. Aeff → Fit energy range: [70 GeV, 1 TeV]
A Crab Nebula analysis sanity check for the same nights as OP 313 can be checked here: https://www.lst1.iac.es/wiki/index.php/Crab_Sanity_check_JorgeOtero.
- Theta2 plots
Energy dependent cuts: 0.7gh efficiency, 0.7theta efficiency (maximum theta cut restricted to 0.2 to avoid ON-OFF overlapping in post-DL3 analysis)
Different energy bins:
Nightwise theta2 plots:
- SED
Assuming a PWL+EBL spectral model (Saldana-Lopez 2021 EBL model):
Spectral parameters for the high state (December 2023): F0=(1.15 土 0.13)・10e-11 TeV-1 cm-2 s-1 alpha=2.07 土 0.51 E0=100 GeV
Spectral parameters for the low state (January 2024): F0=(2.6 土 2.5)・10e-12 TeV-1 cm-2 s-1 alpha=2.4 土 1.4 E0=450 GeV
- Lightcurve
Complete lightcurve above 100 GeV assuming the SED obtained for the high state in December 2023:
Average flux: F=(4.69 土 0.10)・10e-11 cm-2 s-1 = 10% Crab Nebula flux
Separating high and low state:
Average flux for the high state (December 2023): F=(1.35 土 0.19)・10e-10 cm-2 s-1 = 28.3% Crab Nebula flux
Average flux for the low state (January 2024) compatible with 0 flux.
Analysis 3 (Seiya Nozaki)[edit]
Standard MC dark analysis
Source dependent analysis
-Zd~15-54 deg
-MC prod: 20230901_v0.10.4_allsky_base_prod/…/dec_3476
- 2023-12-09
- run 15942, 15943, 15944, 15945
- standard MC (dark)
- 2023-12-10
- run 15973, 15974, 15975, 15976, 15977, 15978
- standard MC (dark)
- 2023-12-11
- run 16007, 16008, 16009, 16010, 16011
- standard MC (dark)
- 2023-12-12
- run 16042, 16043, 16044, 16045, 16046
- standard MC (dark)
- 2023-12-13
- run 16079, 16080, 16081, 16082, 16083
- standard MC (dark)
- 2023-12-14
- run 16114, 16115, 16116, 16117, 16118, 16119
- standard MC (dark)
- 2023-12-15
- run 16150, 16151, 16152, 16153, 16154, 16155, 16156, 16157
- standard MC (dark)
- 2023-12-16
- run 16188, 16189, 16190, 16191, 16192, 16193
- standard MC (dark)
- 2023-12-17
- run 16220, 16221, 16222, 16223, 16224, 16225
- standard MC (dark)
- 2024-01-09
- run 16293, 16294, 16295
- standard MC (dark)
Crosschecks of plots to be approved[edit]
ADD HERE SEDs AND LCs OF DANIEL, SEIYA AND JORGE PLOTTED ALL TOGETHER
Fast Analysis (Joshua Baxter)[edit]
- 2023-12-14
- Run Number: [16114, 16115, 16116, 16117, 16118, 16119]
gamma efficiency = 50 %, theta = 68 %
- 2023-12-15
- Run Number: [16150, 16151, 16152, 16153, 16154, 16155, 16156, 16157]
gamma efficiency = 50 %, theta = 68 %
- 2023-12-16
- Run Number: [16187, 16188, 16189, 16190, 16191, 16192, 16193]
gamma efficiency = 50 %, theta = 68 %
Fast Analysis (Elisa Visentin)[edit]
Joint runs: 2023/12/09, 2023/12/10, 2023/12/11, 2023/12/14, 2023/12/15, 2023/12/16, 2023/12/17, 2023/12/18, 2024/01/09
gamma efficiency 90%, theta cut 0.2 deg, MAGIC intensity > 30 pe, LST intensity > 50 pe Data selection: LST runs where > 3 minutes of stereo joint data taking
Fermi LAT[edit]
Mireia[edit]
David[edit]
Analysis time range
tmin = 723081605. tmax = 725846405.
Older analysis
3 hours bin LC from tmin = 723956405.0 (2023-12-11 03:00:00.000) to tmax = 724550405.000 (2023-12-18 00:00:00.000)
Corresponding SED