Fast Offline Analysis gw-alerts241125 November 2024 P071

From my_wiki
Revision as of 12:57, 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

S241125n FOA theta square excess over0.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