Changeset 3792 in Sophya for trunk/Cosmo/RadioBeam/applobe.cc

Timestamp:

Jun 28, 2010, 6:06:00 PM (15 years ago)

Author:

ansari

Message:

Ajout du programme de calcul et sauvegarde reponse interfero, Reza 28/06/2010

File:

• trunk/Cosmo/RadioBeam/applobe.cc (modified) (4 diffs)

trunk/Cosmo/RadioBeam/applobe.cc

@@ -3 +3 @@
     R. Ansari , C. Magneville - Juin 2010 
 
-  Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamArcmin]
+  Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamDoL]
 ---------------------------------------------------------------  */
 
@@ -50 +50 @@
 
   if (narg < 3) {
-    cout << "Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamArcmin]\n" << endl;
+    cout << "Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamDoL]\n" << endl;
     return 1;
   }
@@ -64 +64 @@
   }
   bool fgcorrbeam=false;
-  double tbeamarcmin=15.;
+  double tbeamDoL=135;
   if (narg>4) {
-    tbeamarcmin=atof(arg[4]);
-    if (tbeamarcmin>1.e-3)  fgcorrbeam=true;
+    tbeamDoL=atof(arg[4]);
+    if (tbeamDoL>1)  fgcorrbeam=true;
   }
 
   int rc = 91;
 
-  cout << " ====== applobe :   Input NVSS catalog name= " << inname << " OutName=" << outname;
+  cout << " ====== applobe :   Input skycube  name= " << inname << " OutName=" << outname;
   bool fginmap=true;
   try {
@@ -105 +105 @@
 
     if (fgcorrbeam) {
-      double DoL = 1.22/ArcminToRadian(tbeamarcmin);
+      double DoL = tbeamDoL;
+      double tbeamarcmin = RadianToDegree(1.22/DoL)*60.;
       Four2DResponse tbeam(2, DoL, DoL );
-      cout << "applobe[3]: calling Correct2RefLobe() with target beam:" << tbeamarcmin 
-           << " arcmin -> D/Lambda=" << DoL << endl;
+      cout << "applobe[3]: calling Correct2RefLobe() with target beam D/Lambda=" << DoL  
+           << " -> arcmin " << tbeamarcmin << endl;
       beam.Correct2RefLobe(tbeam, incube, dx, dy, Freq0MHz, dfreq);
     }