AGN paper LST1

People involved

• Chaitanya Priyadarshi
• Daniel Morcuende
• David Sanchez
• Estelle Pons
• Andres Baquero
• Ryuji Takeishi
• Joshua Baxter
• Gaetano Di Marco
• Mireia Nievas
• Lea Heckmann
• Noah Biederbeck
• Monica Vazquez Acosta
• José Luis Contreras

Slack channel

https://cta-north.slack.com/archives/C031QEQ3YDP

Meetings

2022-02-23: Indico[1], Zoom[2]

List of AGNs

Mrk 421

• Name of the source: Mrk 421
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.0308
  • HBL
  • RA: 11 04 27.314 (hh mm ss), Dec: +38 12 31.8 (dd mm ss)
  • RA, Dec in deg (ICRS): 166.113808

Wobble observations, offset = 0.4 deg

Mrk421
Date	Runs (Zd < 30, 17.23 hrs)	Runs (Zd > 30, 11 hrs)
20201212	3222-3227	3216-3221
20201213	3238, 3239
20201214	3247-3250
20210214	3661-3665	3666-3668
20210215	3687, 3688
20210216	3715, 3716	3717-3721
20210217	3732-3736
20210219	3800
20210313	4016	4017-4019
20210316	4097-4100
20210317	4131-4134	4135-4138
20210318		4479
20210320	4184, 4185	4186-4188
20210321		4208
20210412		4386-4388
20210416	4441
20210417	4457-4459
20210418		4479
20210430	4568-4572
20210505	4612-4614
20210508	4612-4614
20210509	4667-4670	4671
20210531		4852, 4853
20210601		48668, 4869
20210606		4939
20210608		4962

Wobble observations, offset = 0.8 or 0.5 deg

Mrk421
Date	Runs, Zd < 30, 1.76 hrs	Runs, Zd > 30, 2.18 hrs
20210212	3617-3619	3620-3624
20210213	3641	3642-3646
20210416	4442, 4443

Data analysis

Mrk 421 Preliminary Analysis with lstchain v0.7.x Priyadarshi

Main Results

Mrk 421 Preliminary Analysis with lstchain v0.7.x Priyadarshi#Analysis Results

Mrk 501

Mrk501 coordinates:

• RA = 16h 53m 52.21s
• DEC = 39deg 45m 37.6s
• z = 0.033640

Wobble observations

Mrk501
Date	Runs, Zd < 30, 19.8 hrs	Runs (Zd > 30, 5.24 hrs)
20200710	2177-2180
20200909		2606-2607
20210317	4139, 4140
20210319		4167-4170
20210320	4191-4193
20210321	4212-4215
20210416	4452, 4455
20210417	4473-4476
20210418	4480	4483
20210419	4496-4498
20210503	4586-4590
20210504	4600-4602	4597-4599
20210506	4625-4628
20210507	4643-4645
20210508	4656, 4657
20210509	4673-4677	4672
20210605	4921-4933	4934
20210606	4944-4947
20210608	4977
20210612	5011
20210613	5039, 5040	5038
20210702		5095-5096
20210708	5202-5205
20210710	5227, 5228
20210801	5438	5439
20210805	5489, 5490
20210808	5550	5551
20210906	6076	6077, 6078

Data analysis

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

A total of 24.5 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

• OFF regions selected - 1 Wobble
• Safe Energy mask for systematics - No
• Energy axes -
  • True: 10 GeV to 100 TeV with 10 bins per decade
  • Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

• Model used - Log Parabola
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~490 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Main Results

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

Preliminary LC and SED, for the full dataset-

1ES 1959+650

1ES1959+650 coordinates:

• RA = 19h 59m 59.8s = 299.999 deg
• DEC = +65deg 08m 55s = 65.149 deg
• z = 0.048

Wobble observations, offset ~ 0.4 deg, Zd > 30

1ES1959+650 8.62 hrs
Date	Runs
20200713	2202
20200809	2370, 2372
20200811	2411-2413
20200812	2420-2422
20210319	4173, 4174, 4176
20210530	4840, 4841, 4844
20210807	5520-5535, 5540, 5542

Wobble observations, offset = 0.2 - 0.6 deg, Zd > 30

1ES1959+650 3.17 hrs
Date	Runs
20200711	2193-2196
20200713	2206, 2207
20200726	2348, 2350-2352
20200809	2371, 2374
20200811	2414

Data analysis

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

A total of 8.69 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

• OFF regions selected - 1 Wobble
• Safe Energy mask for systematics - No
• Energy axes -
  • True: 10 GeV to 100 TeV with 10 bins per decade
  • Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

• Model used - Log Parabola
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~370 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Main Results

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

Preliminary LC and SED, for the full dataset-

1ES 0647+250

1ES0647+250 coordinates:

• RA = 06h 50m 46.5s
• DEC = +25deg 03m 00s
• z = 0.41 [3]

Wobble Observations

1ES0647+250
Date	Runs, Zd < 30, 6.8 hrs	Runs, Zd > 30, 1.64 hrs
20201215	3281-3286	3287, 3288
20201218	3331-3335
20201219	3341-3344	3346-3348
20201220	3357-3364

Data analysis

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

A total of 8.16 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

• OFF regions selected - 1 Wobble
• Safe Energy mask for systematics - No
• Energy axes -
  • True: 10 GeV to 100 TeV with 10 bins per decade
  • Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

• Model used - Log Parabola
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~56 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Main Results

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

Preliminary LC and SED, for the full dataset-

PG 1553+113

PG1553+113 coordinates:

• RA = 15h 55m 44.7s = 238.936 deg
• DEC = +11deg 11m 41s = 11.195 deg
• z = 0.43-0.58
• z = 0.433 Dorigo Jones et al

Data available

PG 1553+113
Date	Runs, Zd < 30, 6.48 hrs	Runs, Zd > 30, 2.76 hrs
20210321		4207
20210407	4333-4336
20210408	4343-4344
20210410	4366-4370	4363-4365
20210420	4514-4516
20210508		4652-4654
20210604	4902-4905
20210605		4930-4932
20210606	4940-4943
20210609	4987-4989

2 runs might be removed at the end

Data analysis

We used the V0.7.3 and V0.7.1 data and used the corresponding IRFs (there is a wide spread of zenith angle) for each run. The background model was computed using the tool of Mathieu de Bony presented at the analysis call

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

A total of 8.78 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

• OFF regions selected - 1 Wobble
• Safe Energy mask for systematics - No
• Energy axes -
  • True: 10 GeV to 100 TeV with 10 bins per decade
  • Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

• Model used - Log Parabola
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~86 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Main Results

Livetime : 8.07 h

Significance : 13.5

Excess : 1624.9

S/B : 0.145

Rate : 3.36 per min

Count map of PG 1553

Theta2 plot

Excess and Significance map

Zd<30	Zd>30
Livetime : 6.34 h Significance : 13.8 Excess : 1462.83 S/B : 0.18	Livetime : 2.33 h Significance : 2.68 Excess : 181.6 S/B : 1.30

Analysis by Chaitanya Priyadarshi using v0.7.x data and MC

Preliminary LC and SED, for the full dataset-

• With Log Parabola model

• With Power Law model

M87

M87 coordinates:

• RA = 12h 30m 49s = 187.7059 deg
• DEC = +12deg 23m 28s = 12.3911 deg
• z = 0.0042

Wobble observations, offset = 0.4 deg

M87
Date	Runs (Zd < 30, 9.26 hrs)	Runs (Zd > 30, 1.94 hrs)
20201217	3738,3739, 3741-3745	3746-3748
20210320	4181, 4182
20210321	4198-4206	4197
20210413	4396-4398
20210417	4460-4465
20210419	4490-4493	4494, 4495
20210420		4506-4509
20210504	4594, 4595
20210507	4636, 4637
20210609	4985, 4986