Difference between revisions of "AGN paper LST1"

From my_wiki
Jump to: navigation, search
(→‎1ES 1959+650)
(→‎1ES 0647+250)
Line 490: Line 490:
  
 
== 1ES 0647+250 ==  
 
== 1ES 0647+250 ==  
 +
'''1ES0647+250 coordinates''':
 +
* RA = 06h 50m 46.5s
 +
* DEC = +25deg 03m 00s
 +
* z = 0.41 [https://pos.sissa.it/395/792]
 +
* z (recently fixed) = 0.433 [https://ui.adsabs.harvard.edu/link_gateway/2022MNRAS.509.4330D/doi:10.1093/mnras/stab3331]
  
 
=== Data available ===
 
=== Data available ===
 +
=== Wobble Observations ===
 +
 +
{| class="wikitable"
 +
|-
 +
!colspan="3"| 1ES0647+250 <br>
 +
|-
 +
| style="background-color: lightblue; text-align:center"|
 +
'''Date'''
 +
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
 +
'''Runs, Zd < 30, 6.8 hrs'''
 +
| style="background-color: lightblue;text-align:center"; style=max-width:8em|
 +
'''Runs, Zd > 30, 1.64 hrs'''
 +
|-
 +
! scope="row"| 20201215
 +
| 3281-3286
 +
| 3287, 3288
 +
|-
 +
! scope="row"| 20201218
 +
| 3331-3335
 +
|
 +
|-
 +
! scope="row"| 20201219
 +
| 3341-3344
 +
| 3346-3348
 +
|-
 +
! scope="row"| 20201220
 +
| 3357-3364
 +
|
 +
|-
 +
|}
 +
  
 
=== Data analysis ===  
 
=== 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 ===
 
=== Main Results ===
 +
 +
==== Analysis by Chaitanya Priyadarshi using v0.7.x data and MC  ====
 +
Preliminary LC and SED, for the full dataset-
 +
 +
[[File:LC_Run_3281_3364_lp_model_e_ref_56GeV_fit_30GeV_100TeV.png]]
 +
[[File:SED_Run_3281_3364_lp_model_e_ref_56GeV_fit_range_w_5_bin_p_decade_10_40.png]]
  
 
== PG 1553+113 ==
 
== PG 1553+113 ==

Revision as of 14:03, 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]
  • z (recently fixed) = 0.433 [2]

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

Data available

In total there is 35 runs analysed

Date

Runs

20210322 4207
20210408 4333-4334-4335-4336
20210409 4343-4344
20210411 4363-4364-4365-4366-4367-4368-4369-4370
20210421 4514-4515-4516
20210508 4652-4653-4654
20210604 4902-4903-4904-4905
20210606 4930-4931-4932-4940-4941-4942-4943
20210609 4987-4988-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

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

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