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

Timestamp:

Sep 8, 2009, 11:50:59 PM (16 years ago)

Author:

ansari

Message:

Ajout calcul spectre moyenne et sigma das BRProcARaw2C, Reza 08/09/2009

File:

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

trunk/AddOn/TAcq/brproc.cc

@@ -13 +13 @@
 #include "histos.h"
 #include "fioarr.h"
+#include "matharr.h"
 #include "timestamp.h"
 #include "ctimer.h"
@@ -116 +117 @@
     TVector<r_4> spectreV1(cfour1.Size());
     TVector<r_4> spectreV2(cfour1.Size());
+    TVector<r_4> moyspecV1(cfour1.Size());   // Moyenne des Spectres 
+    TVector<r_4> moyspecV2(cfour1.Size());
+    TVector<r_4> sigspecV1(cfour1.Size());   // Sigma des Spectres
+    TVector<r_4> sigspecV2(cfour1.Size());
     TVector< complex<r_4> > visiV12( cfour1.Size() );
 
@@ -130 +135 @@
 
     uint_4 ifile = 0;                           
-    uint_4 nzm = 0;  
+    uint_4 nzm = 0;  // Nb de paquets moyennes pour le calcul de chaque spectre
+    uint_4 nmoyspec = 0;  // Nb de spectres moyennes
 
     uint_4 curfc=0;
@@ -236 +242 @@
         spectreV2 /= (r_4)(nzm);
 
+        // pour le calcul des moyennes et sigmas de ces spectres 
+        moyspecV1 += spectreV1;
+        moyspecV2 += spectreV2;
+        sigspecV1 += (spectreV1 && spectreV1);
+        sigspecV2 += (spectreV2 && spectreV2);
+        nmoyspec++;
+
         visiV12 /= complex<r_4>((r_4)nzm, 0.);
 
@@ -297 +310 @@
   }
 ---- */ 
+  if (nmoyspec>0) {  // Calcul des moyennes et sigmas des spectres 
+    r_4 fnms = nmoyspec;
+    moyspecV1 /= fnms;
+    moyspecV2 /= fnms;
+    sigspecV1 /= fnms;
+    sigspecV2 /= fnms;
+    sigspecV1 -= (moyspecV1 && moyspecV1);
+    sigspecV2 -= (moyspecV2 && moyspecV2);
+    sigspecV1 = Sqrt(sigspecV1);
+    sigspecV2 = Sqrt(sigspecV2);
+    TVector<r_4> rsbV1, rsbV2;   // Rapport signal/bruit
+    moyspecV1.DivElt(sigspecV1, rsbV1, false, true);
+    moyspecV2.DivElt(sigspecV2, rsbV2, false, true);
+    sprintf(fname,"%s_ms.ppf",path_.c_str());
+    POutPersist po(fname);
+    po << PPFNameTag("moyspecV1") << moyspecV1; 
+    po << PPFNameTag("moyspecV2") << moyspecV2; 
+    po << PPFNameTag("sigspecV1") << sigspecV1; 
+    po << PPFNameTag("sigspecV2") << sigspecV2; 
+    po << PPFNameTag("rsbV1") << rsbV1; 
+    po << PPFNameTag("rsbV2") << rsbV2; 
+    cout << " BRProcARaw2C::run() created moysigspec file " << fname << card2name_(card_) << endl;
+  }
+
   ts.SetNow();
   tm.SplitQ();