Difference between revisions of "Mrk 501 Xcheck lstchain v0.9.x Heckmann"

From my_wiki
Jump to: navigation, search
(→‎DL2 data)
 
(16 intermediate revisions by the same user not shown)
Line 13: Line 13:
  
  
 +
 +
== DL1 data and models ==
 +
 +
The DL1 files of the real data are stored here -
 +
/remote/ceph/user/h/heckmann/Mrk501_LST/data/DL1/v0.9.3/tailcut84/
 +
 +
Models used:
 +
/fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_3476/
  
 
== DL2 data ==
 
== DL2 data ==
  
 
The DL2 files of the real data are stored here -
 
The DL2 files of the real data are stored here -
−
/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/Mrk501/v09x/tailcut84/DL2/real/src_indep/20220531_dec3476/{Date}
+
/fefs/aswg/workspace/lea.heckmann/Data/Mrk501/v09x/tailcut84/DL2/real/src_indep/20230127_dec3476/
  
 
Link to MC DL2 files used:
 
Link to MC DL2 files used:
−
/fefs/aswg/data/mc/DL2/AllSky/20220531_dec3476_std/TestingDataset/dec_3476/
+
/fefs/aswg/data/mc/DL2/AllSky/20230127_v0.9.12_base_prod/TestingDataset/dec_3476/
 +
 
 +
IRFs:
 +
/fefs/aswg/workspace/lst1-agn-paper/LST1/Mrk501/Analysis2/source-independent/irfs/irf_20230127_v0.9.12_base_prod_az_tel_node_theta*
  
 
==DL3 data selection==
 
==DL3 data selection==
Line 40: Line 51:
 
==High-level analysis==
 
==High-level analysis==
  
−
Runs included:  
+
Data selection:
 +
 
 +
max_zen = 52
 +
min_livetime = 60s
 +
wobble_offset_diff = 0.1 (from 0.4)
 +
max_ped_std = 1.8 phe
 +
pixel rate of >10 and >30 phe pulses in cosmics
 +
 
 +
 
 +
153 Runs included:  
  
−
[2606 2607 2608 2610 4167 4168 4169 4452 4453 4454 4455 4473 4474 4475
+
[2177 2178 2179 2606 2607 2608 2610 4139 4140 4167 4168 4169 4170 4191
−
  4476 4480 4483 4484 4496 4497 4498 4586 4587 4588 4597 4598 4599 4600
+
4192 4193 4194 4195 4212 4213 4214 4215 4452 4455 4473 4474 4475 4476
−
  4601 4602 4625 4626 4627 4628 4642 4643 4644 4645 4656 4657 4672 4673
+
  4480 4483 4484 4496 4497 4498 4586 4587 4588 4597 4598 4599 4600 4601
−
  4674 4675 4676 4677 4922 4933 4934 4944 4945 4946 4947 5011 5094 5095
+
  4602 4625 4626 4627 4628 4642 4643 4644 4645 4656 4657 4672 4673 4674
−
5134 5202 5203 5204 5205 5206 5438 5439 5489 5490 5550 5551 6076 6077
+
  4675 4676 4677 4922 4933 4934 4944 4945 4946 4947 5011 5037 5038 5039
−
6078 6863 6864 6865 6880 6881 6882 6928 6930 6931 6959 6978 6979 6980
+
5040 5094 5095 5134 5202 5203 5204 5205 5206 5227 5228 5438 5439 5489
−
6981 7155 7156 7158 7353 7354 7356 7357 7358 7359 7360 7361 7362 7363
+
5490 5550 5551 6076 6077 6078 6862 6863 6864 6865 6879 6880 6881 6882
−
7364 7381 7382 7383 7384 7385 7401 7402 7403 7482 7483 7485 7564 7565
+
6928 6929 6930 6931 6958 6959 6978 6979 6980 6981 7155 7156 7158 7351
−
7635 7636 7637 7718 7719 7720 7721 7722 7723 7724 7725 7775 7776 7777
+
7353 7354 7355 7356 7357 7358 7359 7360 7361 7362 7363 7364 7381 7382
−
7778 7860 7861 7862 7863 7875 7876 7877 8274 8275 8276 8325 8326 8381
+
7383 7384 7385 7401 7402 7403 7482 7483 7564 7565 7635 7636 7637 7718
−
8382 8383 8416 8417 8418 8419]
+
7719 7720 7721 7722 7723 7724 7775 7777 7778 7875 7876 8276 8416]  
  
 
* lstchain Tools to generate source-independent IRF and DL3
 
* lstchain Tools to generate source-independent IRF and DL3
Line 67: Line 87:
 
For spectral analysis,  
 
For spectral analysis,  
 
* Model used -  
 
* Model used -  
−
** Log Parabola + EBL (Dominguez) (6.5 sigmas preference over Power Law + EBL)
+
** Power law or Log Parabola + EBL (Dominguez) - depending on if log parabola is preferred at >3sigma
−
* Reference energy used by calculating the decorrelation energy with a Power Law fit, 561.5 GeV (Full dataset)
+
* Reference energy used by calculating the decorrelation energy with a Power Law fit, 574 GeV (Full dataset)
 
* Spectral Fit energy range - 10 GeV to 100 TeV with 5 bins per decade
 
* Spectral Fit energy range - 10 GeV to 100 TeV with 5 bins per decade
 
* Spectral fit with gammapy default optimization - minuit
 
* Spectral fit with gammapy default optimization - minuit
Line 82: Line 102:
 
* Full Dataset
 
* Full Dataset
  
−
[[File:COMP_LC.png|500 px]]
+
[[File:COMP_LC.png|1800 px]]
  
 
[[File:COMP_SED.png|500 px]]
 
[[File:COMP_SED.png|500 px]]

Latest revision as of 12:38, 15 June 2023

Xcheck following the analysis in: Mrk 501 Preliminary Analysis with lstchain v0.9.x Priyadarshi w IRF interpolation

General information[edit]

  • Name of the source: Mrk 501
  • Brief description of the source:
    • Object type : AGN
    • Redshift (z) : 0.034
    • RA: 16 53 52.216 (hh mm ss), Dec: +39 45 36.613 (dd mm ss)
    • RA, Dec in deg (ICRS): 253.46756699, 39.76017034
  • Analysis by Heckmann Lea (MPP - heckmann@mpp.mpg.de)


DL1 data and models[edit]

The DL1 files of the real data are stored here - /remote/ceph/user/h/heckmann/Mrk501_LST/data/DL1/v0.9.3/tailcut84/

Models used: /fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_3476/

DL2 data[edit]

The DL2 files of the real data are stored here - /fefs/aswg/workspace/lea.heckmann/Data/Mrk501/v09x/tailcut84/DL2/real/src_indep/20230127_dec3476/

Link to MC DL2 files used: /fefs/aswg/data/mc/DL2/AllSky/20230127_v0.9.12_base_prod/TestingDataset/dec_3476/

IRFs: /fefs/aswg/workspace/lst1-agn-paper/LST1/Mrk501/Analysis2/source-independent/irfs/irf_20230127_v0.9.12_base_prod_az_tel_node_theta*

DL3 data selection[edit]

For energy-dependent cuts[edit]

Information about your DL3 data selection.

  • intensity > 50
  • r: [0, 1]
  • wl: [0.01, 1]
  • leakage_intensity_width_2: [0, 1]
  • source-independent
    • gh_efficiency: 0.8
    • theta_containment: 0.8


High-level analysis[edit]

Data selection: 

max_zen = 52 
min_livetime = 60s 
wobble_offset_diff = 0.1 (from 0.4)
max_ped_std = 1.8 phe
pixel rate of >10 and >30 phe pulses in cosmics


153 Runs included: 

[2177 2178 2179 2606 2607 2608 2610 4139 4140 4167 4168 4169 4170 4191
4192 4193 4194 4195 4212 4213 4214 4215 4452 4455 4473 4474 4475 4476
4480 4483 4484 4496 4497 4498 4586 4587 4588 4597 4598 4599 4600 4601
4602 4625 4626 4627 4628 4642 4643 4644 4645 4656 4657 4672 4673 4674
4675 4676 4677 4922 4933 4934 4944 4945 4946 4947 5011 5037 5038 5039
5040 5094 5095 5134 5202 5203 5204 5205 5206 5227 5228 5438 5439 5489
5490 5550 5551 6076 6077 6078 6862 6863 6864 6865 6879 6880 6881 6882
6928 6929 6930 6931 6958 6959 6978 6979 6980 6981 7155 7156 7158 7351
7353 7354 7355 7356 7357 7358 7359 7360 7361 7362 7363 7364 7381 7382
7383 7384 7385 7401 7402 7403 7482 7483 7564 7565 7635 7636 7637 7718
7719 7720 7721 7722 7723 7724 7775 7777 7778 7875 7876 8276 8416]
  • lstchain Tools to generate source-independent IRF and DL3
  • Science Tool: gammapy 1.0
  • point-like IRF, 1D analysis

For DL3 data reduction to DL4 dataset,

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

For spectral analysis,

  • Model used -
    • Power law or Log Parabola + EBL (Dominguez) - depending on if log parabola is preferred at >3sigma
  • Reference energy used by calculating the decorrelation energy with a Power Law fit, 574 GeV (Full dataset)
  • Spectral Fit energy range - 10 GeV to 100 TeV with 5 bins per decade
  • Spectral fit with gammapy default optimization - minuit


For LC,

  • E_min = 100 Gev

Analysis Results[edit]

Preliminary LC and SEDs in comparison

  • Full Dataset

COMP LC.png

COMP SED.png

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Mrk_501_Xcheck_lstchain_v0.9.x_Heckmann&oldid=7799"