Changeset 3358 in Sophya for trunk/Cosmo/SimLSS/cmvobserv3d.cc
- Timestamp:
- Oct 23, 2007, 12:19:27 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/SimLSS/cmvobserv3d.cc
r3354 r3358 539 539 540 540 //----------------------------------------------------------------- 541 cout<<endl<<"\n--- Converting fluctuations into mass"<<endl;542 fluct3d.TurnFluct2Mass();543 nm = fluct3d.MeanSigma2(rm,rs2);544 PrtTim(">>>> End Converting fluctuations into mass");545 546 541 if(no_poisson) { 547 542 cout<<"\n--- Converting !!!DIRECTLY!!! mass into HI mass: mass per pixel =" 548 543 <<mass_by_pixel<<endl; 549 rm = fluct3d.Turn Mass2MeanNumber(mass_by_mpc3); // ici on doit donner Msol/Mpc^3544 rm = fluct3d.TurnFluct2MeanNumber(mass_by_mpc3); // ici on doit donner Msol/Mpc^3 550 545 } else { 551 546 cout<<"\n--- Converting mass into galaxy number: gal per pixel =" 552 547 <<ngal_by_mpc3*fluct3d.GetDVol()<<endl; 553 rm = fluct3d.Turn Mass2MeanNumber(ngal_by_mpc3);548 rm = fluct3d.TurnFluct2MeanNumber(ngal_by_mpc3); // ici on doit donner Ngal/Mpc^3 554 549 } 555 550 nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200); 556 551 nm = fluct3d.MeanSigma2(rm,rs2); 552 if(rm>0.) cout<<"normalised sigma(dM/M) is "<<sqrt(rs2)/rm<<endl; 557 553 PrtTim(">>>> End Converting mass into galaxy number or mass"); 558 554
Note:
See TracChangeset
for help on using the changeset viewer.