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

Timestamp:

Aug 1, 2007, 11:04:27 AM (18 years ago)

Author:

cmv

Message:

variation de la taille du lobe avec la frequence cmv 01/08/2007

File:

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

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -15 +15 @@
 
 /* --- Check Peterson at al. astro-ph/0606104 v1
-cmvdefsurv -z 0.0025 -x 1 -U 0.75,0.3,0.7,-1,1 -V 300 -O 400000,6000 -N 75 -M 6.156e9 -F 3 -2 1.5
+cmvdefsurv -U 0.75,0.3,0.7,-1,1 -V 300 -z 0.0025,0.2,Z -x 1,90,A -O 400000,6000 -N 75 -M 6.156e9 -F 3 -2 1.5
  --- */
 
@@ -38 +38 @@
       <<"   \'M\' : en distance (pour X-Y-Z) : resolution et largeur Mpc"<<endl
       <<"----------------"<<endl
-      <<" -L lobewidth : taille du lobe d\'observation (FWHM) en arcmin (def= 1\')"<<endl
-      <<"                Si <=0 alors angle solide du lobe = celui du pixel"<<endl
       <<" -O surf,tobs : surface effective (m^2) et temps d\'observation (s)"<<endl
       <<" -N Tsys : temperature du system (K)"<<endl
+      <<" -L lobewidth,freqlob : taille du lobe d\'observation (FWHM) en arcmin (def= 1\')"<<endl
+      <<"                        pour la frequence freqlob en MHz"<<endl
+      <<"            Si lobewidth<=0 : l'angle solide du lobe = celui du pixel"<<endl
+      <<"            Si freqlob<=0 : la frequence de reference est celle du redshift etudie"<<endl
       <<" -2 : two polarisations measured"<<endl
       <<" -M  : masse de HI de reference (MSol), si <=0 mean schechter in pixel"<<endl
@@ -79 +81 @@
  double redshift = 0.;
  double tobs = 6000., surfeff = 400000.;
- double lobewidth = 1.;  // taille du lobe d'observation en arcmin
+ // taille du lobe d'observation en arcmin pour la frequence
+ double lobewidth0 = -1., lobefreq0 = Fr_HyperFin_Par*1.e3;
  double Tsys=75.;
  // a 408 MHz (Haslam) + evol index a -2.6
@@ -106 +109 @@
     break;
   case 'L' :
-    sscanf(optarg,"%lf",&lobewidth);
+    sscanf(optarg,"%lf,%lf",&lobewidth0,&lobefreq0);
     break;
   case 'N' :
@@ -316 +319 @@
  cout<<"dang lumi = "<<dlum<<" in ["<<dlumlim[0]<<","<<dlumlim[1]<<"] Mpc"<<endl;
 
+ double lobewidth = lobewidth0;   // ArcMin
+ if(lobefreq0<=0.) lobefreq0 = nuhiz*1.e3; // MHz
+ // La taille angulaire du lobe change avec la frequence donc avec le redshift
+ lobewidth *= lobefreq0/(nuhiz*1.e3);
+ cout<<"\n--- Lobe: width="<<lobewidth0<<" pour "<<lobefreq0<<" MHz"<<endl
+     <<"          changed to "<<lobewidth<<" pour "<<nuhiz*1.e3<<" MHz"<<endl;
  double slobe = lobewidth/2.35482; // sigma du lobe en arcmin
  double lobecyl = sqrt(8.)*slobe; // diametre du lobe cylindrique equiv en arcmin
@@ -321 +330 @@
  double nlobes = rad2min(adtx)*rad2min(adty)/lobearea;
  if(lobewidth<=0.) nlobes = 1.;
- cout<<"\nBeam FWHM = "<<lobewidth<<"\' -> sigma = "<<slobe<<"\' -> "
+ cout<<"Beam FWHM = "<<lobewidth<<"\' -> sigma = "<<slobe<<"\' -> "
      <<" Dcyl = "<<lobecyl<<"\' -> area = "<<lobearea<<" arcmin^2"<<endl;
  cout<<"Number of beams in one transversal pixel = "<<nlobes<<endl;