Build_plots_from_galvotable00 Input votable that contains information from ellipse and galfit executions 2012-09-07 07:32:21.998 UTC <?xml version='1.0'?> <VOTABLE version="1.1" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.ivoa.net/xml/VOTable/v1.1 http://www.ivoa.net/xml/VOTable/v1.1" xmlns="http://www.ivoa.net/xml/VOTable/v1.1"> <!-- ! VOTable written by STIL version 3.0-3 (uk.ac.starlink.votable.VOTableWriter) ! at 2012-08-09T13:16:24 !--> <RESOURCE> <TABLE nrows="3"> <DESCRIPTION> Table of sources detected in image </DESCRIPTION> <PARAM arraysize="2" datatype="char" name="QUERY_STATUS" value="OK"/> <FIELD arraysize="88" datatype="char" name="imagen"/> <FIELD arraysize="1" datatype="char" name="band"/> <FIELD arraysize="85" datatype="char" name="cat"/> <FIELD arraysize="85" datatype="char" name="sex"/> <FIELD arraysize="85" datatype="char" name="sexparam"/> <FIELD arraysize="91" datatype="char" name="checkimage_name"/> <FIELD arraysize="88" datatype="char" name="galfit"/> <FIELD arraysize="92" datatype="char" name="galfit_exp"/> <FIELD arraysize="92" datatype="char" name="galfit2nd"/> <FIELD arraysize="92" datatype="char" name="galfitgal"/> <FIELD arraysize="45" datatype="char" name="psf"/> <FIELD arraysize="92" datatype="char" name="model"/> <FIELD arraysize="94" datatype="char" name="ellip_tab"/> <FIELD arraysize="94" datatype="char" name="ellip_resultheader"/> <FIELD arraysize="94" datatype="char" name="ellip_resultdata"/> <FIELD arraysize="94" datatype="char" name="ellip_tab2"/> <FIELD arraysize="94" datatype="char" name="ellip_resultheader2"/> <FIELD arraysize="94" datatype="char" name="ellip_resultdata2"/> <FIELD arraysize="94" datatype="char" name="ellip_script1"/> <FIELD arraysize="94" datatype="char" name="ellip_script2"/> <FIELD arraysize="94" datatype="char" name="ellip_plot"/> <FIELD datatype="short" name="CIG"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="float" name="ra"/> <FIELD datatype="float" name="dec"/> <FIELD datatype="float" name="mag_zeropoint"/> <FIELD datatype="float" name="satur_level"/> <FIELD datatype="float" name="gain"/> <FIELD datatype="float" name="mB"/> <FIELD datatype="float" name="mBcorr"/> <FIELD datatype="short" name="nx"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="short" name="ny"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="float" name="zp"/> <FIELD datatype="short" name="xconvsize"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="short" name="yconvsize"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="short" name="bar"> <VALUES null='-32768'/> </FIELD> <FIELD datatype="int" name="NUMBER" ucd="meta.record"> <DESCRIPTION>Running object number</DESCRIPTION> <VALUES null='-2147483648'/> </FIELD> <FIELD datatype="float" name="X_IMAGE" ucd="pos.cartesian.x;pos.barycenter;instr.det;meta.main" unit="pix"> <DESCRIPTION>Object position along x</DESCRIPTION> </FIELD> <FIELD datatype="float" name="Y_IMAGE" ucd="pos.cartesian.y;pos.barycenter;instr.det;meta.main" unit="pix"> <DESCRIPTION>Object position along y</DESCRIPTION> </FIELD> <FIELD datatype="double" name="ALPHA_SKY" ucd="pos.eq.ra;pos.barycenter" unit="deg"> <DESCRIPTION>Right ascension of barycenter (native)</DESCRIPTION> </FIELD> <FIELD datatype="double" name="DELTA_SKY" ucd="pos.eq.dec;pos.barycenter" unit="deg"> <DESCRIPTION>Declination of barycenter (native)</DESCRIPTION> </FIELD> <FIELD datatype="float" name="MAG_AUTO" ucd="phot.mag;meta.main" unit="mag"> <DESCRIPTION>Kron-like elliptical aperture magnitude</DESCRIPTION> </FIELD> <FIELD datatype="float" name="MAGERR_AUTO" ucd="stat.stdev;phot.mag;meta.main" unit="mag"> <DESCRIPTION>RMS error for AUTO magnitude</DESCRIPTION> </FIELD> <FIELD datatype="int" name="ISOAREA_IMAGE" ucd="phys.area" unit="pix2"> <DESCRIPTION>Isophotal area above Analysis threshold</DESCRIPTION> <VALUES null='-2147483648'/> </FIELD> <FIELD datatype="float" name="FWHM_IMAGE" ucd="phys.size.diameter;instr.det.psf" unit="pix"> <DESCRIPTION>FWHM assuming a gaussian core</DESCRIPTION> </FIELD> <FIELD datatype="float" name="BACKGROUND" ucd="instr.skyLevel" unit="ct"> <DESCRIPTION>Background at centroid position</DESCRIPTION> </FIELD> <FIELD datatype="float" name="CLASS_STAR" ucd="src.class.starGalaxy"> <DESCRIPTION>S/G classifier output</DESCRIPTION> </FIELD> <FIELD datatype="float" name="FLUX_BEST" ucd="phot.flux" unit="ct"> <DESCRIPTION>Best of FLUX_AUTO and FLUX_ISOCOR</DESCRIPTION> </FIELD> <FIELD datatype="float" name="ELLIPTICITY" ucd="src.ellipticity;instr.det"> <DESCRIPTION>1 - B_IMAGE/A_IMAGE</DESCRIPTION> </FIELD> <FIELD datatype="float" name="THETA_IMAGE" ucd="pos.posAng;instr.det;meta.main" unit="deg"> <DESCRIPTION>Position angle (CCW/x)</DESCRIPTION> </FIELD> <FIELD datatype="float" name="FLUX_RADIUS" ucd="phys.size.radius;instr.det" unit="pix"> <DESCRIPTION>Fraction-of-light radii</DESCRIPTION> </FIELD> <FIELD datatype="float" name="KRON_RADIUS" ucd="arith.factor;arith.ratio"> <DESCRIPTION>Kron apertures in units of A or B</DESCRIPTION> </FIELD> <FIELD datatype="float" name="A_IMAGE" ucd="phys.size.smajAxis;instr.det;meta.main" unit="pix"> <DESCRIPTION>Profile RMS along major axis</DESCRIPTION> </FIELD> <FIELD datatype="float" name="B_IMAGE" ucd="phys.size.sminAxis;instr.det;meta.main" unit="pix"> <DESCRIPTION>Profile RMS along minor axis</DESCRIPTION> </FIELD> <FIELD datatype="float" name="ERRA_IMAGE" ucd="stat.stdev;stat.max;pos.errorEllipse;instr.det;meta.main" unit="pix"> <DESCRIPTION>RMS position error along major axis</DESCRIPTION> </FIELD> <FIELD datatype="float" name="ERRB_IMAGE" ucd="stat.stdev;stat.min;pos.errorEllipse;instr.det;meta.main" unit="pix"> <DESCRIPTION>RMS position error along minor axis</DESCRIPTION> </FIELD> <FIELD datatype="float" name="axes_rate"/> <FIELD datatype="double" name="dist"/> <FIELD datatype="float" name="rotated_theta"/> <FIELD datatype="double" name="a_pix"/> <FIELD datatype="double" name="sma0"/> <FIELD datatype="double" name="minsma"/> <FIELD datatype="double" name="maxsma"/> <FIELD datatype="double" name="step"/> <FIELD datatype="double" name="ellipticity_not0"/> <FIELD datatype="double" name="sma02"/> <FIELD datatype="double" name="minsma2"/> <FIELD datatype="double" name="maxsma2"/> <FIELD datatype="double" name="step2"/> <FIELD datatype="float" name="x_2"> <DESCRIPTION>Position x</DESCRIPTION> </FIELD> <FIELD datatype="float" name="y_2"> <DESCRIPTION>Position y</DESCRIPTION> </FIELD> <FIELD datatype="float" name="mauto_2"> <DESCRIPTION>Integrated magnitude</DESCRIPTION> </FIELD> <FIELD datatype="float" name="flux_rad0_5_2_big"> <DESCRIPTION>R_s (disk scale-length) * 1.5 [pix]</DESCRIPTION> </FIELD> <FIELD datatype="float" name="flux_rad0_5_2_small"> <DESCRIPTION>R_s (disk scale-length) * 0.25 [pix]</DESCRIPTION> </FIELD> <FIELD datatype="float" name="flux_rad0_5_2_medium"> <DESCRIPTION>R_s (disk scale-length) * 0.5 [pix]</DESCRIPTION> </FIELD> <FIELD datatype="float" name="axes_rate_2"> <DESCRIPTION>Axis ratio (b/a)</DESCRIPTION> </FIELD> <FIELD datatype="float" name="theta_2"> <DESCRIPTION>Position angle</DESCRIPTION> </FIELD> <FIELD datatype="float" name="sky_2"> <DESCRIPTION>Sky background at center of fitting region</DESCRIPTION> </FIELD> <FIELD datatype="double" name="x_disk"/> <FIELD datatype="double" name="y_disk"/> <FIELD datatype="float" name="Mag_disk"/> <FIELD datatype="float" name="hR_disk"/> <FIELD datatype="float" name="ba_disk"/> <FIELD datatype="float" name="PA_disk"/> <FIELD datatype="double" name="x_bul"/> <FIELD datatype="double" name="y_bul"/> <FIELD datatype="float" name="Mag_bul"/> <FIELD datatype="float" name="re_bul"/> <FIELD datatype="float" name="n_bul"/> <FIELD datatype="float" name="ba_bul"/> <FIELD datatype="float" name="PA_bul"/> <FIELD datatype="double" name="x_bar"/> <FIELD datatype="double" name="y_bar"/> <FIELD datatype="float" name="Mag_bar"/> <FIELD datatype="float" name="re_bar"/> <FIELD datatype="float" name="n_bar"/> <FIELD datatype="float" name="ba_bar"/> <FIELD datatype="float" name="PA_bar"/> <FIELD datatype="double" name="hr_disk_arcsec"/> <FIELD datatype="double" name="re_bul_arcsec"/> <FIELD datatype="double" name="re_bar_arcsec"/> <DATA> <TABLEDATA> <TR> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/fbCIG0011_g.fit</TD> <TD>g</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g.cat</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g.sex</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0000.param</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_seg.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_exp.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_BuD.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_gal.galfit</TD> <TD>/home/julian/Documents/psf/CIG0011_g_psf.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_model.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_cd.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_data.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_2.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_cd_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_data_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_ellip1.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g_ellip2.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_g.png</TD> <TD>11</TD> <TD>3.6340444</TD> <TD>-0.7356951</TD> <TD>28.337</TD> <TD>62160.0</TD> <TD>3.855</TD> <TD>14.8</TD> <TD>14.11</TD> <TD>950</TD> <TD>951</TD> <TD>24.0077</TD> <TD>95</TD> <TD>95</TD> <TD>0</TD> <TD>1</TD> <TD>461.176</TD> <TD>464.708</TD> <TD>3.6327212</TD> <TD>-0.7373948</TD> <TD>13.9152</TD> <TD>0.0018</TD> <TD>12425</TD> <TD>192.7</TD> <TD>6.947886</TD> <TD>0.0</TD> <TD>587137.6</TD> <TD>0.604</TD> <TD>75.3</TD> <TD>43.97</TD> <TD>3.87</TD> <TD>50.991</TD> <TD>20.205</TD> <TD>0.0824</TD> <TD>0.0333</TD> <TD>0.39624637</TD> <TD>7.754454436015714</TD> <TD>-14.699997</TD> <TD>78.14472912597657</TD> <TD>60.56216507263185</TD> <TD>1.3024121520996095</TD> <TD>78.14472912597655</TD> <TD>7.814472912597656</TD> <TD>0.6039999723434448</TD> <TD>50.79407393188477</TD> <TD>1.953618228149414</TD> <TD>119.17071191711428</TD> <TD>1.953618228149414</TD> <TD>460.2907</TD> <TD>463.3397</TD> <TD>14.01</TD> <TD>50.40435</TD> <TD>8.400725</TD> <TD>16.80145</TD> <TD>0.4413</TD> <TD>-11.961</TD> <TD>9.32</TD> <TD>462.5996</TD> <TD>463.047</TD> <TD>14.069</TD> <TD>48.0003</TD> <TD>0.429</TD> <TD>-12.126</TD> <TD>458.3548</TD> <TD>463.5709</TD> <TD>15.7355</TD> <TD>25.0611</TD> <TD>25.0611</TD> <TD>0.4507</TD> <TD>-12.3401</TD> <TD>-999.0</TD> <TD>-999.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>19.00811933898926</TD> <TD>9.924195602416992</TD> <TD>-39.204</TD> </TR> <TR> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/fbCIG0011_i.fit</TD> <TD>i</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i.cat</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i.sex</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0000.param</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_seg.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_exp.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_BDB.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_gal.galfit</TD> <TD>/home/julian/Documents/psf/CIG0011_i_psf.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_model.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_cd.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_data.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_2.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_cd_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_data_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_ellip1.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i_ellip2.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_i.png</TD> <TD>11</TD> <TD>3.6340444</TD> <TD>-0.7356951</TD> <TD>27.731</TD> <TD>51453.0</TD> <TD>6.565</TD> <TD>14.8</TD> <TD>14.11</TD> <TD>950</TD> <TD>950</TD> <TD>23.402</TD> <TD>95</TD> <TD>95</TD> <TD>1</TD> <TD>1</TD> <TD>460.526</TD> <TD>464.399</TD> <TD>3.6326604</TD> <TD>-0.737462</TD> <TD>13.2804</TD> <TD>0.0023</TD> <TD>7375</TD> <TD>35.28</TD> <TD>7.460361</TD> <TD>0.0</TD> <TD>602903.1</TD> <TD>0.513</TD> <TD>77.2</TD> <TD>34.322</TD> <TD>4.75</TD> <TD>33.114</TD> <TD>16.112</TD> <TD>0.0672</TD> <TD>0.0343</TD> <TD>0.4865616</TD> <TD>8.079840762117804</TD> <TD>-12.800003</TD> <TD>62.28742950439453</TD> <TD>48.27275786590577</TD> <TD>1.038123825073242</TD> <TD>62.28742950439453</TD> <TD>6.228742950439453</TD> <TD>0.5130000114440918</TD> <TD>40.48682917785644</TD> <TD>1.5571857376098632</TD> <TD>94.98832999420166</TD> <TD>1.5571857376098632</TD> <TD>459.9873</TD> <TD>464.5281</TD> <TD>13.4899</TD> <TD>34.0089</TD> <TD>5.66815</TD> <TD>11.3363</TD> <TD>0.4445</TD> <TD>-11.4531</TD> <TD>0.815</TD> <TD>462.4898</TD> <TD>464.4001</TD> <TD>13.5804</TD> <TD>43.8435</TD> <TD>0.4418</TD> <TD>-11.9599</TD> <TD>459.9158</TD> <TD>465.435</TD> <TD>17.1831</TD> <TD>5.2484</TD> <TD>5.2484</TD> <TD>0.5734</TD> <TD>-2.6001</TD> <TD>459.0488</TD> <TD>464.0722</TD> <TD>14.6718</TD> <TD>24.6591</TD> <TD>24.6591</TD> <TD>0.4292</TD> <TD>-11.4535</TD> <TD>17.362025299072265</TD> <TD>2.0783664836883546</TD> <TD>9.76500343322754</TD> </TR> <TR> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/fbCIG0011_r.fit</TD> <TD>r</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r.cat</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r.sex</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0000.param</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_seg.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_exp.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_BuD.galfit</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_gal.galfit</TD> <TD>/home/julian/Documents/psf/CIG0011_r_psf.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_model.fits</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_cd.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_data.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_2.tab</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_cd_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_data_2.ellip</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_ellip1.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r_ellip2.py</TD> <TD>/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/workflows/CIG0011_r.png</TD> <TD>11</TD> <TD>3.6340444</TD> <TD>-0.7356951</TD> <TD>28.089</TD> <TD>50631.0</TD> <TD>4.6</TD> <TD>14.8</TD> <TD>14.11</TD> <TD>950</TD> <TD>950</TD> <TD>23.7596</TD> <TD>95</TD> <TD>95</TD> <TD>0</TD> <TD>1</TD> <TD>461.293</TD> <TD>464.531</TD> <TD>3.6327027</TD> <TD>-0.7373858</TD> <TD>13.4906</TD> <TD>0.0018</TD> <TD>10750</TD> <TD>172.15</TD> <TD>7.017813</TD> <TD>0.0</TD> <TD>690821.7</TD> <TD>0.606</TD> <TD>74.6</TD> <TD>38.147</TD> <TD>4.19</TD> <TD>45.671</TD> <TD>18.004</TD> <TD>0.0755</TD> <TD>0.0294</TD> <TD>0.39421076</TD> <TD>7.770138577238733</TD> <TD>-15.400002</TD> <TD>75.77915240478517</TD> <TD>58.7288431137085</TD> <TD>1.262985873413086</TD> <TD>75.77915240478517</TD> <TD>7.577915240478516</TD> <TD>0.6060000061988831</TD> <TD>49.25644906311036</TD> <TD>1.894478810119629</TD> <TD>115.56320741729736</TD> <TD>1.894478810119629</TD> <TD>460.1561</TD> <TD>464.0741</TD> <TD>13.6804</TD> <TD>39.4131</TD> <TD>6.56885</TD> <TD>13.1377</TD> <TD>0.4399</TD> <TD>-11.5987</TD> <TD>0.5109</TD> <TD>473.6168</TD> <TD>475.5044</TD> <TD>15.5186</TD> <TD>52.1401</TD> <TD>0.2826</TD> <TD>-16.5911</TD> <TD>459.4644</TD> <TD>464.0984</TD> <TD>13.6261</TD> <TD>52.741</TD> <TD>52.741</TD> <TD>0.447</TD> <TD>-10.1518</TD> <TD>-999.0</TD> <TD>-999.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>-99.0</TD> <TD>20.647479034423828</TD> <TD>20.885436447143555</TD> <TD>-39.204</TD> </TR> </TABLEDATA> </DATA> </TABLE> </RESOURCE> </VOTABLE> 2012-09-07 07:32:00.953 UTC makePlots_STDOUTmakePlots_STDERRFormat_conversionformatTableIn0formatTableOut0voTable0outputTable00org.purl.wf4ever.astrotavernastilts-activity0.9org.purl.wf4ever.astrotaverna.tpipe.FormatConversionActivity String net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Parallelize 1 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.ErrorBouncenet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Failovernet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Retry 1.0 1000 5000 0 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.InvokeformatTableIn_valuevalue00net.sf.taverna.t2.activitiesstringconstant-activity1.4net.sf.taverna.t2.activities.stringconstant.StringConstantActivity votable net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Parallelize 1 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.ErrorBouncenet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Failovernet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Retry 1.0 1000 5000 0 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.InvokeformatTableOut_valuevalue00net.sf.taverna.t2.activitiesstringconstant-activity1.4net.sf.taverna.t2.activities.stringconstant.StringConstantActivity csv net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Parallelize 1 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.ErrorBouncenet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Failovernet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Retry 1.0 1000 5000 0 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.InvokemakePlotscsvData0STDOUT00STDERR00 It uses makeplotlib (python) to build a plot for every galaxy. If NaN, Inf or Divisions by 0 appears during the calculations, some of the graphics are not represented. 2012-09-07 13:32:57.920 UTC net.sf.taverna.t2.activitiesexternal-tool-activity1.4net.sf.taverna.t2.activities.externaltool.ExternalToolActivity 789663B8-DA91-428A-9F7D-B3F3DA185FD4 default local <?xml version="1.0" encoding="UTF-8"?> <localInvocation><shellPrefix>/bin/sh -c</shellPrefix><linkCommand>/bin/ln -s %%PATH_TO_ORIGINAL%% %%TARGET_NAME%%</linkCommand></localInvocation> 3c9ba421-1a47-486b-a1b5-ccc3461eba81 python make_plots.py csvData 1200 1800 make_plots.py true false false UTF-8 false import sys from numpy import array, arange, sin, cos, zeros, nonzero import numpy from math import pi import pylab as P from scipy.special import gamma from numpy import exp from matplotlib.mlab import csv2rec #read data from file #header: number of lines to skip def read_ellip_data_table(lines, header=11): table=[] for line in lines[header:]: data=line.split() try: #table.append(map(lambda x:float(x),data)) table.append(data) except: print data #m=numpy.matrix(table, dtype=float) m=numpy.matrix(table) listColumns=m.H.tolist() return listColumns # READ GALFIT DATA filename=sys.argv[1] #filename='/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/testEllipse/galfit.csv' #filename='/home/julian/Documents/wf4ever/showcase61.2goldenExampler/2GE/testEllipse/galfit_RO2ndGE_fake.csv' galaxies = csv2rec(filename) #data from galfit mag_disk = galaxies['mag_disk'] # [14.069, 13.5804, 15.5186] mag_disk = array(mag_disk) hr_disk = galaxies['hr_disk'] #[48.0003, 43.8435, 52.1401] hr_disk = array(hr_disk) ba_disk = galaxies['ba_disk'] #[0.429, 0.4418, 0.2826] ba_disk = array(ba_disk) mag_bul = galaxies['mag_bul'] # [15.7355, 17.1831, 13.6261] mag_bul = array(mag_bul) re_bul = galaxies['re_bul'] # [25.0611, 5.2484, 52.741] re_bul = array(re_bul) n_bul = galaxies['n_bul'] #[25.0611, 5.2484, 52.741] n_bul = array(n_bul) ba_bul = galaxies['ba_bul'] #[0.4507, 0.5734, 0.447] ba_bul = array(ba_bul) mag_bar = galaxies['mag_bar'] #[-99.0, 14.6718, -99.0] mag_bar = array(mag_bar) re_bar = galaxies['re_bar'] #[-99.0, 24.6591, -99.0] re_bar = array(re_bar) n_bar = galaxies['n_bar'] #[-99.0, 24.6591, -99.0] n_bar = array(n_bar) ba_bar = galaxies['ba_bar'] #[-99.0, 0.4292, -99.0] ba_bar = array(ba_bar) #READ ELLIPSE RESULTS AND MAKE PLOT FOR EVERY GALAXY for pos_galaxy in range(0, galaxies.size, 1): filename = galaxies['ellip_resultdata2'][pos_galaxy] f=open(filename, 'r') lines=f.readlines() f.close() listColumnsEllip=read_ellip_data_table(lines, 0) #(0:39) #listColumnsEllip=read_ellip_data_table(lines, 41) #(41:100) #ELLIPSE RESULTS #data from ellipse try: sma = map(lambda x:float(x),listColumnsEllip[0]) #sma = [14.8005 , 34.534, 54.2675] sma = array(sma) mag = map(lambda x:float(x),listColumnsEllip[17])#mag = [-5.100928, -4.544101, -3.811355] mag= array(mag) errM = map(lambda x:float(x),listColumnsEllip[18]) #errM = [0.01270742, 0.02156172, 0.01376063] errM = array(errM) except ValueError: #Remove 'INDEF ' value from data file (columns 0, 17 and 18) while ('INDEF' in listColumnsEllip[0]): pos = listColumnsEllip[0].index('INDEF') del listColumnsEllip[0][pos] del listColumnsEllip[17][pos] del listColumnsEllip[18][pos] while ('INDEF' in listColumnsEllip[17]): pos = listColumnsEllip[17].index('INDEF') del listColumnsEllip[0][pos] del listColumnsEllip[17][pos] del listColumnsEllip[18][pos] while ('INDEF' in listColumnsEllip[18]): pos = listColumnsEllip[18].index('INDEF') del listColumnsEllip[0][pos] del listColumnsEllip[17][pos] del listColumnsEllip[18][pos] finally: sma = map(lambda x:float(x),listColumnsEllip[0]) #sma = [14.8005 , 34.534, 54.2675] sma = array(sma) mag = map(lambda x:float(x),listColumnsEllip[17])#mag = [-5.100928, -4.544101, -3.811355] mag= array(mag) errM = map(lambda x:float(x),listColumnsEllip[18]) #errM = [0.01270742, 0.02156172, 0.01376063] errM = array(errM) #data operations from ellipse a_arc = sma * 0.396 flux = 10**(mag/-2.5) zp = galaxies['zp'] # [24.007745571418038, 23.402005887495726, 23.75962950480297] zp = array(zp) mag_i = (-2.5) * numpy.log10(flux / 53.907456) + zp[pos_galaxy] mu = mag_i + 2.5 * numpy.log10(0.396 ** 2) #pasamos a mag/arcsec2 #data operations from galfit #sobreescribimos el resultado de galfit #sersic: m_ser = mag_bul[pos_galaxy] flux_ser = 10. ** (m_ser / (-2.5)) re_ser = re_bul[pos_galaxy] n_ser = n_bul[pos_galaxy] ba_ser = ba_bul[pos_galaxy] k_ser = 1.9992 * n_ser - 0.3271 #expdisk: m_exp = mag_disk[pos_galaxy] flux_exp = 10. ** (m_exp / (-2.5)) rs_exp = hr_disk[pos_galaxy] ba_exp = ba_disk[pos_galaxy] #bar m_bb = mag_bar[pos_galaxy] flux_bb = 10. ** (m_bb / (-2.5)) re_bb = re_bar[pos_galaxy] n_bb = n_bar[pos_galaxy] ba_bb = ba_bar[pos_galaxy] k_bb = 1.9992 * n_bb - 0.3271 xx = arange(1000.) * (max(a_arc) + 5.) / 1000. xxres = a_arc #funcion sersic mue_ser = flux_ser / (2. * pi * re_ser ** 2 * exp(k_ser) * n_ser * k_ser ** (-2. * n_ser) * gamma(2. * n_ser) * ba_ser) mu_ser = mue_ser * exp(-k_ser * ((xx / re_ser) ** (1. / n_ser) - 1.)) mu_ser_res = mue_ser * exp(-k_ser * ((xxres / re_ser) ** (1. / n_ser) - 1.)) #funcion expdisk mu0_exp = flux_exp / (2 * pi * rs_exp ** 2 * ba_exp) mu_exp = mu0_exp * exp(-xx / rs_exp) mu_exp_res = mu0_exp * exp(-xxres / rs_exp) #funcion bar num = sma.size if n_bb != -99.: mue_bb = flux_bb / (2. * pi * re_bb ** 2 * exp(k_bb) * n_bb * k_bb ** (-2. * n_bb) * gamma(2. * n_bb) * ba_bb) mu_bb = mue_bb * exp(-k_bb * ((xx / re_bb) ** (1. / n_bb) - 1.)) mu_bb_res = mue_bb * exp(-k_bb * ((xxres / re_bb) ** (1. / n_bb) - 1.)) #else: # mu_bb = zeros(1000., float) # mu_bb_res = zeros(num, float) if n_bb != -99.: difmu = -2.5 * numpy.log10(mu_exp_res + mu_ser_res + mu_bb_res) - mu else: difmu = -2.5 * numpy.log10(mu_exp_res + mu_ser_res) - mu #ellipse plot fig = P.figure() #figure(num=None, figsize=(8.27, 11.69), dpi=100) #ax1 = P.subplot(211) ax1 = P.subplot2grid((3,1), (0, 0), rowspan=2) P.errorbar(a_arc,mu,yerr=errM,capsize=3, label='ellipse', fmt='ro', markersize=2, color='black') #P.plot(x,numpy.cos(1.3*x),linewidth=3, color='red', label='zorder=1',zorder = 1) # bottom P.plot(xx,-2.5*numpy.log10(mu_ser),label='Bulge', color='red') # bottom P.plot(xx,-2.5*numpy.log10(mu_exp),label='Disk', color='green') # bottom if n_bb != -99.: P.plot(xx,-2.5*numpy.log10(mu_bb),label='Bar', color='yellow') P.plot(xx,-2.5*numpy.log10(mu_exp + mu_ser + mu_bb),label='Disk+bulge+bar', color='blue') else: P.plot(xx,-2.5*numpy.log10(mu_exp + mu_ser),label='Disk+bulge', color='blue') P.xlim(0., a_arc.max()+5.) #This 2 lines are commented because some plots are out of this range for some galaxies P.ylim(28., mu.min()-1.5) P.ylabel(r'$\mu(i -mag arsec^{-2})$') P.legend(fancybox=True, shadow=True, prop=dict(size=8)) ax2= P.subplot2grid((3,1), (2, 0), sharex=ax1) P.ylim(-2., 2) if numpy.isnan(difmu).sum() !=len(difmu) : P.plot(a_arc, difmu, 'o-', markersize=2) P.ylabel(r'$\Delta \mu(i -mag arsec^{-2})$') P.xlabel('r(arcsec)') # #P.show() P.savefig(galaxies['ellip_plot'][pos_galaxy]) csvData csvData true false false UTF-8 false false false false true true 0 false net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Parallelize 1 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.ErrorBouncenet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Failovernet.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.Retry 1.0 1000 5000 0 net.sf.taverna.t2.coreworkflowmodel-impl1.4net.sf.taverna.t2.workflowmodel.processor.dispatch.layers.InvokeFormat_conversionformatTableInformatTableIn_valuevalueFormat_conversionformatTableOutformatTableOut_valuevalueFormat_conversionvoTablevotablemakePlotscsvDataFormat_conversionoutputTablemakePlots_STDOUTmakePlotsSTDOUTmakePlots_STDERRmakePlotsSTDERR a45b5cd7-a7fa-4f1f-b792-ab9621b4af12 2012-08-22 13:00:22.916 UTC 1e4f46e2-14d5-4d60-81ae-c4c05f93e97c 2012-09-07 09:06:57.588 UTC 0fb71bc8-96fc-4056-a580-83be534b9578 2012-08-24 06:37:10.334 UTC Build plots from galfit and ellipse results 2012-08-30 09:50:17.730 UTC 49fd9419-883b-49a9-b63e-f8121acd815b 2012-08-24 06:48:48.690 UTC Julian Garrido 2012-08-24 06:48:45.384 UTC 51984655-6907-42f6-83fe-25b214a0d5e1 2012-08-22 12:49:51.737 UTC bff728ca-8d83-4ac7-aeb6-8f63467a9183 2012-08-22 12:46:51.798 UTC 5edd5e43-1a48-4502-a0bf-c1b95becc10d 2012-08-10 06:42:48.42 UTC 04950e69-3c7a-4b42-ad7b-0ab66d524143 2012-08-10 06:45:35.160 UTC 871d5b74-ed2e-4607-a0bd-75ee95e782ae 2012-08-09 11:19:58.595 UTC 14d1cacb-f3af-4e6a-bf29-4c4bfd476cc2 2012-08-09 12:52:54.336 UTC ef2224e8-e0d3-425d-a4da-104b0d30282a 2012-09-07 13:32:59.347 UTC 3f6b815b-91c6-4b72-a1d1-095c8589fb03 2012-08-30 09:50:40.172 UTC b2cd94c3-9e9f-4111-be25-bfe482117026 2012-08-22 12:49:12.753 UTC d4c562b2-2452-4dff-a70e-4aa82ee31a82 2012-08-22 12:48:41.150 UTC ee3e8f23-f378-4947-8a6e-f84f585575c6 2012-09-07 07:32:23.274 UTC 482ec305-de8d-4e5c-85f0-1a215817e20c 2012-08-09 12:59:09.553 UTC 89045f2e-fd0d-474a-9ac8-6e8c92368198 2012-08-22 12:54:43.817 UTC cb1d7196-261e-4e57-923c-3d01e9475341 2012-08-30 09:51:05.89 UTC c2df3dc5-0dfa-4b12-b72f-1500a2f470c9 2012-08-10 06:44:59.361 UTC Build plots from galfit and ellipse results (http://www.myexperiment.org/workflows/3068.html, http://www.myexperiment.org/workflows/3052.html). 2012-08-30 09:50:45.977 UTC 1efe009f-959c-4077-b132-bcea384d43a6 2012-09-07 07:26:00.509 UTC 677ebe6f-9358-4262-8216-ff4b27baf694 2012-08-09 11:20:52.391 UTC