Changeset 3331 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.cc

Timestamp:

Oct 2, 2007, 3:11:22 PM (18 years ago)

Author:

cmv

Message:

mise en place evolution seulement sur distance au plan Z , cmv 02/10/2007

File:

• trunk/Cosmo/SimLSS/genefluct3d.cc (modified) (6 diffs)

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -744 +744 @@
 
 //-------------------------------------------------------------------
-void GeneFluct3D::ApplyGrowthFactor(void)
+void GeneFluct3D::ApplyGrowthFactor(int type_evol)
 // Apply Growth to real space
 // Using the correspondance between redshift and los comoving distance
 // describe in vector "zred_" "loscom_"
-{
- if(lp_>0) cout<<"--- ApplyGrowthFactor ---"<<endl;
+// type_evol = 1 : evolution de la puissance du bruit avec la distance a l'observateur
+//             2 : evolution de la puissance du bruit avec la distance du plan Z
+//             (tous les pixels d'un plan Z sont mis au meme redshift z que celui du milieu)
+{
+ if(lp_>0) cout<<"--- ApplyGrowthFactor: evol="<<type_evol<<endl;
  check_array_alloc();
 
@@ -756 +759 @@
    cout<<bla<<endl; throw ParmError(bla);
  }
+ if(type_evol<1 || type_evol>2) {
+   char *bla = "GeneFluct3D::ApplyGrowthFactor_Error: bad type_evol value";
+   cout<<bla<<endl; throw ParmError(bla);
+ }
 
  InterpFunc interpinv(loscom2zred_min_,loscom2zred_max_,loscom2zred_);
@@ -762 +769 @@
  //CHECK: Histo hgr(0.9*zred_[0],1.1*zred_[n-1],1000);
  for(long i=0;i<Nx_;i++) {
-   double dx2 = xobs_[0] - i*Dx_; dx2 *= dx2;
+   double dx2 = DXcom(i); dx2 *= dx2;
    for(long j=0;j<Ny_;j++) {
-     double dy2 = xobs_[1] - j*Dy_; dy2 *= dy2;
+     double dy2 = DYcom(j); dy2 *= dy2;
      for(long l=0;l<Nz_;l++) {
-       double dz2 = xobs_[2] - l*Dz_; dz2 *= dz2;
-       dz2 = sqrt(dx2+dy2+dz2);
-       double z = interpinv(dz2);
+       double dz = DZcom(l);
+       if(type_evol==1) dz = sqrt(dx2+dy2+dz*dz);
+         else dz = fabs(dz); // tous les plans Z au meme redshift
+       double z = interpinv(dz);
        //CHECK: hgr.Add(z);
        double dzgr = (*growth_)(z);   // interpolation par morceau
@@ -1416 +1424 @@
 }
 
-void GeneFluct3D::AddNoise2Real(double snoise,bool with_evol)
+void GeneFluct3D::AddNoise2Real(double snoise,int type_evol)
 // add noise to every pixels (meme les <=0 !)
-{
- if(lp_>0) cout<<"--- AddNoise2Real: snoise = "<<snoise<<" evol="<<with_evol<<endl;
- check_array_alloc();
+// type_evol = 0 : pas d'evolution de la puissance du bruit
+//             1 : evolution de la puissance du bruit avec la distance a l'observateur
+//             2 : evolution de la puissance du bruit avec la distance du plan Z
+//                 (tous les plans Z sont mis au meme redshift z de leur milieu)
+{
+ if(lp_>0) cout<<"--- AddNoise2Real: snoise = "<<snoise<<" evol="<<type_evol<<endl;
+ check_array_alloc();
+
+ if(type_evol<0) type_evol = 0;
+ if(type_evol>2) {
+   char *bla = "GeneFluct3D::AddNoise2Real_Error: bad type_evol value";
+   cout<<bla<<endl; throw ParmError(bla);
+ }
 
  vector<double> correction;
  InterpFunc *intercor = NULL;
 
- if(with_evol) {
+ if(type_evol>0) {
    // Sigma_Noise(en mass) :
    //      Slim ~ 1/sqrt(DNu) * sqrt(nlobe)   en W/m^2Hz
@@ -1454 +1472 @@
  } 
 
- double dx2=0., dy2=0., dz2=0.;
+ double corrlim[2] = {1.,1.};
  for(long i=0;i<Nx_;i++) {
-   dx2 = DXcom(i); dx2 *= dx2;
+   double dx2 = DXcom(i); dx2 *= dx2;
    for(long j=0;j<Ny_;j++) {
-     dy2 = DYcom(j); dy2 *= dy2;
+     double dy2 = DYcom(j); dy2 *= dy2;
      for(long l=0;l<Nz_;l++) {
        double corr = 1.;
-       if(with_evol) {
-         dz2 = DZcom(l); dz2 *= dz2; dz2 = sqrt(dx2+dy2+dz2);
-         corr = (*intercor)(dz2);
+       if(type_evol>0) {
+         double dz = DZcom(l);
+         if(type_evol==1) dz = sqrt(dx2+dy2+dz*dz);
+           else dz = fabs(dz); // tous les plans Z au meme redshift
+         corr = (*intercor)(dz);
+         if(corr<corrlim[0]) corrlim[0]=corr; else if(corr>corrlim[1]) corrlim[1]=corr;
        }
        int_8 ip = IndexR(i,j,l);
@@ -1470 +1491 @@
    }
  }
+ if(type_evol>0)
+   cout<<"correction factor range: ["<<corrlim[0]<<","<<corrlim[1]<<"]"<<endl;
 
  if(intercor!=NULL) delete intercor;