<!-- 
RSS generated by JIRA (8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b) at Sat Feb 10 15:28:19 JST 2024

It is possible to restrict the fields that are returned in this document by specifying the 'field' parameter in your request.
For example, to request only the issue key and summary append 'field=key&field=summary' to the URL of your request.
-->
<rss version="0.92" >
<channel>
    <title>PFS-JIRA</title>
    <link>https://pfspipe.ipmu.jp/jira</link>
    <description>This file is an XML representation of an issue</description>
    <language>en-us</language>    <build-info>
        <version>8.3.4</version>
        <build-number>803005</build-number>
        <build-date>13-09-2019</build-date>
    </build-info>


<item>
            <title>[REDMINE1D-5] [RM-6508]  adaptative 2nd pass redshift window size </title>
                <link>https://pfspipe.ipmu.jp/jira/browse/REDMINE1D-5</link>
                <project id="11002" key="REDMINE1D">1D Redmine </project>
                    <description>&lt;p&gt;&lt;em&gt;&lt;font color=&quot;#505f79&quot;&gt; Created on 2021-05-21 08:24:00 by Didier Vibert. % Done: 50&lt;/font&gt;&lt;/em&gt;&lt;/p&gt;


&lt;p&gt;given the behavior of the distance shift from 1st pass candidate to 2nd pass solution (see  #6051, and attached plots: &lt;a href=&quot;https://projets.lam.fr/attachments/7362&quot; class=&quot;external-link&quot; rel=&quot;nofollow&quot;&gt;https://projets.lam.fr/attachments/7362&lt;/a&gt; &amp;amp; &lt;a href=&quot;https://projets.lam.fr/attachments/7363&quot; class=&quot;external-link&quot; rel=&quot;nofollow&quot;&gt;https://projets.lam.fr/attachments/7363&lt;/a&gt;),&lt;/p&gt;

&lt;p&gt;if we want to be conservative to keep all candidates  future solution inside the second pass window, we need a big enough size, which will waste some time for the majority of spectra, thus an alternative is to handle specifically the small number of spectra for which the solutions moves away. &lt;/p&gt;

&lt;p&gt;Morevover , from #6360, at least some of them are explained by  large velocity dispersion leading to a larger pdf peak width.&lt;br/&gt;
Several possibility (not exclusive):&lt;/p&gt;

&lt;ul&gt;
	&lt;li&gt;if no maxima is found in the second pass window (max is at border) enlarge the window size (only on the size where the max is at border ? ) and recompute 2nd pass&lt;/li&gt;
	&lt;li&gt;since 2nd pass window is used after fitting the velocity dispersion, one could redefine the window size given the fitted velocity, for large velocity we could enlarge the window size before even starting to use the 2nd pass window (@COperatorLineModel::RecomputeAroundCandidates@)&lt;/li&gt;
&lt;/ul&gt;
</description>
                <environment></environment>
        <key id="16709">REDMINE1D-5</key>
            <summary>[RM-6508]  adaptative 2nd pass redshift window size </summary>
                <type id="3" iconUrl="https://pfspipe.ipmu.jp/jira/secure/viewavatar?size=xsmall&amp;avatarId=10518&amp;avatarType=issuetype">Task</type>
                                            <priority id="10000" iconUrl="https://pfspipe.ipmu.jp/jira/images/icons/priorities/medium.svg">Normal</priority>
                        <status id="3" iconUrl="https://pfspipe.ipmu.jp/jira/images/icons/statuses/inprogress.png" description="This issue is being actively worked on at the moment by the assignee.">In Progress</status>
                    <statusCategory id="4" key="indeterminate" colorName="yellow"/>
                                    <resolution id="-1">Unresolved</resolution>
                                        <assignee username="r2j.migrate">Redmine-Jira Migtation</assignee>
                                    <reporter username="r2j.migrate">Redmine-Jira Migtation</reporter>
                        <labels>
                    </labels>
                <created>Fri, 4 Jun 2021 01:19:06 +0000</created>
                <updated>Fri, 12 Jan 2024 18:44:24 +0000</updated>
                                                                                <due></due>
                            <votes>0</votes>
                                    <watches>1</watches>
                                                                <comments>
                            <comment id="36624" author="r2j.migrate" created="Fri, 12 Jan 2024 18:38:46 +0000"  >&lt;p&gt;Comment by Mira Sarkis on 2021-07-06 08:21:37:&lt;br/&gt;
Note:&lt;/p&gt;
&lt;ul&gt;
	&lt;li&gt;Given that we precompute the continuum in the secondpass prior to @estimateSecondPass@, and&lt;/li&gt;
	&lt;li&gt;Given that @estimateSecondPass@ uses the @precomputedContinuum@ results in 2 cases: refitfirstpass or retryall.&lt;br/&gt;
if we decide to &quot;correct/weigh&quot; the halfwindowSize using the velocity (i.e., absoption or emission) computed in @estimateSecondPass@, then we need to recompute the continuum for the newly added redshifts. &lt;/li&gt;
&lt;/ul&gt;


&lt;p&gt;As a first version, I call @::precomputecontinuum@ after @estimatesecondpass@, only on candidates with estimated velocities&amp;gt;300. &lt;/p&gt;

&lt;p&gt;Concerning the windowsize correction:&lt;/p&gt;
&lt;ul&gt;
	&lt;li&gt;correction heurstically starts from velocity = 250 by a step of 100. Velocity ranges are:
	&lt;ul&gt;
		&lt;li&gt;&lt;span class=&quot;error&quot;&gt;&amp;#91;20:249&amp;#93;&lt;/span&gt; - no correction (N=0),&lt;/li&gt;
		&lt;li&gt;&lt;span class=&quot;error&quot;&gt;&amp;#91;250:349&amp;#93;&lt;/span&gt; - simple correction (N=1),&lt;/li&gt;
		&lt;li&gt;&lt;span class=&quot;error&quot;&gt;&amp;#91;350:449&amp;#93;&lt;/span&gt; - double correction(N=2),&lt;/li&gt;
		&lt;li&gt;&lt;span class=&quot;error&quot;&gt;&amp;#91;450:500&amp;#93;&lt;/span&gt; - triple correction (N=3)&lt;/li&gt;
	&lt;/ul&gt;
	&lt;/li&gt;
	&lt;li&gt;correction is set heuristically to a multiple (N) of 0.001, as follows
	&lt;ul&gt;
		&lt;li&gt;corrected_wdw = (m_halfwdwSize + N*0.001)*(1+zcand)&lt;/li&gt;
	&lt;/ul&gt;
	&lt;/li&gt;
&lt;/ul&gt;

</comment>
                            <comment id="36625" author="r2j.migrate" created="Fri, 12 Jan 2024 18:38:52 +0000"  >&lt;p&gt;Comment by Vincent Le Brun on 2021-07-06 09:25:34:&lt;br/&gt;
Is it necessary to use the value of the velocity? we could use the width and the location of the peak : if the peak is on the edge of the window, or if its width is more than xx% of the window, we go back to the original size&lt;/p&gt;</comment>
                            <comment id="36626" author="r2j.migrate" created="Fri, 12 Jan 2024 18:39:05 +0000"  >&lt;p&gt;Comment by Mira Sarkis on 2021-07-06 11:37:22:&lt;br/&gt;
Vincent Le Brun wrote in #note-4:&lt;br/&gt;
&amp;gt; Is it necessary to use the value of the velocity? &lt;br/&gt;
To my understanding, we have observed a link between the estimated velocity and the width of the peak thus the idea of &quot;pre-estimating or &quot;customizing&quot; a window size for each candidate. &lt;br/&gt;
&amp;gt; we could use the width and the location of the peak : if the peak is on the edge of the window, or if its width is more than xx% of the window, we go back to the original size&lt;/p&gt;

&lt;p&gt;I&apos;m not sure what do you mean by &quot;original size&quot;?!&lt;br/&gt;
We can imagine fitting around candidates first within the &lt;span class=&quot;error&quot;&gt;&amp;#91;zc-halfwdwsize; zc+halfwdwsize&amp;#93;&lt;/span&gt;, and if, as you say, the peak is on the border or its width is much bigger than the secondpass window (passed in param.json), then we decide to extend the fitting by including redshifts from &lt;span class=&quot;error&quot;&gt;&amp;#91;zc - customHalfwdwsize; zc+customHalfwdwsize&amp;#93;&lt;/span&gt;.&lt;/p&gt;



</comment>
                            <comment id="36627" author="r2j.migrate" created="Fri, 12 Jan 2024 18:39:21 +0000"  >&lt;p&gt;Comment by Vincent Le Brun on 2021-07-09 14:09:23:&lt;br/&gt;
Mira Sarkis wrote in #note-5:&lt;br/&gt;
&amp;gt; Vincent Le Brun wrote in #note-4:&lt;br/&gt;
&amp;gt; &amp;gt; Is it necessary to use the value of the velocity? &lt;br/&gt;
&amp;gt; To my understanding, we have observed a link between the estimated velocity and the width of the peak thus the idea of &quot;pre-estimating or &quot;customizing&quot; a window size for each candidate. &lt;br/&gt;
&amp;gt; &amp;gt; we could use the width and the location of the peak : if the peak is on the edge of the window, or if its width is more than xx% of the window, we go back to the original size&lt;br/&gt;
&amp;gt; &lt;br/&gt;
&amp;gt; I&apos;m not sure what do you mean by &quot;original size&quot;?!&lt;br/&gt;
the original window size that allowed to keep all solutions (0.01?)&lt;br/&gt;
&amp;gt; We can imagine fitting around candidates first within the &lt;span class=&quot;error&quot;&gt;&amp;#91;zc-halfwdwsize; zc+halfwdwsize&amp;#93;&lt;/span&gt;, and if, as you say, the peak is on the border or its width is much bigger than the secondpass window (passed in param.json), then we decide to extend the fitting by including redshifts from &lt;span class=&quot;error&quot;&gt;&amp;#91;zc - customHalfwdwsize; zc+customHalfwdwsize&amp;#93;&lt;/span&gt;.&lt;br/&gt;
let&apos;s try this&lt;/p&gt;
</comment>
                            <comment id="36628" author="r2j.migrate" created="Fri, 12 Jan 2024 18:39:40 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2021-07-12 14:28:40:&lt;br/&gt;
I agree with the Vincent proposition. Instead of building several window sizes: ( m_halfwdwSize + N*0.001)*(1+zcand) ) &lt;/p&gt;

&lt;p&gt;it is simpler to use only two: a fine one (using m_halfwdwSize) for the majority of cases, and then if it fails (large velocity, peak at border) go directly to a large one (as we had before 1e-2)&lt;/p&gt;

&lt;p&gt;(note: this change will render #6052 usefull). &lt;/p&gt;</comment>
                            <comment id="36629" author="r2j.migrate" created="Fri, 12 Jan 2024 18:39:47 +0000"  >&lt;p&gt;Comment by Mira Sarkis on 2021-07-13 10:01:28:&lt;br/&gt;
Didier Vibert wrote in #note-7:&lt;br/&gt;
&amp;gt; I agree with the Vincent proposition. Instead of building several window sizes: ( m_halfwdwSize + N*0.001)*(1+zcand) ) &lt;br/&gt;
&amp;gt; &lt;br/&gt;
&amp;gt; it is simpler to use only two: a fine one (using m_halfwdwSize) for the majority of cases, and then if it fails (large velocity, peak at border) go directly to a large one (as we had before 1e-2)&lt;br/&gt;
&amp;gt; &lt;/p&gt;

&lt;p&gt;So you are saying to not customize the window size based on the estimated velocity but simply enlarging the window size if no peaks appear in the PDF slots?&lt;/p&gt;

&lt;p&gt;I am not sure this is a good idea given that we can identify border-peaks &lt;b&gt;only&lt;/b&gt; after marginalizing the PDF in CPDFz::Compute and calling ::searchMaxPDFcandidates (especially that we dont care about Chi2 peaks) which happens at a relatively advanced stage in CMethodLinemodel. &lt;/p&gt;

</comment>
                            <comment id="36630" author="r2j.migrate" created="Fri, 12 Jan 2024 18:39:56 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2021-07-13 14:46:48:&lt;br/&gt;
Mira Sarkis wrote in #note-8:&lt;br/&gt;
&amp;gt; Didier Vibert wrote in #note-7:&lt;br/&gt;
&amp;gt; &amp;gt; I agree with the Vincent proposition. Instead of building several window sizes: ( m_halfwdwSize + N*0.001)*(1+zcand) ) &lt;br/&gt;
&amp;gt; &amp;gt; &lt;br/&gt;
&amp;gt; &amp;gt; it is simpler to use only two: a fine one (using m_halfwdwSize) for the majority of cases, and then if it fails (large velocity, peak at border) go directly to a large one (as we had before 1e-2)&lt;br/&gt;
&amp;gt; &amp;gt; &lt;br/&gt;
&amp;gt; &lt;br/&gt;
&amp;gt; So you are saying to not customize the window size based on the estimated velocity but simply enlarging the window size if no peaks appear in the PDF slots?&lt;br/&gt;
&amp;gt; &lt;br/&gt;
not exactly: I still propose to customize the window size based on the estimated velocity  , but not to &quot;fine&quot; tune: ie only 2 sizes: small, default one, and if velocity is larger than threshold, big size. there no big gain at increasing by a few steps the size: go directly to the larger size.&lt;/p&gt;

&lt;p&gt;&amp;gt; I am not sure this is a good idea given that we can identify border-peaks &lt;b&gt;only&lt;/b&gt; after marginalizing the PDF in CPDFz::Compute and calling ::searchMaxPDFcandidates (especially that we dont care about Chi2 peaks) which happens at a relatively advanced stage in CMethodLinemodel.&lt;/p&gt;

&lt;p&gt;yes that&apos;s true: that&apos;s why I proposed at some point to split the work in two, knowing that handling the recompute of the second pass after pdf computation needs some more changes. Moreover, if the job  enlarging with velocity is sufficient to remove all the cases where the window is to small, then no need to implement the second step which is harder. &lt;/p&gt;</comment>
                            <comment id="36631" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:05 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2021-07-23 13:48:27:&lt;br/&gt;
je viens de penser, peut &#234;tre aussi une autre piste &#224; creuser: calculer la largeur du pic pour chaque candidat 1st pass (on ne le fait pas actuellement) et, en fonction, adapter la taille de la fen&#234;tre 2nd passe  ? &lt;/p&gt;

&lt;p&gt;Ceci pourrait &#234;tre fait en plus de l&apos;adaptation avec la vitesse.&lt;/p&gt;

&lt;p&gt;ce qui m&apos;y a fait penser, c&apos;est que l&apos;adaptation &#224;  partir de la vitesse ne se fait qu&apos;avec des raies en &#233;mission, mais quid des spectres dont le pic est essentiellement produit par les raies d&apos;absorption ou du continu tout simplement. On sait qu&apos;on peut obtenir des bons z &#224; bas z, sans raies d&#8217;&#233;mission, gr&#226;ce au continu. Les pics sont probablement larges dans ces cas-l&#224; .&lt;/p&gt;</comment>
                            <comment id="36632" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:14 +0000"  >&lt;p&gt;Comment by Vincent Le Brun on 2021-07-23 14:20:59:&lt;br/&gt;
je confirme &lt;/p&gt;</comment>
                            <comment id="36633" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:26 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2021-07-30 09:48:01:&lt;/p&gt;
&lt;h2&gt;&lt;a name=&quot;recherched%27unproxypourdeltazfinal&quot;&gt;&lt;/a&gt;recherche d&apos;un proxy pour delta_z final&lt;/h2&gt;

&lt;p&gt;apr&#232;s quelques investigations par Mira, sur quelques spectres &#224; raie d&#8217;&#233;mission large, il s&apos;av&#232;re que:&lt;/p&gt;

&lt;ul&gt;
	&lt;li&gt;le dz 1st pass n&apos;est pas un bon guess du dz deuxi&#232;me passe (ce qui est attendu), il n&apos;en demeure pas moins qu&apos;il faut le conserver comme minorant en g&#233;n&#233;ral pour prendre en compte les spectres sans raies d&#8217;&#233;mission&lt;/li&gt;
&lt;/ul&gt;


&lt;ul&gt;
	&lt;li&gt;l&apos;estimation de la largeur du pic de la pdf (deltaz/(1+z) &#224; partir de la largeur de la raie sigma_lambda,  en consid&#233;rant que la pdf est l&apos;auto-correlation du profil est fausse (et ne se justifie pas th&#233;oriquement...)&lt;/li&gt;
&lt;/ul&gt;


&lt;div class=&apos;table-wrap&apos;&gt;
&lt;table class=&apos;confluenceTable&apos;&gt;&lt;tbody&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.spectrum    &lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.HalfwdwSize_parameter   &lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.1st pass candidate	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.Redshift values	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.Deltaz	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.sigma_lambda/lambda*(1+z)	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.Deltaz/(1+Z)	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.sigma_lambda/lambda	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;_.sigma/delatZ&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6.spc_Euclid_GC_E2E_200k		  &lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6. 0.002	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP0	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.07145	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000476565	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00679762	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000230063482102	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.003281575707838	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;14.2637835342503&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP1	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.0600863	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000261602	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00347875	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000246774248474	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.003281572453111	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;13.2978723404255&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP2	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.05989	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000426925	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00675966	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00020725621271	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.003281563578638	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;15.833366516367&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP3	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.71516	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000742926	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00890999	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00027362144404	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.003281570883484	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;11.9931056390542&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP4	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.02415	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000513293	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00439567	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000253584467554	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.002171612775733	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;8.56366636599369&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;\9.	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6.spc_SDSS_Highz_QSO	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6. 0.001	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP0	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;3.36839	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00017212	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00540008	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;3.94012439365533E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001236171678811	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;31.3739251684871&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP1	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.85996	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000433905	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00756644	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;6.32518265412626E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001102986023242	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;17.4380106244454&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP2	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.23523	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;9.49E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00276312	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;4.24564809885336E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001236168090085	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;29.1124716710971&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP3	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.79706	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00131643	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00679706	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000193676383613	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.16325212886367&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP4	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.73475	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.0021114	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00673475	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000313508296522	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;3.18970825044994&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;\9.	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6. se8_unit	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6. 0.002	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP0	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.38202	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000720764	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00673397	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000302585200796	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.002826999773302	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;9.34282233851857&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP1	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;3.19272	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00145552	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.0118528	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000347154114751	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.002826995363392	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;8.1433439595471&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP2	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;2.12224	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00089721	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00882655	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000287360997233	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.002826992800041	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;9.83777487990549&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP3	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.32507	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00063692	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00403199	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000273935838491	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001734137036734	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;6.33044966400804&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP4	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;2.13946	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.0011213	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00887523	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000357163333822	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.002826992540118	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;7.91512530098992&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;\9.	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6.	sp8_unit	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;/6. 0.001	&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP0	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;2.20769	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000173049	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00320769	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.39481683080348E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;18.5363105247647&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP1	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;1.38918	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000127554	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00238918	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.33881917645385E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;18.7307336500619&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP2	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.822757	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;9.35E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00182276	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;5.12959215079136E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001000001645858	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;19.5019969978163&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP3	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;2.74923	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000253563	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00374923	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;6.76306868343633E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;14.7861872591821&lt;/td&gt;
&lt;/tr&gt;
&lt;tr&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;FP4	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;3.27571	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.000372548	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.00427571	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;8.71312600714267E-05	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;0.001	&lt;/td&gt;
&lt;td class=&apos;confluenceTd&apos;&gt;11.4769372000387&lt;/td&gt;
&lt;/tr&gt;
&lt;/tbody&gt;&lt;/table&gt;
&lt;/div&gt;


&lt;h2&gt;&lt;a name=&quot;estimationanalytiquededeltaz%3A&quot;&gt;&lt;/a&gt;estimation analytique de delta_z: &lt;/h2&gt;

&lt;p&gt;Finalement pour estimer la largeur du pic de la pdf, il vaut mieux utiliser l&apos;expression analytique donn&#233;e par Hagen et.al., 2007 (attachment:ao-46-22-5374.pdf), dans le cas de  fit d&apos;une Gaussienne 1D, en fonction des param&#232;tres estim&#233;s de la Gaussienne :&lt;/p&gt;

&lt;p&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/attach/noimage.png&quot; imagetext=&quot;Capture%20d%E2%80%99%C3%A9cran%20de%202021-07-27%2010-59-43.png&quot; align=&quot;absmiddle&quot; border=&quot;0&quot; /&gt; (eq 16)&lt;/p&gt;

&lt;p&gt;o&#249;&lt;/p&gt;
&lt;ul&gt;
	&lt;li&gt;x est la position estim&#233;e (lambda pour nous, et donc x = lambda_rest *(z+1))&lt;/li&gt;
	&lt;li&gt;var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt; est la variance sur la position de la Gaussienne fitt&#233;e (x=lambda) et donc, delta_z =  sqrt( var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt;)/lambda_rest = sqrt( var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt;) * (1+z) / x&lt;/li&gt;
	&lt;li&gt;sigma est l&apos;&#233;cart-type du bruit par pixel (m&#234;me unit&#233; que le spectre)&lt;/li&gt;
	&lt;li&gt;delta_x la taille du pixel (m&#234;me unit&#233; que x, donc en Angtrom)&lt;/li&gt;
	&lt;li&gt;w est la largeur de la fit&#233;e de la  gaussienne (donc la largeur de la raie en incluant LSF et dispersion de vitesse)&lt;/li&gt;
	&lt;li&gt;A est l&apos;amplitude fit&#233;e de la Gaussienne (m&#234;me unit&#233; que le spectre)&lt;/li&gt;
&lt;/ul&gt;


&lt;p&gt;Pour utiliser cette expression il faut disposer de tous les param&#232;tres fitt&#233;s et principalement la largeur. La boucle de fit de la vitesse se fait avant le calcul de la pdf sur la fen&#234;tre seconde passe (sur une petite plage z autour de la solution 1stpass et en utilisant les r&#233;sultats du fit 1stpass), on peut donc r&#233;cup&#233;rer la valeur de l&apos;amplitude fit&#233;e de la raie la plus intense correspondant, au meilleur fit de la vitesse, ainsi que le niveau de bruit &#224; la longueur d&apos;onde de la raie (il faudrait r&#233;cup&#233;rer ces valeurs &#224; cet endroit dans le code: &lt;a href=&quot;https://gitlab.lam.fr/CPF/cpf-redshift/-/blob/release-0.22/RedshiftLibrary/src/lib/operator/linemodel.cpp#L1806&quot; class=&quot;external-link&quot; rel=&quot;nofollow&quot;&gt;https://gitlab.lam.fr/CPF/cpf-redshift/-/blob/release-0.22/RedshiftLibrary/src/lib/operator/linemodel.cpp#L1806&lt;/a&gt;) pour estimer delta_z et modifier le param&#232;tre HalfwindowSize.&lt;/p&gt;

&lt;p&gt;Une autre solution &#224; essayer, pour ne tenir compte que de la largeur est d&apos;&#233;liminer l&apos;amplitude dans l&apos;expression de var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt; et de le remplacer par un SNR minimal (eg 3.5) de la raie afin d&apos;obtenir un majorant de delta_z. &lt;/p&gt;

&lt;p&gt;Pour cel&#224;, en utilisant l&apos;equation (17) du flux U:&lt;br/&gt;
&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/attach/noimage.png&quot; imagetext=&quot;Capture%20d%E2%80%99%C3%A9cran%20de%202021-07-30%2011-26-32.png&quot; align=&quot;absmiddle&quot; border=&quot;0&quot; /&gt;&lt;/p&gt;

&lt;p&gt;et l&apos;expression de la variance ,var(U), du FLUX estim&#233; (equation 19):&lt;br/&gt;
&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/attach/noimage.png&quot; imagetext=&quot;Capture%20d%E2%80%99%C3%A9cran%20de%202021-07-30%2011-28-20.png&quot; align=&quot;absmiddle&quot; border=&quot;0&quot; /&gt; &lt;img src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/attach/noimage.png&quot; imagetext=&quot;Capture%20d%E2%80%99%C3%A9cran%20de%202021-07-30%2011-28-44.png&quot; align=&quot;absmiddle&quot; border=&quot;0&quot; /&gt;&lt;/p&gt;

&lt;p&gt;et en posant SNR=U/sqrt(var(U)),&lt;/p&gt;

&lt;p&gt;on obtient (calcul &#224; v&#233;rifier...):&lt;/p&gt;

&lt;p&gt;&lt;del&gt;delta_z = (1+z)*sqrt(var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt;)/x =(1+z)/x*sqrt(4/3*pi) sqrt(w)/SNR&lt;/del&gt;&lt;/p&gt;

&lt;p&gt;&lt;b&gt;apr&#232;s correction:  delta_z = (1+z)*sqrt(var&lt;img class=&quot;emoticon&quot; src=&quot;https://pfspipe.ipmu.jp/jira/images/icons/emoticons/error.png&quot; height=&quot;16&quot; width=&quot;16&quot; align=&quot;absmiddle&quot; alt=&quot;&quot; border=&quot;0&quot;/&gt;)/x =(1+z)/x*sqrt(4/3) w/SNR&lt;/b&gt; (eq. 1)&lt;/p&gt;

&lt;h2&gt;&lt;a name=&quot;TODO&quot;&gt;&lt;/a&gt;TODO&lt;/h2&gt;

&lt;ul&gt;
	&lt;li&gt;v&#233;rifier le calcul de l&apos;expression de delta_z en fonction de (w, SNR)  &lt;b&gt;DONE&lt;/b&gt;&lt;/li&gt;
	&lt;li&gt;faire un run sur un dataset statistiquement repr&#233;sentatif, eg Euclid_GC_E2E_200k, en prenant une taille de fen&#234;tre 2nd pass conservative (eg la valeur qu&apos;on utilisait avant: 1e-2 par exemple).&lt;br/&gt;
Avec les sorties de ce run, r&#233;cup&#233;rer pour chaque spectre, l&apos;amplitude, la largeur et le niveau de bruit sous la raie la plus intense et comparer avec l&apos;expression (delta_x~13 Angstrom pour Euclid), comparer aussi avec l&apos;expression en fonction du SNR en v&#233;rifiant qu&apos;avec SNR=3.5 on a bien une borne sup et que celle-ci est l&#224; la limite des plus grands delta_z (ie qu&apos;elle ne soit pas trop grande).&lt;/li&gt;
&lt;/ul&gt;



</comment>
                            <comment id="36634" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:32 +0000"  >&lt;p&gt;Comment by Mira Sarkis on 2021-10-04 13:03:59:&lt;br/&gt;
Quelques r&#233;sultats:&lt;br/&gt;
En se basant sur les r&#233;sultats de 5000 spectres Euclid_GC_E2E_200k (cosmos) dont les valeurs de redshift couvre la globalit&#233; du range &lt;span class=&quot;error&quot;&gt;&amp;#91;0:5&amp;#93;&lt;/span&gt;. &lt;br/&gt;
Pour chaque spectre et pour chacun de ces 5 candidats: &lt;br/&gt;
1) je calcule les valeurs de analyticalDeltaz selon les deux formules analytiques: eq. 16 (Meth1)  et eq. 1 (Meth2)  de la note &lt;a href=&quot;#note-14]&quot;&gt;[https://projets.lam.fr/issues/6508#note-14]&lt;/a&gt;&lt;br/&gt;
2)  je compare ces valeurs entre elles memes et puis &#224; estimatedDeltaz (ou RedshiftUncertainty) obtenu  par amazed.&lt;/p&gt;

&lt;ul&gt;
	&lt;li&gt;Le ratio entre analyticalDeltaz_Meth2/analyticalDeltaz_Meth1 varie entre 2.13 et 4.2, avec &lt;b&gt;une moyenne de ratio = 2.5&lt;/b&gt;  (70% des donn&#233;espr&#233;sente un ratio &amp;lt; 2.5).&lt;/li&gt;
&lt;/ul&gt;


&lt;ul&gt;
	&lt;li&gt;estimatedDeltaz varie principalement entre &lt;span class=&quot;error&quot;&gt;&amp;#91;10-3; 0.003&amp;#93;&lt;/span&gt; comme le montre la distribution des valeurs dans le plot suivant : &lt;span class=&quot;image-wrap&quot; style=&quot;&quot;&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/secure/attachment/17025/17025_EstimatedDeltazDistribution.png&quot; style=&quot;border: 0px solid black&quot; /&gt;&lt;/span&gt;&lt;br/&gt;
Les plots suivants montrent la distribution des differences entre estimatedDeltaz et les analyticalDeltaz:&lt;br/&gt;
Le premier plot correspond &#224; la distribution cumul&#233;e et dont l&apos;axe X est en log10. 50% de la population pr&#233;sente une diff&#233;rence aux alentours de 10^-3, ie., de meme ordre que estimatedDeltaz .&lt;br/&gt;
&lt;span class=&quot;image-wrap&quot; style=&quot;&quot;&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/secure/attachment/17024/17024_cumultaiveErrDist_Meth1_Meth2_Deltaz_AbsDifference.png&quot; style=&quot;border: 0px solid black&quot; /&gt;&lt;/span&gt; &lt;span class=&quot;image-wrap&quot; style=&quot;&quot;&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/secure/attachment/17026/17026_noncumulativeDist_Meth1_Meth2_deltaz.png&quot; style=&quot;border: 0px solid black&quot; /&gt;&lt;/span&gt;&lt;/li&gt;
&lt;/ul&gt;
</comment>
                            <comment id="36635" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:40 +0000"  >&lt;p&gt;Comment by Vincent Le Brun on 2021-10-04 16:40:55:&lt;br/&gt;
du coup la il manque la comparaison entre les m&#233;thodes analytiques et l&apos;estimation par Amazed ?&lt;/p&gt;</comment>
                            <comment id="36636" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:47 +0000"  >&lt;p&gt;Comment by Mira Sarkis on 2021-10-08 08:19:25:&lt;/p&gt;

&lt;p&gt;Vincent Le Brun wrote in #note-17:&lt;br/&gt;
&amp;gt; du coup la il manque la comparaison entre les m&#233;thodes analytiques et l&apos;estimation par Amazed ?&lt;/p&gt;

&lt;p&gt;Dans la #note-16, j&apos;essayais de dire que la &lt;b&gt;diff&#233;rence&lt;/b&gt; entre l&apos;estimation et les m&#233;thodes analytiques * est de l&apos;ordre de la valeur de deltaz estim&#233;*.  Mais il nous semble possible d&apos;utiliser le calcul analytique comme un majorant pour deltaz estim&#233;: on a une premi&#232;re conclusion encore &#224; v&#233;rifier que deltazEstim&#233; = R * computedDeltaz avec R = 10 pour eq et R = 4 pour eq.&lt;/p&gt;

&lt;p&gt;Les plots suivants correspondent &#224; des histogramme 2D de la variation des valeurs de deltaz estim&#233; (yAxis) par amazed en fonction des valeurs calcul&#233;es (xAxis) par les deux m&#233;thodes. &lt;br/&gt;
&lt;span class=&quot;image-wrap&quot; style=&quot;&quot;&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/secure/attachment/17028/17028_2DHist_Hagen2007.png&quot; style=&quot;border: 0px solid black&quot; /&gt;&lt;/span&gt;&lt;br/&gt;
&lt;span class=&quot;image-wrap&quot; style=&quot;&quot;&gt;&lt;img src=&quot;https://pfspipe.ipmu.jp/jira/secure/attachment/17027/17027_2DHist_SNRformula.png&quot; style=&quot;border: 0px solid black&quot; /&gt;&lt;/span&gt;&lt;/p&gt;

&lt;p&gt;TODO:&lt;br/&gt;
Didier propose de refaire ce meme travail en utilisant les valeurs provenant de la firstpass (&lt;a href=&quot;https://projets.lam.fr/issues/6683&quot; class=&quot;external-link&quot; rel=&quot;nofollow&quot;&gt;https://projets.lam.fr/issues/6683&lt;/a&gt;) plutot que ceux de la secondpass, sur les validations-tests pfs7 et Euclid noiseless, avec secondwdwsize large de = 5E-2&lt;/p&gt;</comment>
                            <comment id="36637" author="r2j.migrate" created="Fri, 12 Jan 2024 18:40:56 +0000"  >&lt;p&gt;Comment by Vincent Le Brun on 2023-04-26 12:13:21:&lt;br/&gt;
following the discussion on Mattermost, and given that for Euclid the LSF can be as high as 150A, which gives a relative redshift interval of 1e-2, we should maybe include the LSF in the calculation of the width. The method is to be determined&lt;/p&gt;</comment>
                            <comment id="36638" author="r2j.migrate" created="Fri, 12 Jan 2024 18:41:13 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2023-05-25 15:16:46:&lt;br/&gt;
Finalement, il s&apos;av&#232;re difficile de trouver un bon proxy du delta_z (et donc de la largeur minimale de la fen&#234;tre 2nd passe) &#224; partir des donn&#233;es first-pass ou dispersion de vitesse.&lt;/p&gt;

&lt;p&gt;Du coup, apr&#232;s le merge de l&apos;issue #7529, l&apos;id&#233;e est d&apos;aller plus loin que le warning et de reprendre le calcul de la pdf avec une fen&#234;tre &#233;largie. Je rebaisse le %r&#233;alis&#233; pour correspondre &#224; l&apos;ajout de cette impl&#233;mentation.&lt;/p&gt;</comment>
                            <comment id="36639" author="r2j.migrate" created="Fri, 12 Jan 2024 18:41:19 +0000"  >&lt;p&gt;Comment by Didier Vibert on 2024-01-12 16:47:57:&lt;br/&gt;
should wait multiobs #6028&lt;/p&gt;</comment>
                    </comments>
                    <attachments>
                            <attachment id="17028" name="2DHist_Hagen2007.png" size="57798" author="r2j.migrate" created="Fri, 12 Jan 2024 18:44:08 +0000"/>
                            <attachment id="17027" name="2DHist_SNRformula.png" size="54232" author="r2j.migrate" created="Fri, 12 Jan 2024 18:43:57 +0000"/>
                            <attachment id="17019" name="ao-46-22-5374.pdf" size="701296" author="r2j.migrate" created="Fri, 12 Jan 2024 18:41:35 +0000"/>
                            <attachment id="17020" name="Capture d&#8217;&#233;cran de 2021-07-27 10-59-43.png" size="9071" author="r2j.migrate" created="Fri, 12 Jan 2024 18:42:01 +0000"/>
                            <attachment id="17021" name="Capture d&#8217;&#233;cran de 2021-07-30 11-26-32.png" size="7688" author="r2j.migrate" created="Fri, 12 Jan 2024 18:42:38 +0000"/>
                            <attachment id="17022" name="Capture d&#8217;&#233;cran de 2021-07-30 11-28-20.png" size="2312" author="r2j.migrate" created="Fri, 12 Jan 2024 18:42:53 +0000"/>
                            <attachment id="17023" name="Capture d&#8217;&#233;cran de 2021-07-30 11-28-44.png" size="3380" author="r2j.migrate" created="Fri, 12 Jan 2024 18:43:09 +0000"/>
                            <attachment id="17024" name="cumultaiveErrDist_Meth1_Meth2_Deltaz_AbsDifference.png" size="28966" author="r2j.migrate" created="Fri, 12 Jan 2024 18:43:19 +0000"/>
                            <attachment id="17025" name="EstimatedDeltazDistribution.png" size="36452" author="r2j.migrate" created="Fri, 12 Jan 2024 18:43:34 +0000"/>
                            <attachment id="17026" name="noncumulativeDist_Meth1_Meth2_deltaz.png" size="43216" author="r2j.migrate" created="Fri, 12 Jan 2024 18:43:46 +0000"/>
                    </attachments>
                <subtasks>
                    </subtasks>
                <customfields>
                                                <customfield id="customfield_10500" key="com.atlassian.jira.plugins.jira-development-integration-plugin:devsummary">
                        <customfieldname>Development</customfieldname>
                        <customfieldvalues>
                            
                        </customfieldvalues>
                    </customfield>
                                                                                                                                                                                                            <customfield id="customfield_10010" key="com.pyxis.greenhopper.jira:gh-lexo-rank">
                        <customfieldname>Rank</customfieldname>
                        <customfieldvalues>
                            <customfieldvalue>0|zzs780:</customfieldvalue>

                        </customfieldvalues>
                    </customfield>
                                                                                                                                                                                                                                                                    </customfields>
    </item>
</channel>
</rss>