add prepare_header step into quantization_and_output

ObservationSim/sim_steps/readout_output.py

@@ -35,6 +35,9 @@ def apply_gain(self, chip, filt, tel, pointing, catalog, obs_param):
 
 def quantization_and_output(self, chip, filt, tel, pointing, catalog, obs_param):
 
+    if not hasattr(self, 'h_ext'):
+        _, _ = self.prepare_headers(chip=chip, pointing=pointing)
+
     if obs_param["format_output"] == True:
         self.chip_output.Log_info("  Apply 1*16 format")
         chip.img = chip_utils.formatOutput(GSImage=chip.img)