Flat Field coefficients computation + the rise and decay times

src/nectarchain/data/container/gain_container.py

@@ -41,6 +41,10 @@ class GainContainer(NectarCAMContainer):
     )
     low_gain = Field(type=np.ndarray, dtype=np.float64, ndim=2, description="low gain")
     pixels_id = Field(type=np.ndarray, dtype=np.uint16, ndim=1, description="pixel ids")
+    charge = Field(type=np.ndarray, dtype=np.float64, ndim=1, description="charge")
+    charge_std = Field(
+        type=np.ndarray, dtype=np.float64, ndim=1, description="charge_std"
+    )
 
     @classmethod
     def from_hdf5(cls, path, group_name="data"):

src/nectarchain/makers/component/photostatistic_algorithm.py

@@ -56,6 +56,8 @@ def __init__(
             high_gain=np.zeros((self.npixels, 3)),
             low_gain=np.zeros((self.npixels, 3)),
             pixels_id=self._pixels_id,
+            charge=np.zeros((self.npixels, 1)),
+            charge_std=np.zeros((self.npixels, 1)),
         )
 
     @classmethod
@@ -141,6 +143,9 @@ def run(self, pixels_id: np.ndarray = None, **kwargs) -> None:
             (self.gainLG * mask, gainLG_err, gainLG_err)
         ).T
         self._results.is_valid = mask
+        log.info("Trying to write charges")
+        self._results.charge = self.meanChargeHG
+        self._results.charge_std = self.sigmaChargeHG
 
         figpath = kwargs.get("figpath", False)
         if figpath:

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)