Difference between revisions of "NGC 1275"

From my_wiki
Jump to: navigation, search
(→‎Cycle I (including major flare Dec 2022))
(→‎IC 310 dec line comparison)
 
(174 intermediate revisions by 3 users not shown)
Line 3: Line 3:
 
* '''Name of the source''': NGC 1275
 
* '''Name of the source''': NGC 1275
 
* '''Brief description of the source''':  
 
* '''Brief description of the source''':  
−
:- Object type : Seyfert 2 Galaxy
+
:- Object type : FR-I Radio Galaxy + Spiral galaxy (in interaction)
 
:- RA: 3 19 48.16 (hh mm ss), Dec: +41 30 42.1088 (dd mm ss)
 
:- RA: 3 19 48.16 (hh mm ss), Dec: +41 30 42.1088 (dd mm ss)
  
Line 16: Line 16:
 
= Cycle I (including major flare Dec 2022) =
 
= Cycle I (including major flare Dec 2022) =
  
−
== Mireia Nievas ==
+
== Mireia ==
 +
 
 +
* Producer: Mireia Nievas Rosillo <mnievas@iac.es>
  
 
=== Data acquisition ===
 
=== Data acquisition ===
 +
 +
=== old datacheck notebook ===
 +
 +
(below some key points from Abelardo's old datacheck selection notebook, only showing those relevant for data selection)
  
 
* ZD and AZ distributions
 
* ZD and AZ distributions
Line 30: Line 36:
 
[[File:NGC1275_CycleI_MoonConditions2_mnievas.png|720px]]
 
[[File:NGC1275_CycleI_MoonConditions2_mnievas.png|720px]]
  
−
== Preliminary results ==
+
* Rate of cosmics
 +
 
 +
[[File:NGC1275_CycleI_RateOfCosmics_mnievas.png|720px]]
 +
 
 +
* Rate of pulses above a given pulse size
 +
 
 +
[[File:NGC1275_CycleI_RateAboveCharge_mnievas.png|720px]]
 +
 
 +
* NSB after all cuts
 +
 
 +
[[File:NGC1275_CycleI_NSB_aftercuts_mnievas.png|500px]]
 +
 
 +
=== new data_select notebook ===
 +
 
 +
[[File:NGC1275_CycleI_ds1_mnievas.png|720px]]
 +
 
 +
[[File:NGC1275_CycleI_ds2_mnievas.png|720px]]
 +
 
 +
[[File:NGC1275_CycleI_ds4_mnievas.png|720px]]
 +
 
 +
[[File:NGC1275_CycleI_ds5_mnievas.png|720px]]
 +
 
 +
=== Data selection ===
 +
 
 +
(old selection using Abelardo's old datacheck notebook)
 +
 
 +
* 204 wobble runs for the selected source (out of 246)
 +
* Observation time: 53.56 hours (out of 65.53 hours)
 +
* (most data in Perseus-MA mode, shared with IC310, some runs in NGC1275 pointing mode)
 +
* (moon data is removed)
 +
 
 +
<nowiki>
 +
1 : 2022-12-20 : [11522, 11528, 11529, 11530, 11531, 11532, 11522, 11528, 11529, 11530, 11531, 11532]
 +
2 : 2022-12-21 : [11546, 11547, 11548, 11549, 11550]
 +
3 : 2022-12-22 : [11566, 11568, 11569, 11570, 11571, 11572, 11573, 11574, 11566, 11568, 11569, 11570, 11571, 11572, 11573, 11574]
 +
4 : 2022-12-23 : [11595, 11596, 11597, 11598, 11599, 11600]
 +
5 : 2023-01-13 : [11618, 11619]
 +
6 : 2023-01-14 : [11641, 11642, 11643, 11644, 11641, 11642, 11643, 11644]
 +
7 : 2023-01-15 : [11660, 11661, 11662, 11663, 11660, 11661, 11662, 11663]
 +
8 : 2023-08-15 : [13935, 13936, 13937, 13938, 13939, 13940, 13941, 13942, 13943, 13944]
 +
9 : 2023-08-17 : [13987, 13988, 13989, 13990, 13991, 13992, 13993, 13994, 13995, 13996]
 +
10 : 2023-08-18 : [14020, 14021, 14022, 14023]
 +
11 : 2023-08-20 : [14069, 14070, 14071, 14072, 14073]
 +
12 : 2023-08-21 : [14093, 14094, 14095, 14096, 14097, 14098, 14099, 14100]
 +
13 : 2023-09-11 : [14258]
 +
14 : 2023-09-16 : [14456]
 +
15 : 2023-09-18 : [14495]
 +
16 : 2023-09-19 : [14590]
 +
17 : 2023-09-20 : [14623, 14624, 14625, 14626]
 +
18 : 2023-09-21 : [14668, 14669]
 +
19 : 2023-09-22 : [14700, 14701, 14702, 14703, 14704]
 +
20 : 2023-09-23 : [14712, 14713, 14715, 14718, 14719, 14720, 14721]
 +
21 : 2023-09-24 : [14731, 14732, 14733, 14736, 14737, 14738]
 +
22 : 2023-10-09 : [15010, 15011, 15012, 15013, 15014, 15015, 15016, 15017]
 +
23 : 2023-10-10 : [15039, 15040]
 +
24 : 2023-10-11 : [15070, 15071, 15072, 15073, 15074]
 +
25 : 2023-10-12 : [15117, 15118, 15119, 15120]
 +
26 : 2023-11-08 : [15321, 15322, 15323, 15324, 15325, 15326, 15327]
 +
27 : 2023-11-10 : [15401, 15402, 15403, 15404, 15405, 15406, 15407]
 +
28 : 2023-11-11 : [15436, 15437, 15438]
 +
29 : 2023-11-12 : [15470, 15471, 15472, 15473]
 +
30 : 2023-11-16 : [15627, 15628, 15629, 15630, 15631]
 +
31 : 2023-11-17 : [15660, 15661, 15662, 15663, 15664]
 +
32 : 2023-11-19 : [15739, 15740, 15741]
 +
33 : 2023-12-01 : [15786, 15787, 15788]
 +
34 : 2023-12-08 : [15889, 15890, 15891, 15892, 15893]
 +
35 : 2024-01-29 : [16423, 16424, 16425, 16426]
 +
36 : 2024-01-30 : [16438, 16439, 16440, 16441, 16442]
 +
37 : 2024-02-02 : [16492, 16493, 16494, 16495, 16496, 16497, 16498, 16499, 16500, 16501, 16502, 16503]
 +
</nowiki>
 +
 
 +
 
 +
(new selection using news data_select notebook, 07.02.2024)
 +
 
 +
Total number of runs: 142 (in 34 nights) out of 246
 +
Total observation time: 38.47 hours out of 65.53 h
 +
 
 +
<nowiki>
 +
20221220 [11522 11528 11529 11530 11531 11532]
 +
20221221 [11546 11547 11548 11549 11550]
 +
20221222 [11566 11568 11569 11570 11571 11572 11573 11574]
 +
20221223 [11595 11596 11597 11598 11599 11600]
 +
20230113 [11618 11619]
 +
20230114 [11641 11642 11643 11644]
 +
20230115 [11660 11661 11662 11663]
 +
20230817 [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
 +
20230911 [14258]
 +
20230916 [14456]
 +
20230918 [14495]
 +
20230919 [14590]
 +
20230920 [14623 14624 14625]
 +
20230921 [14668]
 +
20230922 [14700 14701 14702 14703 14704]
 +
20230923 [14712 14713 14715 14718 14719 14720 14721]
 +
20230924 [14730 14731 14732 14733 14736 14737 14738]
 +
20231006 [14930 14931]
 +
20231007 [14954 14955 14956 14957 14958]
 +
20231010 [15040]
 +
20231012 [15117 15118 15119]
 +
20231108 [15321 15322 15323 15324 15325 15326 15327]
 +
20231110 [15401 15402 15403 15404 15405 15406 15407]
 +
20231111 [15436 15437 15438]
 +
20231112 [15470 15471 15472 15473]
 +
20231113 [15517]
 +
20231116 [15627 15628 15629 15630 15631]
 +
20231117 [15660 15661 15662 15663 15664]
 +
20231119 [15739 15740]
 +
20231201 [15786 15787 15788]
 +
20231208 [15889 15890 15891 15892 15893]
 +
20240129 [16423 16424 16425 16426]
 +
20240130 [16438 16439 16440 16441 16442]
 +
20240202 [16492 16493 16495 16496 16497 16498 16500 16503]
 +
 
 +
</nowiki>
 +
 
 +
=== Analysis details ===
 +
 
 +
Software and analysis notebook:
 +
 
 +
* Files (inc. provence): /fefs/aswg/workspace/mireia.nievas/NGC1275
 +
* Pipeline: /fefs/aswg/workspace/mireia.nievas/NGC1275/NGC1275_to_DL3.ipynb
 +
* lstchain: v0.10.5 originally, v0.10.7 now.
 +
* gammapy: 1.1
 +
* ctapipe: 0.19.2
 +
* astropy: 5.3.4
 +
 
 +
Source:
 +
 
 +
* Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
 +
* Closest (chosen) dec track: 48.22 deg.
 +
* Cleaning: default, tail-cuts 8-4.
 +
* Analysis type: source-independent
 +
* IRFs: point-like, energy-dependent gammaness/theta2 cuts
 +
 
 +
Starting points:
 +
 
 +
* MC: DL2
 +
* Models
 +
* Data: DL1b (merged files with no images)
 +
 
 +
Results
 +
 
 +
==== Analysis details ====
 +
 
 +
* Total significance (per energy bin)
 +
 
 +
[[File:NGC1275_CycleI_AveragePerEnergyBin_mnievas.png|720px]]
 +
 
 +
 
 +
* December flare
 +
 
 +
[[File:NGC1275_CycleI_FlareDecTheta2_mnievas.png|720px]]
 +
 
 +
 
 +
* January flare
 +
 
 +
(note that the highest flux night was with Moon, so not included for now)
 +
 
 +
[[File:NGC1275_CycleI_FlareJanTheta2_mnievas.png|720px]]
 +
 
 +
* Post-flare
 +
 
 +
[[File:NGC1275_CycleI_PostFlareTheta2_mnievas.png|720px]]
 +
 
 +
* Per bayesian block bin (from LC analysis, see below)
 +
 
 +
{| class="wikitable"
 +
|+ Theta2 plots for different bayesian blocks
 +
|-
 +
! BB0 !! BB1
 +
|-
 +
| [[File:NGC1275_CycleI_Theta2_bb0_mnievas.png|320px]]
 +
|| [[File:NGC1275_CycleI_Theta2_bb1_mnievas.png|320px]]
 +
|-
 +
! BB2 !! BB3
 +
|-
 +
| [[File:NGC1275_CycleI_Theta2_bb3_mnievas.png|320px]]
 +
|| [[File:NGC1275_CycleI_Theta2_bb4_mnievas.png|320px]]
 +
|-
 +
! BB4 !! BB5
 +
|-
 +
| [[File:NGC1275_CycleI_Theta2_bb5_mnievas.png|320px]]
 +
|| [[File:NGC1275_CycleI_Theta2_bb6_mnievas.png|320px]]
 +
|-
 +
! BB6 !!
 +
|-
 +
| [[File:NGC1275_CycleI_Theta2_bb7_mnievas.png|320px]]
 +
||
 +
|-
 +
|}
 +
 
 +
==== TS SkyMap ====
 +
 
 +
* Average (all selected runs)
 +
 
 +
[[File:NGC1275_CycleI_TSSkymapAverage_mnievas.png|720px]]
 +
 
 +
* Dec flare
 +
 
 +
[[File:NGC1275_CycleI_TSSkymapDecFlare_mnievas.png|720px]]
 +
 
 +
* Jan flare
 +
 
 +
[[File:NGC1275_CycleI_TSSkymapJanFlare_mnievas.png|720px]]
 +
 
 +
 
 +
* Night of 2023-12-08 to 2023-12-09 (~ 1.72h)
 +
 
 +
Both NGC1275 and IC310 detectable at the same time (!!). This mini-flare from IC310 will be covered in a dedicated wiki.
 +
 
 +
[[File:NGC1275_CycleI_TSSkymapWithIC310_20231208_mnievas.png|720px]]
 +
 
 +
==== Light Curve ====
 +
 
 +
[[File:NGC1275_CycleI_LightCurveWithBB_mnievas.png|720px]]
 +
 
 +
==== SED ====
 +
 
 +
* Average
 +
 
 +
[[File:NGC1275_CycleI_AverageSED_mnievas.png|600px]]
 +
 
 +
* Flux-Binned SED
 +
 
 +
(from the 2nd plot, I show also the previous state for comparison)
 +
 
 +
 
 +
[[File:NGC1275_CycleI_SEDWithBB_mnievas.png|720px]]
 +
 
 +
* December
 +
 
 +
[[File:NGC1275_CycleI_DecemberFlare_mnievas.png|720px]]
 +
 
 +
* January
 +
 
 +
[[File:NGC1275_CycleI_JanuaryFlare_mnievas.png|720px]]
 +
 
 +
Bayesian blocks from the flux in the run-wise Light Curve to define flux bins, then computing one SED per block.
 +
 
 +
* Fermi+LST
 +
 
 +
==== MWL Campaign ====
 +
 
 +
== Jorge ==
 +
 
 +
=== Analysis details ===
 +
 
 +
*Software and analysis notebook:
 +
 
 +
*lstchain: v0.10.5
 +
 
 +
*Closest (chosen) dec track: 48.22 deg.
 +
 
 +
*Cleaning: default, tail-cuts 8-4.
 +
 
 +
*Analysis type: source-independent
 +
 
 +
*IRFs: point-like, interpolated, energy-dependent 70% efficiency gammaness/theta2 cuts
 +
 
 +
*Starting points:
 +
 
 +
MC: /fefs/aswg/data/mc/AllSky/20230901_v0.10.4_allsky_base_prod/TestingDataset/dec_4822/
 +
 
 +
Models: /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_4822/
 +
 
 +
=== Data selection ===
 +
 
 +
New data quality selection notebook used, default settings + min_drdi_at_422pe = 1.2 instead of the default 1.5
 +
 
 +
Total selected time: 36.2 h (for the run list and plots with details on the data quality selection, see Mireia's analysis)
 +
 
 +
=== NGC 1275 ===
 +
 
 +
Coordinates: RA=49.95066663, DEC=+41.5116969
 +
 
 +
==== Theta2 plots ====
 +
 
 +
All data stacked
 +
 
 +
[[File:NGC1275_20221220-20240202_endep_Otero.png|thumb|centre|400px|All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
All data stacked in different energy ranges
 +
 
 +
[[File:NGC1275_20221220-20240202_endep_energybins_Otero.png|thumb|centre|1000px|All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
Theta2 plots in the bayesian blocks defined by Mireia:
 +
{| align='center'
 +
| [[File:NGC1275_MJD_59933.830_59933.834_endep_Otero.png|thumb|centre|400px|BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:NGC1275_MJD_59934.069_59934.125_endep_Otero.png|thumb|centre|400px|BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:NGC1275_MJD_59934.916_59936.991_endep_Otero.png|thumb|centre|400px|BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
{| align='center'
 +
| [[File:NGC1275_MJD_59936.991_59937.015_endep_Otero.png|thumb|centre|400px|BB4, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:NGC1275_MJD_59937.032_59959.916_endep_Otero.png|thumb|centre|400px|BB5, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:NGC1275_MJD_59959.917_60268.153_endep_Otero.png|thumb|centre|400px|BB6, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
{| align='center'
 +
| [[File:NGC1275_MJD_60279.871_60286.988_endep_Otero.png|thumb|centre|400px|BB7, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:NGC1275_MJD_60338.857_60343.023_endep_Otero.png|thumb|centre|400px|BB8, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
 
 +
==== Skymaps ====
 +
 
 +
All data stacked
 +
Skymap centered at the position of NGC 1275
 +
Skymap between 200 GeV and 700 GeV
 +
 
 +
[[File:NGC1275_skymap_20221220-20240202_200-700GeV_endep_Otero.png|thumb|centre|700px|Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
 
 +
Night of 2023-12-08 to 2023-12-09 (~ 1.72h)
 +
Both NGC1275 and IC310 detectable at the same time (!!).
 +
Skymap centered at the position of NGC 1275, the off-centered source is IC 310.
 +
Skymap between 200 GeV and 700 GeV
 +
 
 +
[[File:NGC1275_skymap_MJD_2023-12-08_200-700GeV_endep_Otero.png|thumb|centre|700px|Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
 
 +
==== SEDs ====
 +
 
 +
SED configuration:
 +
*E_reco=[31.6227766,31622.7766] GeV
 +
*E_true=[10,50000] GeV
 +
*n_bins_reco=18
 +
*n_bins_true=50
 +
*E_min_fit=31.6227766 GeV
 +
*E_max_fit=31622.7766 GeV
 +
*Aeff cut > 5%
 +
*1 OFF region
 +
 
 +
SED averaging all data
 +
 
 +
[[File:SED_NGC1275_average_20221220-20240202_endep_Otero.png|thumb|centre|400px|All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
Log-parabola:
 +
-F0=(6.41 土 0.54)・10e-11 TeV-1 cm-2 s-1
 +
-alpha=2.91 土 0.09
 +
-beta=0.21 土 0.07
 +
-E0=300 GeV
 +
 
 +
SED during each bayesian blocks defined by Mireia.
 +
There is clear spectral variability over the different bayesian blocks:
 +
 
 +
BB SED configuration:
 +
*E_reco=[31.6227766,31622.7766] GeV
 +
*E_true=[10,50000] GeV
 +
*n_bins_reco=18
 +
*n_bins_true=50
 +
*E_min_fit=31.6227766 GeV
 +
*E_max_fit=31622.7766 GeV
 +
*Aeff cut > 5%
 +
*1 OFF region
 +
 
 +
These SEDs are observed, using a log-parabola spectral model.
 +
 
 +
{| align='center'
 +
| [[File:SED_NGC1275_MJD_59933.830_59933.834_endep_Otero.png|thumb|centre|400px|BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_NGC1275_MJD_59934.069_59934.125_endep_Otero.png|thumb|centre|400px|BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_NGC1275_MJD_59934.916_59936.991_endep_Otero.png|thumb|centre|400px|BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
{| align='center'
 +
| [[File:SED_NGC1275_MJD_59936.991_59937.015_endep_Otero.png|thumb|centre|400px|BB4, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_NGC1275_MJD_59937.032_59959.916_endep_Otero.png|thumb|centre|400px|BB5, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_NGC1275_MJD_59959.917_60268.153_endep_Otero.png|thumb|centre|400px|BB6, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
{| align='center'
 +
| [[File:SED_NGC1275_MJD_60279.871_60286.988_endep_Otero.png|thumb|centre|400px|BB7, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_NGC1275_MJD_60338.857_60343.023_endep_Otero.png|thumb|centre|400px|BB8, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
 
 +
==== Lightcurve ====
 +
 
 +
NGC 1275 lightcurve above 100 GeV
 +
This lightcurve was obtained using as reference the average spectrum over all the studied period, presented above in this wiki
 +
 
 +
{| align='center'
 +
| [[File:LC_NGC1275_complete_nightwise_endep_Otero.png|thumb|centre|400px|Nightwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:LC_NGC1275_complete_runwise_endep_Otero.png|thumb|centre|400px|Runwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
 
 +
=== IC 310 ===
 +
 
 +
Coordinates: RA=49.179078 deg, DEC=+41.32489 deg
 +
 
 +
The data taken centered at NGC 1275 are not used for the analysis of IC 310 due to some OFF regions being located on top of NGC 1275. Therefore, only those with the Perseus-MA pointing are used.
 +
 
 +
==== Theta2 plots ====
 +
 
 +
All data stacked
 +
 
 +
[[File:IC310_20230817-20240202_endep_Otero.png|thumb|centre|400px|All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
All data stacked in different energy ranges
 +
 
 +
[[File:IC310_20230817-20240202_endep_energybins_Otero.png|thumb|centre|1000px|All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
Most of the signal is concentrated below ~800 GeV.
 +
There is a hint of a signal stacking all data, but we do not reach detection. However, we do see a clear detection on the night of 2023/12/08:
 +
 
 +
[[File:IC310_20231208_energybins_Otero.png|thumb|centre|400px|All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
[[File:IC310_20231208_endep_energybins_Otero.png|thumb|centre|1000px|All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
 
 +
The signal extends until ~3-4 TeV. Probably optimizing a bit the energy range of the theta2 plot, the significance is higher, but the
 +
detection is clear and solid.
 +
 
 +
==== Skymaps ====
 +
 
 +
All data stacked
 +
Skymap centered at the position of IC 310
 +
Skymap between 200 GeV and 700 GeV
 +
No detection or hint of IC 310, but NGC 1275 is visible away from the center, as seen from the previous skymaps centered at that position.
 +
 
 +
[[File:IC310_skymap_20221220-20240202_200-700GeV_endep_Otero.png|thumb|centre|700px|Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
 
 +
Night of 2023-12-08 to 2023-12-09 (~ 1.72h)
 +
Both NGC1275 and IC310 detectable at the same time (!!).
 +
Skymap centered at the position of IC 310, the off-centered source is NGC 1275.
 +
Skymap between 200 GeV and 700 GeV
 +
 
 +
[[File:IC310_skymap_MJD_2023-12-08_200-700GeV_endep_Otero.png|thumb|centre|700px|Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
 
 +
==== SEDs ====
 +
 
 +
SED configuration:
 +
*E_reco=[31.6227766,31622.7766] GeV
 +
*E_true=[10,50000] GeV
 +
*n_bins_reco=18
 +
*n_bins_true=50
 +
*E_min_fit=31.6227766 GeV
 +
*E_max_fit=31622.7766 GeV
 +
*Aeff cut > 5%
 +
*1 OFF region
 +
 
 +
SED averaging all data
 +
 
 +
[[File:SED_IC310_average_20221220-20240202_endep_Otero.png|thumb|centre|400px|All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]] 
 +
 
 +
Log-parabola:
 +
-F0=(2.61 土 0.57)・10e-11 TeV-1 cm-2 s-1
 +
-alpha=2.44 土 0.49
 +
-beta=0.44 土 0.31
 +
-E0=300 GeV
 +
 
 +
SED during each bayesian blocks defined by Mireia.
 +
There is clear spectral variability over the different bayesian blocks:
 +
 
 +
BB SED configuration:
 +
*E_reco=[31.6227766,31622.7766] GeV
 +
*E_true=[10,50000] GeV
 +
*n_bins_reco=18
 +
*n_bins_true=50
 +
*E_min_fit=31.6227766 GeV
 +
*E_max_fit=31622.7766 GeV
 +
*Aeff cut > 5%
 +
*1 OFF region
 +
 
 +
These SEDs are observed, using a log-parabola spectral model.
 +
 
 +
{| align='center'
 +
| [[File:SED_IC310_BB1_endep_Otero.png|thumb|centre|400px|BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_IC310_BB2_endep_Otero.png|thumb|centre|400px|BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:SED_IC310_BB3_endep_Otero.png|thumb|centre|400px|BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
 
 +
==== Lightcurve ====
 +
 
 +
IC310 lightcurve above 100 GeV
 +
This lightcurve was obtained using as reference the average spectrum over all the studied period, presented above in this wiki
 +
 
 +
{| align='center'
 +
| [[File:LC_IC310_complete_nightwise_endep_Otero.png|thumb|centre|400px|Nightwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
| [[File:LC_IC310_complete_runwise_endep_Otero.png|thumb|centre|400px|Runwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs]]
 +
|}
 +
 
 +
== Tora ==
 +
 
 +
* Producer: Tora Therese Høiland Arnesen <tora.arnesen@iac.es>
 +
 
 +
=== Tora: NGC 1275 ===
 +
 
 +
==== Tora NGC1275: Data acquisition ====
 +
 
 +
Dates up to: 2024-06-01 (last obs 2024-02-10)
 +
 
 +
Selection based on Abelardo's data_quality notebook
 +
 
 +
* File: /fefs/aswg/workspace/tora.arnesen/perseus/NGC1275/NGC1275_data_quality.ipynb
 +
* astropy: 5.3.4
 +
 
 +
[[File:Parameters_NGC1275_selection Large.png|720px]]
 +
 
 +
[[File:RelLightYield_NGC1275_selection.png|720px]]
 +
 
 +
[[File:Nsubruns_NGC1275_selection.png|720px]]
 +
 
 +
[[File:CosmicsTriggerRate_NGC1275_selection.png|720px]]
 +
 
 +
[[File:NGC1275SelectionRuns.png|720px]]
 +
 
 +
==== Tora NGC1275: Data selection ====
 +
 
 +
Total number of runs: 153 out of 238
 +
 
 +
Total number of nights: 37 out of 50
 +
 
 +
Total observation time (hours): 41.02 out of 63.88
 +
 
 +
<nowiki>
 +
20221220 : [11522 11528 11529 11530 11531 11532]
 +
20221221 : [11546 11547 11548 11549 11550]
 +
20221222 : [11566 11568 11569 11570 11571 11572 11573 11574]
 +
20221223 : [11595 11596 11597 11598 11599 11600]
 +
20230113 : [11618 11619]
 +
20230114 : [11641 11642 11643 11644]
 +
20230115 : [11660 11661 11662 11663]
 +
20230817 : [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
 +
20230911 : [14258]
 +
20230916 : [14456]
 +
20230918 : [14495]
 +
20230919 : [14590]
 +
20230920 : [14623 14624 14625]
 +
20230921 : [14668]
 +
20230922 : [14700 14701 14702 14703 14704]
 +
20230923 : [14712 14713 14715 14718 14719 14720 14721]
 +
20230924 : [14730 14731 14732 14733 14736 14737 14738]
 +
20231006 : [14930 14931]
 +
20231007 : [14954 14955 14956 14957 14958]
 +
20231010 : [15040]
 +
20231012 : [15117 15118 15119]
 +
20231108 : [15321 15322 15323 15324 15325 15326 15327]
 +
20231110 : [15401 15402 15403 15404 15405 15406 15407]
 +
20231111 : [15436 15437 15438]
 +
20231112 : [15470 15471 15472 15473]
 +
20231113 : [15517]
 +
20231116 : [15627 15628 15629 15630 15631]
 +
20231117 : [15660 15661 15662 15663 15664]
 +
20231119 : [15739 15740]
 +
20231201 : [15786 15787 15788]
 +
20231208 : [15889 15890 15891 15892 15893]
 +
20240129 : [16423 16424 16425 16426]
 +
20240130 : [16438 16439 16440 16441 16442]
 +
20240202 : [16492 16493 16495 16496 16497 16498 16500 16503]
 +
20240207 : [16618]
 +
20240208 : [16650 16651 16652 16653 16654 16655]
 +
20240210 : [16664 16665 16666 16667]
 +
</nowiki>
 +
 
 +
==== Tora NGC1275: Analysis details ====
 +
 
 +
* Directory: /fefs/aswg/workspace/tora.arnesen/perseus/NGC1275
 +
* Pipeline:
 +
** ./launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822/lstchain_config_2024-01-31.json
 +
** ./launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_4822 ./irf_dl3_tool_config.json gh_eff_0.7_th_cont_0.7_20240417_TT
 +
** ./launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_4822/gh_eff_0.7_th_cont_0.7_20240417_TT NGC\ 1275
 +
** lstchain_create_dl3_index_files -d ./DL3 -p dl3\*.fits.gz
 +
* lstchain: v0.10.9
 +
 
 +
Source:
 +
 
 +
* Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
 +
* Closest (chosen) dec track: 48.22 deg.
 +
* Cleaning: default, tail-cuts 8-4.
 +
* Analysis type: source-independent
 +
* IRFs: point-like, energy-dependent gammaness/theta2 cuts (0.7/0.7)
 +
 
 +
Starting points:
 +
 
 +
* MC: DL2
 +
* Models
 +
* Data: DL1b (merged files with no images)
 +
 
 +
 
 +
==== Tora NGC1275: Theta2 plots ====
 +
 
 +
* Total significance (per energy bin)
 +
 
 +
[[File:Theta2_total_TT_all_NGC1275.png|thumb|centre|1000px|All data, energy binned]]
 +
 
 +
* Per bayesian block bin (from LC analysis, see below)
 +
 
 +
{| class="wikitable"
 +
|+ Theta2 plots for different bayesian blocks
 +
|-
 +
! BB0 !! BB1
 +
|-
 +
| [[File:TT_all_NGC1275_theta2_excess_plots_BB0.png|600px]]
 +
|| [[File:TT_all_NGC1275_theta2_excess_plots_BB1.png|600px]]
 +
|-
 +
! BB2 !! BB3
 +
|-
 +
| [[File:TT_all_NGC1275_theta2_excess_plots_BB2.png|600px]]
 +
|| [[File:TT_all_NGC1275_theta2_excess_plots_BB3.png|600px]]
 +
|-
 +
! BB4 !! BB5
 +
|-
 +
| [[File:TT_all_NGC1275_theta2_excess_plots_BB4.png|600px]]
 +
|| [[File:TT_all_NGC1275_theta2_excess_plots_BB5.png|600px]]
 +
|-
 +
! BB6 !! BB7
 +
|-
 +
| [[File:TT_all_NGC1275_theta2_excess_plots_BB6.png|600px]]
 +
|| [[File:TT_all_NGC1275_theta2_excess_plots_BB7.png|600px]]
 +
|-
 +
|}
 +
 
 +
==== Tora NGC1275: Lightcurves ====
 +
 
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
 
 +
[[File:LCs_NGC1275_all_TT_LST1 Large.png|thumb|centre|700px|Lighcurve LP runwise (top) and nightwise (bottom), compared to MAGIC average flux for Crab and NGC1275]]
 +
 
 +
[[File:LC_BB_LST_NGC1275_all_TT.png|thumb|centre|720px|Lighcurve LP with runwise, night-wise, and bayesian block-wise flux points]]
 +
 
 +
[[File:LCs_all_models_NGC1275_all_TT_LST1 Medium.png|thumb|centre|700px|Lighcurves runwise (top) and nightwise (bottom), for Power-Law (PL), Log-Parabola (LP), and Exponential cutoff Power-Law (EPWL). Compared to MAGIC average flux for Crab and NGC1275]]
 +
 
 +
==== Tora NGC1275: SEDs ====
 +
 
 +
 
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
 
 +
[[File:SED_LP_NGC1275_all_TT_LST1.png|thumb|centre|720px|All data, LP model, compared with Magic 2018 LP model]] 
 +
 
 +
Log-parabola spectral model (no EBL correction):
 +
-F0 = (5.3 土 0.5)・10e-11 TeV-1 cm-2 s-1
 +
-alpha = 2.8 土 0.1
 +
-beta = 0.13 土 0.07
 +
-E0 = 0.3 TeV
 +
 
 +
[[File:SEDs_allBB_LST_NGC1275_all_TT.png|thumb|centre|720px|SED per bayesian block]]
 +
 
 +
[[File:SED_hatch_all_models_NGC1275_all_TT_LST1.png|thumb|centre|700px| SED all data, for Power-Law (PL), Log-Parabola (LP), and Exponential cutoff Power-Law (EPWL). Compared with Magic 2018 LP model and EPWL model.]]
 +
 
 +
Poer-Law spectral model (no EBL correction):
 +
-F0 = (4.6 土 0.3)・10e-11 TeV-1 cm-2 s-1
 +
-index = 2.8 土 0.1
 +
-E0 = 0.3 TeV
 +
 
 +
Exponential cutoff Power-Law spectral model (no EBL correction):
 +
-F0 = (6 土 1)・10e-11 TeV-1 cm-2 s-1
 +
-alpha = 1 (frozen)
 +
-index = 2.6 土 0.2
 +
-lambda = 0.5 土 0.4  TeV-1
 +
-E0 = 0.3 TeV
 +
 
 +
==== Tora NGC1275: Skymaps ====
 +
 
 +
All data stacked
 +
Skymaps centered at the position of NGC 1275
 +
Skymaps between 0.2 TeV and 7 TeV
 +
 
 +
[[File:Significance_Off_and_Excess_maps_all_NGC1275_subIC310.png|thumb|centre|700px|| Significance, off significance, and excess maps]]
 +
 
 +
[[File:Background_residuals_map_all_NGC1275_subIC310.png|thumb|centre|700px| Background residuals map]]
 +
 
 +
=== NGC 1275 Cross-check SED ===
 +
 
 +
Dates up to: 2024-02-02
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
[[File:SED_comparison_NGC1275_TT_LST1.png|thumb|centre|720px|All data, LP model, compared with Magic 2018 LP models]]
 +
 
 +
=== Tora: IC 310 ===
 +
 
 +
==== Tora IC310: Data acquisition ====
 +
 
 +
Dates up to: 2024-05-28 (last obs 2024-02-10)
 +
 
 +
Selection based on Abelardo's data_quality notebook
 +
 
 +
* File: /fefs/aswg/workspace/tora.arnesen/perseus/IC310/IC310_data_quality.ipynb
 +
* astropy: 5.3.4
 +
 
 +
[[File:Parameters_IC310_selection Large.png|720px]]
 +
 
 +
[[File:RelLightYield_IC310_selection.png|720px]]
 +
 
 +
[[File:Nsubruns_IC310_selection.png|720px]]
 +
 
 +
[[File:CosmicsTriggerRate_IC310_selection.png|720px]]
 +
 
 +
[[File:IC310SelectionRuns.png|720px]]
 +
 
 +
==== Tora IC310: Data selection ====
 +
 
 +
Total number of runs: 133 out of 217
 +
 
 +
Total number of nights: 37 out of 49
 +
 
 +
Total observation time (hours): 36.95 out of 59.57
 +
 
 +
<nowiki>
 +
20221220 : [11529 11530]
 +
20221221 : [11548 11549]
 +
20221222 : [11568 11569 11573 11574]
 +
20221223 : [11598 11599]
 +
20230113 : [11619]
 +
20230114 : [11641 11644]
 +
20230115 : [11660 11663]
 +
20230817 : [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
 +
20230911 : [14258]
 +
20230916 : [14456]
 +
20230918 : [14495]
 +
20230919 : [14590]
 +
20230920 : [14623 14624 14625]
 +
20230921 : [14668]
 +
20230922 : [14700 14701 14702 14703 14704]
 +
20230923 : [14712 14713 14715 14718 14719 14720 14721]
 +
20230924 : [14730 14731 14732 14733 14736 14737 14738]
 +
20231006 : [14930 14931]
 +
20231007 : [14954 14955 14956 14957 14958]
 +
20231010 : [15040]
 +
20231012 : [15117 15118 15119]
 +
20231108 : [15321 15322 15323 15324 15325 15326 15327]
 +
20231110 : [15401 15402 15403 15404 15405 15406 15407]
 +
20231111 : [15436 15437 15438]
 +
20231112 : [15470 15471 15472 15473]
 +
20231113 : [15517]
 +
20231116 : [15627 15628 15629 15630 15631]
 +
20231117 : [15660 15661 15662 15663 15664]
 +
20231119 : [15739 15740]
 +
20231201 : [15786 15787 15788]
 +
20231208 : [15889 15890 15891 15892 15893]
 +
20240129 : [16423 16424 16425 16426]
 +
20240130 : [16438 16439 16440 16441 16442]
 +
20240202 : [16492 16493 16495 16496 16497 16498 16500 16503]
 +
20240207 : [16618]
 +
20240208 : [16650 16651 16652 16653 16654 16655]
 +
20240210 : [16664 16665 16666 16667]
 +
</nowiki>
 +
 
 +
==== Tora IC310: Analysis details ====
 +
 
 +
* Directory: /fefs/aswg/workspace/tora.arnesen/perseus/IC310
 +
* Pipeline:
 +
** ./IC310_launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822/lstchain_config_2024-01-31.json
 +
** ./IC310_launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_4822 ./irf_dl3_tool_config.json gh_eff_0.7_th_cont_0.7_20240510_TT
 +
** ./IC310_launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_4822/gh_eff_0.7_th_cont_0.7_20240510_TT IC\ 310
 +
** lstchain_create_dl3_index_files -d ./DL3 -p dl3\*.fits.gz
 +
* lstchain: v0.10.9
 +
* gammapy: 1.2 (SED & LCs) / 1.1 (Theta2 & Skymaps)
 +
 
 +
Source:
 +
 
 +
* Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
 +
* Closest (chosen) dec track: 48.22 deg.
 +
* Cleaning: default, tail-cuts 8-4.
 +
* Analysis type: source-independent
 +
* IRFs: point-like, energy-dependent gammaness/theta2 cuts (0.7/0.7)
 +
 
 +
Starting points:
 +
 
 +
* MC: DL2
 +
* Models
 +
* Data: DL1b (merged files with no images)
 +
 
 +
 
 +
==== Tora IC310: Theta2 plots ====
 +
 
 +
* Total significance (per energy bin)
 +
 
 +
[[File:Theta2_total_TT_all_IC310 Medium.png|thumb|centre|1000px|All data, energy binned]]
 +
 
 +
* Per bayesian block bin (from LC analysis, see below)
 +
 
 +
{| class="wikitable"
 +
|+ Theta2 plots for different bayesian blocks
 +
|-
 +
! BB0 !! BB1
 +
|-
 +
| [[File:TT_all_IC310_theta2_excess_plots_BB0.png|600px]]
 +
|| [[File:TT_all_IC310_theta2_excess_plots_BB1.png|600px]]
 +
|-
 +
! BB2 !!
 +
|-
 +
| [[File:TT_all_IC310_theta2_excess_plots_BB2.png|600px]]
 +
||
 +
|-
 +
|}
 +
 
 +
==== Tora IC310: Lightcurves ====
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
[[File:LCs_IC310_all_TT_LST1 Large.png|thumb|centre|700px|Lighcurve LP runwise (top) and nightwise (bottom), compared to MAGIC average flux for Crab and NGC1275]]
 +
 
 +
[[File:LC_BB_LST_IC310_all_TT.png|thumb|centre|720px|Lighcurve LP with runwise, night-wise, and bayesian block-wise flux points]]
 +
 
 +
==== Tora IC310: SEDs ====
 +
 
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
[[File:SED_LP_IC310_all_TT_LST1.png|thumb|centre|720px|All data, LP model, compared with Magic 2018 LP model]] 
 +
 
 +
Log-parabola spectral model (no EBL correction):
 +
-F0 = (2.9 土 0.6)・10e-11 TeV-1 cm-2 s-1
 +
-alpha = 2.6 土 0.3
 +
-beta = 0.5 土 0.3
 +
-E0 = 0.3 TeV
 +
 
 +
[[File:SEDs_allBB_LST_IC310_all_TT.png|thumb|centre|720px|SED per bayesian block]]
 +
 
 +
==== Tora IC310: Skymaps ====
 +
 
 +
All data stacked
 +
Skymaps centered at the position of IC 310
 +
Skymaps between 0.2 TeV and 7 TeV
 +
 
 +
[[File:Significance_Off_and_Excess_maps_all_IC310_subNGC1275.png|thumb|centre|700px| Significance, off significance, and excess maps]]
 +
 
 +
[[File:Background_residuals_map_all_IC310_subNGC1275.png|thumb|centre|700px| Background residuals map]]
 +
 
 +
=== IC 310 Cross-check SED ===
 +
 
 +
Last date: 2024-02-02
 +
 
 +
SED & LC configuration:
 +
* e_reco_min = 0.0316227766 TeV
 +
* e_reco_max = 31.6227766 TeV
 +
* e_true_min = 0.01 TeV
 +
* e_true_max = 50 TeV
 +
* e_reco_bin = 18
 +
* e_true_bin = 50
 +
* per_decade = False
 +
* e_fit_min = 0.0316227766 TeV
 +
* e_fit_max = 31.6227766 TeV
 +
* e_lc_min = 0.09999999999467535 TeV
 +
* e_lc_max = 49.99999999999999 TeV
 +
* Aeff cut > 5%
 +
* safe_min_energy = 0.04 TeV
 +
* safe_max_energy = 7 TeV
 +
* 1 OFF region
 +
 
 +
[[File:SED_comparison_IC310_TT_LST1.png|thumb|centre|720px|All data up to 2024-02-02, LP models, compared with Tora's NGC1275 LP model]]
 +
 
 +
=== IC 310 dec line comparison ===
 +
 
 +
Comparison analysis using dec 34.76 and dec 48.22
 +
 
 +
[[File:LCs_dec_comparison_IC310_3476_all_TT_LST1 Large.png|thumb|centre|720px| Light curves for LP models using dec 34.76 and dec 48.22, compared with Tora's NGC1275 LP model]]
 +
 
 +
[[File:SED_dec_comparison_IC310_3476_all_TT_LST1.png|thumb|centre|720px| SED for LP models using dec 34.76 and dec 48.22, compared with Tora's NGC1275 LP model]]

Latest revision as of 15:41, 11 June 2024

General information[edit]

  • Name of the source: NGC 1275
  • Brief description of the source:
- Object type : FR-I Radio Galaxy + Spiral galaxy (in interaction)
- RA: 3 19 48.16 (hh mm ss), Dec: +41 30 42.1088 (dd mm ss)

NGC 1275 data analysis[edit]

Related analysis pages[edit]

Analysis_NGC1275_Dec_2022_w_IRF_interpolation_Priyadarshi (Dec 2022, source-independent, IRF interpolation, energy-dependent cuts. Dec 20 flare for ATel https://www.astronomerstelegram.org/?read=15819 )

Analysis_NGC1275_Dec_2022_source_dep_Nozaki (Dec 20 Flare, source-dependent, IRF interpolation, energy-dependent cuts for ATel https://www.astronomerstelegram.org/?read=15819 )

Cycle I (including major flare Dec 2022)[edit]

Mireia[edit]

  • Producer: Mireia Nievas Rosillo <mnievas@iac.es>

Data acquisition[edit]

old datacheck notebook[edit]

(below some key points from Abelardo's old datacheck selection notebook, only showing those relevant for data selection)

  • ZD and AZ distributions

NGC1275 CycleI ZD AZ hist mnievas.png

  • NSB vs moon altitude and vs run number

NGC1275 CycleI MoonConditions mnievas.png

NGC1275 CycleI MoonConditions2 mnievas.png

  • Rate of cosmics

NGC1275 CycleI RateOfCosmics mnievas.png

  • Rate of pulses above a given pulse size

NGC1275 CycleI RateAboveCharge mnievas.png

  • NSB after all cuts

NGC1275 CycleI NSB aftercuts mnievas.png

new data_select notebook[edit]

NGC1275 CycleI ds1 mnievas.png

NGC1275 CycleI ds2 mnievas.png

NGC1275 CycleI ds4 mnievas.png

NGC1275 CycleI ds5 mnievas.png

Data selection[edit]

(old selection using Abelardo's old datacheck notebook)

  • 204 wobble runs for the selected source (out of 246)
  • Observation time: 53.56 hours (out of 65.53 hours)
  • (most data in Perseus-MA mode, shared with IC310, some runs in NGC1275 pointing mode)
  • (moon data is removed)
1 : 2022-12-20 : [11522, 11528, 11529, 11530, 11531, 11532, 11522, 11528, 11529, 11530, 11531, 11532]
2 : 2022-12-21 : [11546, 11547, 11548, 11549, 11550]
3 : 2022-12-22 : [11566, 11568, 11569, 11570, 11571, 11572, 11573, 11574, 11566, 11568, 11569, 11570, 11571, 11572, 11573, 11574]
4 : 2022-12-23 : [11595, 11596, 11597, 11598, 11599, 11600]
5 : 2023-01-13 : [11618, 11619]
6 : 2023-01-14 : [11641, 11642, 11643, 11644, 11641, 11642, 11643, 11644]
7 : 2023-01-15 : [11660, 11661, 11662, 11663, 11660, 11661, 11662, 11663]
8 : 2023-08-15 : [13935, 13936, 13937, 13938, 13939, 13940, 13941, 13942, 13943, 13944]
9 : 2023-08-17 : [13987, 13988, 13989, 13990, 13991, 13992, 13993, 13994, 13995, 13996]
10 : 2023-08-18 : [14020, 14021, 14022, 14023]
11 : 2023-08-20 : [14069, 14070, 14071, 14072, 14073]
12 : 2023-08-21 : [14093, 14094, 14095, 14096, 14097, 14098, 14099, 14100]
13 : 2023-09-11 : [14258]
14 : 2023-09-16 : [14456]
15 : 2023-09-18 : [14495]
16 : 2023-09-19 : [14590]
17 : 2023-09-20 : [14623, 14624, 14625, 14626]
18 : 2023-09-21 : [14668, 14669]
19 : 2023-09-22 : [14700, 14701, 14702, 14703, 14704]
20 : 2023-09-23 : [14712, 14713, 14715, 14718, 14719, 14720, 14721]
21 : 2023-09-24 : [14731, 14732, 14733, 14736, 14737, 14738]
22 : 2023-10-09 : [15010, 15011, 15012, 15013, 15014, 15015, 15016, 15017]
23 : 2023-10-10 : [15039, 15040]
24 : 2023-10-11 : [15070, 15071, 15072, 15073, 15074]
25 : 2023-10-12 : [15117, 15118, 15119, 15120]
26 : 2023-11-08 : [15321, 15322, 15323, 15324, 15325, 15326, 15327]
27 : 2023-11-10 : [15401, 15402, 15403, 15404, 15405, 15406, 15407]
28 : 2023-11-11 : [15436, 15437, 15438]
29 : 2023-11-12 : [15470, 15471, 15472, 15473]
30 : 2023-11-16 : [15627, 15628, 15629, 15630, 15631]
31 : 2023-11-17 : [15660, 15661, 15662, 15663, 15664]
32 : 2023-11-19 : [15739, 15740, 15741]
33 : 2023-12-01 : [15786, 15787, 15788]
34 : 2023-12-08 : [15889, 15890, 15891, 15892, 15893]
35 : 2024-01-29 : [16423, 16424, 16425, 16426]
36 : 2024-01-30 : [16438, 16439, 16440, 16441, 16442]
37 : 2024-02-02 : [16492, 16493, 16494, 16495, 16496, 16497, 16498, 16499, 16500, 16501, 16502, 16503]


(new selection using news data_select notebook, 07.02.2024)

Total number of runs: 142 (in 34 nights) out of 246 Total observation time: 38.47 hours out of 65.53 h 

20221220 [11522 11528 11529 11530 11531 11532]
20221221 [11546 11547 11548 11549 11550]
20221222 [11566 11568 11569 11570 11571 11572 11573 11574]
20221223 [11595 11596 11597 11598 11599 11600]
20230113 [11618 11619]
20230114 [11641 11642 11643 11644]
20230115 [11660 11661 11662 11663]
20230817 [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
20230911 [14258]
20230916 [14456]
20230918 [14495]
20230919 [14590]
20230920 [14623 14624 14625]
20230921 [14668]
20230922 [14700 14701 14702 14703 14704]
20230923 [14712 14713 14715 14718 14719 14720 14721]
20230924 [14730 14731 14732 14733 14736 14737 14738]
20231006 [14930 14931]
20231007 [14954 14955 14956 14957 14958]
20231010 [15040]
20231012 [15117 15118 15119]
20231108 [15321 15322 15323 15324 15325 15326 15327]
20231110 [15401 15402 15403 15404 15405 15406 15407]
20231111 [15436 15437 15438]
20231112 [15470 15471 15472 15473]
20231113 [15517]
20231116 [15627 15628 15629 15630 15631]
20231117 [15660 15661 15662 15663 15664]
20231119 [15739 15740]
20231201 [15786 15787 15788]
20231208 [15889 15890 15891 15892 15893]
20240129 [16423 16424 16425 16426]
20240130 [16438 16439 16440 16441 16442]
20240202 [16492 16493 16495 16496 16497 16498 16500 16503]

Analysis details[edit]

Software and analysis notebook:

  • Files (inc. provence): /fefs/aswg/workspace/mireia.nievas/NGC1275
  • Pipeline: /fefs/aswg/workspace/mireia.nievas/NGC1275/NGC1275_to_DL3.ipynb
  • lstchain: v0.10.5 originally, v0.10.7 now.
  • gammapy: 1.1
  • ctapipe: 0.19.2
  • astropy: 5.3.4

Source:

  • Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
  • Closest (chosen) dec track: 48.22 deg.
  • Cleaning: default, tail-cuts 8-4.
  • Analysis type: source-independent
  • IRFs: point-like, energy-dependent gammaness/theta2 cuts

Starting points:

  • MC: DL2
  • Models
  • Data: DL1b (merged files with no images)

Results

Analysis details[edit]

  • Total significance (per energy bin)

NGC1275 CycleI AveragePerEnergyBin mnievas.png


  • December flare

720px


  • January flare

(note that the highest flux night was with Moon, so not included for now)

720px

  • Post-flare

720px

  • Per bayesian block bin (from LC analysis, see below)
Theta2 plots for different bayesian blocks
BB0 BB1
NGC1275 CycleI Theta2 bb0 mnievas.png NGC1275 CycleI Theta2 bb1 mnievas.png
BB2 BB3
NGC1275 CycleI Theta2 bb3 mnievas.png NGC1275 CycleI Theta2 bb4 mnievas.png
BB4 BB5
NGC1275 CycleI Theta2 bb5 mnievas.png NGC1275 CycleI Theta2 bb6 mnievas.png
BB6
NGC1275 CycleI Theta2 bb7 mnievas.png

TS SkyMap[edit]

  • Average (all selected runs)

NGC1275 CycleI TSSkymapAverage mnievas.png

  • Dec flare

NGC1275 CycleI TSSkymapDecFlare mnievas.png

  • Jan flare

NGC1275 CycleI TSSkymapJanFlare mnievas.png


  • Night of 2023-12-08 to 2023-12-09 (~ 1.72h)

Both NGC1275 and IC310 detectable at the same time (!!). This mini-flare from IC310 will be covered in a dedicated wiki.

NGC1275 CycleI TSSkymapWithIC310 20231208 mnievas.png

Light Curve[edit]

NGC1275 CycleI LightCurveWithBB mnievas.png

SED[edit]

  • Average

NGC1275 CycleI AverageSED mnievas.png

  • Flux-Binned SED

(from the 2nd plot, I show also the previous state for comparison)


NGC1275 CycleI SEDWithBB mnievas.png

  • December

720px

  • January

720px

Bayesian blocks from the flux in the run-wise Light Curve to define flux bins, then computing one SED per block.

  • Fermi+LST

MWL Campaign[edit]

Jorge[edit]

Analysis details[edit]

  • Software and analysis notebook:
  • lstchain: v0.10.5
  • Closest (chosen) dec track: 48.22 deg.
  • Cleaning: default, tail-cuts 8-4.
  • Analysis type: source-independent
  • IRFs: point-like, interpolated, energy-dependent 70% efficiency gammaness/theta2 cuts
  • Starting points:
MC: /fefs/aswg/data/mc/AllSky/20230901_v0.10.4_allsky_base_prod/TestingDataset/dec_4822/

Models: /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_4822/

Data selection[edit]

New data quality selection notebook used, default settings + min_drdi_at_422pe = 1.2 instead of the default 1.5

Total selected time: 36.2 h (for the run list and plots with details on the data quality selection, see Mireia's analysis)

NGC 1275[edit]

Coordinates: RA=49.95066663, DEC=+41.5116969

Theta2 plots[edit]

All data stacked
All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
All data stacked in different energy ranges
All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Theta2 plots in the bayesian blocks defined by Mireia:
BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB4, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB5, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB6, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB7, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB8, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

Skymaps[edit]

All data stacked
Skymap centered at the position of NGC 1275
Skymap between 200 GeV and 700 GeV
Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Night of 2023-12-08 to 2023-12-09 (~ 1.72h)
Both NGC1275 and IC310 detectable at the same time (!!).
Skymap centered at the position of NGC 1275, the off-centered source is IC 310.
Skymap between 200 GeV and 700 GeV
Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

SEDs[edit]

SED configuration:

  • E_reco=[31.6227766,31622.7766] GeV
  • E_true=[10,50000] GeV
  • n_bins_reco=18
  • n_bins_true=50
  • E_min_fit=31.6227766 GeV
  • E_max_fit=31622.7766 GeV
  • Aeff cut > 5%
  • 1 OFF region
SED averaging all data
All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Log-parabola:
-F0=(6.41 土 0.54)・10e-11 TeV-1 cm-2 s-1
-alpha=2.91 土 0.09
-beta=0.21 土 0.07
-E0=300 GeV

SED during each bayesian blocks defined by Mireia.
There is clear spectral variability over the different bayesian blocks:

BB SED configuration:

  • E_reco=[31.6227766,31622.7766] GeV
  • E_true=[10,50000] GeV
  • n_bins_reco=18
  • n_bins_true=50
  • E_min_fit=31.6227766 GeV
  • E_max_fit=31622.7766 GeV
  • Aeff cut > 5%
  • 1 OFF region
These SEDs are observed, using a log-parabola spectral model.
BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB4, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB5, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB6, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB7, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB8, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

Lightcurve[edit]

NGC 1275 lightcurve above 100 GeV
This lightcurve was obtained using as reference the average spectrum over all the studied period, presented above in this wiki
Nightwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Runwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

IC 310[edit]

Coordinates: RA=49.179078 deg, DEC=+41.32489 deg

The data taken centered at NGC 1275 are not used for the analysis of IC 310 due to some OFF regions being located on top of NGC 1275. Therefore, only those with the Perseus-MA pointing are used.

Theta2 plots[edit]

All data stacked
All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
All data stacked in different energy ranges
All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Most of the signal is concentrated below ~800 GeV.
There is a hint of a signal stacking all data, but we do not reach detection. However, we do see a clear detection on the night of 2023/12/08:
All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
All data, energy binned, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
The signal extends until ~3-4 TeV. Probably optimizing a bit the energy range of the theta2 plot, the significance is higher, but the 
detection is clear and solid.

Skymaps[edit]

All data stacked
Skymap centered at the position of IC 310
Skymap between 200 GeV and 700 GeV
No detection or hint of IC 310, but NGC 1275 is visible away from the center, as seen from the previous skymaps centered at that position.
Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Night of 2023-12-08 to 2023-12-09 (~ 1.72h)
Both NGC1275 and IC310 detectable at the same time (!!).
Skymap centered at the position of IC 310, the off-centered source is NGC 1275.
Skymap between 200 GeV and 700 GeV
Skymap between 200 GeV and 700 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

SEDs[edit]

SED configuration:

  • E_reco=[31.6227766,31622.7766] GeV
  • E_true=[10,50000] GeV
  • n_bins_reco=18
  • n_bins_true=50
  • E_min_fit=31.6227766 GeV
  • E_max_fit=31622.7766 GeV
  • Aeff cut > 5%
  • 1 OFF region
SED averaging all data
All data, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Log-parabola:
-F0=(2.61 土 0.57)・10e-11 TeV-1 cm-2 s-1
-alpha=2.44 土 0.49
-beta=0.44 土 0.31
-E0=300 GeV

SED during each bayesian blocks defined by Mireia.
There is clear spectral variability over the different bayesian blocks:

BB SED configuration:

  • E_reco=[31.6227766,31622.7766] GeV
  • E_true=[10,50000] GeV
  • n_bins_reco=18
  • n_bins_true=50
  • E_min_fit=31.6227766 GeV
  • E_max_fit=31622.7766 GeV
  • Aeff cut > 5%
  • 1 OFF region
These SEDs are observed, using a log-parabola spectral model.
BB1, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB2, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
BB3, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

Lightcurve[edit]

IC310 lightcurve above 100 GeV
This lightcurve was obtained using as reference the average spectrum over all the studied period, presented above in this wiki
Nightwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs
Runwise lightcurve above 100 GeV, energy dependent cuts (0.7 gammaness, 0.7 theta), interpolated IRFs

Tora[edit]

  • Producer: Tora Therese Høiland Arnesen <tora.arnesen@iac.es>

Tora: NGC 1275[edit]

Tora NGC1275: Data acquisition[edit]

Dates up to: 2024-06-01 (last obs 2024-02-10)

Selection based on Abelardo's data_quality notebook

  • File: /fefs/aswg/workspace/tora.arnesen/perseus/NGC1275/NGC1275_data_quality.ipynb
  • astropy: 5.3.4

Parameters NGC1275 selection Large.png

RelLightYield NGC1275 selection.png

Nsubruns NGC1275 selection.png

CosmicsTriggerRate NGC1275 selection.png

NGC1275SelectionRuns.png

Tora NGC1275: Data selection[edit]

Total number of runs: 153 out of 238

Total number of nights: 37 out of 50

Total observation time (hours): 41.02 out of 63.88 

20221220 : [11522 11528 11529 11530 11531 11532]
20221221 : [11546 11547 11548 11549 11550]
20221222 : [11566 11568 11569 11570 11571 11572 11573 11574]
20221223 : [11595 11596 11597 11598 11599 11600]
20230113 : [11618 11619]
20230114 : [11641 11642 11643 11644]
20230115 : [11660 11661 11662 11663]
20230817 : [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
20230911 : [14258]
20230916 : [14456]
20230918 : [14495]
20230919 : [14590]
20230920 : [14623 14624 14625]
20230921 : [14668]
20230922 : [14700 14701 14702 14703 14704]
20230923 : [14712 14713 14715 14718 14719 14720 14721]
20230924 : [14730 14731 14732 14733 14736 14737 14738]
20231006 : [14930 14931]
20231007 : [14954 14955 14956 14957 14958]
20231010 : [15040]
20231012 : [15117 15118 15119]
20231108 : [15321 15322 15323 15324 15325 15326 15327]
20231110 : [15401 15402 15403 15404 15405 15406 15407]
20231111 : [15436 15437 15438]
20231112 : [15470 15471 15472 15473]
20231113 : [15517]
20231116 : [15627 15628 15629 15630 15631]
20231117 : [15660 15661 15662 15663 15664]
20231119 : [15739 15740]
20231201 : [15786 15787 15788]
20231208 : [15889 15890 15891 15892 15893]
20240129 : [16423 16424 16425 16426]
20240130 : [16438 16439 16440 16441 16442]
20240202 : [16492 16493 16495 16496 16497 16498 16500 16503]
20240207 : [16618]
20240208 : [16650 16651 16652 16653 16654 16655]
20240210 : [16664 16665 16666 16667]

Tora NGC1275: Analysis details[edit]

  • Directory: /fefs/aswg/workspace/tora.arnesen/perseus/NGC1275
  • Pipeline:
    • ./launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822/lstchain_config_2024-01-31.json
    • ./launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_4822 ./irf_dl3_tool_config.json gh_eff_0.7_th_cont_0.7_20240417_TT
    • ./launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_4822/gh_eff_0.7_th_cont_0.7_20240417_TT NGC\ 1275
    • lstchain_create_dl3_index_files -d ./DL3 -p dl3\*.fits.gz
  • lstchain: v0.10.9

Source:

  • Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
  • Closest (chosen) dec track: 48.22 deg.
  • Cleaning: default, tail-cuts 8-4.
  • Analysis type: source-independent
  • IRFs: point-like, energy-dependent gammaness/theta2 cuts (0.7/0.7)

Starting points:

  • MC: DL2
  • Models
  • Data: DL1b (merged files with no images)


Tora NGC1275: Theta2 plots[edit]

  • Total significance (per energy bin)
All data, energy binned
  • Per bayesian block bin (from LC analysis, see below)
Theta2 plots for different bayesian blocks
BB0 BB1
TT all NGC1275 theta2 excess plots BB0.png TT all NGC1275 theta2 excess plots BB1.png
BB2 BB3
TT all NGC1275 theta2 excess plots BB2.png TT all NGC1275 theta2 excess plots BB3.png
BB4 BB5
TT all NGC1275 theta2 excess plots BB4.png TT all NGC1275 theta2 excess plots BB5.png
BB6 BB7
TT all NGC1275 theta2 excess plots BB6.png TT all NGC1275 theta2 excess plots BB7.png

Tora NGC1275: Lightcurves[edit]

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region


Lighcurve LP runwise (top) and nightwise (bottom), compared to MAGIC average flux for Crab and NGC1275
Lighcurve LP with runwise, night-wise, and bayesian block-wise flux points
Lighcurves runwise (top) and nightwise (bottom), for Power-Law (PL), Log-Parabola (LP), and Exponential cutoff Power-Law (EPWL). Compared to MAGIC average flux for Crab and NGC1275

Tora NGC1275: SEDs[edit]

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region


All data, LP model, compared with Magic 2018 LP model

Log-parabola spectral model (no EBL correction): 

-F0 = (5.3 土 0.5)・10e-11 TeV-1 cm-2 s-1
-alpha = 2.8 土 0.1
-beta = 0.13 土 0.07
-E0 = 0.3 TeV
SED per bayesian block
SED all data, for Power-Law (PL), Log-Parabola (LP), and Exponential cutoff Power-Law (EPWL). Compared with Magic 2018 LP model and EPWL model.

Poer-Law spectral model (no EBL correction): 

-F0 = (4.6 土 0.3)・10e-11 TeV-1 cm-2 s-1
-index = 2.8 土 0.1
-E0 = 0.3 TeV

Exponential cutoff Power-Law spectral model (no EBL correction): 

-F0 = (6 土 1)・10e-11 TeV-1 cm-2 s-1
-alpha = 1 (frozen)
-index = 2.6 土 0.2
-lambda = 0.5 土 0.4  TeV-1
-E0 = 0.3 TeV

Tora NGC1275: Skymaps[edit]

All data stacked Skymaps centered at the position of NGC 1275 Skymaps between 0.2 TeV and 7 TeV

Significance, off significance, and excess maps
Background residuals map

NGC 1275 Cross-check SED[edit]

Dates up to: 2024-02-02

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region
All data, LP model, compared with Magic 2018 LP models

Tora: IC 310[edit]

Tora IC310: Data acquisition[edit]

Dates up to: 2024-05-28 (last obs 2024-02-10)

Selection based on Abelardo's data_quality notebook

  • File: /fefs/aswg/workspace/tora.arnesen/perseus/IC310/IC310_data_quality.ipynb
  • astropy: 5.3.4

Parameters IC310 selection Large.png

RelLightYield IC310 selection.png

Nsubruns IC310 selection.png

CosmicsTriggerRate IC310 selection.png

IC310SelectionRuns.png

Tora IC310: Data selection[edit]

Total number of runs: 133 out of 217

Total number of nights: 37 out of 49

Total observation time (hours): 36.95 out of 59.57 

20221220 : [11529 11530]
20221221 : [11548 11549]
20221222 : [11568 11569 11573 11574]
20221223 : [11598 11599]
20230113 : [11619]
20230114 : [11641 11644]
20230115 : [11660 11663]
20230817 : [13987 13988 13989 13990 13991 13992 13993 13994 13995 13996]
20230911 : [14258]
20230916 : [14456]
20230918 : [14495]
20230919 : [14590]
20230920 : [14623 14624 14625]
20230921 : [14668]
20230922 : [14700 14701 14702 14703 14704]
20230923 : [14712 14713 14715 14718 14719 14720 14721]
20230924 : [14730 14731 14732 14733 14736 14737 14738]
20231006 : [14930 14931]
20231007 : [14954 14955 14956 14957 14958]
20231010 : [15040]
20231012 : [15117 15118 15119]
20231108 : [15321 15322 15323 15324 15325 15326 15327]
20231110 : [15401 15402 15403 15404 15405 15406 15407]
20231111 : [15436 15437 15438]
20231112 : [15470 15471 15472 15473]
20231113 : [15517]
20231116 : [15627 15628 15629 15630 15631]
20231117 : [15660 15661 15662 15663 15664]
20231119 : [15739 15740]
20231201 : [15786 15787 15788]
20231208 : [15889 15890 15891 15892 15893]
20240129 : [16423 16424 16425 16426]
20240130 : [16438 16439 16440 16441 16442]
20240202 : [16492 16493 16495 16496 16497 16498 16500 16503]
20240207 : [16618]
20240208 : [16650 16651 16652 16653 16654 16655]
20240210 : [16664 16665 16666 16667]

Tora IC310: Analysis details[edit]

  • Directory: /fefs/aswg/workspace/tora.arnesen/perseus/IC310
  • Pipeline:
    • ./IC310_launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_4822/lstchain_config_2024-01-31.json
    • ./IC310_launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_4822 ./irf_dl3_tool_config.json gh_eff_0.7_th_cont_0.7_20240510_TT
    • ./IC310_launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_4822/gh_eff_0.7_th_cont_0.7_20240510_TT IC\ 310
    • lstchain_create_dl3_index_files -d ./DL3 -p dl3\*.fits.gz
  • lstchain: v0.10.9
  • gammapy: 1.2 (SED & LCs) / 1.1 (Theta2 & Skymaps)

Source:

  • Coordinates: RA=49.95066663, DEC=41.5116969 (solved with astropy's SkyCoord / CDS name resolver)
  • Closest (chosen) dec track: 48.22 deg.
  • Cleaning: default, tail-cuts 8-4.
  • Analysis type: source-independent
  • IRFs: point-like, energy-dependent gammaness/theta2 cuts (0.7/0.7)

Starting points:

  • MC: DL2
  • Models
  • Data: DL1b (merged files with no images)


Tora IC310: Theta2 plots[edit]

  • Total significance (per energy bin)
All data, energy binned
  • Per bayesian block bin (from LC analysis, see below)
Theta2 plots for different bayesian blocks
BB0 BB1
TT all IC310 theta2 excess plots BB0.png TT all IC310 theta2 excess plots BB1.png
BB2
TT all IC310 theta2 excess plots BB2.png

Tora IC310: Lightcurves[edit]

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region
Lighcurve LP runwise (top) and nightwise (bottom), compared to MAGIC average flux for Crab and NGC1275
Lighcurve LP with runwise, night-wise, and bayesian block-wise flux points

Tora IC310: SEDs[edit]

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region
All data, LP model, compared with Magic 2018 LP model

Log-parabola spectral model (no EBL correction): 

-F0 = (2.9 土 0.6)・10e-11 TeV-1 cm-2 s-1
-alpha = 2.6 土 0.3
-beta = 0.5 土 0.3
-E0 = 0.3 TeV
SED per bayesian block

Tora IC310: Skymaps[edit]

All data stacked Skymaps centered at the position of IC 310 Skymaps between 0.2 TeV and 7 TeV

Significance, off significance, and excess maps
Background residuals map

IC 310 Cross-check SED[edit]

Last date: 2024-02-02

SED & LC configuration:

  • e_reco_min = 0.0316227766 TeV
  • e_reco_max = 31.6227766 TeV
  • e_true_min = 0.01 TeV
  • e_true_max = 50 TeV
  • e_reco_bin = 18
  • e_true_bin = 50
  • per_decade = False
  • e_fit_min = 0.0316227766 TeV
  • e_fit_max = 31.6227766 TeV
  • e_lc_min = 0.09999999999467535 TeV
  • e_lc_max = 49.99999999999999 TeV
  • Aeff cut > 5%
  • safe_min_energy = 0.04 TeV
  • safe_max_energy = 7 TeV
  • 1 OFF region
All data up to 2024-02-02, LP models, compared with Tora's NGC1275 LP model

IC 310 dec line comparison[edit]

Comparison analysis using dec 34.76 and dec 48.22

Light curves for LP models using dec 34.76 and dec 48.22, compared with Tora's NGC1275 LP model
SED for LP models using dec 34.76 and dec 48.22, compared with Tora's NGC1275 LP model
Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=NGC_1275&oldid=10509"