diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py
index c62edf0cc356d61090df7e4f1a080416d8bb8bf9..3ed95dc5440ff5c9b810f0d0e8049b55d957fef2 100644
--- a/csst_mci_sim/csst_mci_sim.py
+++ b/csst_mci_sim/csst_mci_sim.py
@@ -821,7 +821,7 @@ def psfMaker_IDW(px, py, PSFMat, cen_col, cen_row, dn=5, IDWindex=3, OnlyNeighbo
         dist = np.sqrt((ref_col - cen_col[idx])
                        ** 2 + (ref_row - cen_row[idx])**2)
 
-        #if IDWindex == 1:
+        # if IDWindex == 1:
         psfWeight[ipsf] = dist**IDWindex
         # if IDWindex == 2:
         #     psfWeight[ipsf] = dist**2
@@ -1051,8 +1051,8 @@ class MCIsimulator():
             (self.information['ysize'], self.information['xsize']), dtype=float)
         return
 
-##########################################IDWindex#####################################
-##############################################################################
+# ###############################IDWindex#####################################
+# ############################################################################
     def _loadGhostModel(self):
         """
         Reads in a ghost model from a FITS file and stores the data to self.ghostModel.
@@ -1285,7 +1285,12 @@ class MCIsimulator():
             # SED flux*filter_flux*tel_eff
             eff_arr = Nphoton*ft_effnew*tel_effnew
 
-            Sumphoton[ch] = np.sum(eff_arr)
+            S    checkInputList(input_ra_list, input_nstars)
+    checkInputList(input_dec_list, input_nstars)
+    checkInputList(input_pmra_list, input_nstars)
+    checkInputList(input_pmdec_list, input_nstars)
+    checkInputList(input_rv_list, input_nstars)
+    checkInputList(input_parallax_list, input_nstars)umphoton[ch] = np.sum(eff_arr)
 
             # seven wavelength for the seleced filter
             iwave = self.filter_psf[ch]['psf_iwave'][:]
@@ -1818,16 +1823,12 @@ class MCIsimulator():
                     star_output[ch][int(nlayccd)-1, 0] = fpx  # ra
                     star_output[ch][int(nlayccd)-1, 1] = fpy  # dec
 
-            #######################################################################
+            # ################################################################
 
             # ###### use  SED_code to generate star SED   #######
             if get_file_extension(starcat) == '.fits':
                 umag = self.star['gmag'][j]; gmag = self.star['gmag'][j]; rmag = self.star['rmag'][j]; imag = self.star['imag'][j]; zmag = self.star['zmag'][j];
 
-
-
-
-
             # SED of j-th star
             else:
                 umag = self.star['umag'][j]