Difference between revisions of "AGN paper LST1"

From my_wiki
Jump to: navigation, search
(→‎1ES 0647+250)
(→‎PG 1553+113)
Line 558: Line 558:
  
 
== PG 1553+113 ==
 
== 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 [https://ui.adsabs.harvard.edu/abs/2022MNRAS.509.4330D/abstract Dorigo Jones et al]
  
 
=== Data available ===
 
=== Data available ===
−
In total there is 35 runs analysed
 
  
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
 
+
!colspan="3"| PG 1553+113 <br>
 +
|-
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
 
'''Date'''
 
'''Date'''
−
| style="background-color: lightblue;text-align:left"|
+
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
−
'''Runs'''
+
'''Runs, Zd < 30, 6.48 hrs'''
 +
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
 +
'''Runs, Zd > 30, 2.76 hrs'''
 
|-
 
|-
−
! scope="row"| 20210322
+
! scope="row"| 20210321
 +
|
 
| 4207
 
| 4207
 +
|-
 +
! scope="row"| 20210407
 +
| 4333-4336
 +
|
 
|-
 
|-
 
! scope="row"| 20210408
 
! scope="row"| 20210408
−
| 4333-4334-4335-4336
 
−
|-
 
−
! scope="row"| 20210409
 
 
| 4343-4344
 
| 4343-4344
 +
|
 
|-
 
|-
−
! scope="row"| 20210411
+
! scope="row"| 20210410
−
| 4363-4364-4365-4366-4367-4368-4369-4370
+
| 4366-4370
 +
| 4363-4365
 
|-
 
|-
−
! scope="row"| 20210421
+
! scope="row"| 20210420
−
| 4514-4515-4516
+
| 4514-4516
 +
|
 
|-
 
|-
 
! scope="row"| 20210508
 
! scope="row"| 20210508
−
| 4652-4653-4654
+
|
 +
| 4652-4654
 
|-
 
|-
 
! scope="row"| 20210604
 
! scope="row"| 20210604
−
| 4902-4903-4904-4905
+
| 4902-4905
 +
|
 +
|-
 +
! scope="row"| 20210605
 +
|
 +
| 4930-4932
 
|-
 
|-
 
! scope="row"| 20210606
 
! scope="row"| 20210606
−
| 4930-4931-4932-4940-4941-4942-4943
+
| 4940-4943
 +
|
 
|-
 
|-
 
! scope="row"| 20210609
 
! scope="row"| 20210609
−
| 4987-4988-4989
+
| 4987-4989
 +
|
 
|-
 
|-
 
|}
 
|}
 +
 
2 runs might be removed at the end
 
2 runs might be removed at the end
  
Line 604: Line 626:
 
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.
 
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  
 
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 ===
 
=== Main Results ===
Line 666: Line 707:
 
[[File:PG1553_Fev2022_significance_mapzdsup30.png|150px]]
 
[[File:PG1553_Fev2022_significance_mapzdsup30.png|150px]]
 
|}
 
|}
 +
 +
==== Analysis by Chaitanya Priyadarshi using v0.7.x data and MC  ====
 +
Preliminary LC and SED, for the full dataset-
 +
 +
[[File:LC_Run_4207_4989_lp_model_e_ref_86GeV_fit_30GeV_100TeV.png]]
 +
[[File:SED_Run_4207_4989_lp_model_e_ref_86GeV_fit_range_w_5_bin_p_decade_10_40.png]]

Revision as of 14:19, 18 February 2022

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

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)

Gammaness cuts distribution 80% efficiency different MC pointing.png Theta cuts distribution 80% efficiency different MC pointing.png

Effective area 80% efficiency different MC pointing.png Energy dispersion 80% efficiency different MC pointing.png

Mrk 421

Mrk421 coordinates:

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

Data available

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-

LC Run 3216 4962 lp model e ref 342GeV fit 300GeV 100TeV.png

SED Run 3216 4962 lp model e ref 342GeV fit range w 5 bin p decade 10 40.png

Mrk 501

Mrk501 coordinates:

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

Data available

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-

LC Run 2177 6078 lp model e ref 487GeV fit 300GeV 100TeV.png

SED Run 2177 6078 lp model e ref 487GeV fit range w 5 bin p decade 10 40.png

1ES 1959+650

1ES1959+650 coordinates:

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

Data available

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-

LC Run 2203 5542 lp model e ref 371GeV fit 30GeV 100TeV.png SED Run 2203 5542 lp model e ref 371GeV fit range w 5 bin p decade 10 40.png

1ES 0647+250

1ES0647+250 coordinates:

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

Data available

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-

LC Run 3281 3364 lp model e ref 56GeV fit 30GeV 100TeV.png SED Run 3281 3364 lp model e ref 56GeV fit range w 5 bin p decade 10 40.png

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

PG1553 Fev2022 countmap.png


Theta2 plot


PG1553 Theta2.png

Excess and Significance map


PG1553 Fev2022 significance map.png


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

PG1553 Fev2022 significance map zdinf30.png

PG1553 Fev2022 significance mapzdsup30.png

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

Preliminary LC and SED, for the full dataset-

LC Run 4207 4989 lp model e ref 86GeV fit 30GeV 100TeV.png SED Run 4207 4989 lp model e ref 86GeV fit range w 5 bin p decade 10 40.png

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=AGN_paper_LST1&oldid=4140"