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

Timestamp:

Sep 27, 2007, 6:21:55 PM (18 years ago)

Author:

cmv

Message:

possibilite de ne pas faire poisson sur Ngal cmv 27/09/2007

File:

• trunk/Cosmo/SimLSS/genefluct3d.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -1054 +1054 @@
 }
 
+double  GeneFluct3D::TurnMass2HIMass(double m_by_mpc3)
+{
+ if(lp_>0) cout<<"--- TurnMass2HIMass : "<<m_by_mpc3<<" Msol/Mpc^3"<<endl;
+
+ double mall=0., mgood=0.;
+ int_8  nall=0,  ngood=0;
+ for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) for(long l=0;l<Nz_;l++) {
+   int_8 ip = IndexR(i,j,l);
+   mall += data_[ip]; nall++;
+   if(data_[ip]>0.) {mgood += data_[ip]; ngood++;}
+ }
+ if(ngood>0) mgood /= (double)ngood;
+ if(nall>0)  mall  /= (double)nall;
+ if(lp_>0) cout<<"...ngood="<<ngood<<" mgood="<<mgood
+               <<",  nall="<<nall<<" mall="<<mall<<endl;
+ if(ngood<=0 || mall<=0.) {
+   cout<<"TurnMass2HIMass_Error: ngood="<<ngood<<" <=0 || mall="<<mall<<" <=0"<<endl;
+   throw RangeCheckError("TurnMass2HIMass_Error: ngood<=0 || mall<=0");
+ }
+
+ // On doit mettre m*Vol masse de HI dans notre survey
+ // On en met uniquement dans les pixels de masse >0.
+ // On MET a zero les pixels <0
+ // On renormalise sur les pixels>0 pour qu'on ait m*Vol masse de HI
+ //   comme on ne prend que les pixels >0, on doit normaliser
+ //   a la moyenne de <1+d_rho/rho> sur ces pixels
+ //   (rappel sur tout les pixels <1+d_rho/rho>=1)
+ // masse de HI a mettre ds 1 px:
+ double dm = m_by_mpc3*Vol_/ (mgood/mall) /(double)ngood;
+ if(lp_>0) cout<<"...HI mass density move from "
+               <<m_by_mpc3*Vol_/double(NRtot_)<<" to "<<dm<<" / pixel"<<endl;
+
+ double sum = 0.;
+ for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) for(long l=0;l<Nz_;l++) {
+   int_8 ip = IndexR(i,j,l);
+   if(data_[ip]<=0.) data_[ip] = 0.;
+   else {
+     data_[ip] *= dm; // par coherence on multiplie aussi les <=0
+     sum += data_[ip];  // mais on ne les compte pas
+   }
+ }
+
+ if(lp_>0) cout<<sum<<"...mass HI put into survey / "<<m_by_mpc3*Vol_<<endl;
+
+ return sum;
+}
+
 double GeneFluct3D::TurnMass2MeanNumber(double n_by_mpc3)
 // do NOT treate negative or nul values
 {
- if(lp_>0) cout<<"--- TurnMass2MeanNumber ---"<<endl;
+ if(lp_>0) cout<<"--- TurnMass2MeanNumber : "<<n_by_mpc3<<" gal/Mpc^3"<<endl;
 
  double mall=0., mgood=0.;