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

Timestamp:

Jan 31, 2007, 7:17:31 PM (19 years ago)

Author:

ansari

Message:

modifs et corrections diverses ds treccyl.cc et recbeam.pic, Reza 31/01/2007

File:

• trunk/Cosmo/RadioBeam/treccyl.cc (modified) (6 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 = 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 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 gmean = 1.;    // gain moyen
 static double gsig = 0.;     // sigma des gains
@@ -149 +149 @@
 
   // BRSourceGen sg;
-  int nsf = 7;
+  int nsf = 6;
   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/150);
+  BRSourceGen sg(frq, nsf, M_PI/3, M_PI/60);
+  /*
   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.;
+  for(is=nsf; is<2*nsf; is++)     sg.angY(is) = M_PI/30.;
+  for(is=2*nsf; is<3*nsf; is++)     sg.angY(is) = -M_PI/30.;
+  */
   // sg.WritePPF(string("brsrcm.ppf"));
   // BRSourceGen  sg(string("brsrcm.ppf"));
@@ -177 +179 @@
   mcyl.SetGains(gmean, gsig, nantgz);
 
+  /*
+  for(int kkc=0; kkc<12; kkc++) 
+    mcyl.AddCylinder(5.*kkc);
+  */
+ 
   mcyl.AddCylinder(0.);
   mcyl.AddCylinder(5.);
-  mcyl.AddCylinder(10.);
-  mcyl.AddCylinder(18.);
-  mcyl.AddCylinder(30.);
-  mcyl.AddCylinder(46.);
-  mcyl.AddCylinder(64.);
-
+  mcyl.AddCylinder(15.);
+  mcyl.AddCylinder(35.);
+  mcyl.AddCylinder(65.);
+ 
   mcyl.SetSources(sg);
   
@@ -190 +195 @@
 
   //  mcyl.ReconstructCylinderPlaneS(true);
-  mcyl.ReconstructSourceBox(15, M_PI/500.);
+  mcyl.ReconstructSourceBox(10, M_PI/600.);
 
   POutPersist po(ppfname);
@@ -219 +224 @@
   int kfmin, kfmax;
   po << PPFNameTag("recsrcbox") << mcyl.getRecSrcBox();
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[0] - 1;  kfmax = kfmin+2;
+  kfmin = mcyl.getRecSrcBox().SizeZ()/0.5*frq[0] - 2;  kfmax = kfmin+2;
+  cout << "testmulticyl/Info: slice0 kfmin=" << kfmin << " kfmax=" << kfmax << endl;
   {
   TMatrix<r_4> slice0 = mcyl.getRecXYSlice(kfmin, kfmax);
   po << PPFNameTag("recXYf0") << slice0;
   }
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[1] - 1;  kfmax = kfmin+2;
+  kfmin = mcyl.getRecSrcBox().SizeZ()/0.5*frq[1] - 2;  kfmax = kfmin+2;
+  cout << "testmulticyl/Info: slice1 kfmin=" << kfmin << " kfmax=" << kfmax << endl;
   {
   TMatrix<r_4> slice1 = mcyl.getRecXYSlice(kfmin, kfmax);
   po << PPFNameTag("recXYf1") << slice1;
   }
-  kfmin = mcyl.getRecSrcBox().SizeZ()*frq[2] - 1;  kfmax = kfmin+2;
+  kfmin = mcyl.getRecSrcBox().SizeZ()/0.5*frq[2] - 2;  kfmax = kfmin+2;
+  cout << "testmulticyl/Info: slice2 kfmin=" << kfmin << " kfmax=" << kfmax << endl;
   {
   TMatrix<r_4> slice2 = mcyl.getRecXYSlice(kfmin, kfmax);