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

Timestamp:

Apr 18, 2011, 5:30:44 PM (14 years ago)

Author:

ansari

Message:

Corrections diverses: choix lobe gaussien/triangle et specif DishDiameter au lieu de DoL ds applobe/calcpk2, possibilite application lobe freq.independante ds applobe, Reza 18/04/2011

File:

• trunk/Cosmo/RadioBeam/lobe.cc (modified) (1 diff)

trunk/Cosmo/RadioBeam/lobe.cc

@@ -18 +18 @@
 }
 
+/* --Methode-- */
+void BeamEffect::ApplyLobe(TArray< TF >& a, double dx, double dy, double f0)
+// dx, dy en radioans, f0, df en MHz
+{
+  Timer tm("BeamEffect::ApplyLobe");
+  FFTWServer ffts(true);                     
+  ffts.setNormalize(true); 
+  
+  H21Conversions conv;
+  conv.setFrequency(f0);
+  fresp_.setLambda(conv.getLambda());
+
+  TArray< complex<TF> > fourAmp;
+  double dkx = DeuxPI/(double)a.SizeX()/dx;
+  double dky = DeuxPI/(double)a.SizeY()/dy;
+  
+  for(sa_size_t kz=0; kz<a.SizeZ(); kz++) {
+    TArray< TF > slice( a(Range::all(), Range::all(), kz) );
+    ffts.FFTForward(slice, fourAmp);
+    ApplyLobeK2D(fresp_, fourAmp, dkx, dky);
+    ffts.FFTBackward(fourAmp, slice, true);
+    if (kz%20==0)  cout << "BeamEffect::ApplyLobe() done kz=" << kz << " / a.SizeZ()=" << a.SizeZ() << endl; 
+  }
+  double mean, sigma;
+  TF min, max;
+  a.MinMax(min, max);
+  MeanSigma(a, mean, sigma);
+  cout << " BeamEffect::ApplyLobe() - Result Min=" << min << " Max=" << max 
+       << " Mean=" << mean << " Sigma=" << sigma << endl; 
+  return;
+}