Changeset 3658 in Sophya for trunk/AddOn/TAcq/brproc.cc

Timestamp:

Oct 18, 2009, 11:10:33 AM (16 years ago)

Author:

ansari

Message:

1/ Gestion multi-fibre ds RAcqMemZoneMgr et les DMAReader/DiskWriter (brproc.cc)
2/ Possibilite d'ajout de mot cle ds l'entete FITS par la classe MiniFITSFile

Reza , 18/10/2009

File:

• trunk/AddOn/TAcq/brproc.cc (modified) (5 diffs)

trunk/AddOn/TAcq/brproc.cc

@@ -151 +151 @@
     uint_8 firsttt=0;
     bool fgfirst=true;
+    double moysig[2]={0.,0.};
+    double sigsig[2]={0.,0.};
+    uint_8 nbsig[2]={0,0};
+
     for (uint_4 kmz=0; kmz<nmax_; kmz++) {
       if (stop_) break;
@@ -176 +180 @@
             r_4 vts=(r_4)(*(paq.Data1()+j));
             h1.Add((r_8)vts);
+            moysig[0] += (double)vts;
+            sigsig[0] += ((double)vts)*((double)vts);
+            nbsig[0]++;
             vx(j) = vts-127.5;
           }
@@ -201 +208 @@
             r_4 vts=(r_4)(*(paq.Data2()+j));
             h2.Add((r_8)vts);
+            moysig[1] += (double)vts;
+            sigsig[1] += ((double)vts)*((double)vts);
+            nbsig[1]++;
             vx(j) = vts-127.5;
           }
@@ -306 +316 @@
           po << PPFNameTag(tagh1) << h1; 
           po << PPFNameTag(tagh2) << h2; 
-        }
+
+          double sspvmax[3] = {0.,0.,0.};
+          int_4 sspvmaxidx[3] = {-1,-1,-1};
+          for(int jji=1;jji<visiV12.Size()-1;jji++) {
+            r_4 zmv2 = Zmod2(visiV12(jji)); 
+            if (zmv2>sspvmax[2]) { sspvmax[2]=zmv2; sspvmaxidx[2]=jji; }
+          }
+          TVector<r_4>& sspv = spectreV1;
+          for(int ic=0; ic<2; ic++) {
+            if (ic==1) sspv = spectreV2;
+            for(int jji=1;jji<sspv.Size()-1;jji++) 
+              if (sspv(jji)>sspvmax[ic]) { sspvmax[ic]=sspv(jji); sspvmaxidx[ic]=jji; }
+            if (nbsig[ic] < 1) { moysig[ic]=sigsig[ic]=-1.;  }
+            else {
+              moysig[ic] /= (double)nbsig[ic];
+              sigsig[ic] /= (double)nbsig[ic];
+              sigsig[ic] -= (moysig[ic]*moysig[ic]);
+              sigsig[ic] = sqrt(sigsig[ic]);
+              cout << "===Voie " << ic << " Moy=" << moysig[ic] << " Sig=" << sigsig[ic]
+                   << " MaxSpec Amp= " << sqrt(sspvmax[ic])/double(pq.DataSize()/2/2) 
+                   << " Pos=" << sspvmaxidx[ic] << "  (NPts=" << nbsig[ic] << ")" << endl;
+            }
+          }
+          cout << "=== Voie1x2 MaxSpec Amp= " << sqrt(sqrt(sspvmax[2])/double(pq.DataSize()/2/2))
+               << " Pos=" << sspvmaxidx[2] << endl;
+        }  // fin if (fghist_)
+
         if (fgtimfreq) {
           timfreqV1 /= (r_4)nzm;
@@ -320 +356 @@
           h1.Zero();
           h2.Zero();
+          moysig[0]=moysig[1]=0.;
+          sigsig[0]=sigsig[1]=0.;
+          nbsig[0]=nbsig[1]=0;
         }
         if (fgtimfreq) {