Changeset 408 in Sophya for trunk/Poubelle/archTOI.old/archeopsfile.cc

Timestamp:

Sep 22, 1999, 10:57:36 AM (26 years ago)

Author:

ansari

Message:

nouvelles toi

File:

• trunk/Poubelle/archTOI.old/archeopsfile.cc (modified) (1 diff)

trunk/Poubelle/archTOI.old/archeopsfile.cc

@@ -881 +881 @@
 }
 
+double ArcheopsFile::getMuVBolo2(int ibolo, int imesure) { // microvolts, filtre 
+  // On commence par trouver la valeur d'offset, en fittant une droite sur les soustractions
+  if (!lastBolo() || !llastBolo()) return 0;
+  block_type_reglage* reglage = lastReglage();
+  block_type_param*   param = lastParam();
+  if (!reglage) return 0;
+  if (!param) return 0;
+  double sumx=0, sumy=0, sumxy=0, sumx2=0;
+  int n=20;
+  int s = imesure % 2 ? 1 : -1;
+  for (int i=1; i<=n; i++) {
+    double x = -i;
+    double y = s*(getRawBolo(ibolo, imesure-i+1) - getRawBolo(ibolo, imesure-i))/2;
+    s = -s;
+    sumx += x;
+    sumy += y;
+    sumxy += x*y;
+    sumx2 += x*x;
+  }
+  double a = (sumxy/n - (sumx/n)*(sumy/n)) / (sumx2/n - (sumx/n)*(sumx/n));
+  double b = (sumy/n) - a*(sumx/n);
+  s = imesure % 2 ? 1 : -1;
+  double y = getRawBolo(ibolo, imesure);
+  y = y-s*b;
+  return bol_micro_volt(y,(double)param->param.bolo[ibolo].bolo_gain*gain_ampli(reglage->reglage.bolo[ibolo]));
+}
+
 double ArcheopsFile::getGainAmpli(int ibolo) {
      return gain_ampli(lastReglage()->reglage.bolo[ibolo]);
 }
+
+#define capa(uu)        ((parametr.bolo[uu].bolo_bebo==10)?0.000868 * (double)parametr.bolo[uu].bolo_capa:0.001 * (double)parametr.bolo[uu].bolo_capa)          // capa en pF
+#define pt_micA(uu)     ((4096. * 22. * 20.) / capa(uu) )                       //  pts / micro Amperes
+
+
+#define pt_micV(uu)     (((double)parametr.bolo[uu].bolo_diviseur) / 2441.)     // pts / microvolt
+
+
+double ArcheopsFile::getDACV(int ibolo) {
+  double dacV = dac_V(lastReglage()->reglage.bolo[ibolo]);
+  param_bolo& parametr = lastParam()->param;
+  return dacV / pt_micV(ibolo);
+}
+
+double ArcheopsFile::getDACI(int ibolo) {
+  double dacI = dac_I(lastReglage()->reglage.bolo[ibolo]);
+  param_bolo& parametr = lastParam()->param;
+  return dacI / pt_micA(ibolo);
+}
+
 
 // SST, gyros...