Changeset 3154 in Sophya for trunk/Cosmo/SimLSS/cmvconcherr.cc

Timestamp:

Jan 19, 2007, 5:19:49 PM (19 years ago)

Author:

cmv

Message:

modifs d'options de cmvobserv3d.cc
creation de init_fftw() et deplacement a l'initialisation
creation de SetObservator pour positionner l'observateur + ecriture FITS et PPF

(mais le codage des distance/redshift reste a coder)

cmv 19/01/2007

File:

• trunk/Cosmo/SimLSS/cmvconcherr.cc (modified) (3 diffs)

trunk/Cosmo/SimLSS/cmvconcherr.cc

@@ -119 +119 @@
    nread++;
    for(int i=0;i<nherr;i++) herrconc->AddBin(i,herr(i));
-   sum+=herr(itest); sum2+=herr(itest)*herr(itest), nsum++;
+   sum+=herr(itest); sum2+=herr(itest)*herr(itest); nsum++;
  }
  herrconc->ToVariance();
@@ -182 +182 @@
    nread++;
    for(int i=0;i<nherrx;i++)
-     for(int j=0;j<nherrx;j++) herrconc->AddBin(i,j,herr(i,j));
-   sum+=herr(itestx,itesty); sum2+=herr(itestx,itesty)*herr(itestx,itesty), nsum++;
+     for(int j=0;j<nherry;j++) herrconc->AddBin(i,j,herr(i,j));
+   sum+=herr(itestx,itesty); sum2+=herr(itestx,itesty)*herr(itestx,itesty); nsum++;
  }
  herrconc->ToVariance();
+ if(nsum>0) {sum/=nsum; sum2/=nsum; sum2-=sum*sum;}
  cout<<"Test bin "<<itestx<<","<<itesty
      <<" mean="<<sum<<" sigma^2="<<sum2<<"  nsum="<<nsum<<endl;
@@ -214 +215 @@
 openppf cmvconcherr.ppf
 
+zone 2 2
 disp herrconc "hbincont err"
 disp herrconc "hbinerr"
 disp herrconc "hbinent"
 
+zone
 n/plot herrconc.val%log10(x) x>0 ! "connectpoints"
-n/plot herrconc.sqrt(err2)%log10(x) x>0&&err2>0. ! "connectpoints same red"
-
-n/plot herrconc.sqrt(err2)/val%log10(x) x>0&&err2>0.&&val>0 ! "connectpoints"
+n/plot herrconc.sqrt(err2)%log10(x) x>0&&err2>0 ! "connectpoints same red"
+
+n/plot herrconc.sqrt(err2)/val%log10(x) x>0&&err2>0&&val>0 ! "connectpoints"
 
 #### Histo2DErr 2D
 openppf cmvconcherr2.ppf
 
-disp herrconc2 "hbincont"
-disp herrconc2 "hbinerr"
-disp herrconc2 "hbinent"
+zone 2 2
+imag herrconc2 "hbincont"
+imag herrconc2 "hbinerr"
+imag herrconc2 "hbinent"
 */