Difference between revisions of "AGN paper LST1"

From my_wiki
Jump to: navigation, search
(→‎People involved)
(→‎List of AGNs)
Line 52: Line 52:
 
=== Data available ===  
 
=== Data available ===  
  
−
=== Wobble (0.4 deg), Zd < 30 ===
+
=== Wobble observations, offset = 0.4 deg ===
  
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
!colspan="3"| Mrk421 <br> Wobble (0.4 deg)
+
!colspan="3"| Mrk421 <br>
 
|-  
 
|-  
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
Line 172: Line 172:
  
  
−
=== Wobble (0.8 or 0.5 deg) ===
+
=== Wobble observations, offset = 0.8 or 0.5 deg ===
  
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
!colspan="3"| Mrk421 <br> Wobble (0.8 or 0.5 deg)
+
!colspan="3"| Mrk421 <br>
 
|-  
 
|-  
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
Line 242: Line 242:
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
!colspan="3"| Mrk501 <br> Wobble observations
+
!colspan="3"| Mrk501 <br>
 
|-  
 
|-  
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
Line 407: Line 407:
 
=== Data available ===  
 
=== Data available ===  
  
−
=== Wobble observations, offset~0.4,  Zd > 30 ===
+
=== Wobble observations, offset ~ 0.4 deg,  Zd > 30 ===
  
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
!colspan="2"| 1ES1959+650 <br> Wobble observations, offset~0.4, Zd > 30, 8.62 hrs
+
!colspan="2"| 1ES1959+650 <br> 8.62 hrs
 
|-  
 
|-  
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
Line 442: Line 442:
 
|}
 
|}
  
−
=== Wobble observations, offset 0.2-0.6,  Zd > 30 ===
+
=== Wobble observations, offset = 0.2 - 0.6 deg,  Zd > 30 ===
  
 
{| class="wikitable"  
 
{| class="wikitable"  
 
|-  
 
|-  
−
!colspan="2"| 1ES1959+650 <br> Wobble observations, offset 0.2-0.6 deg, Zd > 30, 3.17 hrs
+
!colspan="2"| 1ES1959+650 <br> 3.17 hrs
 
|-  
 
|-  
 
| style="background-color: lightblue; text-align:center"|
 
| style="background-color: lightblue; text-align:center"|
Line 727: Line 727:
 
[[File:LC_Run_4207_4989_pl_model_e_ref_86GeV_fit_30GeV_100TeV.png]]
 
[[File:LC_Run_4207_4989_pl_model_e_ref_86GeV_fit_30GeV_100TeV.png]]
 
[[File:SED_Run_4207_4989_pl_model_e_ref_86GeV_fit_range_w_5_bin_p_decade_10_40.png]]
 
[[File:SED_Run_4207_4989_pl_model_e_ref_86GeV_fit_range_w_5_bin_p_decade_10_40.png]]
 +
 +
 +
== M87 ==
 +
'''M87 coordinates''':
 +
* RA = 12h 30m 49s = 187.7059 deg
 +
* DEC = +12deg 23m 28s = 12.3911 deg
 +
* z = 0.0042
 +
 +
=== Data available ===
 +
 +
=== Wobble observations, offset = 0.4 deg ===
 +
 +
{| class="wikitable"
 +
|-
 +
!colspan="3"| M87 <br>
 +
|-
 +
| style="background-color: lightblue; text-align:center"|
 +
'''Date'''
 +
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
 +
'''Runs (Zd < 30, 9.26 hrs)'''
 +
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
 +
'''Runs (Zd > 30, 1.94 hrs)'''
 +
|-
 +
! scope="row"| 20201217
 +
| 3738,3739, 3741-3745
 +
| 3746-3748
 +
|-
 +
! scope="row"| 20210320
 +
| 4181, 4182
 +
|
 +
|-
 +
! scope="row"| 20210321
 +
| 4198-4206
 +
| 4197
 +
|-
 +
! scope="row"| 20210413
 +
| 4396-4398
 +
|
 +
|-
 +
! scope="row"| 20210417
 +
| 4460-4465
 +
|
 +
|-
 +
! scope="row"| 20210419
 +
| 4490-4493
 +
| 4494, 4495
 +
|-
 +
! scope="row"| 20210420
 +
|
 +
| 4506-4509
 +
|-
 +
! scope="row"| 20210504
 +
| 4594, 4595
 +
|
 +
|-
 +
! scope="row"| 20210507
 +
| 4636, 4637
 +
|
 +
|-
 +
! scope="row"| 20210609
 +
| 4985, 4986
 +
|
 +
|-
 +
|}
 +
 +
=== Data analysis ===
 +
 +
=== Main Results ===

Revision as of 09:43, 23 February 2022

Contents

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

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 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

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

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 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-

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 [3]

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-

  • With Log Parabola model

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

  • With Power Law model

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


M87

M87 coordinates:

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

Data available

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

Data analysis

Main Results

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