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

Timestamp:

Feb 1, 2007, 3:15:16 PM (19 years ago)

Author:

ansari

Message:

petites ameliorations + reglage parametres -> version +/- utilisable , Reza 01/02/2007

File:

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

trunk/Cosmo/RadioBeam/treccyl.cc

@@ -21 +21 @@
 #include "mbeamcyl.h"
 
+/* 
+   Projet BAORadio / HSHS
+   Programme de simulation pour reconstruction de lobe radio.
+   programme principal de test
+
+   R. Ansari - LAL      Jan 2007
+
+*/
 
 // Declaration des fonctions de ce fichier
 static int test1cyl(string& ppfname);
-static int testmulticyl(string& ppfname);
+static int testmulticyl(string& ppfname, int ncyl=5);
 
 //-----------------------------------------------------------
@@ -33 +41 @@
 static double freq0 = 2.;  // frequence de base
 static double da = 0.25;     // pas des antennes le long du cylindre
+// ATTENTION : les parametres suivants sont relies a MR/da
+static double maxangX = M_PI/3.; // angle max en X ( +/- )
+static double maxangY = M_PI/60.; // angle max en Y ( +/- )
+static int nsrcmax = 50;  // Nb total de sources - en un plan
+
 static double snoise = 1.0;  // sigma du bruit 
 static double tjit = 0.05;   // sigma du jitter en temps
@@ -56 +69 @@
 string ppfname = "treccyl.ppf";
 int act = 1;
+int ncyl = 5; 
 if (narg < 3) {
   cout << "Usage: treccyl act ppfname [PrtLev=0] \n"
-       << "  act= X ou XY  , ppfname=  treccyl.ppf par defaut" << endl;
+       << " -act= X ou XY5 ou XY12  (5 ou 12 cylindres) \n"
+       << " -ppfname=  treccyl.ppf par defaut" << endl;
   return 1;
 }
-if (strcmp(arg[1],"XY") == 0) act = 2;
+if (strcmp(arg[1],"XY5") == 0) { act = 2;  ncyl = 5; }
+else if (strcmp(arg[1],"XY12") == 0) { act = 2;  ncyl = 12; }
 if (narg > 2)  ppfname = arg[2];
 if (narg > 3)  prtlevel = atoi(arg[3]);
 
 int rc = 0;
+cout << ">>>> treccyl : " << arg[1] << " PPFName=" << ppfname << endl;
 try {
-  if (act == 2) rc = testmulticyl(ppfname);
+  if (act == 2) rc = testmulticyl(ppfname, ncyl);
   else rc = test1cyl(ppfname);
-  cout << "treccy/Info - FIN  " << endl;
 }
   catch (PThrowable& exc) {
@@ -85 +101 @@
   } 
 
- return rc;
+cout << ">>>> treccyl ------- FIN ----------- Rc=" << rc << endl;
+return rc;
 }
 
@@ -95 +112 @@
   // BRSourceGen sg;
   int nsrc = 60; 
-  BRSourceGen sg(nsrc);
+  BRSourceGen sg(nsrcmax, maxangX, 0.);
   //  sg.WritePPF(string("brsrc1.ppf"));
 
@@ -119 +136 @@
   cout << "treccy/test1cyl:  signal vectors OK " << endl;
   PrtTim("test1cyl:[1] ");
+  
+  cout << "--- treccy/test1cyl: Saving to PPF file " << ppfname << endl;
 
   POutPersist po(ppfname);
@@ -145 +164 @@
 
 //--- Fonction de test : reconstruction cube AngX-AngY-Frequence (multi-cylindre)
-int testmulticyl(string& ppfname)
+int testmulticyl(string& ppfname, int ncyl)
 {
+
+  if ((ncyl != 5) && (ncyl != 12)) ncyl = 5;
 
   // BRSourceGen sg;
@@ -156 +177 @@
   
   
-  cout << "treccy/testmulticyl: BRSourceGen sg([frq=0.1,0.27,0.38], " << nsf << ")" << endl;
-  BRSourceGen sg(frq, nsf, M_PI/3, M_PI/60);
-  /*
+  cout << "testmulticyl: BRSourceGen sg([frq=0.1,0.27,0.38], " << nsf 
+       << "," << maxangX << "," << maxangY << ")" << endl;
+  BRSourceGen sg(frq, nsf, maxangX, maxangY);
+  
   int is;
-  for(is=0; is<nsf; is++)     sg.angY(is) = 0.;
-  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.;
-  */
+  double fay[6] = {-0.7,-0.5,0.,0.,0.5,0.7};
+  for(is=0; is<3*nsf; is++) {
+    int ism = is%nsf;
+    sg.angX(is) = maxangX*(ism-2.5)/3.;
+    sg.angY(is) = maxangY*fay[ism];
+  }
   // sg.WritePPF(string("brsrcm.ppf"));
   // BRSourceGen  sg(string("brsrcm.ppf"));
-  cout << "=== testmulticyl: BRSourceGen NbSrc= " << sg.NbSources() 
-       << " NbRecep=" << MR << " NSamples=" << NE << endl;
+  cout << "=== testmulticyl: NbSrc= " << sg.NbSources() 
+       << " NbRecep=" << MR << " NSamples=" << NE << " NCyl=" << ncyl << endl;
   if (prtlevel > 1)  sg.Print(cout);
   
@@ -179 +203 @@
   mcyl.SetGains(gmean, gsig, nantgz);
 
-  /*
-  for(int kkc=0; kkc<12; kkc++) 
-    mcyl.AddCylinder(5.*kkc);
-  */
- 
-  mcyl.AddCylinder(0.);
-  mcyl.AddCylinder(5.);
-  mcyl.AddCylinder(15.);
-  mcyl.AddCylinder(35.);
-  mcyl.AddCylinder(65.);
- 
+  if (ncyl == 12) {
+    cout << " --- testmulticyl/ NCyl= " << ncyl << " posY=0 ... " << (ncyl-1)*5. 
+         << " (EqualSpacing)" << endl;
+    for(int kkc=0; kkc<12; kkc++) {
+      cout << "..." << kkc << " - mcyl.AddCylinder(posY= " << 5.*kkc << " )" << endl;
+      mcyl.AddCylinder(5.*kkc);
+    }
+  }
+  else {
+    cout << " --- testmulticyl/ NCyl= " << ncyl << " posY=0 ... 55 (IrregularSpacing)" << endl;
+    double posyc[5] = {0.,5.,15.,35.,55.};
+    for(int kkc=0; kkc<5; kkc++) {
+      cout << "..." << kkc << " - mcyl.AddCylinder(posY= " << posyc[kkc] << " )" << endl;
+      mcyl.AddCylinder(posyc[kkc]);
+    }
+  }
+
   mcyl.SetSources(sg);
   
@@ -195 +225 @@
 
   //  mcyl.ReconstructCylinderPlaneS(true);
-  mcyl.ReconstructSourceBox(10, M_PI/600.);
-
+  mcyl.ReconstructSourceBox(10, maxangY/10.);
+
+  cout << "--- treccy/testmulticyl: Saving to PPF file " << ppfname << endl;
   POutPersist po(ppfname);