Difference between revisions of "M87 TTHA"

From my_wiki
Jump to: navigation, search
(→‎M87 analysis (TTHA))
Line 88: Line 88:
 
20240607 : [17676 17677 17678 17679 17680 17681]
 
20240607 : [17676 17677 17678 17679 17680 17681]
 
  <nowiki>
 
  <nowiki>
 +
  
 
==Monte Carlo information==
 
==Monte Carlo information==
Line 101: Line 102:
  
  
−
==DL1 data==
+
==DL2 Data==
−
 
 
−
Please include your settings to produce your DL1a and DL1b files.  This includes whether you use LSTOSA, specific versions of lstchain, dllab scripts, cleaning levels, and calibration information.
 
−
 
 
−
'''Example''' :
 
−
 
 
−
Reprocessed using DL1a files produced by LSTOSA (lstchain v0.7.3) and dl1ab script (v0.7.5)
 
−
 
 
−
* original DL1a files
 
−
 
 
−
/fefs/aswg/data/real/DL1/20210808/v0.7.3/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
 
−
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/XXX/zenith_20deg/south_pointing/20210416_v0.7.3_prod5_trans_80_local_taicut_8_4/
 
−
 
 
−
* lstchain v0.7.5
 
−
* real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning
 
−
 
 
−
  /fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/lstchain_standard_config.json
 
−
 
 
−
* MC: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF+NSB tuning
 
−
  "image_modifier": {
 
−
    "increase_nsb": true,
 
−
    "extra_noise_in_dim_pixels": 0.74,
 
−
    "extra_bias_in_dim_pixels": 0.38,
 
−
    "transition_charge": 8,
 
−
    "extra_noise_in_bright_pixels": 0.052,
 
−
    "increase_psf": true,
 
−
    "smeared_light_fraction": 0.125
 
−
  },
 
−
 
 
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/lstchain_dl1ab_tune_MC_to_BLLac_config.json
 
  
−
* Produced DL1b files
+
Create DL2 files by running
  
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/data/
+
  ./M87_launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931/lstchain_config_2024-01-31.json
−
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/Run0XXXX
 
  
 +
DL2 files stored in
  
−
==Random forest==
+
/fefs/aswg/workspace/tora.arnesen/M87/DL2
  
−
Please include your specific .json files used in producing your RFs and for which source analysis (source dependent or independent), and variabules used in the RF. 
 
  
−
'''Example'''
+
==IRF Files==
  
−
* lstchain-0.7.6.dev242+g2086cb9
+
Use copy of school config file
 +
* 70 % efficiency for both gammaness and theta gamma-ray selection cuts
 +
* Point-like IRF
  
−
* source-indep (diffuse gamma, src_r<3deg, disp_norm)
+
  /fefs/aswg/workspace/tora.arnesen/M87/irf_dl3_tool_config.json
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcindep/RF/lstchain_standard_config.json
 
  
 +
Create IRF files by running
  
−
* source-dep
+
  ./M87_launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_931 ./irf_dl3_tool_config.json 931_gh_eff_0.7_th_cont_0.7_20240609_TT
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcdep/RF/lstchain_src_dep_config.json
 
  
−
==DL2 data==
+
IRF files stored in
  
−
Information about your DL2 data and settings such as: specific versions of lstchain, specific .json version used.
+
/fefs/aswg/workspace/tora.arnesen/M87/IRF
  
−
'''Example''':
 
  
−
* lstchain-0.7.6.dev242+g2086cb9
+
==DL3 Data==
  
−
* source-indep
+
Create DL3 files by running
−
/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcindep/DL2/data/
 
−
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/srcindep/DL2/data/dl2_LST-1.Run0XXXX_merged.h5
 
  
−
* source-dep
+
  ./M87_launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_931/931_gh_eff_0.7_th_cont_0.7_20240609_TT M87
−
  /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcdep/DL2/data/
 
−
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/srcdep/DL2/data/dl2_LST-1.Run0XXXX_merged.h5
 
  
−
==DL3 data selection==
+
DL3 files stored in
  
−
Information about your DL3 data selection.
+
/fefs/aswg/workspace/tora.arnesen/M87/DL3
  
−
'''Example'''
 
−
*intensity > 50
 
−
*r: [0,  1 ]
 
−
*wl: [0.1, 1 ]
 
−
*leakage_intensity_width_2: [0, 0.2 ]
 
−
 
−
* source-indep
 
−
**fixed_gh_cut: 0.3
 
−
**fixed_theta_cut: 0.2
 
−
 
−
* source-dep
 
−
**fixed_gh_cut: 0.7
 
−
**fixed_alpha_cut: 10
 
  
 
==High-level analysis==
 
==High-level analysis==
  
−
Please put any information about the production of higher level analysis here.
+
* lstchain v0.10.9 (v0.10.7 for Theta2 and Skymaps)
−
 
+
* Gammapy 1.2 (1.1 for Theta2 and Skymaps)
−
'''Example'''
+
* 1D analysis (point-like IRF)
−
 
 
−
*lstchain to generate source-dep IRF, DL3
 
−
*Science Tool: gammapy 0.18.2
 
−
* point-like IRF, 1D analysis
 
  
 
==Analysis Results==
 
==Analysis Results==

Revision as of 13:27, 10 June 2024

M87 analysis (TTHA)

General information

  • Name of the source: M87
  • Brief description of the source:
    • Object type : AGN / Radiogal / FRI
    • Redshift (z) : z=0.00420 (16.52 Mpc)
    • Other relevant information: [TeVCat] [Simbad]
    • RA: 12 30 49.42338414 (hh mm ss), Dec: +12 23 28.0436859 (dd mm ss)
    • RA, Dec in deg (ICRS): 187.70593075958334, 12.391123293916666
  • Analysis by Tora T. H. Arnesen (Instituto de Astrofísica de Canarias - <tora.arnesen@iac.es>)

Data-selection information

Data selected with data quality notebook 

/fefs/aswg/workspace/tora.arnesen/M87/M87_data_quality.ipynb
  • Path to datacheck files
    • files = glob.glob("/fefs/aswg/data/real/OSA/DL1DataCheck_LongTerm/night_wise/all/DL1_datacheck_20*")
  • Dates of data-selection:
    • first_date = 20191123
    • last_date = 20240609
  • ZD range (deg):
    • min_zenith = 0 * u.deg
    • max_zenith = 60 * u.deg
  • maximum diffuse NSB
    • max_diffuse_nsb_std = 2.3
  • Pointing stability
    • max_pointing_dec_std = 0.011
  • Cherenkov Transparency parameters
    • min_mean_fit_p = -3.
    • max_LS_periodogram_maxamplitude = 1e-2
    • min_drdi_index = -2.35
    • max_drdi_index = -2.1
    • min_drdi_at_422pe = 1.25
    • min_fraction_around_mode = 0.8
    • max_intensity_at_half_peak_rate = 70

Selected Data

Total number of runs: 130 out of 234

Total number of nights: 36 out of 49

Total observation time (hours): 33.70 out of 60.15 

<nowiki>

20210413 : [4396 4397 4398] 20210417 : [4460 4461 4462 4463 4464 4465] 20210419 : [4494 4495] 20210420 : [4508 4509] 20210504 : [4594 4595] 20210507 : [4636 4637] 20220301 : [7185 7186 7187] 20220323 : [7444] 20220407 : [7712 7713 7714 7715 7716] 20220408 : [7746] 20220409 : [7773] 20220503 : [8092] 20220506 : [8209 8211 8212 8213 8214] 20230412 : [12662] 20230418 : [12744 12745 12746 12751 12752 12753 12754] 20230419 : [12775 12776 12778] 20230420 : [12796 12797 12798 12799] 20230614 : [13300 13301 13302 13303] 20240110 : [16315 16316 16317 16318 16319] 20240210 : [16692 16693 16694 16695] 20240212 : [16751 16752 16753 16754 16755] 20240304 : [16960 16961 16962 16963] 20240310 : [17034 17035 17037 17045 17047 17048] 20240311 : [17072 17075 17076 17083 17084 17085 17086] 20240313 : [17103 17104 17105] 20240315 : [17160 17161] 20240502 : [17244 17245 17246] 20240503 : [17253 17254 17255 17256 17258 17259 17261 17262 17263 17264 17265 17266] 20240504 : [17284 17289 17292] 20240507 : [17359] 20240514 : [17403 17404] 20240527 : [17449 17450 17451 17452 17453] 20240601 : [17570] 20240603 : [17612 17613 17614 17615] 20240604 : [17641 17642 17643 17644] 20240607 : [17676 17677 17678 17679 17680 17681] 

<nowiki>


Monte Carlo information

  • All-sky MC production used: 20240131_allsky_v0.10.5_all_dec_base
    • Declination line (training): dec_931

Random forest

Source declination is 9.31 deg. The two closest MC declination tracks are 9.31 deg and 18.02 deg. Chose dec 9.31 deg as it is closer to the dec of the source at 12.32 deg. 

/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931


DL2 Data

Create DL2 files by running 

./M87_launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931/lstchain_config_2024-01-31.json

DL2 files stored in 

/fefs/aswg/workspace/tora.arnesen/M87/DL2


IRF Files

Use copy of school config file

  • 70 % efficiency for both gammaness and theta gamma-ray selection cuts
  • Point-like IRF
/fefs/aswg/workspace/tora.arnesen/M87/irf_dl3_tool_config.json

Create IRF files by running 

./M87_launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_931 ./irf_dl3_tool_config.json 931_gh_eff_0.7_th_cont_0.7_20240609_TT

IRF files stored in 

/fefs/aswg/workspace/tora.arnesen/M87/IRF


DL3 Data

Create DL3 files by running 

./M87_launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_931/931_gh_eff_0.7_th_cont_0.7_20240609_TT M87

DL3 files stored in 

/fefs/aswg/workspace/tora.arnesen/M87/DL3


High-level analysis

  • lstchain v0.10.9 (v0.10.7 for Theta2 and Skymaps)
  • Gammapy 1.2 (1.1 for Theta2 and Skymaps)
  • 1D analysis (point-like IRF)

Analysis Results

Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here.

Theta2 plot

Significance map

Excess map

Spectral results

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