Flat Field coefficients computation using the model-based approach + requered changes in containers files. And the script to compute raise and decay times.

Cleaned branch.

Author: Anastasiia Mikhno

src/nectarchain/user_scripts/ajardinb/flatfield_example.py

@@ -104,6 +104,7 @@ def get_bad_pixels(output_from_FlatFieldComponent):
     hi_lo = get_hi_lo_ratio(output_from_FlatFieldComponent)
 
     amp_int_per_pix_per_event = FlatFieldOutput.amp_int_per_pix_per_event[:, :, :]
+
     mean_amp_int_per_pix = np.mean(amp_int_per_pix_per_event, axis=0)
     mean_amp = np.mean(mean_amp_int_per_pix, axis=1)
     std_amp = np.std(mean_amp_int_per_pix, axis=1)

src/nectarchain/user_scripts/amikhno/analysis_photostat_simple.py

@@ -24,6 +24,8 @@
 
 from nectarchain.data.container import GainContainer
 
+plt.style.use("/home/amikhno/Downloads/plot_style.mpltstyle")
+
 parser = argparse.ArgumentParser(description="Run NectarCAM photostatistics analysis")
 
 parser.add_argument("-r", "--run-number", required=True, help="Run number")
@@ -58,9 +60,10 @@
 run_path = args.run_path + f"/runs/NectarCAM.Run{run_number}.0000.fits.fz"
 filename_ps = (
     args.analysis_file + f"/PhotoStat/PhotoStatisticNectarCAM_FFrun{run_number}"
-    f"_GlobalPeakWindowSum_window_width_8_Pedrun{run_number}_FullWaveformSum.h5"
+    f"_LocalPeakWindowSum_window_shift_4_window_width_16_Pedrun{run_number}_FullWaveformSum.h5"
 )
 
+
 # --- Example usage ---
 if not os.path.exists(filename_ps):
     print(f"[INFO] {filename_ps} not found, running gain_PhotoStat_computation.py...")
@@ -80,10 +83,12 @@
         "--SPE_run_number",
         run_spe_number,
         "--method",
-        "GlobalPeakWindowSum",
+        "LocalPeakWindowSum",
         "--extractor_kwargs",
-        '{"window_width":8}',
+        '{"window_shift": 4, "window_width":16}',
         "--overwrite",
+        "--SPE_config",
+        "nominal",
         "-v",
         "INFO",
         "--reload_events",

src/nectarchain/user_scripts/amikhno/compute_rise_decay_times.py

@@ -12,6 +12,8 @@
 from ctapipe_io_nectarcam import NectarCAMEventSource
 from matplotlib import pyplot as plt
 
+plt.style.use("/home/amikhno/Downloads/plot_style.mpltstyle")
+
 parser = argparse.ArgumentParser(
     description="Run computation of mean parameters of waveform"
 )