Fast Offline Analysis gw-alerts241125 November 2024 P071

From my_wiki
Revision as of 12:54, 27 November 2024 by Sweta.menon (talk | contribs)
Jump to: navigation, search

Go back to Logbook_of_ToO

General information

  • Name of the source: gw-alerts241125 or S241125n
  • Brief description of the source: Swift- BAT position(GCN Circular 38308) of a gravitational wave alert(GCN Circular 38305)
    • Object type : Binary Black hole merger
    • Redshift (z) : NA
    • Other relevant information
    • RA: 03h 52m 18.96s (hh mm ss), Dec: +69d 41’ 20.4″ (dd mm ss)
    • RA, Dec in deg (ICRS): 58.079, +69.689
  • Analysis by:
    • Analyzer 1 Sweta Menon (INAF-Rome - sweta.menon@inaf.it)

Data-taking information

  • Dates of data-taking: 2024-11-25
    • Date (time of observation h)" 3.4
      • ZD range (deg): Zenith Range
      • Observations: (see example below)
20241125:   [19799 19801 19802 19804 19805 19806 19807 19808 19809 19810 19811]
  • Data quality selection: O hours of good data it is not mandatory to apply data selection for these preliminary analyses. But if you do it and select good-quality runs, please also add the plot in which the distribution of the different quality parameters are shown together with the limit values.
  • Joint observations with MAGIC?: Yes except run 19879,19880
  • Information on data-taking (observations summaries): Bad weather with poor transmission rates Write here if there was any problem with data-taking. (check ELOG for weather info, rate values and e.g. add link https://www.lst1.iac.es/datacheck/observation_overview/2024/03/run_summary_2024_03_17.pdf)

General analysis information

  • source-independent analysis
  • lstchain version: (e.g. v0.10.10)
  • DL1b standard production (v0.10/tailcut84), run-wise merged files in:

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19799.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19800.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19801.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19802.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19803.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19804.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19805.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19806.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19807.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19808.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19809.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19810.h5

/fefs/aswg/data/real/DL1/20241125/v0.10/tailcut84/dl1_LST-1.Run19811.h5


Monte Carlo information

Write below the MC production ID among the existing MC

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


Random forest

Path to the models used: (example below)

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

DL2 data

The config file to produce the DL2 files must be the one used to produce the RF models. Should be located together with model's files.


Files in:

/fefs/aswg/workspace/sweta.menon/FOA/S241125n/DL2/*

DL3 data & IRF

Files in:

/fefs/aswg/workspace/sweta.menon/FOA/S241125n/DL3/*

Configuration used for producing IRF/DL3 files (you can leave standard settings):

 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.01, 1],
     "leakage_intensity_width_2": [0, 1],
     "event_type": [32, 32]
   }
 },
 "DL3Cuts": {
   "min_event_p_en_bin": 100,
   "min_gh_cut": 0.1,
   "max_gh_cut": 0.95,
   "min_theta_cut": 0.0,
   "max_theta_cut": 0.2,
   "fill_theta_cut": 0.2,
   "allowed_tels": [1]
 },
  • Standard cuts: Intensity > 50 p.e., 70 % efficiency for both gammaness and theta gamma-ray selection energy-dependent cuts.
  • Recommended to produce unzipped DL3 files since the Gammapy analysis Gammapy is faster.
  • You can set the max theta cut (and filling value to the same value) to e.g. 0.2 deg. In this way, you can use more off-regions in Gammapy.
  • Point-like IRF.

High-level analysis and results

Analysis settings

  • Gammapy version: 1.1
  • 1D analysis (point-like IRF)

Theta2 plot

Settings:

  • Theta2 cut for signal extraction: X.XX deg2
  • Normalization range: X.X - Y.Y deg2
  • Number of off positions:
  • Energy range:

Daily results

  • YYYY-MM-DD

S241125n FOA heta square excess under 0.2TeV.jpg

File:S241125n FOA heta square excess over 0.2TeV.jpg

Stacked observations

Significance and Excess sky maps

Spectral results

Include the spectral parameters (see example below) and the SED plot (best-fit model + spectral flux points)

For example:

  type      name      value         unit        error      min    max frozen is_norm link
-------- ---------- ---------- -------------- --------- --------- --- ------ ------- ----
spectral      index 2.8255e+00                3.298e-01 1.000e+00 nan  False   False     
spectral  amplitude 1.7730e-09 TeV-1 s-1 cm-2 4.988e-10       nan nan  False    True     
spectral  reference 1.0000e+02            GeV 0.000e+00       nan nan   True   False       


Light curve