Latest revision as of 11:28, 29 April 2022

Overview[edit]

Source-dependent/independent analysis using LSTCHAIN v0.9.4 (MC production processed with LSTCHAIN v0.9.3, R0 to DL1 process by LSTOSA using LSTCHAIN v0.9.2)
Cleaning method tail cut 84 and dynamic cleaning. Also, the MC files are tuned to match the NSB in the FoV and the PSF.
Observations used in the analysis:

Monte Carlo information[edit]

Link to MC files used:

- Particle types: gamma off0.4deg, proton

- ZD (deg): 40

- AZ (deg): 180

- Other information: MC production through PR in lstMCpipe Github webpage (https://github.com/cta-observatory/lstmcpipe/pull/185)

- MC files are tuned to match the NSB in the FoV and the PSF.

DL1 data[edit]

Processed with LSTOSA using LSTCHAIN v0.9.2.

- Dynamic cleaning and tail cut cleaning with pedestal threshold applied by the the standard parameters using pipeline LSTOSA.

 "dynamic_cleaning": {
   "apply": true,
   "threshold": 267,
   "fraction_cleaning_intensity": 0.03
 }

Real data[edit]

Parameters for the tail cut cleaning with pedestal threshold

 "tailcuts_clean_with_pedestal_threshold": {
   "picture_thresh":8,
   "boundary_thresh":4,
   "sigma":2.5,
   "keep_isolated_pixels":false,
   "min_number_picture_neighbors":2,
   "use_only_main_island":false,
   "delta_time": 2
 }

Original DL1a files (processed up to DL1 by LSTOSA using LSTCHAIN v0.9.2)

/fefs/aswg/data/real/DL1/{}/v0.9/tailcut84/dl1_LST-1.RunXXXXX.XXXX.h5

MC data[edit]

Tunned NSB

   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 0.6,
   "extra_bias_in_dim_pixels": 0.269,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 0.566,

Tunned PSF

   "increase_psf": true,
   "smeared_light_fraction": 0.175

Original DL1a files

/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/{}/

Random forest[edit]

Dispnorm parameterisation

Source-dependent

- LSTCHAIN v0.9.4

- Config file:

/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json

- path:

/fefs/aswg/workspace/arnau.aguasca/Analysis/models/20200609_XXXX/srcdep/v0.9.4/

Source-independent

- Processed with lstMCpipe using LSTCHAIN v0.9.3.

- Path:

/fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/

DL1 to DL2 data[edit]

Real data[edit]

Source-dependent

- LSTCHAIN v0.9.4

- Config file:

//fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json

Source-independent

- LSTCHAIN v0.9.4

- Config file:

/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_standard_config_tailcut84_dispnorm_v0.9.3.json

MC data[edit]

Source-dependent

- LSTCHAIN v0.9.4

- Config file:

/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json

Source-independent

- Processed with lstMCpipe using LSTCHAIN v0.9.3.

- Config file:

/fefs/aswg/data/mc/analysis_logs/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/off0.4deg/lstchain_config_NSBtuning.json

DL3 data selection[edit]

Source-dependent

- Quality cuts

 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2]
   }
 },
 "DL3FixedCuts": {
   "fixed_gh_cut": 0.8,
   "fixed_alpha_cut": 15,
   "allowed_tels": [1]
 }

Source-independent

- Quality cuts

High-level analysis[edit]

Source-dependent

- IRF: point-like, single-offset

/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/XXXX

- Produced DL3 files

/fefs/aswg/workspace/arnau.aguasca/XXXX

Source-independent

Analysis results[edit]

High level analysis performed with gammapy-v0.19

Alpha plot[edit]

Filters in EventSelector applied:

 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2]
   }
 }

See section Crab Check to know more about the justification of the applied cuts. Stacked alpha plots:

Alpha plot using runs in 20210809, 20210810 and 20210812.
Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812

Daily alpha plots:

Alpha plot using runs in 20210809.
Alpha plot using runs in 20210810.
Alpha plot using runs in 20210812.

Energy threshold for the analysis

Theta plot[edit]

Stacked theta plots:

Theta2 plot using runs in 20210809, 20210810 and 20210812.
Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812

Daily theta plots:

Theta2 plot using runs in 20210809.
Theta2 plot using runs in 20210810.
Theta2 plot using runs in 20210812.

Energy threshold for the analysis

Spectral results[edit]

Source-dependent

Stacked SED (Aug 9-12)

Daily SEDs

20210809 SED
20210810 SED
20210812 SED

HESS and MAGIC comparison

HESS/Fermi and LST-1 SED.
MAGIC and LST-1 SED (Aug 9-12).

Daily MAGIC/LST-1 SEDs:

MAGIC and LST-1 SED (Aug 9).
MAGIC and LST-1 SED (Aug 10).
MAGIC and LST-1 SED (Aug 12).

Source-independent

Stacked SED (Aug 9-12)

Daily SEDs

20210809 SED
20210810 SED
20210812 SED

HESS and MAGIC comparison

HESS/Fermi and LST-1 SED.
MAGIC and LST-1 SED (Aug 9-12).

No daily SED can be plotted because the SED from LST-1 using gc=0.6 and tc=0.2 deg are not reliable.

Light curve results[edit]

Source-dependent

HESS/Fermi and LST-1 light curve.
MAGIC and LST-1 light curve.

Source-independent

- Notice that since the daily SEDs with quality cuts (gc=0.6 and tc=0.2 deg) are not reliable, the daily flux are computed with the stacked SED using data between 20210809-12!!!
- No comparison with MAGIC light curve is done because the energy threshold using (gc=0.6 and tc=0.2 deg) is higher than 100 GeV.

HESS/Fermi and LST-1 light curve.

Crab Check[edit]

Run selection using Crab_data_selection_notebook: https://indico.cta-observatory.org/event/3984/contributions/32878/attachments/20998/29572/Crab_data_selection_notebook.pdf

Threshold for the rate30_mask is set to 4.3. I exclude run 6016 because I obtain a strange SED at low energies.

Data reduction Specifications[edit]

Source-dependent

Cuts at alpha cut = 15 deg are studied.

Energy threshold for the analysis

gammaness cut>0.6 and alpha cut<15 deg, int>50 p.e.
gammaness cut>0.7 and alpha cut<15 deg, int>50 p.e.
gammaness cut>0.8 and alpha cut<15 deg, int>50 p.e.
gammaness cut>0.9 and alpha cut<15 deg, int>100 p.e.

Source-independent

Cuts at theta cut = 0.2 deg are studied.

Energy threshold for the analysis

gammaness cut>0.6 and theta cut< 0.2 deg, int>50 p.e.
gammaness cut>0.7 and theta cut< 0.2 deg, int>50 p.e.
gammaness cut>0.8 and theta cut< 0.2 deg, int>50 p.e.
gammaness cut>0.9 and theta cut<0.2 deg, int>100 p.e.

I decide to use only the loose cuts in order to have a low energy threshold.

Spectral results[edit]

Source-dependent

Loose cuts, 20 November 2020. 2 bins per decade
Hard cuts, 20 November 2020. 2 bins per decade

Loose cuts, September 2021. 2 bins per decade
Hard cuts, September 2021. 2 bins per decade

Loose cuts, March 2022. 2 bins per decade
Hard cuts, March 2022. 2 bins per decade

All runs. 4 bins per decade

Conclusions: I decide to apply the following cuts: gammaness cut = 0.8, alpha cut = 15 deg and intensity = 50 p.e.

Source-independent

File:Crab srcind 20201120 loose.jpg
Loose cuts, 20 November 2020
File:Crab srcind 202109 loose.jpg
Loose cuts, September 2021
File:Crab srcind 202203 loose.jpg
Loose cuts, March 2022

File:Crab srcind all loose.jpg
All runs

@@ Line 146: / Line 146: @@
    }
 See section [[#Crab Check|Crab Check]] to know more about the justification of the applied cuts.
+Stacked alpha plots:
 <div><ul>
 <li style="display: inline-block;"> [[File:Alpha plot 09to12.png|thumb|none|400px|Alpha plot using runs in 20210809, 20210810 and 20210812.]] </li>
 <li style="display: inline-block;"> [[File:Cumulative SNr sig gc08 ac15.png|thumb|none|450px|Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812]] </li>
+</ul></div>
+Daily alpha plots:
+<div><ul>
+<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210809 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210809.]] </li>
+<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210810 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210810.]] </li>
+<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210812 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210812.]] </li>
 </ul></div>
 * Energy threshold for the analysis
 [[File:Energy thd srcdep RSOph int50 gcut08 acut15.jpg|400px]]
 === Theta plot ===
+Stacked theta plots:
+<div><ul>
+<li style="display: inline-block;"> [[File:Theta2 plot 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210809, 20210810 and 20210812.]] </li>
+<li style="display: inline-block;"> [[File:Cumulative SNr sqrt-ts 3off gc06 tc02.png|thumb|none|450px|Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812]] </li>
+</ul></div>
+Daily theta plots:
+<div><ul>
+<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210809 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210809.]] </li>
+<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210810 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210810.]] </li>
+<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210812 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210812.]] </li>
+</ul></div>
+* Energy threshold for the analysis
+[[File:Energy thd srcind RSOph int50 gcut06 tcut02.jpg|400px]]
 === Spectral results ===
 *; Source-dependent
+<div><ul>
+<li style="display: inline-block;"> [[File:RSOph SED 09to12 int50 gc08 ac15.png|thumb|none|400px|Stacked SED (Aug 9-12)]] </li>
+</ul></div>
+:Daily SEDs
+<div><ul>
+<li style="display: inline-block;"> [[File:Srcdep spec 20210809 gc08 ac15.png|thumb|none|400px|20210809 SED]] </li>
+<li style="display: inline-block;"> [[File:Srcdep spec 20210810 gc08 ac15.png|thumb|none|400px|20210810 SED]] </li>
+<li style="display: inline-block;"> [[File:Srcdep spec 20210812 gc08 ac15.png|thumb|none|400px|20210812 SED]]</li>
+</ul></div>
+:*; HESS and MAGIC comparison
+<div><ul>
+<li style="display: inline-block;"> [[File:LST-1HESSLAT SED RSOph srcdep 09to12 int50 gc08 ac15.png|thumb|none|400px|HESS/Fermi and LST-1 SED.]] </li>
+<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 09to12 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9-12).]] </li>
+</ul></div>
+::Daily MAGIC/LST-1 SEDs:
+<div><ul>
+<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210809 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9).]] </li>
+<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210810 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 10).]] </li>
+<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210812 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 12).]] </li>
+</ul></div>
 *; Source-independent
+<div><ul>
+<li style="display: inline-block;"> [[File:Srcind spec 20210809to12 gc06 tc02.png|thumb|none|400px|Stacked SED (Aug 9-12)]] </li>
+</ul></div>
+:Daily SEDs
+<div><ul>
+<li style="display: inline-block;"> [[File:Srcind spec 20210809 gc06 tc02.png|thumb|none|400px|20210809 SED]] </li>
+<li style="display: inline-block;"> [[File:Srcind spec 20210810 gc06 tc02.png|thumb|none|400px|20210810 SED]] </li>
+<li style="display: inline-block;"> [[File:Srcind spec 20210812 gc06 tc02.png|thumb|none|400px|20210812 SED]]</li>
+</ul></div>
+:*; HESS and MAGIC comparison
+<div><ul>
+<li style="display: inline-block;"> [[File:LST-1HESSLAT SED RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|400px|HESS/Fermi and LST-1 SED.]] </li>
+<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9-12).]] </li>
+</ul></div>
+No daily SED can be plotted because the SED from LST-1 using gc=0.6 and tc=0.2 deg are not reliable.
 === Light curve results ===
 *; Source-dependent
+<div><ul>
+<li style="display: inline-block;"> [[File:LST-1HESSLAT lightcurve RSOph srcind 09to12 int50 gc08 ac15.png|thumb|none|500px|HESS/Fermi and LST-1 light curve.]] </li>
+<li style="display: inline-block;"> [[File:LST-1MAGIC lightcurve RSOph srcind 09to12 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 light curve.]] </li>
+</ul></div>
 *; Source-independent
+:;- Notice that since the daily SEDs with quality cuts (gc=0.6 and tc=0.2 deg) are not reliable, the daily flux are computed with the stacked SED using data between 20210809-12!!!
+:;- No comparison with MAGIC light curve is done because the energy threshold using (gc=0.6 and tc=0.2 deg) is higher than 100 GeV.
+<div><ul>
+<li style="display: inline-block;"> [[File:LST-1HESSLAT lightcurve RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|500px|HESS/Fermi and LST-1 light curve.]] </li>
+</ul></div>

Difference between revisions of "Analysis A. Aguasca-Cabot lstchain v0.9.X"

Latest revision as of 11:28, 29 April 2022

Contents

Overview[edit]

Monte Carlo information[edit]

DL1 data[edit]

Real data[edit]

MC data[edit]

Random forest[edit]

DL1 to DL2 data[edit]

Real data[edit]

MC data[edit]

DL3 data selection[edit]

High-level analysis[edit]

Analysis results[edit]

Alpha plot[edit]

Theta plot[edit]

Spectral results[edit]

Light curve results[edit]

Crab Check[edit]

Data reduction Specifications[edit]

Spectral results[edit]

Navigation menu

Personal tools

Namespaces

Variants

Views

More

Search

Navigation

Tools