| | 1 | = 2013-02-27 = |
| | 2 | |
| | 3 | Example: |
| | 4 | Original: |
| | 5 | {{{ |
| | 6 | |
| | 7 | float size = (axes.major > source->moments->Mrf) ? axes.major : source->moments->Mrf; |
| | 8 | PAR[PM_PAR_LENGTH] = 1.5*2.35*size; // a tophat of length L has L = 1.5 * 2.35 * sigma |
| | 9 | |
| | 10 | #X Y ID Mxx Mxx Myy Mr1 EXT_F EXT_L EXT_W EXT_Q |
| | 11 | 77.558716 2633.799072 -98 : 11.551100 1.738350 10.857200 6.066510 : 16.974348 32.950000 1.984060 0.725825 |
| | 12 | #IT F L W Q oChi nChi delt dlin rho lambda |
| | 13 | 0 62.666386 21.384449 1.984059 0.686898 : 5817.757812 3362.864990 nan 1238.471436 1.982196 0.001000 |
| | 14 | 1 16.869297 22.949951 1.984059 0.686511 : 3362.864990 3327.216797 0.849738 22.782276 1.564734 0.000250 |
| | 15 | 2 16.829279 24.949951 1.984059 0.689350 : 3327.216797 3292.971436 0.012339 20.549458 1.666485 0.000063 |
| | 16 | 3 16.904570 26.949951 1.984059 0.697310 : 3292.971436 3260.270020 0.011854 18.246912 1.792162 0.000016 |
| | 17 | 4 16.923841 28.949951 1.984059 0.710121 : 3260.270020 3231.121338 0.011319 16.398548 1.777516 0.000004 |
| | 18 | 5 16.916101 30.949951 1.984059 0.721472 : 3231.121338 3207.896484 0.010090 12.160172 1.909912 0.000001 |
| | 19 | XX 16.974512 32.949951 1.984059 0.725825 : 3207.896484 |
| | 20 | }}} |
| | 21 | |
| | 22 | Sqrt: |
| | 23 | {{{ |
| | 24 | float size = (axes.major > sqrt(source->moments->Mrf)) ? axes.major : sqrt(source->moments->Mrf); |
| | 25 | PAR[PM_PAR_LENGTH] = 1.5*2.35*size; // a tophat of length L has L = 1.5 * 2.35 * sigma |
| | 26 | |
| | 27 | #X Y ID Mxx Mxx Myy Mr1 EXT_F EXT_L EXT_W EXT_Q |
| | 28 | 77.559525 2633.799561 -103 : 11.602900 1.715950 10.896300 6.099250 : 17.214936 29.877100 1.992040 0.654534 |
| | 29 | #IT F L W Q oChi nChi delt dlin rho lambda |
| | 30 | 0 62.998722 12.710316 1.992041 0.683864 : 4785.610352 3625.804443 nan 586.573059 1.977257 0.001000 |
| | 31 | 1 19.728992 13.877067 1.992041 0.690744 : 3625.804443 3593.194092 0.398012 14.094786 2.313647 0.000250 |
| | 32 | 2 19.357067 15.877067 1.992041 0.688212 : 3593.194092 3559.262207 0.011191 14.204083 2.388882 0.000063 |
| | 33 | 3 18.790325 17.877068 1.992041 0.654087 : 3559.262207 3522.943359 0.011644 17.820257 2.038065 0.000016 |
| | 34 | 4 18.215334 19.877068 1.992041 0.623698 : 3522.943359 3488.975830 0.012464 18.180664 1.868333 0.000004 |
| | 35 | 5 18.069387 21.877068 1.992041 0.611393 : 3488.975830 3458.826904 0.011657 14.959837 2.015325 0.000001 |
| | 36 | 6 18.013340 23.877068 1.992041 0.616009 : 3458.826904 3429.208984 0.010346 14.725613 2.011320 0.000000 |
| | 37 | 7 17.703453 25.877068 1.992041 0.629182 : 3429.208984 3399.503418 0.010164 16.940977 1.753474 0.000000 |
| | 38 | 8 17.366875 27.877068 1.992041 0.643639 : 3399.503418 3371.803467 0.010194 17.004745 1.628954 0.000000 |
| | 39 | XX 17.215097 29.877068 1.992041 0.654534 : 3371.803467 |
| | 40 | }}} |
| | 41 | |
| | 42 | Adjustment: |
| | 43 | {{{ |
| | 44 | float size = (axes.major > sqrt(source->moments->Mrf)) ? axes.major : sqrt(source->moments->Mrf); |
| | 45 | if (size > psfAxes.major) { size -= psfAxes.major; } |
| | 46 | else { size = psfAxes.major; } |
| | 47 | PAR[PM_PAR_LENGTH] = 1.5*2.35*size; // a tophat of length L has L = 1.5 * 2.35 * sigma |
| | 48 | #X Y ID Mxx Mxx Myy Mr1 EXT_F EXT_L EXT_W EXT_Q |
| | 49 | 77.559189 2633.799561 -98 : 11.558200 1.736740 10.870700 6.072910 : 5.257746 7.728160 1.978440 0.787694 |
| | 50 | #IT F L W Q oChi nChi delt dlin rho lambda |
| | 51 | 0 62.719723 5.728158 1.978445 0.687694 : 3845.604736 3837.298584 nan 124.292549 0.066827 0.001000 |
| | 52 | XX 5.257778 7.728158 1.978445 0.787694 : 3837.298584 |
| | 53 | }}} |
