AGN paper LST1

People involved

• Chaitanya Priyadarshi
• Daniel Morcuende
• David Sanchez
• Estelle Pons
• Andres Baquero
• Ryuji Takeishi
• Joshua Baxter
• Gaetano Di Marco
• Mireia Nievas
• Monica Vazquez Acosta

Slack channel

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

Meetings

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

List of AGNs

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

Using MC with 2 different sky pointing, where the RF model was trained using the corresponding MC DL1 files

• For zenith < 30 deg, MC (point-gamma, proton and electron) and RF models used -

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/gamma/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/ /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/

• For zenith > 30 deg, MC (point-gamma, proton and electron) and RF models used -

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/ /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/

For DL3 cuts, energy-dependent cuts with an efficiency of 80% were used for gammaness and theta(only for point-like IRFs) Using the DL3 Tools of lstchain v0.9.2 and with gammapy developer build after the stable release of v0.19 (to extract signal from energy-dependent theta regions)

Mrk 421

Mrk421 coordinates:

• RA = 11h 04m 27.314s
• DEC = 38deg 12m 31.80s
• z = 0.0308

Wobble (0.4 deg), Zd < 30

Mrk421 Wobble (0.4 deg)
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 (0.8 or 0.5 deg)

Mrk421 Wobble (0.8 or 0.5 deg)
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

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

A total of 27 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, ~340 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-

Mrk 501

Mrk501 coordinates:

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

Wobble observations

Mrk501 Wobble observations
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, Zd > 30

1ES1959+650 Wobble observations, offset~0.4, Zd > 30, 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, Zd > 30

1ES1959+650 Wobble observations, offset 0.2-0.6 deg, Zd > 30, 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