Changeset 3163 in Sophya for trunk/Cosmo/RadioBeam/treccyl.cc

Timestamp:

Jan 30, 2007, 7:05:40 PM (19 years ago)

Author:

ansari

Message:

debug divers ds RadioBeam - Reza 30/01/2007

File:

• trunk/Cosmo/RadioBeam/treccyl.cc (modified) (7 diffs)

trunk/Cosmo/RadioBeam/treccyl.cc

@@ -33 +33 @@
 static double freq0 = 2.;  // frequence de base
 static double da = 0.25;     // pas des antennes le long du cylindre
-static double snoise = 1.0;  // sigma du bruit 
-static double tjit = 0.05;   // sigma du jitter en temps
-static double tos = 0.02;    // sigma des offsets en temps
+static double snoise = 0.; //1.0;   sigma du bruit 
+static double tjit = 0.; //0.05;    sigma du jitter en temps
+static double tos = 0.; //0.02;     sigma des offsets en temps
 static double gmean = 1.;    // gain moyen
 static double gsig = 0.;     // sigma des gains
@@ -149 +149 @@
 
   // BRSourceGen sg;
-  int nsf = 10;
+  int nsf = 7;
   vector<double> frq;
   frq.push_back(0.1);
@@ -157 +157 @@
   
   cout << "treccy/testmulticyl: BRSourceGen sg([frq=0.1,0.27,0.38], " << nsf << ")" << endl;
-  BRSourceGen sg(frq, nsf, M_PI/3, M_PI/1000);
+  BRSourceGen sg(frq, nsf, M_PI/3, M_PI/150);
+  int is;
+  for(is=0; is<nsf; is++)     sg.angY(is) = 0.;
+  for(is=nsf; is<2*nsf; is++)     sg.angY(is) = M_PI/150.;
+  for(is=2*nsf; is<3*nsf; is++)     sg.angY(is) = -M_PI/150.;
   // sg.WritePPF(string("brsrcm.ppf"));
   // BRSourceGen  sg(string("brsrcm.ppf"));
@@ -174 +178 @@
 
   mcyl.AddCylinder(0.);
-  mcyl.AddCylinder(100.);
-  mcyl.AddCylinder(300.);
-  mcyl.AddCylinder(700.);
-  mcyl.AddCylinder(1500.);
+  mcyl.AddCylinder(5.);
+  mcyl.AddCylinder(10.);
+  mcyl.AddCylinder(18.);
+  mcyl.AddCylinder(30.);
+  mcyl.AddCylinder(46.);
+  mcyl.AddCylinder(64.);
 
   mcyl.SetSources(sg);
@@ -184 +190 @@
 
   //  mcyl.ReconstructCylinderPlaneS(true);
-  mcyl.ReconstructSourceBox();
+  mcyl.ReconstructSourceBox(15, M_PI/500.);
 
   POutPersist po(ppfname);
@@ -192 +198 @@
 
   {
-  TMatrix<r_4> srcplane0 = module(mcyl.GetCylinder(0).getRecSrcPlane());
+    //  TMatrix<r_4> srcplane0 = module(mcyl.GetCylinder(0).getRecSrcPlane());
+  TMatrix< complex<r_4> > srcplane0 = mcyl.GetCylinder(0).getRecSrcPlane();
   po << PPFNameTag("recsrcplane0") << srcplane0;
   }
 
   {
-  TMatrix<r_4> srcplane3 = module(mcyl.GetCylinder(3).getRecSrcPlane());
+    //  TMatrix<r_4> srcplane2 = module(mcyl.GetCylinder(3).getRecSrcPlane());
+  TMatrix< complex<r_4> > srcplane2 = mcyl.GetCylinder(2).getRecSrcPlane();
+  po << PPFNameTag("recsrcplane2") << srcplane2;
+  }
+
+  {
+    //  TMatrix<r_4> srcplane3 = module(mcyl.GetCylinder(3).getRecSrcPlane());
+  TMatrix< complex<r_4> > srcplane3 = mcyl.GetCylinder(0).getRecSrcPlane();
   po << PPFNameTag("recsrcplane3") << srcplane3;
   }
@@ -205 +219 @@
   int kfmin, kfmax;
   po << PPFNameTag("recsrcbox") << mcyl.getRecSrcBox();
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[0] - 3;  kfmax = kfmin+6;
+  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[0] - 1;  kfmax = kfmin+2;
   {
   TMatrix<r_4> slice0 = mcyl.getRecXYSlice(kfmin, kfmax);
   po << PPFNameTag("recXYf0") << slice0;
   }
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[1] - 3;  kfmax = kfmin+6;
+  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[1] - 1;  kfmax = kfmin+2;
   {
   TMatrix<r_4> slice1 = mcyl.getRecXYSlice(kfmin, kfmax);
   po << PPFNameTag("recXYf1") << slice1;
   }
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[2] - 3;  kfmax = kfmin+6;
+  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[2] - 1;  kfmax = kfmin+2;
   {
   TMatrix<r_4> slice2 = mcyl.getRecXYSlice(kfmin, kfmax);