diff --git a/scripts/makecds.py b/scripts/makecds.py
new file mode 100644
index 0000000..5aaa4db
--- /dev/null
+++ b/scripts/makecds.py
@@ -0,0 +1,683 @@
+import cdspyreadme
+from astropy.io import fits
+from astropy.table import Table
+
+def make_euds_cds(infile1=None,outfile1=None,infile2=None,outfile2=None,fluxlim=None):
+
+    print("Reading {}".format(infile1))
+
+    #########################
+    """ Read eUDS catalog """
+    #########################
+    
+    tablemaker = cdspyreadme.CDSTablesMaker()
+
+    r = Table.read(infile1)
+    tab=Table([
+        r['ID_SRC'],           # 0
+        r['NAME'],             # 1
+        r['RA'],               # 2
+        r['DEC'],              # 3
+        r['RADEC_ERR'],        # 4 Positional error (68% confidence)
+        r['EXT'],              # 5 Source extent 
+        r['EXT_ERR'],          # 6 Extent uncertainty
+        r['EXT_LIKE'],         # 7
+        r['DET_LIKE'],         # 8 
+        r['ML_RATE'],          # 9
+        r['ML_RATE_ERR'],      # 10 
+        r['ML_CTS'],           # 11
+        r['ML_CTS_ERR'],       # 12 
+        r['ML_FLUX'],          # 13
+        r['ML_FLUX_ERR'],      # 14 
+        r['ML_EXP'],           # 15
+        r['ML_BKG'],           # 16
+        r['DR12_IAU_NAME'] ,   # 17 4XMM J021738.8-051257
+        r['DET_LIKE_1']    ,   # 18
+        r['ML_RATE_1']     ,   # 19
+        r['ML_RATE_ERR_1'] ,   # 20
+        r['ML_CTS_1']      ,   # 21
+        r['ML_CTS_ERR_1']  ,   # 22
+        r['ML_FLUX_1']     ,   # 23
+        r['ML_FLUX_ERR_1'] ,   # 24
+        r['ML_EXP_1']      ,   # 25
+        r['ML_BKG_1']      ,   # 26
+        r['DET_LIKE_2']    ,   # 27
+        r['ML_RATE_2']     ,   # 28
+        r['ML_RATE_ERR_2'] ,   # 29
+        r['ML_CTS_2']      ,   # 30
+        r['ML_CTS_ERR_2']  ,   # 31
+        r['ML_FLUX_2']     ,   # 32
+        r['ML_FLUX_ERR_2'] ,   # 33
+        r['ML_EXP_2']      ,   # 34
+        r['ML_BKG_2']      ,   # 35
+        r['DET_LIKE_3']    ,   # 36
+        r['ML_RATE_3']     ,   # 37
+        r['ML_RATE_ERR_3'] ,   # 38
+        r['ML_CTS_3']      ,   # 39
+        r['ML_CTS_ERR_3']  ,   # 40
+        r['ML_FLUX_3']     ,   # 41
+        r['ML_FLUX_ERR_3'] ,   # 42
+        r['ML_EXP_3']      ,   # 43
+        r['ML_BKG_3']      ,   # 44
+        r['DET_LIKE_4']    ,   # 45
+        r['ML_RATE_4']     ,   # 46
+        r['ML_RATE_ERR_4'] ,   # 47
+        r['ML_CTS_4']      ,   # 48
+        r['ML_CTS_ERR_4']  ,   # 49
+        r['ML_FLUX_4']     ,   # 50
+        r['ML_FLUX_ERR_4'] ,   # 51
+        r['ML_EXP_4']      ,   # 52
+        r['ML_BKG_4']      ,   # 53
+    ])
+
+    table = tablemaker.addTable(tab, name=outfile1, description="eUDS X-ray catalog, converted from eUDS.fits")
+
+    col = table.get_column('ID_SRC')           # 0
+    col.name="SrcID"
+    col.set_format("I3")
+    col.description="Source ID (ID_SRC)"
+
+    col = table.get_column('NAME')             # 1
+    col.name="Name"
+    col.set_format("A21")
+    col.description="Source name, eUDS JHHMMSS.s+DDMMSS"
+
+    col = table.get_column('RA')               # 2
+    col.set_format("F10.6")
+    col.name="RAdeg"
+    col.description="RA (J2000) (RA)"
+
+    col = table.get_column('DEC')              # 3
+    col.name="DEdeg"
+    col.set_format("F10.6")
+    col.description="Dec (J2000) (RA)"
+
+    col = table.get_column('RADEC_ERR')        # 4 Positional error (68% confidence)
+    col.set_format("F6.3")
+    col.name="ePos"
+    col.description="Positional uncertainty (RADEC_ERR)"
+
+    col = table.get_column('EXT')              # 5 Source extent 
+    col.set_format("F5.2")
+    col.name="Ext"
+    col.description="Source extent (EXT)"
+
+    col = table.get_column('EXT_ERR')          # 6 Extent uncertainty
+    col.set_format("F5.2")
+    col.name="e_Ext"
+    col.description="Extent error (EXT_ERR)"
+
+    col = table.get_column('EXT_LIKE')         # 7
+    col.set_format("F10.6")
+    col.name="ExtLike"
+    col.description="Extent likelihood (EXT_LIKE)"
+
+    col = table.get_column('DET_LIKE')         # 8 
+    col.set_format("F12.6")
+    col.name="DetLike"
+    col.description="Detection likelihood measured by PSF-fitting; 0.3-2.3 keV (DET_LIKE)"
+
+    col = table.get_column('ML_RATE')          # 9
+    col.name="MLRate"
+    col.set_format("F8.6")
+    col.description="Source count rate measured by PSF-fitting; 0.3-2.3 keV (ML_RATE)"
+
+    col = table.get_column('ML_RATE_ERR')      # 10 
+    col.name="e_MLRate"
+    col.set_format("F8.6")
+    col.description="1 sigma count rate error; 0.3-2.3 keV (ML_RATE_ERR)"
+
+    col = table.get_column('ML_CTS')           # 11
+    col.name="MLCts"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.3-2.3 keV (ML_CTS)"
+
+    col = table.get_column('ML_CTS_ERR')       # 12 
+    col.name="e_MLCts"
+    col.set_format("F11.6")
+    col.description="1 sigma counts error; 0.3-2.3 keV (ML_CTS_ERR)"
+
+    col = table.get_column('ML_FLUX')          # 13
+    col.name="MLFlux"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.3-2.3 keV (ML_FLUX)"
+
+    col = table.get_column('ML_FLUX_ERR')      # 14 
+    col.name="e_MLFlux"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.3-2.3 keV (ML_FLUX_ERR)"
+
+    col = table.get_column('ML_EXP')           # 15
+    col.name="MLExp"
+    col.set_format("E12.6")
+    col.description="Vignetted exposure value; 0.3-2.3 keV (ML_EXP)"
+
+    col = table.get_column('ML_BKG')           # 16
+    col.name="MLBkg"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.3-2.3 keV (ML_BKG)"
+
+    col = table.get_column('DR12_IAU_NAME')    # 17 4XMM J021738.8-051257
+    col.name="DR12Name"
+    col.set_format("A21")
+    col.description="4XMM-DR12 source name used for forced photometry (DR12_IAU_NAME)"
+
+    col = table.get_column('DET_LIKE_1')       # 18
+    col.name="DetLike1"
+    col.set_format("F12.6")
+    col.description="Detection likelihood measured by PSF-fitting; 0.3-0.6 keV (DET_LIKE_1)"
+
+    col = table.get_column('ML_RATE_1')        # 19
+    col.name="MLRate1"
+    col.set_format("F8.6")
+    col.description="Source count rate measured by PSF-fitting; 0.3-0.6 keV (ML_RATE_1)"
+
+    col = table.get_column('ML_RATE_ERR_1')    # 20
+    col.name="e_MLRate1"
+    col.set_format("F8.6")
+    col.description="1 sigma count rate error; 0.3-0.6 keV (ML_RATE_ERR_1)"
+
+    col = table.get_column('ML_CTS_1')         # 21
+    col.name="MLCts1"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.3-0.6 keV (ML_CTS_1)"
+
+    col = table.get_column('ML_CTS_ERR_1')     # 22
+    col.name="e_MLCts1"
+    col.set_format("F11.6")
+    col.description="1 sigma counts error; 0.3-0.6 keV (ML_CTS_ERR_1)"
+
+    col = table.get_column('ML_FLUX_1')        # 23
+    col.name="MLFlux1"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.3-0.6 keV (ML_FLUX_1)"
+
+    col = table.get_column('ML_FLUX_ERR_1')    # 24
+    col.name="e_MLFlux1"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.3-0.6 keV (ML_FLUX_ERR_1)"
+
+    col = table.get_column('ML_EXP_1')         # 25
+    col.name="MLExp1"
+    col.set_format("E12.6")
+    col.description="Vignetted exposure value; 0.3-0.6 keV (ML_EXP_1)"
+
+    col = table.get_column('ML_BKG_1')         # 26
+    col.name="MLBkg1"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.3-0.6 keV (ML_BKG_1)"
+
+    col = table.get_column('DET_LIKE_2')       # 27
+    col.name="DetLike2"
+    col.set_format("F12.6")
+    col.description="Detection likelihood measured by PSF-fitting; 0.6-2.3 keV (DET_LIKE_2)"
+
+    col = table.get_column('ML_RATE_2')        # 28
+    col.name="MLRate2"
+    col.set_format("F8.6")
+    col.description="Source count rate measured by PSF-fitting; 0.6-2.3 keV (ML_RATE_2)"
+
+    col = table.get_column('ML_RATE_ERR_2')    # 29
+    col.name="e_MLRate2"
+    col.set_format("F8.6")
+    col.description="1 sigma count rate error; 0.6-2.3 keV (ML_RATE_ERR_2)"
+
+    col = table.get_column('ML_CTS_2')         # 30
+    col.name="MLCts2"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.6-2.3 keV (ML_CTS_2)"
+
+    col = table.get_column('ML_CTS_ERR_2')     # 31
+    col.name="e_MLCts2"
+    col.set_format("F11.6")
+    col.description="1 sigma counts error; 0.6-2.3 keV (ML_CTS_ERR_2)"
+
+    col = table.get_column('ML_FLUX_2')        # 32
+    col.name="MLFlux2"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.6-2.3 keV (ML_FLUX_2)"
+
+    col = table.get_column('ML_FLUX_ERR_2')    # 33
+    col.name="e_MLFlux2"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.6-2.3 keV (ML_FLUX_ERR_2)"
+
+    col = table.get_column('ML_EXP_2')         # 34
+    col.name="MLExp2"
+    col.set_format("E12.6")
+    col.description="Vignetted exposure value; 0.6-2.3 keV (ML_EXP_2)"
+
+    col = table.get_column('ML_BKG_2')         # 35
+    col.name="MLBkg2"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.6-2.3 keV (ML_BKG_2)"
+
+    col = table.get_column('DET_LIKE_3')       # 36
+    col.name="DetLike3"
+    col.set_format("F12.6")
+    col.description="Detection likelihood measured by PSF-fitting; 2.3-5 keV (DET_LIKE_3)"
+
+    col = table.get_column('ML_RATE_3')        # 37
+    col.name="MLRate3"
+    col.set_format("F8.6")
+    col.description="Source count rate measured by PSF-fitting; 2.3-5 keV (ML_RATE_3)"
+
+    col = table.get_column('ML_RATE_ERR_3')    # 38
+    col.name="e_MLRate3"
+    col.set_format("F8.6")
+    col.description="1 sigma count rate error; 2.3-5 keV (ML_RATE_ERR_3)"
+
+    col = table.get_column('ML_CTS_3')         # 39
+    col.name="MLCts3"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 2.3-5 keV (ML_CTS_3)"
+
+    col = table.get_column('ML_CTS_ERR_3')     # 40
+    col.name="e_MLCts3"
+    col.set_format("F11.6")
+    col.description="1 sigma counts error; 2.3-5 keV (ML_CTS_ERR_3)"
+
+    col = table.get_column('ML_FLUX_3')        # 41
+    col.name="MLFlux3"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 2.3-5 keV (ML_FLUX_3)"
+
+    col = table.get_column('ML_FLUX_ERR_3')    # 42
+    col.name="e_MLFlux3"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 2.3-5 keV (ML_FLUX_ERR_3)"
+
+    col = table.get_column('ML_EXP_3')         # 43
+    col.name="MLExp3"
+    col.set_format("E12.6")
+    col.description="Vignetted exposure value; 2.3-5 keV (ML_EXP_3)"
+
+    col = table.get_column('ML_BKG_3')         # 44
+    col.name="MLBkg3"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 2.3-5 keV (ML_BKG_3)"
+
+    col = table.get_column('DET_LIKE_4')       # 45
+    col.name="DetLike4"
+    col.set_format("F12.6")
+    col.description="Detection likelihood measured by PSF-fitting; 5-8 keV (DET_LIKE_4)"
+
+    col = table.get_column('ML_RATE_4')        # 46
+    col.name="MLRate4"
+    col.set_format("F8.6")
+    col.description="Source count rate measured by PSF-fitting; 5-8 keV (ML_RATE_4)"
+
+    col = table.get_column('ML_RATE_ERR_4')    # 47
+    col.name="e_MLRate4"
+    col.set_format("F8.6")
+    col.description="1 sigma count rate error; 5-8 keV (ML_RATE_ERR_4)"
+
+    col = table.get_column('ML_CTS_4')         # 48
+    col.name="MLCts4"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 5-8 keV (ML_CTS_4)"
+
+    col = table.get_column('ML_CTS_ERR_4')     # 49
+    col.name="e_MLCts4"
+    col.set_format("F11.6")
+    col.description="1 sigma counts error; 5-8 keV (ML_CTS_ERR_4)"
+
+    col = table.get_column('ML_FLUX_4')        # 50
+    col.name="MLFLux4"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 5-8 keV (ML_FLUX_4)"
+
+    col = table.get_column('ML_FLUX_ERR_4')    # 51
+    col.name="e_MLFlux4"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 5-8 keV (ML_FLUX_ERR_4)"
+
+    col = table.get_column('ML_EXP_4')         # 52
+    col.name="MLExp4"
+    col.set_format("E12.6")
+    col.description="Vignetted exposure value; 5-8 keV (ML_EXP_4)"
+
+    col = table.get_column('ML_BKG_4')         # 53
+    col.name="MLBkg4"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 5-8 keV (ML_BKG_4)"
+
+    #####################################
+    """ Read 4XMM-DR12 forced catalog """
+    #####################################
+    print("Reading {}".format(infile2))
+
+    r = Table.read(infile2)
+    tab2=Table([
+            r['DR12_IAU_NAME'],                                         # 0
+            r['DR12_SRCID'],                                            # 1
+            r['DR12_RA'],                                               # 2
+            r['DR12_DEC'],                                              # 3
+            r['DR12_RADEC_ERR'],                                        # 4
+            r['DET_LIKE_0']   ,   # 5
+            r['DET_LIKE_1']    ,   # 6
+            r['DET_LIKE_2']    ,   # 7
+            r['DET_LIKE_3']    ,   # 8
+            r['DET_LIKE_4']    ,   # 9
+            r['ML_RATE_0']     ,   # 10
+            r['ML_RATE_1']     ,   # 11
+            r['ML_RATE_2']     ,   # 12
+            r['ML_RATE_3']     ,   # 13
+            r['ML_RATE_4']     ,   # 14
+            r['ML_RATE_ERR_0']     ,   # 15
+            r['ML_RATE_ERR_1']     ,   # 16
+            r['ML_RATE_ERR_2']     ,   # 17
+            r['ML_RATE_ERR_3']     ,   # 18
+            r['ML_RATE_ERR_4']     ,   # 19
+            r['ML_CTS_0']     ,   # 20 
+            r['ML_CTS_1']     ,   # 21
+            r['ML_CTS_2']     ,   # 22
+            r['ML_CTS_3']     ,   # 23
+            r['ML_CTS_4']     ,   # 24
+            r['ML_CTS_ERR_0']     ,   # 25 
+            r['ML_CTS_ERR_1']     ,   # 26 
+            r['ML_CTS_ERR_2']     ,   # 27 
+            r['ML_CTS_ERR_3']     ,   # 28 
+            r['ML_CTS_ERR_4']     ,   # 29 
+            r['ML_FLUX_0']     ,   # 30
+            r['ML_FLUX_1']     ,   # 31
+            r['ML_FLUX_2']     ,   # 32
+            r['ML_FLUX_3']     ,   # 33
+            r['ML_FLUX_4']     ,   # 34
+            r['ML_FLUX_ERR_0']     ,   # 35
+            r['ML_FLUX_ERR_1']     ,   # 36
+            r['ML_FLUX_ERR_2']     ,   # 37
+            r['ML_FLUX_ERR_3']     ,   # 38
+            r['ML_FLUX_ERR_4']     ,   # 39
+            r['ML_EXP_0']     ,   # 40
+            r['ML_EXP_1']     ,   # 41
+            r['ML_EXP_2']     ,   # 42
+            r['ML_EXP_3']     ,   # 43
+            r['ML_EXP_4']     ,   # 44
+            r['ML_BKG_0']     ,   # 45
+            r['ML_BKG_1']     ,   # 46
+            r['ML_BKG_2']     ,   # 47
+            r['ML_BKG_3']     ,   # 48
+            r['ML_BKG_4']     ,   # 49
+            r['CONF'],   # 50
+    ])
+
+    table2 = tablemaker.addTable(tab2, name=outfile2, description="eUDS forced photometry of 4XMM-DR12 sources, converted from eUDS_4XMM-DR12.fits")
+
+    col = table2.get_column('DR12_IAU_NAME')                                         # 0
+    col.name="DR12Name"
+    col.set_format("A21")
+    col.description="4XMM-DR12 source IAU name (DR12_IAU_NAME)"
+
+    col = table2.get_column('DR12_SRCID')                                            # 1
+    col.name="DR12SrcID"
+    col.set_format("I15")
+    col.description="Source ID (64-bit integer, DR12_SRCID)"
+
+    col = table2.get_column('DR12_RA')                                               # 2
+    col.name="RAdeg"
+    col.set_format("F10.6")
+    col.description="Right ascension (J2000, DR12_RA)"
+
+    col = table2.get_column('DR12_DEC')                                              # 3
+    col.name="DEdeg"
+    col.set_format("F10.6")
+    col.description="Declination (J2000, DR12_DEC)"
+
+    col = table2.get_column('DR12_RADEC_ERR')                                        # 4
+    col.name="ePos"
+    col.set_format("F6.3")
+    col.description="Positional error (DR12_RADEC_ERR)"
+
+    col = table2.get_column('DET_LIKE_0')      # 5
+    col.name="DetLike0"
+    col.set_format("F12.6")
+    col.description="Detection likelihood; 0.3-2.3 keV (DET_LIKE_0)"
+
+    col = table2.get_column('DET_LIKE_1')       # 6
+    col.name="DetLike1"
+    col.set_format("F12.6")
+    col.description="Detection likelihood; 0.3-0.6 keV (DET_LIKE_1)"
+
+    col = table2.get_column('DET_LIKE_2')       # 7
+    col.name="DetLike2"
+    col.set_format("F12.6")
+    col.description="Detection likelihood; 0.6-2.3 keV (DET_LIKE_2)"
+
+    col = table2.get_column('DET_LIKE_3')       # 8
+    col.name="DetLike3"
+    col.set_format("F12.6")
+    col.description="Detection likelihood; 2.3-5 keV (DET_LIKE_3)"
+
+    col = table2.get_column('DET_LIKE_4')       # 9
+    col.name="DetLike4"
+    col.set_format("F12.6")
+    col.description="Detection likelihood; 5-8 keV (DET_LIKE_4)"
+
+    col = table2.get_column('ML_RATE_0')        # 10
+    col.name="MLRate0"
+    col.set_format("F10.6")
+    col.description="Source count rate measured by PSF-fitting; 0.3-2.3 keV (ML_RATE_0)"
+
+    col = table2.get_column('ML_RATE_1')        # 11
+    col.name="MLRate1"
+    col.set_format("F10.6")
+    col.description="Source count rate measured by PSF-fitting; 0.3-0.6 keV (ML_RATE_1)"
+
+    col = table2.get_column('ML_RATE_2')        # 12
+    col.name="MLRate2"
+    col.set_format("F10.6")
+    col.description="Source count rate measured by PSF-fitting; 0.6-2.3 keV (ML_RATE_2)"
+
+    col = table2.get_column('ML_RATE_3')        # 13
+    col.name="MLRate3"
+    col.set_format("F10.6")
+    col.description="Source count rate measured by PSF-fitting; 2.3-5 keV (ML_RATE_3)"
+
+    col = table2.get_column('ML_RATE_4')        # 14
+    col.name="MLRate4"
+    col.set_format("F10.6")
+    col.description="Source count rate measured by PSF-fitting; 5-8 keV (ML_RATE_4)"
+
+    col = table2.get_column('ML_RATE_ERR_0')        # 15
+    col.name="e_MLRate0"
+    col.set_format("F10.6")
+    col.description="1 sigma count rate error; 0.3-2.3 keV (ML_RATE_ERR_0)"
+
+    col = table2.get_column('ML_RATE_ERR_1')        # 16
+    col.name="e_MLRate1"
+    col.set_format("F10.6")
+    col.description="1 sigma count rate error; 0.3-0.6 keV (ML_RATE_ERR_1)"
+
+    col = table2.get_column('ML_RATE_ERR_2')        # 17
+    col.name="e_MLRate2"
+    col.set_format("F10.6")
+    col.description="1 sigma count rate error; 0.6-2.3 keV (ML_RATE_ERR_2)"
+
+    col = table2.get_column('ML_RATE_ERR_3')        # 18
+    col.name="e_MLRate3"
+    col.set_format("F10.6")
+    col.description="1 sigma count rate error; 2.3-5 keV (ML_RATE_ERR_3)"
+
+    col = table2.get_column('ML_RATE_ERR_4')        # 19
+    col.name="e_MLRate4"
+    col.set_format("F10.6")
+    col.description="1 sigma count rate error; 5-8 keV (ML_RATE_ERR_4)"
+
+    col = table2.get_column('ML_CTS_0')        # 20 
+    col.name="MLCts0"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.3-2.3 keV (ML_CTS_0)"
+
+    col = table2.get_column('ML_CTS_1')        # 21
+    col.name="MLCts1"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.3-0.6 keV (ML_CTS_1)"
+
+    col = table2.get_column('ML_CTS_2')        # 22
+    col.name="MLCts2"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 0.6-2.3 keV (ML_CTS_2)"
+
+    col = table2.get_column('ML_CTS_3')        # 23
+    col.name="MLCts3"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 2.3-5 keV (ML_CTS_3)"
+
+    col = table2.get_column('ML_CTS_4')        # 24
+    col.name="MLCts4"
+    col.set_format("F11.6")
+    col.description="Source net counts measured from count rate; 5-8 keV (ML_CTS_4)"
+
+    col = table2.get_column('ML_CTS_ERR_0')        # 25 
+    col.name="e_MLCts0"
+    col.set_format("F12.6")
+    col.description="1 sigma counts error; 0.3-2.3 keV (ML_CTS_ERR_0)"
+
+    col = table2.get_column('ML_CTS_ERR_1')        # 26 
+    col.name="e_MLCts1"
+    col.set_format("F12.6")
+    col.description="1 sigma counts error; 0.3-0.6 keV (ML_CTS_ERR_1)"
+
+    col = table2.get_column('ML_CTS_ERR_2')        # 27 
+    col.name="e_MLCts2"
+    col.set_format("F12.6")
+    col.description="1 sigma counts error; 0.6-2.3 keV (ML_CTS_ERR_2)"
+
+    col = table2.get_column('ML_CTS_ERR_3')        # 28 
+    col.name="e_MLCts3"
+    col.set_format("F12.6")
+    col.description="1 sigma counts error; 2.3-5 keV (ML_CTS_ERR_3)"
+
+    col = table2.get_column('ML_CTS_ERR_4')        # 29 
+    col.name="e_MLCts4"
+    col.set_format("F12.6")
+    col.description="1 sigma counts error; 5-8 keV (ML_CTS_ERR_4)"
+
+    col = table2.get_column('ML_FLUX_0')        # 30
+    col.name="MLFlux0"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.3-2.3 keV (ML_FLUX_0)"
+
+    col = table2.get_column('ML_FLUX_1')        # 31
+    col.name="MLFlux1"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.3-0.6 keV (ML_FLUX_1)"
+
+    col = table2.get_column('ML_FLUX_2')        # 32
+    col.name="MLFlux2"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 0.6-2.3 keV (ML_FLUX_2)"
+
+    col = table2.get_column('ML_FLUX_3')        # 33
+    col.name="MLFlux3"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 2.3-5 keV (ML_FLUX_3)"
+
+    col = table2.get_column('ML_FLUX_4')        # 34
+    col.name="MLFlux4"
+    col.set_format("E12.6")
+    col.description="Source flux converted from count rate; 5-8 keV (ML_FLUX_4)"
+
+    col = table2.get_column('ML_FLUX_ERR_0')        # 35
+    col.name="e_MLFlux0"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.3-2.3 keV (ML_FLUX_ERR_0)"
+
+    col = table2.get_column('ML_FLUX_ERR_1')        # 36
+    col.name="e_MLFlux1"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.3-0.6 keV (ML_FLUX_ERR_1)"
+
+    col = table2.get_column('ML_FLUX_ERR_2')        # 37
+    col.name="e_MLFlux2"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 0.6-2.3 keV (ML_FLUX_ERR_2)"
+
+    col = table2.get_column('ML_FLUX_ERR_3')        # 38
+    col.name="e_MLFlux3"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 2.3-5 keV (ML_FLUX_ERR_3)"
+
+    col = table2.get_column('ML_FLUX_ERR_4')        # 39
+    col.name="e_MLFlux4"
+    col.set_format("E12.6")
+    col.description="1 sigma flux error; 5-8 keV (ML_FLUX_ERR_4)"
+
+    col = table2.get_column('ML_EXP_0')        # 40
+    col.name="MLExp0"
+    col.set_format("E10.6")
+    col.description="Vignetted exposure value; 0.3-2.3 keV (ML_EXP_0)"
+
+    col = table2.get_column('ML_EXP_1')        # 41
+    col.name="MLExp1"
+    col.set_format("E10.6")
+    col.description="Vignetted exposure value; 0.3-0.6 keV (ML_EXP_1)"
+
+    col = table2.get_column('ML_EXP_2')        # 42
+    col.name="MLExp2"
+    col.set_format("E10.6")
+    col.description="Vignetted exposure value; 0.6-2.3 keV (ML_EXP_2)"
+
+    col = table2.get_column('ML_EXP_3')        # 43
+    col.name="MLExp3"
+    col.set_format("E10.6")
+    col.description="Vignetted exposure value; 2.3-5 keV (ML_EXP_3)"
+
+    col = table2.get_column('ML_EXP_4')        # 44
+    col.name="MLExp4"
+    col.set_format("E10.6")
+    col.description="Vignetted exposure value; 5-8 keV (ML_EXP_4)"
+
+    col = table2.get_column('ML_BKG_0')        # 45
+    col.name="MLBkg0"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.3-2.3 keV (ML_BKG_0)"
+
+    col = table2.get_column('ML_BKG_1')        # 46
+    col.name="MLBkg1"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.3-0.6 keV (ML_BKG_1)"
+
+    col = table2.get_column('ML_BKG_2')        # 47
+    col.name="MLBkg2"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 0.6-2.3 keV (ML_BKG_2)"
+
+    col = table2.get_column('ML_BKG_3')        # 48
+    col.name="MLBkg3"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 2.3-5 keV (ML_BKG_3)"
+
+    col = table2.get_column('ML_BKG_4')        # 49
+    col.name="MLBkg4"
+    col.set_format("F12.6")
+    col.description="Background at the source position; 5-8 keV (ML_BKG_4)"
+
+    col = table2.get_column('CONF')   # 50
+    col.name="Conf"
+    col.set_format("I1")
+    col.description="True if source is located within 60 arcsec from another source with 0.3-2.3 keV flux > 1E−14 erg/s/cm**2 (CONF)"
+
+    tablemaker.writeCDSTables()
+
+    # Customize ReadMe output
+    
+    tablemaker.title = "eUDS: The SRG/eROSITA X-ray Survey of the UKIDSS Ultra Deep Survey Field. Catalogue of Sources."
+    tablemaker.author = 'Krivonos R.'
+    tablemaker.authors="Gilfanov M., Medvedev P., Sazonov S., Sunyaev R."
+    tablemaker.date = 2024
+    tablemaker.bibcode = "2024MNRAS.528.1264K"
+    tablemaker.keywords = "catalogues - surveys - X-ray: general"
+    tablemaker.abstract = "The eROSITA X-ray telescope on board the Spectrum-Roentgen-Gamma (SRG) spacecraft observed the field of the UKIDSS Ultra-Deep Survey (UDS) in August-September 2019, during its flight to Sun-Earth L2 point. The resulting eROSITA UDS (or eUDS) survey was thus the first eROSITA X-ray imaging survey, which demonstrated the capability of the telescope to perform uniform observations of large sky areas. With a moderate single-camera exposure of 150 ks, eUDS covered ∼5 deg^2^ with the limiting flux ranging between 4*10^-15^ and 5*10^-14^erg/cm^2^/s in the 0.3-2.3 keV band. We present a catalogue of 647 sources detected at likelihood >10 (∼4sigma) during the eUDS. The catalogue provides information on the source fluxes in the main energy band 0.3-2.3 keV and forced photometry in a number of bands between 0.3 and 8 keV. Using the deeper 4XMM-DR12 catalogue, we have identified 22 strongly variable objects that have brightened or faded by at least a factor of ten during the eROSITA observations compared to previous observations by XMM-Newton. We also provide a catalogue of 22 sources detected by eROSITA in the hard energy band of 2.3-5 keV."
+    tablemaker.more_description = "The catalogue (fits/eUDS.fits) contains 647 sources detected in the 0.3-2.3 keV band using eSASS adopting a detection likelihood threshold of 10. The catalogue provides information on the source fluxes in the main energy band 0: 0.3-2.3 keV and forced photometry in 1: 0.3-0.6, 2: 0.6-2.3, 3: 2.3-5.0 and 4: 5.0-8.0 keV bands. The catalogue of forced PSF-fitting fluxes (fits/eUDS_4XMM-DR12.fits) contains 3603 XMM-Newton sources from 4XMM-DR12 catalogue with the eUDS exposure of more than 100 s. The catalogue provides information on the source forced photometry in the energy bands 0: 0.3-2.3, 1: 0.3-0.6, 2: 0.6-2.3, 3: 2.3-5.0 and 4: 5.0-8.0 keV. To warn about the potential impact of spatial confusion with a nearby bright source, we added a confusion flag (see original paper for details)."
+    
+    tablemaker.putRef("J/A+A/661/A1","The eFEDS X-ray catalogs (Brunner+, 2022)")
+    tablemaker.putRef("J/A+A/661/A38","SRG/ART-XC 1st year all-sky X-ray survey (Pavlinsky+, 2022)")
+    tablemaker.putRef("https://www.srg.cosmos.ru/eUDS","eUDS catalog")
+    #tablemaker.makeReadMe()
+    with open("ReadMe.tmp","w") as fd:
+        tablemaker.makeReadMe(out=fd)
+
+make_euds_cds(infile1="../uds/products/eUDS.fits",outfile1='table1.dat',
+              infile2="../uds/products/eUDS_4XMM-DR12.fits",outfile2='table2.dat')