Changeset 3909 in Sophya for trunk/AddOn/TAcq/brviscalc.cc

Timestamp:

Nov 21, 2010, 2:26:12 PM (15 years ago)

Author:

ansari

Message:

---

1/ Mise en place d'informations auxiliaires au niveau de RAcqMemZoneMgr ,

FillTime (DATEOBS) pour chaque zone et identificateur des fibres

2/ remplissage de FillTime(DATEOBS) Et FiberId's de MemZoneMgr par

les classes PCIEMultiReader et EthernetReader

3/ Ajout de mots-cle (ds les datacards) pour numeros d'identification des

fibres et parametres de controles du DMA (@fiberids @waitenddma)
et corrections/ameliorations des modes acq (raw1c ordfft2c ...)

4/ Parametrage des WaitEndDMA ds PCIEWrapperInterface et implementation

au niveau de DMAMgrV6 (dmamgrv6.h .cc)

5/ Gestion d'echec DMA -> arret acquisition ds PCIEMultiReader,EthernetReader

modification TimeOut DMA entre le premier DMA et les suivants pour permettre
un arret propre rapide

6/ Adaptation de la classe BRVisibilityCalculator pour utilisation

FillTime (DATEOBS) et FiberIds de MemZoneMgr - amelioration de
l'identification des numeros de canaux pour les matrices de visibilites

7/ Petites ameliorations (ecriture mot-cles) pour MultiDataSaver et

amelioration classe BRMultiFitsReader pour lecture FillTime (DATEOBS)
et FiberIds depuis les fichiers FITS et remplissage de RAcqMemZoneMgr

Reza, 21/11/2010

File:

• trunk/AddOn/TAcq/brviscalc.cc (modified) (9 diffs)

trunk/AddOn/TAcq/brviscalc.cc

@@ -32 +32 @@
 
   uint_4 maxnpairs = (2*memgr_.NbFibres()+1)*memgr_.NbFibres();
-  chanum_.SetSize(maxnpairs);
-  sa_size_t k=0;
-  for(size_t i=0; i<2*memgr_.NbFibres(); i++)  vpdata_[i]=NULL;
-  for(size_t i=0; i<2*memgr_.NbFibres(); i++) {
-    for(size_t j=i; j<2*memgr_.NbFibres(); j++) {
-      chanum_(k) = (i+1)*100+(j+1);  k++;
-    }
-  }
+  chanids_.SetSize(2*memgr_.NbFibres());
+  chanpairnumall_.SetSize(maxnpairs);
+  chanpairsall_.SetSize(maxnpairs,2);
+  for(size_t i=0; i<2*memgr_.NbFibres(); i++)   vpdata_[i]=NULL;
   SelectPairs();
 
@@ -82 +78 @@
   po << visdt_;
   if (calcid_ == 0) {
-    POutPersist poc(outpath_+"chanum.ppf");
-    poc << chanum_;
     
     if (fgcktt_) {
@@ -130 +124 @@
   nbpairs_=nbpairs;
   vismtx_.SetSize(nbpairs_, paq.DataSize()/4); 
+
+  chanpairnum_.SetSize(nbpairs_);
+  chanpairs_.SetSize(nbpairs_,2);
+
   return nbpairs_;
 }
@@ -147 +145 @@
 }
 
+/* --Methode-- */
+void BRVisibilityCalculator::UpdateChanIds()
+{
+
+  for(size_t i=0; i<memgr_.NbFibres(); i++) { 
+    chanids_(2*i)=memgr_.FiberId(i)*2-1; 
+    chanids_(2*i+1)=memgr_.FiberId(i)*2; 
+  } 
+  sa_size_t k=0;     // numero de ligne dans la matrice des visibilites 
+  for(size_t i=0; i<vpdata_.size(); i++) {
+    for(size_t j=i; j<vpdata_.size(); j++) {
+      chanpairnumall_(k)=chanids_(i)*CHANPAIRCONVFAC+chanids_(j);
+      chanpairsall_(k,0)=chanids_(i);  chanpairsall_(k,1)=chanids_(j);     k++;
+    }
+  }
+  sa_size_t kpair=0;  
+  k=0;  // numero de ligne dans la matrice des visibilites 
+  for(size_t i=0; i<vpdata_.size(); i++) {
+    for(size_t j=i; j<vpdata_.size(); j++) {
+      kpair++;
+      if (kpair<(pairst_+1))  continue;
+      if (kpair>=(pairst_+nbpairs_+1))  break;
+      chanpairnum_(k)=chanids_(i)*CHANPAIRCONVFAC+chanids_(j);
+      chanpairs_(k,0)=chanids_(i);  chanpairs_(k,1)=chanids_(j);     k++;
+    }
+  }
+
+  string filename;
+  filename = outpath_+"chanum.ppf";
+  if (nbcalc_>1) {
+    char sbuff[32];
+    sprintf(sbuff,"chanum_%d.ppf",(int)calcid_);
+    filename = outpath_+sbuff;
+  }
+  POutPersist poc(outpath_+"chanum.ppf");
+  poc << PPFNameTag("chanids") << chanids_; 
+  poc << PPFNameTag("chanpairs") << chanpairs_; 
+  poc << PPFNameTag("chanpairnum") << chanpairnum_; 
+  poc << PPFNameTag("chanpairsall") << chanpairsall_; 
+  poc << PPFNameTag("chanpairnumal") << chanpairnumall_; 
+  cout << "BRVisibilityCalculator::UpdateChanIds() Channel Ids/Pairs saved to PPF file " << filename << endl; 
+  return;
+}
 
 /* --Methode-- */
@@ -207 +248 @@
         break;  
       }
+      cts_=memgr_.GetAuxData(mid)->FillTime();   // get associated date/time (DATEOBS)
+
       for(size_t fib=0; fib<(size_t)memgr_.NbFibres(); fib++) {
         fbuff_[fib] = memgr_.GetMemZone(mid,fib);
@@ -238 +281 @@
         vismtx_ = complex<r_4>((r_4)0.,(r_4)0.);
         moyfc_=moytt_=0.;
-        first_tmstamp_.SetNow();  // Current date and time 
+        //      first_tmstamp_.SetNow();  // Current date and time 
+        first_tmstamp_=cts_;  // Current date and time 
       }
 
@@ -326 +370 @@
     first_fc_=curfc_[0];
     first_tt_= (vpaq_[0].TimeTag()-ttfirst_[0]);
-    first_tmstamp_.SetNow();  // Current date and time 
+    // first_tmstamp_.SetNow();  // Current date and time 
+    first_tmstamp_=cts_;  // Current date and time 
   }
 
@@ -439 +484 @@
       for(sa_size_t jf=jf1_; jf<jf2_; jf++) {
         xnt[2]=jf;
-        xnt[3]=chanum_(rv+pairst_);
+        xnt[3]=chanpairnumall_(rv+pairst_);
         xnt[4]=vismtx_(rv,jf).real()/(r_4)(nmean_);
         xnt[5]=vismtx_(rv,jf).imag()/(r_4)(nmean_);
@@ -458 +503 @@
         }
         xnt[2]=jf+djf_/2;
-        xnt[3]=chanum_(rv+pairst_);
+        xnt[3]=chanpairnumall_(rv+pairst_);
         xnt[4]=moyreal/(r_4)(nmean_*djf_);
         xnt[5]=moyimag/(r_4)(nmean_*djf_);