Changeset 3341 in Sophya for trunk/Cosmo/SimLSS/cmvdefsurv.cc

Timestamp:

Oct 5, 2007, 7:23:27 PM (18 years ago)

Author:

cmv

Message:

encore mise au point du calcul du bruit cmv 05/10/2007

File:

• trunk/Cosmo/SimLSS/cmvdefsurv.cc (modified) (4 diffs)

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -445 +445 @@
  // - Interferometre a N antennes (N grand):
  // Slim ->  2 * k * Tsys / [ Atot * Sqrt(dnu*Tobs) ]
- //      C'est aussi la formule pour un telescope unique
+ //      C'est aussi la formule pour un telescope unique de surface Atot
  // --- On ne mesure qu'une seule polarisation
  // Ces formules sont valables si on mesure 1 polarisation:
@@ -503 +503 @@
      fac = 1.;
      cout<<"\n...Observation limits for a N (large) telescope interferometer (with complex correlator)"<<endl
-         <<"     (weak source limit sensitivity in a synthetised image)"<<endl;
+         <<"     (weak source limit sensitivity in a synthetised image)"<<endl
+         <<"     Also valid for a single dish telescope."<<endl;
    } else continue;
 
@@ -509 +510 @@
                / sqrt(fac*(dnuhiz*1.e+9)*tobs) /Jansky2Watt_cst;
    if(ya2polar) slim /= sqrt(2.);
-   SsN = ssig_2polar/slim;
+   SsN = ssig_2polar / slim;
    smass = mhiref / ssig_2polar * slim;
    cout<<"for 1 lobe:"<<endl
@@ -515 +516 @@
        <<"   S/N     = "<<SsN<<endl
        <<"   Mass HI = "<<smass<<" Msol"<<endl;
+
    slim_nl = slim * sqrt(nlobes);
-   SsN_nl = ssig/slim_nl;
-   smass_nl = mhiref/ssig*slim_nl;
+   SsN_nl = ssig_2polar / slim_nl;
+   smass_nl = mhiref / ssig_2polar * slim_nl;
    cout<<"for "<<nlobes<<" lobes:"<<endl
        <<"   Flux    = "<<slim_nl<<" Jy"<<endl