photometry_evaluation.ipynb

   "source": [
    "hdu = fits.open(cat_path_list_l1[0])\n",
    "print(\"Column names of L1 processed photometry catalog\", hdu[1].data.names)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a15e7b01",
   "metadata": {},
   "source": [
    "# 2. Properties of injected sources: (chip_id: 08, g band)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "id": "88efbb11",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Colum names in injected truth catalog:  ['obj_ID', 'ID_chip', 'filter', 'xImage', 'yImage', 'ra', 'dec', 'ra_orig', 'dec_orig', 'z', 'mag', 'obj_type', 'pm_ra', 'pm_dec', 'RV', 'parallax']\n",
      "CSST_MSC_MS_SCIE_20231022045923_20231022050153_10160000000_08_img_L1_injected.cat\n"
     ]
    },
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       "  \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<svg xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"597.215625pt\" height=\"330.391pt\" viewBox=\"0 0 597.215625 330.391\" xmlns=\"http://www.w3.org/2000/svg\" version=\"1.1\">\n",
       " <metadata>\n",
       "  <rdf:RDF xmlns:dc=\"http://purl.org/dc/elements/1.1/\" xmlns:cc=\"http://creativecommons.org/ns#\" xmlns:rdf=\"http://www.w3.org/1999/02/22-rdf-syntax-ns#\">\n",
       "   <cc:Work>\n",
       "    <dc:type rdf:resource=\"http://purl.org/dc/dcmitype/StillImage\"/>\n",
       "    <dc:date>2024-04-23T07:15:22.169281</dc:date>\n",
       "    <dc:format>image/svg+xml</dc:format>\n",
       "    <dc:creator>\n",
       "     <cc:Agent>\n",
       "      <dc:title>Matplotlib v3.8.4, https://matplotlib.org/</dc:title>\n",
       "     </cc:Agent>\n",
       "    </dc:creator>\n",
       "   </cc:Work>\n",
       "  </rdf:RDF>\n",
       " </metadata>\n",
       " <defs>\n",
       "  <style type=\"text/css\">*{stroke-linejoin: round; stroke-linecap: butt}</style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M 0 330.391 \n",
       "L 597.215625 330.391 \n",
       "L 597.215625 0 \n",
       "L 0 0 \n",
       "z\n",
       "\" style=\"fill: #ffffff\"/>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"patch_2\">\n",
       "    <path d=\"M 120.047812 288.430125 \n",
       "L 477.167813 288.430125 \n",
       "L 477.167813 22.318125 \n",
       "L 120.047812 22.318125 \n",
       "z\n",
       "\" style=\"fill: #ffffff\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_3\">\n",
       "    <path d=\"M 136.28054 288.430125 \n",
       "L 142.773631 288.430125 \n",
       "L 142.773631 288.430125 \n",
       "L 136.28054 288.430125 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_4\">\n",
       "    <path d=\"M 142.773631 288.430125 \n",
       "L 149.266722 288.430125 \n",
       "L 149.266722 284.342383 \n",
       "L 142.773631 284.342383 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_5\">\n",
       "    <path d=\"M 149.266722 288.430125 \n",
       "L 155.759813 288.430125 \n",
       "L 155.759813 282.298512 \n",
       "L 149.266722 282.298512 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_6\">\n",
       "    <path d=\"M 155.759813 288.430125 \n",
       "L 162.252903 288.430125 \n",
       "L 162.252903 282.298512 \n",
       "L 155.759813 282.298512 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_7\">\n",
       "    <path d=\"M 162.252903 288.430125 \n",
       "L 168.745994 288.430125 \n",
       "L 168.745994 284.342383 \n",
       "L 162.252903 284.342383 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_8\">\n",
       "    <path d=\"M 168.745994 288.430125 \n",
       "L 175.239085 288.430125 \n",
       "L 175.239085 286.386254 \n",
       "L 168.745994 286.386254 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_9\">\n",
       "    <path d=\"M 175.239085 288.430125 \n",
       "L 181.732176 288.430125 \n",
       "L 181.732176 284.342383 \n",
       "L 175.239085 284.342383 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_10\">\n",
       "    <path d=\"M 181.732176 288.430125 \n",
       "L 188.225267 288.430125 \n",
       "L 188.225267 284.342383 \n",
       "L 181.732176 284.342383 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_11\">\n",
       "    <path d=\"M 188.225267 288.430125 \n",
       "L 194.718358 288.430125 \n",
       "L 194.718358 276.166899 \n",
       "L 188.225267 276.166899 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_12\">\n",
       "    <path d=\"M 194.718358 288.430125 \n",
       "L 201.211449 288.430125 \n",
       "L 201.211449 280.254641 \n",
       "L 194.718358 280.254641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_13\">\n",
       "    <path d=\"M 201.211449 288.430125 \n",
       "L 207.70454 288.430125 \n",
       "L 207.70454 284.342383 \n",
       "L 201.211449 284.342383 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_14\">\n",
       "    <path d=\"M 207.70454 288.430125 \n",
       "L 214.197631 288.430125 \n",
       "L 214.197631 272.079157 \n",
       "L 207.70454 272.079157 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_15\">\n",
       "    <path d=\"M 214.197631 288.430125 \n",
       "L 220.690722 288.430125 \n",
       "L 220.690722 270.035286 \n",
       "L 214.197631 270.035286 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_16\">\n",
       "    <path d=\"M 220.690722 288.430125 \n",
       "L 227.183813 288.430125 \n",
       "L 227.183813 278.21077 \n",
       "L 220.690722 278.21077 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_17\">\n",
       "    <path d=\"M 227.183812 288.430125 \n",
       "L 233.676903 288.430125 \n",
       "L 233.676903 280.254641 \n",
       "L 227.183812 280.254641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_18\">\n",
       "    <path d=\"M 233.676903 288.430125 \n",
       "L 240.169994 288.430125 \n",
       "L 240.169994 276.166899 \n",
       "L 233.676903 276.166899 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_19\">\n",
       "    <path d=\"M 240.169994 288.430125 \n",
       "L 246.663085 288.430125 \n",
       "L 246.663085 267.991415 \n",
       "L 240.169994 267.991415 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_20\">\n",
       "    <path d=\"M 246.663085 288.430125 \n",
       "L 253.156176 288.430125 \n",
       "L 253.156176 272.079157 \n",
       "L 246.663085 272.079157 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_21\">\n",
       "    <path d=\"M 253.156176 288.430125 \n",
       "L 259.649267 288.430125 \n",
       "L 259.649267 276.166899 \n",
       "L 253.156176 276.166899 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_22\">\n",
       "    <path d=\"M 259.649267 288.430125 \n",
       "L 266.142358 288.430125 \n",
       "L 266.142358 286.386254 \n",
       "L 259.649267 286.386254 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_23\">\n",
       "    <path d=\"M 266.142358 288.430125 \n",
       "L 272.635449 288.430125 \n",
       "L 272.635449 259.815931 \n",
       "L 266.142358 259.815931 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_24\">\n",
       "    <path d=\"M 272.635449 288.430125 \n",
       "L 279.12854 288.430125 \n",
       "L 279.12854 265.947544 \n",
       "L 272.635449 265.947544 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_25\">\n",
       "    <path d=\"M 279.12854 288.430125 \n",
       "L 285.621631 288.430125 \n",
       "L 285.621631 263.903673 \n",
       "L 279.12854 263.903673 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_26\">\n",
       "    <path d=\"M 285.621631 288.430125 \n",
       "L 292.114722 288.430125 \n",
       "L 292.114722 270.035286 \n",
       "L 285.621631 270.035286 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_27\">\n",
       "    <path d=\"M 292.114722 288.430125 \n",
       "L 298.607812 288.430125 \n",
       "L 298.607812 265.947544 \n",
       "L 292.114722 265.947544 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_28\">\n",
       "    <path d=\"M 298.607813 288.430125 \n",
       "L 305.100903 288.430125 \n",
       "L 305.100903 259.815931 \n",
       "L 298.607813 259.815931 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_29\">\n",
       "    <path d=\"M 305.100903 288.430125 \n",
       "L 311.593994 288.430125 \n",
       "L 311.593994 249.596577 \n",
       "L 305.100903 249.596577 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_30\">\n",
       "    <path d=\"M 311.593994 288.430125 \n",
       "L 318.087085 288.430125 \n",
       "L 318.087085 239.377222 \n",
       "L 311.593994 239.377222 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_31\">\n",
       "    <path d=\"M 318.087085 288.430125 \n",
       "L 324.580176 288.430125 \n",
       "L 324.580176 227.113996 \n",
       "L 318.087085 227.113996 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_32\">\n",
       "    <path d=\"M 324.580176 288.430125 \n",
       "L 331.073267 288.430125 \n",
       "L 331.073267 229.157867 \n",
       "L 324.580176 229.157867 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_33\">\n",
       "    <path d=\"M 331.073267 288.430125 \n",
       "L 337.566358 288.430125 \n",
       "L 337.566358 245.508835 \n",
       "L 331.073267 245.508835 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_34\">\n",
       "    <path d=\"M 337.566358 288.430125 \n",
       "L 344.059449 288.430125 \n",
       "L 344.059449 194.41206 \n",
       "L 337.566358 194.41206 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_35\">\n",
       "    <path d=\"M 344.059449 288.430125 \n",
       "L 350.55254 288.430125 \n",
       "L 350.55254 208.719157 \n",
       "L 344.059449 208.719157 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_36\">\n",
       "    <path d=\"M 350.55254 288.430125 \n",
       "L 357.045631 288.430125 \n",
       "L 357.045631 163.753996 \n",
       "L 350.55254 163.753996 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_37\">\n",
       "    <path d=\"M 357.045631 288.430125 \n",
       "L 363.538722 288.430125 \n",
       "L 363.538722 165.797867 \n",
       "L 357.045631 165.797867 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_38\">\n",
       "    <path d=\"M 363.538722 288.430125 \n",
       "L 370.031813 288.430125 \n",
       "L 370.031813 145.359157 \n",
       "L 363.538722 145.359157 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_39\">\n",
       "    <path d=\"M 370.031813 288.430125 \n",
       "L 376.524903 288.430125 \n",
       "L 376.524903 110.613351 \n",
       "L 370.031813 110.613351 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_40\">\n",
       "    <path d=\"M 376.524903 288.430125 \n",
       "L 383.017994 288.430125 \n",
       "L 383.017994 141.271415 \n",
       "L 376.524903 141.271415 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_41\">\n",
       "    <path d=\"M 383.017994 288.430125 \n",
       "L 389.511085 288.430125 \n",
       "L 389.511085 90.174641 \n",
       "L 383.017994 90.174641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_42\">\n",
       "    <path d=\"M 389.511085 288.430125 \n",
       "L 396.004176 288.430125 \n",
       "L 396.004176 47.253351 \n",
       "L 389.511085 47.253351 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_43\">\n",
       "    <path d=\"M 396.004176 288.430125 \n",
       "L 402.497267 288.430125 \n",
       "L 402.497267 73.823673 \n",
       "L 396.004176 73.823673 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_44\">\n",
       "    <path d=\"M 402.497267 288.430125 \n",
       "L 408.990358 288.430125 \n",
       "L 408.990358 90.174641 \n",
       "L 402.497267 90.174641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_45\">\n",
       "    <path d=\"M 408.990358 288.430125 \n",
       "L 415.483449 288.430125 \n",
       "L 415.483449 34.990125 \n",
       "L 408.990358 34.990125 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_46\">\n",
       "    <path d=\"M 415.483449 288.430125 \n",
       "L 421.97654 288.430125 \n",
       "L 421.97654 118.788835 \n",
       "L 415.483449 118.788835 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_47\">\n",
       "    <path d=\"M 421.97654 288.430125 \n",
       "L 428.469631 288.430125 \n",
       "L 428.469631 124.920448 \n",
       "L 421.97654 124.920448 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_48\">\n",
       "    <path d=\"M 428.469631 288.430125 \n",
       "L 434.962722 288.430125 \n",
       "L 434.962722 169.885609 \n",
       "L 428.469631 169.885609 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_49\">\n",
       "    <path d=\"M 434.962722 288.430125 \n",
       "L 441.455813 288.430125 \n",
       "L 441.455813 190.324319 \n",
       "L 434.962722 190.324319 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_50\">\n",
       "    <path d=\"M 441.455813 288.430125 \n",
       "L 447.948903 288.430125 \n",
       "L 447.948903 216.894641 \n",
       "L 441.455813 216.894641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_51\">\n",
       "    <path d=\"M 447.948903 288.430125 \n",
       "L 454.441994 288.430125 \n",
       "L 454.441994 216.894641 \n",
       "L 447.948903 216.894641 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_52\">\n",
       "    <path d=\"M 454.441994 288.430125 \n",
       "L 460.935085 288.430125 \n",
       "L 460.935085 251.640448 \n",
       "L 454.441994 251.640448 \n",
       "z\n",
       "\" clip-path=\"url(#p55ff305299)\" style=\"fill: #1f77b4\"/>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_1\">\n",
       "    <g id=\"xtick_1\">\n",
       "     <g id=\"line2d_1\">\n",
       "      <defs>\n",
       "       <path id=\"m790fd9ce94\" d=\"M 0 0 \n",
       "L 0 3.5 \n",
       "\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"136.28054\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_1\">\n",
       "      <!-- 16 -->\n",
       "      <g transform=\"translate(129.91804 303.028562) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-31\" d=\"M 794 531 \n",
       "L 1825 531 \n",
       "L 1825 4091 \n",
       "L 703 3866 \n",
       "L 703 4441 \n",
       "L 1819 4666 \n",
       "L 2450 4666 \n",
       "L 2450 531 \n",
       "L 3481 531 \n",
       "L 3481 0 \n",
       "L 794 0 \n",
       "L 794 531 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "        <path id=\"DejaVuSans-36\" d=\"M 2113 2584 \n",
       "Q 1688 2584 1439 2293 \n",
       "Q 1191 2003 1191 1497 \n",
       "Q 1191 994 1439 701 \n",
       "Q 1688 409 2113 409 \n",
       "Q 2538 409 2786 701 \n",
       "Q 3034 994 3034 1497 \n",
       "Q 3034 2003 2786 2293 \n",
       "Q 2538 2584 2113 2584 \n",
       "z\n",
       "M 3366 4563 \n",
       "L 3366 3988 \n",
       "Q 3128 4100 2886 4159 \n",
       "Q 2644 4219 2406 4219 \n",
       "Q 1781 4219 1451 3797 \n",
       "Q 1122 3375 1075 2522 \n",
       "Q 1259 2794 1537 2939 \n",
       "Q 1816 3084 2150 3084 \n",
       "Q 2853 3084 3261 2657 \n",
       "Q 3669 2231 3669 1497 \n",
       "Q 3669 778 3244 343 \n",
       "Q 2819 -91 2113 -91 \n",
       "Q 1303 -91 875 529 \n",
       "Q 447 1150 447 2328 \n",
       "Q 447 3434 972 4092 \n",
       "Q 1497 4750 2381 4750 \n",
       "Q 2619 4750 2861 4703 \n",
       "Q 3103 4656 3366 4563 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-31\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-36\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_2\">\n",
       "     <g id=\"line2d_2\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"190.389631\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_2\">\n",
       "      <!-- 18 -->\n",
       "      <g transform=\"translate(184.027131 303.028562) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-38\" d=\"M 2034 2216 \n",
       "Q 1584 2216 1326 1975 \n",
       "Q 1069 1734 1069 1313 \n",
       "Q 1069 891 1326 650 \n",
       "Q 1584 409 2034 409 \n",
       "Q 2484 409 2743 651 \n",
       "Q 3003 894 3003 1313 \n",
       "Q 3003 1734 2745 1975 \n",
       "Q 2488 2216 2034 2216 \n",
       "z\n",
       "M 1403 2484 \n",
       "Q 997 2584 770 2862 \n",
       "Q 544 3141 544 3541 \n",
       "Q 544 4100 942 4425 \n",
       "Q 1341 4750 2034 4750 \n",
       "Q 2731 4750 3128 4425 \n",
       "Q 3525 4100 3525 3541 \n",
       "Q 3525 3141 3298 2862 \n",
       "Q 3072 2584 2669 2484 \n",
       "Q 3125 2378 3379 2068 \n",
       "Q 3634 1759 3634 1313 \n",
       "Q 3634 634 3220 271 \n",
       "Q 2806 -91 2034 -91 \n",
       "Q 1263 -91 848 271 \n",
       "Q 434 634 434 1313 \n",
       "Q 434 1759 690 2068 \n",
       "Q 947 2378 1403 2484 \n",
       "z\n",
       "M 1172 3481 \n",
       "Q 1172 3119 1398 2916 \n",
       "Q 1625 2713 2034 2713 \n",
       "Q 2441 2713 2670 2916 \n",
       "Q 2900 3119 2900 3481 \n",
       "Q 2900 3844 2670 4047 \n",
       "Q 2441 4250 2034 4250 \n",
       "Q 1625 4250 1398 4047 \n",
       "Q 1172 3844 1172 3481 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-31\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-38\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_3\">\n",
       "     <g id=\"line2d_3\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"244.498722\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_3\">\n",
       "      <!-- 20 -->\n",
       "      <g transform=\"translate(238.136222 303.028562) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-32\" d=\"M 1228 531 \n",
       "L 3431 531 \n",
       "L 3431 0 \n",
       "L 469 0 \n",
       "L 469 531 \n",
       "Q 828 903 1448 1529 \n",
       "Q 2069 2156 2228 2338 \n",
       "Q 2531 2678 2651 2914 \n",
       "Q 2772 3150 2772 3378 \n",
       "Q 2772 3750 2511 3984 \n",
       "Q 2250 4219 1831 4219 \n",
       "Q 1534 4219 1204 4116 \n",
       "Q 875 4013 500 3803 \n",
       "L 500 4441 \n",
       "Q 881 4594 1212 4672 \n",
       "Q 1544 4750 1819 4750 \n",
       "Q 2544 4750 2975 4387 \n",
       "Q 3406 4025 3406 3419 \n",
       "Q 3406 3131 3298 2873 \n",
       "Q 3191 2616 2906 2266 \n",
       "Q 2828 2175 2409 1742 \n",
       "Q 1991 1309 1228 531 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "        <path id=\"DejaVuSans-30\" d=\"M 2034 4250 \n",
       "Q 1547 4250 1301 3770 \n",
       "Q 1056 3291 1056 2328 \n",
       "Q 1056 1369 1301 889 \n",
       "Q 1547 409 2034 409 \n",
       "Q 2525 409 2770 889 \n",
       "Q 3016 1369 3016 2328 \n",
       "Q 3016 3291 2770 3770 \n",
       "Q 2525 4250 2034 4250 \n",
       "z\n",
       "M 2034 4750 \n",
       "Q 2819 4750 3233 4129 \n",
       "Q 3647 3509 3647 2328 \n",
       "Q 3647 1150 3233 529 \n",
       "Q 2819 -91 2034 -91 \n",
       "Q 1250 -91 836 529 \n",
       "Q 422 1150 422 2328 \n",
       "Q 422 3509 836 4129 \n",
       "Q 1250 4750 2034 4750 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-30\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_4\">\n",
       "     <g id=\"line2d_4\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"298.607813\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_4\">\n",
       "      <!-- 22 -->\n",
       "      <g transform=\"translate(292.245313 303.028562) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-32\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_5\">\n",
       "     <g id=\"line2d_5\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"352.716903\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_5\">\n",
       "      <!-- 24 -->\n",
       "      <g transform=\"translate(346.354403 303.028562) scale(0.1 -0.1)\">\n",
       "       <defs>\n",
       "        <path id=\"DejaVuSans-34\" d=\"M 2419 4116 \n",
       "L 825 1625 \n",
       "L 2419 1625 \n",
       "L 2419 4116 \n",
       "z\n",
       "M 2253 4666 \n",
       "L 3047 4666 \n",
       "L 3047 1625 \n",
       "L 3713 1625 \n",
       "L 3713 1100 \n",
       "L 3047 1100 \n",
       "L 3047 0 \n",
       "L 2419 0 \n",
       "L 2419 1100 \n",
       "L 313 1100 \n",
       "L 313 1709 \n",
       "L 2253 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       </defs>\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-34\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_6\">\n",
       "     <g id=\"line2d_6\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"406.825994\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_6\">\n",
       "      <!-- 26 -->\n",
       "      <g transform=\"translate(400.463494 303.028562) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-36\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_7\">\n",
       "     <g id=\"line2d_7\">\n",
       "      <g>\n",
       "       <use xlink:href=\"#m790fd9ce94\" x=\"460.935085\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_7\">\n",
       "      <!-- 28 -->\n",
       "      <g transform=\"translate(454.572585 303.028562) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-32\"/>\n",
       "       <use xlink:href=\"#DejaVuSans-38\" x=\"63.623047\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_8\">\n",
       "     <!-- MAG_CSST/mag -->\n",
       "     <g transform=\"translate(241.428563 319.79575) scale(0.144 -0.144)\">\n",
       "      <defs>\n",
       "       <path id=\"DejaVuSans-4d\" d=\"M 628 4666 \n",
       "L 1569 4666 \n",
       "L 2759 1491 \n",
       "L 3956 4666 \n",
       "L 4897 4666 \n",
       "L 4897 0 \n",
       "L 4281 0 \n",
       "L 4281 4097 \n",
       "L 3078 897 \n",
       "L 2444 897 \n",
       "L 1241 4097 \n",
       "L 1241 0 \n",
       "L 628 0 \n",
       "L 628 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-41\" d=\"M 2188 4044 \n",
       "L 1331 1722 \n",
       "L 3047 1722 \n",
       "L 2188 4044 \n",
       "z\n",
       "M 1831 4666 \n",
       "L 2547 4666 \n",
       "L 4325 0 \n",
       "L 3669 0 \n",
       "L 3244 1197 \n",
       "L 1141 1197 \n",
       "L 716 0 \n",
       "L 50 0 \n",
       "L 1831 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-47\" d=\"M 3809 666 \n",
       "L 3809 1919 \n",
       "L 2778 1919 \n",
       "L 2778 2438 \n",
       "L 4434 2438 \n",
       "L 4434 434 \n",
       "Q 4069 175 3628 42 \n",
       "Q 3188 -91 2688 -91 \n",
       "Q 1594 -91 976 548 \n",
       "Q 359 1188 359 2328 \n",
       "Q 359 3472 976 4111 \n",
       "Q 1594 4750 2688 4750 \n",
       "Q 3144 4750 3555 4637 \n",
       "Q 3966 4525 4313 4306 \n",
       "L 4313 3634 \n",
       "Q 3963 3931 3569 4081 \n",
       "Q 3175 4231 2741 4231 \n",
       "Q 1884 4231 1454 3753 \n",
       "Q 1025 3275 1025 2328 \n",
       "Q 1025 1384 1454 906 \n",
       "Q 1884 428 2741 428 \n",
       "Q 3075 428 3337 486 \n",
       "Q 3600 544 3809 666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-5f\" d=\"M 3263 -1063 \n",
       "L 3263 -1509 \n",
       "L -63 -1509 \n",
       "L -63 -1063 \n",
       "L 3263 -1063 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-43\" d=\"M 4122 4306 \n",
       "L 4122 3641 \n",
       "Q 3803 3938 3442 4084 \n",
       "Q 3081 4231 2675 4231 \n",
       "Q 1875 4231 1450 3742 \n",
       "Q 1025 3253 1025 2328 \n",
       "Q 1025 1406 1450 917 \n",
       "Q 1875 428 2675 428 \n",
       "Q 3081 428 3442 575 \n",
       "Q 3803 722 4122 1019 \n",
       "L 4122 359 \n",
       "Q 3791 134 3420 21 \n",
       "Q 3050 -91 2638 -91 \n",
       "Q 1578 -91 968 557 \n",
       "Q 359 1206 359 2328 \n",
       "Q 359 3453 968 4101 \n",
       "Q 1578 4750 2638 4750 \n",
       "Q 3056 4750 3426 4639 \n",
       "Q 3797 4528 4122 4306 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-53\" d=\"M 3425 4513 \n",
       "L 3425 3897 \n",
       "Q 3066 4069 2747 4153 \n",
       "Q 2428 4238 2131 4238 \n",
       "Q 1616 4238 1336 4038 \n",
       "Q 1056 3838 1056 3469 \n",
       "Q 1056 3159 1242 3001 \n",
       "Q 1428 2844 1947 2747 \n",
       "L 2328 2669 \n",
       "Q 3034 2534 3370 2195 \n",
       "Q 3706 1856 3706 1288 \n",
       "Q 3706 609 3251 259 \n",
       "Q 2797 -91 1919 -91 \n",
       "Q 1588 -91 1214 -16 \n",
       "Q 841 59 441 206 \n",
       "L 441 856 \n",
       "Q 825 641 1194 531 \n",
       "Q 1563 422 1919 422 \n",
       "Q 2459 422 2753 634 \n",
       "Q 3047 847 3047 1241 \n",
       "Q 3047 1584 2836 1778 \n",
       "Q 2625 1972 2144 2069 \n",
       "L 1759 2144 \n",
       "Q 1053 2284 737 2584 \n",
       "Q 422 2884 422 3419 \n",
       "Q 422 4038 858 4394 \n",
       "Q 1294 4750 2059 4750 \n",
       "Q 2388 4750 2728 4690 \n",
       "Q 3069 4631 3425 4513 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-54\" d=\"M -19 4666 \n",
       "L 3928 4666 \n",
       "L 3928 4134 \n",
       "L 2272 4134 \n",
       "L 2272 0 \n",
       "L 1638 0 \n",
       "L 1638 4134 \n",
       "L -19 4134 \n",
       "L -19 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-2f\" d=\"M 1625 4666 \n",
       "L 2156 4666 \n",
       "L 531 -594 \n",
       "L 0 -594 \n",
       "L 1625 4666 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-6d\" d=\"M 3328 2828 \n",
       "Q 3544 3216 3844 3400 \n",
       "Q 4144 3584 4550 3584 \n",
       "Q 5097 3584 5394 3201 \n",
       "Q 5691 2819 5691 2113 \n",
       "L 5691 0 \n",
       "L 5113 0 \n",
       "L 5113 2094 \n",
       "Q 5113 2597 4934 2840 \n",
       "Q 4756 3084 4391 3084 \n",
       "Q 3944 3084 3684 2787 \n",
       "Q 3425 2491 3425 1978 \n",
       "L 3425 0 \n",
       "L 2847 0 \n",
       "L 2847 2094 \n",
       "Q 2847 2600 2669 2842 \n",
       "Q 2491 3084 2119 3084 \n",
       "Q 1678 3084 1418 2786 \n",
       "Q 1159 2488 1159 1978 \n",
       "L 1159 0 \n",
       "L 581 0 \n",
       "L 581 3500 \n",
       "L 1159 3500 \n",
       "L 1159 2956 \n",
       "Q 1356 3278 1631 3431 \n",
       "Q 1906 3584 2284 3584 \n",
       "Q 2666 3584 2933 3390 \n",
       "Q 3200 3197 3328 2828 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-61\" d=\"M 2194 1759 \n",
       "Q 1497 1759 1228 1600 \n",
       "Q 959 1441 959 1056 \n",
       "Q 959 750 1161 570 \n",
       "Q 1363 391 1709 391 \n",
       "Q 2188 391 2477 730 \n",
       "Q 2766 1069 2766 1631 \n",
       "L 2766 1759 \n",
       "L 2194 1759 \n",
       "z\n",
       "M 3341 1997 \n",
       "L 3341 0 \n",
       "L 2766 0 \n",
       "L 2766 531 \n",
       "Q 2569 213 2275 61 \n",
       "Q 1981 -91 1556 -91 \n",
       "Q 1019 -91 701 211 \n",
       "Q 384 513 384 1019 \n",
       "Q 384 1609 779 1909 \n",
       "Q 1175 2209 1959 2209 \n",
       "L 2766 2209 \n",
       "L 2766 2266 \n",
       "Q 2766 2663 2505 2880 \n",
       "Q 2244 3097 1772 3097 \n",
       "Q 1472 3097 1187 3025 \n",
       "Q 903 2953 641 2809 \n",
       "L 641 3341 \n",
       "Q 956 3463 1253 3523 \n",
       "Q 1550 3584 1831 3584 \n",
       "Q 2591 3584 2966 3190 \n",
       "Q 3341 2797 3341 1997 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "       <path id=\"DejaVuSans-67\" d=\"M 2906 1791 \n",
       "Q 2906 2416 2648 2759 \n",
       "Q 2391 3103 1925 3103 \n",
       "Q 1463 3103 1205 2759 \n",
       "Q 947 2416 947 1791 \n",
       "Q 947 1169 1205 825 \n",
       "Q 1463 481 1925 481 \n",
       "Q 2391 481 2648 825 \n",
       "Q 2906 1169 2906 1791 \n",
       "z\n",
       "M 3481 434 \n",
       "Q 3481 -459 3084 -895 \n",
       "Q 2688 -1331 1869 -1331 \n",
       "Q 1566 -1331 1297 -1286 \n",
       "Q 1028 -1241 775 -1147 \n",
       "L 775 -588 \n",
       "Q 1028 -725 1275 -790 \n",
       "Q 1522 -856 1778 -856 \n",
       "Q 2344 -856 2625 -561 \n",
       "Q 2906 -266 2906 331 \n",
       "L 2906 616 \n",
       "Q 2728 306 2450 153 \n",
       "Q 2172 0 1784 0 \n",
       "Q 1141 0 747 490 \n",
       "Q 353 981 353 1791 \n",
       "Q 353 2603 747 3093 \n",
       "Q 1141 3584 1784 3584 \n",
       "Q 2172 3584 2450 3431 \n",
       "Q 2728 3278 2906 2969 \n",
       "L 2906 3500 \n",
       "L 3481 3500 \n",
       "L 3481 434 \n",
       "z\n",
       "\" transform=\"scale(0.015625)\"/>\n",
       "      </defs>\n",
       "      <use xlink:href=\"#DejaVuSans-4d\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-41\" x=\"86.279297\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-47\" x=\"152.9375\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-5f\" x=\"230.427734\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-43\" x=\"280.427734\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-53\" x=\"350.251953\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-53\" x=\"413.728516\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-54\" x=\"477.205078\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-2f\" x=\"538.289062\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-6d\" x=\"571.980469\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-61\" x=\"669.392578\"/>\n",
       "      <use xlink:href=\"#DejaVuSans-67\" x=\"730.671875\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_2\">\n",
       "    <g id=\"ytick_1\">\n",
       "     <g id=\"line2d_8\">\n",
       "      <defs>\n",
       "       <path id=\"ma58ea94c7b\" d=\"M 0 0 \n",
       "L -3.5 0 \n",
       "\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use xlink:href=\"#ma58ea94c7b\" x=\"120.047812\" y=\"288.430125\" style=\"stroke: #000000; stroke-width: 0.8\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_9\">\n",
       "      <!-- 0 -->\n",
       "      <g transform=\"translate(106.685312 292.229344) scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-30\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",