Changeset 3191 in Sophya for trunk/Cosmo/RadioBeam/multicyl.cc

Timestamp:

Mar 19, 2007, 9:29:05 AM (19 years ago)

Author:

legoff

Message:

Cylinders have different x positions (NS).
This is taken into account in the signal phase shift
and in the exp() term when combining cylinders.
We can have more bins in x than antennas in one cylinder
to see the gain in resolution along x due to several cylinders at different x

File:

• trunk/Cosmo/RadioBeam/multicyl.cc (modified) (2 diffs)

trunk/Cosmo/RadioBeam/multicyl.cc

@@ -84 +84 @@
 void MultiCylinders::ReconstructSourceBox(int halfny, double stepangy)
 {
-
   ReconstructCylinderPlaneS(true);
   TMatrix< complex<r_4> > & mtx = GetCylinder(0).getRecSrcPlane(); 
   // boite 3D X:angX, Y:angY , Z: freq
   sa_size_t sz[5] = {0,0,0,0,0};
-  sz[0] = mtx.NCols();
+//  int nx=mtx.NCols(); // uncomment to go back to nxbin = nantenna
+  int nx=256;
+  sz[0] = nx;
   sz[1] = halfny*2+1;
   sz[2] = mtx.NRows();
@@ -113 +114 @@
 
       double facl = - 2*M_PI*frq*cyl.getCylinderYPos();  // attention au signe -
+//      double facl_x = - 2*M_PI*cyl.getCylinderXPos()/cyl.Da/(double)cyl.NR; 
+      double facl_x = - 2*M_PI*cyl.getCylinderXPos()/cyl.Da; 
       double dphi;
       complex< r_4 > phasefactor;
       int jyy = 0;
       for(int jy=-halfny; jy<=halfny; jy++) {  // Loop over Y angular steps
-        double dphi = facl * sin( (double)jy*stepangy );
-        phasefactor = complex< r_4 > ((r_4) cos(dphi) , (r_4) sin(dphi));
-        for(int ix=0; ix<sz[0]; ix++) {  // Loop over AngX directions
+        for(int ix=0; ix<nx; ix++) {  // Loop over AngX directions
+          double dphi = facl * sin( (double)jy*stepangy )
+                        + facl_x*(double(ix)/double(nx)-1./2.);
+          phasefactor = complex< r_4 > ((r_4) cos(dphi) , (r_4) sin(dphi));
           // sur recp : index ligne -> frequence  ,    index colonne -> angX ,
-          box(ix, jyy, kf) += recp(kf, ix)*phasefactor;
+          int ixx=(int)(ix*(double)cyl.NR/double(nx));
+          box(ix, jyy, kf) += recp(kf, ixx)*phasefactor;
         }  //  
         jyy++;