Changeset 3789 in Sophya for trunk/Cosmo/RadioBeam/lobe.cc

Timestamp:

Jun 27, 2010, 4:57:36 PM (15 years ago)

Author:

ansari

Message:

Corrections diverses, Reza 27/06/2010

File:

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

trunk/Cosmo/RadioBeam/lobe.cc

@@ -6 +6 @@
 
 
-#include "fftwserver.h"    
+#include "fftwserver.h" 
+#include "matharr.h"    
 #include "ctimer.h"    
 
@@ -38 +39 @@
     conv.setFrequency(f0+kz*df);
     fresp_.setLambda(conv.getLambda());
+    //    cout << " DEBUG*" << kz << " lambda=" << conv.getLambda() << " lambda_ratio_=" << fresp_.lambda_ratio_ << endl;
     ApplyLobeK2D(fresp_, fourAmp, dkx, dky);
     ffts.FFTBackward(fourAmp, slice, true);
     if (kz%20==0)  cout << "BeamEffect::ApplyLobe3D() done kz=" << kz << " / a.SizeZ()=" << a.SizeZ() << endl; 
   }
+  double mean, sigma;
+  TF min, max;
+  a.MinMax(min, max);
+  MeanSigma(a, mean, sigma);
+  cout << " BeamEffect::ApplyLobe3D() - Result Min=" << min << " Max=" << max 
+       << " Mean=" << mean << " Sigma=" << sigma << endl; 
   return;
 }
@@ -64 +72 @@
     conv.setFrequency(f0+kz*df);
     fresp_.setLambda(conv.getLambda());
-    Four2DRespRatio rratio(fresp_, rep);
+    Four2DRespRatio rratio(rep, fresp_);
     ApplyLobeK2D(rratio, fourAmp, dkx, dky);
     ffts.FFTBackward(fourAmp, slice, true);
     if (kz%20==0)  cout << "BeamEffect::Correct2RefLobe() done kz=" << kz << " / a.SizeZ()=" << a.SizeZ() << endl; 
   }
+  double mean, sigma;
+  TF min, max;
+  a.MinMax(min, max);
+  MeanSigma(a, mean, sigma);
+  cout << " BeamEffect::Correct2RefLobe() - Result Min=" << min << " Max=" << max 
+       << " Mean=" << mean << " Sigma=" << sigma << endl; 
   return;
 }
@@ -74 +88 @@
 /* --Methode-- */
 void BeamEffect::ApplyLobeK2D(Four2DResponse& rep, TArray< complex<TF> >& fourAmp, double dkx, double dky)
-//  dx, dy en radians, lambda en metres 
 {
   double kxx, kyy;