Changeset 534 for trunk – PI

trunk/Poubelle/archTOI.old/archeopsfile.cc

-              r442
+              r534
 #define utilitaires_de_block_archeops
 #include "archeopsfile.h"
-#include "gyrohandler.h"
 extern "C" {
 …
   block_type_modele         curBlock;
   block_type_modele         peekBlock;
-  SSTHandler   sstHandler;
-  GyroHandler  gyroHandler;
-  GPSParser    gpsParser;
 };
 …
 BlockSet::BlockSet(BlockSet const& x)
-: sstHandler(x.sstHandler), gyroHandler(x.gyroHandler)
+{
   lastParam      = NULL;
 …
+  }
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
+#ifdef __MWERKS__
+#pragma mark -
+#endif
 ArcheopsFile::ArcheopsFile(string const& fname) {
 …
+}
-#define  bitmot                 24              // nb de bit horloge dans un mot ADC
-double ArcheopsFile::perEchant() { // periode d'echantillonage pour le dernier bloc reglage
-  double p,f1,f2,f3;
-  int pp;
-  if (!lastReglage()) return -1;
-  pp=lastReglage()->reglage.horloge.periode;
-  p=pp/5.;
-  f1=1000/p;f2=f1/bitmot;f3=f2*1000./(double)(lastReglage()->reglage.horloge.nb_mesures);
-  return 0.5/f3;                                        //  2 fois la frequence de modulation
+}
-double ArcheopsFile::perBlock() {  // duree (en secondes) correspondant a un bloc bolo
-  return perEchant() * (double)lastParam()->param.n_per_block*2.;
+}
-int ArcheopsFile::nEchBlock() {  // Nb d'echantillons dans un bloc
-  return lastParam()->param.n_per_block*2;
+}
 string  ArcheopsFile::blockKdName() {
 …
+}
-block_type_param* ArcheopsFile::lastParam() {
-  return blockSet->lastParam;
+}
-block_type_journal* ArcheopsFile::lastJournal() {
-  return blockSet->lastJournal;
+}
-block_type_reglage* ArcheopsFile::lastReglage() {
-  return blockSet->lastReglage;
+}
-block_type_dilution* ArcheopsFile::lastDilution() {
-  return blockSet->lastDilution;
+}
-block_type_gps* ArcheopsFile::lastGPS() {
-  return blockSet->lastGPS;
+}
-block_type_une_periode* ArcheopsFile::lastUnePeriode() {
-  return blockSet->lastUnePeriode;
+}
-block_type_synchro_sol* ArcheopsFile::lastSynchroSol() {
-  return blockSet->lastSynchroSol;
+}
-block_type_pointage_sol* ArcheopsFile::lastPointageSol() {
-  return blockSet->lastPointageSol;
+}
-block_type_bolo* ArcheopsFile::lastBolo() {
-  return blockSet->lastBolo;
+}
-block_type_bolo* ArcheopsFile::llastBolo() {
-  return blockSet->llastBolo;
+}
-block_type_gyro* ArcheopsFile::lastGyro() {
-  return blockSet->lastGyro;
+}
-block_type_sst* ArcheopsFile::lastSST() {
-  return blockSet->lastSST;
+}
-block_type_bolo_comprime* ArcheopsFile::lastBoloComp() {
-  return blockSet->lastBoloComp;
+}
-block_type_gyro_comprime* ArcheopsFile::lastGyroComp() {
-  return blockSet->lastGyroComp;
+}
-block_type_sst_comprime* ArcheopsFile::lastSSTComp() {
-  return blockSet->lastSSTComp;
+}
 void ArcheopsFile::forceBlock(block_type_modele* blk) {
   blockSet->setBloc(*blk);
+}
+long ArcheopsFile::searchNextBlock(long pos) {
+  static char* buffer = 0;
+  static int4  debswp = debut_block_mesure;
+  static int4  longmax = taille_maxi_block_archeops*20;
+  if (!buffer) {
+    buffer = new char[longmax];
+#ifdef SWAP
+    bswap4(&debswp);
+#endif
+  }
+  long read = longmax;
+  while (read == longmax) {
+    fseek(f,pos,SEEK_SET);
+    read = fread(buffer,1,longmax,f);
+    //if (read<taille_maxi_block_archeops*2) return -1;
+    // EA 150999 changed i+=4 to i++ -> unaligned, but can lose bytes in flight recorder
+#ifndef __DECCXX
+#define __unaligned
+#endif
+     for (long i=4; i<read; i++) {
+      if (*(__unaligned int4*)(buffer+i) == debswp) {
+        cout << "trying to skip " << i << " bytes to pos="<<pos+i << endl;
+            return pos+i;
+      }
+    }
+    pos += read;
+  }
+  cout << "cannot find block start" << endl;
+  return -1;
+}
 …
        memcpy(&blockSet->curBlock,&blk2,sizeof(blk2));
        curKind = block_bolo;
+       postProcessBlock();
        break;
+     }
      case block_sst_comprime: {
        blockSet->setRawBloc(blockSet->curBlock);
-       if (!blockSet->sstHandler.NeedBlocks()) break; // inutile de decompresser
        block_type_sst blk2;
        block_type_sst_comprime* blk = (block_type_sst_comprime*) &blockSet->curBlock;
 …
        memcpy(&blockSet->curBlock,&blk2,sizeof(blk2));
        curKind = block_sst;
+       //       cout << "process " << numero_block(&blockSet->curBlock) << "\n";
+       blockSet->sstHandler.ProcessBlock((block_type_sst*)&blockSet->curBlock);
+       postProcessBlock();
        break;
+     }
+     case block_sst : {
+       //       cout << "process " << numero_block(&blockSet->curBlock) << "\n";
+       blockSet->sstHandler.ProcessBlock((block_type_sst*)&blockSet->curBlock);
+     }
+     case block_gyro : {
+       blockSet->gyroHandler.ProcessBlock((block_type_gyro*)&blockSet->curBlock);
+     }
+     case block_gps : {
+       blockSet->gpsParser.ProcessBlock((block_type_gps*)&blockSet->curBlock);
+     }
+  }
+}
+// Delegation a tous les producers qui traitent des blocs bruts
+  }
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
+block_type_param* ArcheopsFile::lastParam() {
+  return blockSet->lastParam;
+}
+block_type_journal* ArcheopsFile::lastJournal() {
+  return blockSet->lastJournal;
+}
+block_type_reglage* ArcheopsFile::lastReglage() {
+  return blockSet->lastReglage;
+}
+block_type_dilution* ArcheopsFile::lastDilution() {
+  return blockSet->lastDilution;
+}
+block_type_gps* ArcheopsFile::lastGPS() {
+  return blockSet->lastGPS;
+}
+block_type_une_periode* ArcheopsFile::lastUnePeriode() {
+  return blockSet->lastUnePeriode;
+}
+block_type_synchro_sol* ArcheopsFile::lastSynchroSol() {
+  return blockSet->lastSynchroSol;
+}
+block_type_pointage_sol* ArcheopsFile::lastPointageSol() {
+  return blockSet->lastPointageSol;
+}
+block_type_bolo* ArcheopsFile::lastBolo() {
+  return blockSet->lastBolo;
+}
+block_type_bolo* ArcheopsFile::llastBolo() {
+  return blockSet->llastBolo;
+}
+block_type_gyro* ArcheopsFile::lastGyro() {
+  return blockSet->lastGyro;
+}
+block_type_sst* ArcheopsFile::lastSST() {
+  return blockSet->lastSST;
+}
+block_type_bolo_comprime* ArcheopsFile::lastBoloComp() {
+  return blockSet->lastBoloComp;
+}
+block_type_gyro_comprime* ArcheopsFile::lastGyroComp() {
+  return blockSet->lastGyroComp;
+}
+block_type_sst_comprime* ArcheopsFile::lastSSTComp() {
+  return blockSet->lastSSTComp;
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
 void ArcheopsFile::saveCurBlock() {
 …
 #endif
+long ArcheopsFile::searchNextBlock(long pos) {
+  static char* buffer = 0;
+  static int4  debswp = debut_block_mesure;
+  static int4  longmax = taille_maxi_block_archeops*20;
+  if (!buffer) {
+    buffer = new char[longmax];
+#ifdef SWAP
+    bswap4(&debswp);
+#endif
+  }
+  size_t read = longmax;
+  while (read == longmax) {
+    fseek(f,pos,SEEK_SET);
+    read = fread(buffer,1,longmax,f);
+    //if (read<taille_maxi_block_archeops*2) return -1;
+    // EA 150999 changed i+=4 to i++ -> unaligned, but can lose bytes in flight recorder
+#ifndef __DECCXX
+#define __unaligned
+#endif
+     for (size_t i=4; i<read; i++) {
+      if (*(__unaligned int4*)(buffer+i) == debswp) {
+        cout << "trying to skip " << i << " bytes to pos="<<pos+i << endl;
+        return pos+i;
+      }
+    }
+    pos += read;
+  }
+  cout << "cannot find block start" << endl;
+  return -1;
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
 …
+}
 void ArcheopsFile::setUTCOffset(int UTCOffset) {
+void ArcheopsFile::setUTCOffset(double UTCOffset) {
   utcOffset = UTCOffset;
   startMJD = rawMJD - utcOffset/24.;
 …
+}
-int ArcheopsFile::getBoloBlockNum() {
-  if (!lastBolo()) return -1;
-  return numero_block(lastBolo());
+}
-int ArcheopsFile::getBoloCompBlockNum() {
-  if (!lastBolo()) return -1;
-  return numero_block(lastBoloComp());
+}
-int ArcheopsFile::getReglageBlockNum() {
-  if (!lastReglage()) return -1;
-  return numero_block(lastReglage());
+}
-int ArcheopsFile::getDilutionBlockNum() {
-  if (!lastDilution()) return -1;
-  return numero_block(lastDilution());
+}
-int ArcheopsFile::getSSTBlockNum() {
-  if (!lastSST()) return -1;
-  return numero_block(lastSST());
+}
-int ArcheopsFile::getSSTCompBlockNum() {
-  if (!lastSST()) return -1;
-  return numero_block(lastSSTComp());
+}
-int ArcheopsFile::getGyroBlockNum() {
-  if (!lastGyro()) return -1;
-  return numero_block(lastGyro());
+}
-// GPS
-// $GPGGA,hhmmss.ss,ddmm.mmmm,n,dddmm.mmmm,e,q,ss,y.y,a.a,z,
-int ArcheopsFile::getGPSBlockNum() {
-  if (!lastGPS()) return -1;
-  return numero_block(lastGPS());
+}
-double ArcheopsFile::getGPSUTC() {  // en secondes depuis minuit UTC jour courant
-  if (!lastGPS()) return -9.e99;
-  return blockSet->gpsParser.getUTC();
+}
-double ArcheopsFile::getGPSMJD() {  // modified julian day du dernier bloc GPS, JD - 2450000
-  double t = getGPSUTC()/86400. - 0.5;
-  if (t<0) return t;
-  if (t > (startMJD - int(startMJD))) {
-    return int(startMJD) + t;
-  } else {
-    return int(startMJD) + 1. + t;
+  }
+}
-double ArcheopsFile::getGPSLat() {  // degres, +  = NORD
-  if (!lastGPS()) return -9.e99;
-  return blockSet->gpsParser.getLatitude();
+}
-double ArcheopsFile::getGPSLong() { // degres, +  = EST
-  if (!lastGPS()) return -9.e99;
-  return blockSet->gpsParser.getLongitude();
+}
-double ArcheopsFile::getGPSAlt() { // meters from sea level
-  if (!lastGPS()) return -9.e99;
-  return blockSet->gpsParser.getAltitude();
+}
-bool ArcheopsFile::hasGPSTime() {
-  if (!lastGPS()) return false;
- // return blockSet->gpsParser.hasGPSTime();
-  return blockSet->gpsParser.hasTime();
+}
-bool ArcheopsFile::hasGPSPos() {
-  if (!lastGPS()) return false;
-  return blockSet->gpsParser.hasPosition();
+}
-bool ArcheopsFile::hasGPSAlt() {
-  if (!lastGPS()) return false;
-  return blockSet->gpsParser.hasAltitude();
+}
-// Bolo
-#define  val_DS(j,i)  (blk->data_bolo[j][i]&0x1fffff)
-double ArcheopsFile::getBoloRawMuV(int ibolo, int imesure) {    // microvolts bruts
-  block_type_bolo* blk = imesure >= 0 ? lastBolo() : llastBolo();
-  if (imesure < 0) imesure += nEchBlock();
-  block_type_reglage* reglage = lastReglage();
-  block_type_param*   param = lastParam();
-  if (!blk) return 0;
-  if (!reglage) return 0;
-  if (!param) return 0;
-  int s = imesure % 2 ? 1 : -1;
-  return s*bolo_muV(&param->param, &reglage->reglage, val_DS(ibolo,imesure), ibolo);
+}
-double ArcheopsFile::getBoloRawMuVCN(int ibolo, int imesure) {
-  // Si pas bloc comprime -> 1 bit
-  int lastbit = 1;
-  if (!(lastBoloComp() == NULL) &&
-      (numero_block(lastBoloComp()) != numero_block(lastBolo()))) {
-    unsigned int4* data = lastBoloComp()->data_bolo[ibolo];
-    // Les deux premieres valeurs sont codees directement...
-    if (imesure>=2) {
-      int iExp = (imesure-2)/4 + 1;
-      int expo = data[iExp] & 0xf;
-      lastbit = 1 << expo;
+    }
+  }
-  block_type_reglage* reglage = lastReglage();
-  block_type_param*   param = lastParam();
-  if (!reglage) return 0;
-  if (!param) return 0;
-  double noise = (
-    bolo_muV(&param->param, &reglage->reglage, lastbit, ibolo) -
-    bolo_muV(&param->param, &reglage->reglage, 0, ibolo)) /2.;
-  return noise;
+}
-def_gains
-double ArcheopsFile::getBoloMuV(int ibolo, int imesure) { // microvolts, filtre avec filtre carre
-  double y = (getBoloRawMuV(ibolo, imesure-1) + getBoloRawMuV(ibolo, imesure))/2.;
-  return y;
+}
-double ArcheopsFile::getBoloMuV2(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*(getBoloRawMuV(ibolo, imesure-i+1) - getBoloRawMuV(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 = getBoloRawMuV(ibolo, imesure);
-  y = y-s*b;
-  return y;
+}
-/*
-double ArcheopsFile::getBoloMuV2(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;
-  double sumx=0, sumy=0, sumxy=0, sumx2=0;
-  int n=20;
-  int s = (imesure-n) % 2 ? 1 : -1;
-  static double* data = NULL;
-  if (!data) data = new double[n+1];
-  for (int i=0; i<=n; i++) {
-    data[i] = getBoloRawMuV(ibolo, imesure-n+i);
+  }
-  for (int i=0; i<=n; i++) {
+  }
-  for (int i=1; i<=n; i++) {
-    double x = -i;
-    double y = s*(getBoloRawMuV(ibolo, imesure-i+1) - getBoloRawMuV(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 = getBoloRawMuV(ibolo, imesure);
-  y = y-s*b;
-  return y;
+}
-*/
-double ArcheopsFile::getGainAmpli(int ibolo) {
-     return gain_ampli(lastReglage()->reglage.bolo[ibolo]);
+}
-double ArcheopsFile::getDACV(int ibolo) {
-  block_type_reglage* reglage = lastReglage();
-  block_type_param*   param = lastParam();
-  if (!reglage) return 0;
-  if (!param) return 0;
-  return DAC_muV(&param->param, &reglage->reglage, ibolo);
+}
-double ArcheopsFile::getDACI(int ibolo) {
-  block_type_reglage* reglage = lastReglage();
-  block_type_param*   param = lastParam();
-  if (!reglage) return 0;
-  if (!param) return 0;
-  return DAC_muA(&param->param, &reglage->reglage, ibolo);
+}
-double ArcheopsFile::getBoloMuV2T(int ibolo, int imesure) {
-  return getDACV(ibolo) - getBoloMuV2(ibolo, imesure) ;
+}
-double ArcheopsFile::getBoloRes(int ibolo, int imesure) {
-  double i = getDACI(ibolo); // microAmps
-  double v = getBoloMuV2T(ibolo, imesure); // microVolts
-  double r = v/i; // Ohms
-  return r;
+}
-double ArcheopsFile::getBoloTemp(int ibolo, int imesure) {
-  double r = getBoloRes(ibolo, imesure);
-  block_type_reglage* reglage = lastReglage();
-  block_type_param*   param = lastParam();
-  if (!reglage) return 0;
-  if (!param) return 0;
-  return bolo_temp(&param->param, &reglage->reglage, r, ibolo);
+}
-// Dilution
-int4 ArcheopsFile::getADCDil(int iADC) {
-  block_type_dilution* blk = lastDilution();
-  if (!blk) return 0;
-  if (iADC < 0 || iADC > 47) return 0;
-  return blk->ADC_dil[iADC];
+}
-int4 ArcheopsFile::getSwitchDil() {
-  block_type_dilution* blk = lastDilution();
-  if (!blk) return 0;
-  return blk->switch_dil;
+}
-// SST, gyros...
-void ArcheopsFile::needSSTProcessMask(int mask) {
-  blockSet->sstHandler.NeedProcess(mask);
+}
-long ArcheopsFile::getSSTSignal(int idiode, int imesure) {
-  return blockSet->sstHandler.getSignal(imesure, idiode);
+}
-long ArcheopsFile::getSSTRawSignal(int idiode, int imesure) {
-  return blockSet->sstHandler.getRawSignal(imesure, idiode);
+}
-double ArcheopsFile::getSSTSignalCN(int idiode, int imesure) {
-  return getSSTRawSignalCN(SSTHandler::getDiodPermut(idiode),imesure);
+}
-double ArcheopsFile::getSSTRawSignalCN(int ichannel, int imesure) {
-  // Si pas bloc comprime -> 0.5
-  if (lastSSTComp() == NULL) return .5;
-  if (numero_block(lastSSTComp()) != numero_block(lastSST())) return 0.5;
-  // Attention, on ne transmet pas les canaux 0 et 4....
-  if (ichannel == 0 || ichannel == 4) return 0;
-  int i = ichannel - 1;
-  if (i >= 4) i--;
-  unsigned int4* data = lastSSTComp()->sst[i];
-  // Les deux premieres valeurs sont codees directement...
-  if (imesure<2) return 0.5;
-  int iExp = (imesure-2)/7 + 1;
-  int expo = data[iExp] & 0xf;
-  int noise = 1 << expo;
-  return noise/2.;
+}
-double ArcheopsFile::getSSTStarZ(int istar, int imesure) {
-  return blockSet->sstHandler.getStarZ(imesure+getSSTBlockNum()*nb_per_block*2, istar);
+}
-double ArcheopsFile::getSSTStarF(int istar, int imesure) {
-  return blockSet->sstHandler.getStarF(imesure+getSSTBlockNum()*nb_per_block*2, istar);
+}
-double ArcheopsFile::getSSTStarT(int istar, int imesure) {
-  return blockSet->sstHandler.getStarTime(imesure+getSSTBlockNum()*nb_per_block*2, istar);
+}
-int ArcheopsFile::getNumbStar(int imesure) {
-  return blockSet->sstHandler.getNumbStar(imesure+getSSTBlockNum()*nb_per_block*2);
+}
-long ArcheopsFile::getGyroRaw(int igyro, int imesure) {
-  return blockSet->gyroHandler.getRawSignal(imesure, igyro);
+}
-double ArcheopsFile::getGyroTens(int igyro, int imesure) {
-  return blockSet->gyroHandler.getSignal(imesure, igyro);
+}
-double ArcheopsFile::getGyroSpeed(int igyro, int imesure) {
-  return blockSet->gyroHandler.getSpeed(imesure, igyro);
+}
-// $CHECK$ TBD
-double ArcheopsFile::getAzimut(int /*imesure*/) {return 0;}
-double ArcheopsFile::getPendDirect(int /*imesure*/) {return 0;}
-double ArcheopsFile::getPendOrth(int /*imesure*/) {return 0;}
-double ArcheopsFile::getAlphaAxis(int /*imesure*/) {return 0;}
-double ArcheopsFile::getDeltaAxis(int /*imesure*/) {return 0;}

trunk/Poubelle/archTOI.old/archeopsfile.h

-              r436
+              r534
 #include <stack>
+using namespace std;
+#ifndef nb_max_bolo
+#define _archeops               //  Acquisition Archeops (avec transputer)
+#define programme
+extern "C" {
+#include "archeops.h"
+}
+#endif
+#include "ssthandler.h"
+#include "gpsparser.h"
+#include "ark.h"
+#include "archexc.h"
 enum {
 …
 class BlockSet;
-class ArchExc {
-public:
-  ArchExc(string s) : msg(s) {}
-  string Msg() {return msg;}
-protected:
-  string msg;
-};
 class ArcheopsFile {
 …
   virtual ~ArcheopsFile();
   void setUTCOffset(int UTCOffset=2); // StartMJD : decodage a partir du nom du bloc...
+  void setUTCOffset(double UTCOffset=2); // StartMJD : decodage a partir du nom du bloc...
   double getStartMJD();
 …
   string blockRawKdName();
   string filename() {return fn;}
-  double perEchant(); // periode d'echantillonage pour le dernier bloc reglage
-  double perBlock();  // duree (en secondes) correspondant a un bloc bolo
-                      // utilise infos dernier bloc param et dernier bloc reglage.
-  int nEchBlock(); // Nb d'echantillons dans un bloc bolo
 …
   void forceBlock(block_type_modele* blk);
-  // Decodage des derniers blocs..
-  // GPS
-  // On travaille en MJD = JD - 2450000.0, = 9 octobre 1995 midi UTC
-  // 1er juillet 1999 minuit UTC = JD 2451360.5,   MJD 1360.5
-  // Ces fonctions travaillent sur le dernier bloc GPS
-  int getGPSBlockNum();
-  bool   hasGPSTime();
-  bool   hasGPSPos();
-  bool   hasGPSAlt();
-  double getGPSUTC();  // en secondes depuis minuit UTC
-  double getGPSMJD();  // modified julian day du dernier bloc GPS, JD - 2450000
-  double getGPSLat();  // degres, +  = NORD
-  double getGPSLong(); // degres, +  = EST
-  double getGPSAlt();  // meters from sea level
-  // Bolo
-  int getReglageBlockNum();
-  int getBoloBlockNum();
-  int getBoloCompBlockNum();
-  double getBoloRawMuV(int ibolo, int imesure);   // donnee brute en microvolts et signe change.
-  double getBoloRawMuVCN(int ibolo, int imesure);  // bruit de compression, en microVolts
-  double getBoloMuV(int ibolo, int imesure); // microvolts, filtre avec filtre carre
-  double getBoloMuV2(int ibolo, int imesure); // microvolts, filtre +sioux
-  double getGainAmpli(int ibolo);
-  double getDACV(int ibolo); // offset en microvolt
-  double getDACI(int ibolo); // courant en microAmps
-  double getBoloMuV2T(int ibolo, int imesure); // avec soustraction dacV
-  double getBoloRes(int ibolo, int imesure); // Ohms
-  double getBoloTemp(int ibolo, int imesure); // K
-  // Dilution
-  int getDilutionBlockNum();
-  int4 getADCDil(int iADC);
-  int4 getSwitchDil();
-  // SST, gyros, etc
-  int getSSTBlockNum();
-  int getSSTCompBlockNum();
-  void needSSTProcessMask(int mask);
-  long getSSTRawSignal(int iDiode, int imesure);
-  long getSSTSignal(int iDiode, int imesure);
-  double getSSTRawSignalCN(int iDiode, int imesure);
-  double getSSTSignalCN(int iDiode, int imesure);
-  double getSSTStarZ(int iStar, int imesure);
-  double getSSTStarF(int iStar, int imesure);
-  double getSSTStarT(int iStar, int imesure);
-  int getNumbStar(int imesure);
-  int getGyroBlockNum();
-  long getGyroRaw(int iGyro, int imesure);
-  double getGyroTens(int iGyro, int imesure);
-  double getGyroSpeed(int iGyro, int imesure);
-  double getAzimut(int imesure);
-  double getPendDirect(int imesure);
-  double getPendOrth(int imesure);
-  double getAlphaAxis(int imesure);
-  double getDeltaAxis(int imesure);
   static double decodeMJD(string const& dateString); // sans tenir compte de TU-TLeg
 …
   BlockSet*                 blockSet;
   stack<BlockSet*>          blockStack;
   stack<size_t>             fposStack;
+  stack<long>               fposStack;
   void saveCurBlock();
 …
   long                      curPos;
   long                      peekPos;
   size_t                    fLen;
+  long                      fLen;
   FILE* f;
   string fn;

trunk/Poubelle/archTOI.old/archparam.cc

-              r432
+              r534
   cout << "***END" << endl;
+}
+double ArchParam::AcqParam::SN2MJD(double sample) {
+  return tBlock0 + sample*perEch/86400.;
+}
+double ArchParam::AcqParam::MJD2UTC(double mjd) {
+  return (mjd - utcOrigin)*24. ;
+}
+double ArchParam::AcqParam::UTC2MJD(double utc) {
+  return utc/24 + utcOrigin;
+}
+double ArchParam::AcqParam::SN2UTC(double sample) {
+  return MJD2UTC(SN2MJD(sample));
+}
+double ArchParam::AcqParam::MJD2SN(double mjd) {
+  return (mjd - tBlock0)/perEch*86400.;
+}

trunk/Poubelle/archTOI.old/archparam.h

-              r432
+              r534
     double utcOrigin;
     AcqParam();
+    double SN2MJD(double snum);
+    double SN2UTC(double snum);
+    double MJD2UTC(double mjd);
+    double UTC2MJD(double utc);
+    double MJD2SN(double mjd);
   };

trunk/Poubelle/archTOI.old/archtoi.cc

-              r436
+              r534
 #include "archeopsfile.h"
+#include "toisvr.h"
+#include "toiiter.h"
+#include "toimanager.h"
 #include "archtoi.h"
 #include "archparam.h"
 …
+{
   init();
   svr.readReq(str);
+  iter.readReq(str);
+}
 …
   init();
   ifstream str(filename.c_str());
   svr.readReq(str);
+  iter.readReq(str);
+}
 …
   ostr = NULL;
   requestVersion = "";
+  svr.registerReqHandler(this);
+}
+bool ArchTOI::processTOIReq(string line, string toiname, TOIKind /*kind*/, int index,
+                            bool interp, bool repet, bool flag, bool /*notrig*/)
+  iter.registerReqHandler(this);
+  TOIManager::registerDefaultProducers();
+}
+bool ArchTOI::processTOIReq(TOI const& toi, string line)
+{
   headertoi.push_back(line);
+  if (index>=0) {
+    char idx[10];
+    sprintf(idx,"_%d",index);
+    toiname += idx;
+  }
+  if (flag) {
+    toinames.push_back("flg_"+toiname);
+  }
+  toinames.push_back(toiname);
+  toiflags.push_back(flg((flag?hasflag:0)+((!repet&&!interp)?useNA:0)));
+  tois.push_back(toi);
   return true;
+}
 …
   cout << "starting query" << endl;
   TOIIter iter = svr.doQuery();
+  iter.init();
   (this->*openFile)(outfilename);
   (this->*outHeader)(iter);
   cout << "processing" << endl;
   while (iter.Next()) {
     int nn = iter.getSampleIndex() / iter.getUnderSampling(); // Only for tick mark...
+  while (iter.next()) {
+    int nn = iter.getSampleNum() / iter.getUnderSample(); // Only for tick mark...
     #ifdef __MWERKS__
     yield();
 …
      if (nn%(200*80) == 0) cout << endl;
     // Si rien de dispo parmi les triggering, alors on passe au suivant
+    bool hasValue = false;
+    for (int i=0; i<toiflags.size(); i++) {
+       if (!iter.isTrig(i)) continue;
+       if (iter.canGetValue(i)) {hasValue=true; break;}
+    }
+    if (!hasValue) continue;
+    // normalement, c'est gere par le toiiter...
+    //bool hasValue = false;
+    //for (int i=0; i<toiflags.size(); i++) {
+    //   if (!iter.isTrig(i)) continue;
+    //   if (iter.canGetValue(i)) {hasValue=true; break;}
+    //}
+    //if (!hasValue) continue;
     int icol=0;
     for (int i=0; i<toiflags.size(); i++) {
+    for (int i=0; i<tois.size(); i++) {
       double value = iter.getValue(i);
       bool   ok    = iter.canGetValue(i);
+      bool   isnew = iter.newValue(i);
+      flg    flag  = toiflags[i];
+      if (flag & hasflag) {
+        (this->*outValue)(icol,  (ok && isnew ? 1 : 0));
+        icol++;
+      }
+      if (((flag & useNA)!=0 && !isnew) || !ok)
+     // bool   isnew = iter.newValue(i);
+     // flg    flag  = toiflags[i];
+     // if (flag & hasflag) {
+     //   (this->*outValue)(icol,  (ok && isnew ? 1 : 0));
+     //   icol++;
+     // }
+     //if (((flag & useNA)!=0 && !isnew) || !ok)
+      if (!ok)
         (this->*outValue)(icol,  0, true);
       else
 …
   for (list<string>::iterator i = headeropt.begin(); i != headeropt.end(); i++)
     *ostr << (*i) << '\n';;
   block_type_param* blk = iter.lastParam();
+  block_type_param* blk = iter.getFSet().lastParam();
   if (blk) {
     int nb = blk->param.n_max_bolo;
 …
+}
+list<string> ArchTOI::buildNames() {
+  list<string> toinames;
+  map<string,int> namecount;
+  for (list<TOI>::iterator i = tois.begin(); i != tois.end(); i++) {
+    namecount[(*i).name]++;
+  }
+  map<string,int> curcount;
+  for (list<TOI>::iterator i = tois.begin(); i != tois.end(); i++) {
+    string toiname = (*i).name;
+    if (namecount[(*i).name]>1) {
+      char s[10];
+      sprintf(s,"_%d",curcount[(*i).name]++);
+      toiname += s;
+    }
+    if ((*i).index != TOI::unspec) {
+      char s[10];
+      sprintf(s,"_%d",(*i).index);
+      toiname += s;
+    }
+    toinames.push_back(toiname);
+  }
+  return toinames;
+}
 void ArchTOI::openFile_F(string const& filename) {
 …
+}
 void ArchTOI::outHeader_F(TOIIter& iter) {
+  list<string> toinames = buildNames();
   int ncols=toinames.size();
 …
   char** colunits = new (char*[ncols]);
   int j=0;
+  for (list<string>::iterator i = toinames.begin(); i != toinames.end(); i++,j++) {
+  list<TOI>::iterator kk = tois.begin();
+  for (list<string>::iterator i = toinames.begin(); i != toinames.end(); i++,j++,kk++) {
     colnames[j] = const_cast<char*>((*i).c_str()); // should work for most STL implementations... Check...
     coltypes[j] = "1D";
+    colunits[j] = " ";
+    colunits[j] = const_cast<char*>((*kk).unit.c_str());
+    cout << (*i).c_str() << " " << ((*kk).unit.c_str()) << endl;
+  }
 …
   int ntoireq = headertoi.size();
   fits_write_key(fptr, TINT, "TOIREQ", &ntoireq, NULL, &fitsStatus);
+  for (list<string>::iterator i = headertoi.begin(); i != headertoi.end(); i++,j++) {
+  list<string>::iterator kkk = toinames.begin();
+  for (list<string>::iterator i = headertoi.begin(); i != headertoi.end(); i++,j++,kkk++) {
     char line[80];
     strcpy(line, (*i).c_str());
     char* pline = line;
+    fits_write_keys_str(fptr, "TOIREQ", j, 1, &pline, (char**) NULL, &fitsStatus);
+    char* cmt   = const_cast<char*>((*kkk).c_str());
+    fits_write_keys_str(fptr, "TOIREQ", j, 1, &pline, &cmt, &fitsStatus);
+  }
   j=1;
 …
   // Noms des bolos
   block_type_param* blk = iter.lastParam();
+  block_type_param* blk = iter.getFSet().lastParam();
   if (blk) {
     int nb = blk->param.n_max_bolo;
 …
   fits_report_error(stderr, fitsStatus);  /* print out any error messages */
+}

trunk/Poubelle/archTOI.old/archtoi.h

-              r432
+              r534
 #include <list>
 #include <vector>
 #include "toisvr.h"
+#include "toiiter.h"
 #include "fitsio.h"
 …
             useNA   = 2};
+  virtual bool processTOIReq(string line, string toiname, TOIKind toi, int index,
+                             bool interp, bool repet, bool flag, bool notrig);
+  virtual bool processTOIReq(TOI const& toi, string line);
   virtual bool processOption(string keyw, string args);
 protected:
   void init();
   TOISvr svr;
+  TOIIter iter;
   list<string> headertoi;
   list<string> headeropt;
+  list<string> toinames;
+  vector<flg>  toiflags;
+  list<TOI> tois;
   fmt format;
   string undef;
 …
   void (ArchTOI::*endLine)();
   void (ArchTOI::*closeFile)();
+  list<string> buildNames();
   void openFile_A(string const& filename);

trunk/Poubelle/archTOI.old/archtoimain.cc

-              r436
+              r534
 #include "archeopsfile.h"
 #include "archtoi.h"
 #include "toisvr.h"
+#include "toimanager.h"
 //#include <profiler.h>
 …
     argc = ccommand(&argv);
 #endif
     if (argc != 3) {
       cerr << "usage: archtoi reqfile datafile" << endl;
+      cerr << "or:    archtoi -h reqfile" << endl;
       return(-1);
+    }
 …
 //   ProfilerInit(collectDetailed, bestTimeBase, 10000, 200);
+    if (argv[1] == string("-h")) {
+      ArchTOI toi(argv[2]);
+      TOIManager::dumpAvailTOIs(cout);
+      exit(0);
+    }
+    ArchTOI toi(argv[1]);
+    toi.run(argv[2]);
+    try {
+      ArchTOI toi(argv[1]);
+      toi.run(argv[2]);
+    } catch (ArchExc exc) {
+      cerr << exc.Msg() << endl;
+      exit(-1);
+    }
 //    ProfilerDump("\pprofiler.data");
     return(0);
+}

trunk/Poubelle/archTOI.old/archvers.h

-              r503
+              r534
 #define ARCHVERS_H
 #define ARCHTOI_VERS "1.3"
 #define ARCHTOI_TAG  "V_241099"
+#define ARCHTOI_VERS "2.0"
+#define ARCHTOI_TAG  "V_011199"
 #endif

trunk/Poubelle/archTOI.old/auxgps.h

-              r363
+              r534
   virtual int getLocation(double mjd, double& lat, double& lon, double& alt)=0;
+  virtual bool canGetLocation(double mjd)=0;
+  virtual double getStartMJD()=0;
+  virtual double getEndMJD()=0;
   virtual AuxGPS* clone()=0;

trunk/Poubelle/archTOI.old/auxinterpgps.cc

-              r416
+              r534
   ilast = elems.upper_bound(elem(mjd,0,0,0));
   ilastok = (ilast != elems.end() && ilast != elems.begin());
+}
+bool AuxInterpGPS::canGetLocation(double mjd) {
+  return (mjd > (*elems.begin()).mjd && mjd < (*--elems.end()).mjd);
+}
+double AuxInterpGPS::getStartMJD() {
+  return (*elems.begin()).mjd;
+}
+double AuxInterpGPS::getEndMJD() {
+  return (*--elems.end()).mjd;
+}

trunk/Poubelle/archTOI.old/auxinterpgps.h

-              r363
+              r534
   AuxInterpGPS();
   virtual int getLocation(double mjd, double& lat, double& lon, double& alt);
+  virtual bool canGetLocation(double mjd);
+  virtual double getStartMJD();
+  virtual double getEndMJD();
   void FitsDump(string const& filename);

trunk/Poubelle/archTOI.old/dyffttools.cc

r394	r534
1	1	#include "dyffttools.h"
2	2
	3	// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
3	4
4	5	// Librairie FFTMayer enrobée pour etre utilisable facilement en C++

trunk/Poubelle/archTOI.old/dyffttools.h

-              r394
+              r534
 #ifndef DYFFTTOOLS_SEEN
+#define DYFFTTOOLS
+#define DYFFTTOOLS_SEEN
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 double DYWindow(double index, double indexMax);

trunk/Poubelle/archTOI.old/formepulse.cc

-              r394
+              r534
 #include <fstream.h>
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 #include <math.h>
 …
 #include "ssthardware.h"
 #define LTab (2048)
 #define LSave (128)
+#define LSave (256)
 #define Tpas (5.e-4)
 //float PulseOptOrigin=LTab/16.*Tpas;
 …
 //#define DYTestRun
+/*
+FormePulse::FormePulse() {
+        gainElec=0.;
+        FullwidthHalfMax=0;
+        ShapeArr=new double[LSave];
+        for(int i=0; i<LSave; i++)  ShapeArr[i]=0.;
+}
+*/
 FormePulse::FormePulse(TransFuncElec TF) {
 …
         for(int i=0; i<LSave; i++)  ShapeArr[i]=PulseInReel[i]/Max;
 #ifdef DYTestRun
+//#ifdef DYTestRun
                 //On ecrit le resultat dans un fichier
                 char s2[32]="Resultatfiltrage.txt";
 …
                 dlout<<"Impulsion courant de: "<<SSTFPulseCourant<<" Ampères"<<endl;
+                dlout<<"Maximum d'amplitude: "<<max<< " Volts"<<endl;
+                dlout<<" gain chaine preampli: "<<gainElec<<endl<<endl;
+                dlout<<"Maximum d'amplitude: "<<Max<< " Volts"<<endl;
+                dlout<<" gain chaine preampli: "<<gainElec<<endl;
+                dlout<<"largueur a mi hauteur: "<<FullwidthHalfMax<<endl<<endl;
                 dlout<<"Temps"<<"\tPartie Reelle"<<endl;
 …
                         dlout<<k*Tpas<<"\t"<<ShapeArr[k]<<endl;
+                }
 #endif
+//#endif
 …
         return *this;
+}
 FormePulse::~FormePulse() {
         delete[] ShapeArr;
 …
 double FormePulse::PulseShape(double temps) {
+        double returnValue;
         int ientier=(int)(temps/Tpas);
-        double returnValue;
         if (ientier<0) returnValue=0.;
         else if (ientier>=LSave) returnValue= ShapeArr[127];
+        else if (ientier>=LSave) returnValue= ShapeArr[LSave-1];
         else {
                 double Min=ShapeArr[ientier];
 …
         //en secondes
         return returnVal;
+//      return returnVal=30.e-3;                        //Gros BUGG!!!!!
+}

trunk/Poubelle/archTOI.old/formepulse.h

-              r394
+              r534
 #define FORMEPULSE_SEEN
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 #define  SSTFPulseCourant 1.E-9
 …
 class FormePulse {
 public:
+        FormePulse(TransFuncElec);
+        FormePulse(TransFuncElec TF);
+        FormePulse();                                           // Cas zero lecture
         FormePulse& operator =(const FormePulse&);
         ~FormePulse();

trunk/Poubelle/archTOI.old/pisteetoile.cc

-              r400
+              r534
 FormePulse    PisteEtoile::FPulse(Trapani);
+int SSTnbPasFit= (int) (2*TFExcursion/TPasTFit);
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 PisteEtoile::PisteEtoile(short NoPhotDiod) {
         PasTempsEch=0.005836818076;
+        PasTempsEch=archParam.acq.perEch;
                 //BUGG doit changer     !!!!! Doit etre extrait d'un serveur de param d'acquisition
         NoPhDiod=NoPhotDiod;
         DiodVivante=true;
+        SSTnbPasFit= (int) (2*TFExcursion/TPasTFit);
         AmplNorm=0.;
-//      On echantillonne dix fois plus fin pour avoir une bonne estimation de AmplNorm
-//      Eviter les systematiques
         for(int i=0; i<SSTLongIndexEvt()*10; i++){
                 AmplNorm+=FPulse.PulseShape(i/10.*TPasTFit)*FPulse.PulseShape(i/10.*TPasTFit);
+                AmplNorm+=FPulse.PulseShape(i*PasTempsEch/10.)*FPulse.PulseShape(i*PasTempsEch/10.);
+        }
         AmplNorm/=10.;
 …
 PisteEtoile& PisteEtoile::operator =(PisteEtoile const & x) {
         PasTempsEch=x.PasTempsEch;
-                //BUGG doit changer     !!!!! Doit etre extrait d'un serveur de param d'acquisition
         NoPhDiod=x.NoPhDiod;
         DiodVivante=x.DiodVivante;
 …
 void PisteEtoile::fill(int* pData, int InDebutPiste, int nData) {
 //      if(!isAlive()) return;
         indexDebutPiste=InDebutPiste;
         int* PhDiod=PhDiodArray;
+        indexDebutPiste=InDebutPiste - PhDiodTabLong + nData;
+        int* PhDiod=PhDiodArray + PhDiodTabLong - nData;
         int* pDataC=pData;
         for(int k=0;k<nData;k++,pDataC++, PhDiod++) (*PhDiod)=(*pDataC);
 …
                 if(testEtoile) {
+                // On emballe une etoile.
+                // On emballe une etoile
+                        LastEtoile.InpCurrent= AmplRes/FPulse.GainElectrique();
+                        LastEtoile.InpCurFin= AmplFin/FPulse.GainElectrique();
+                        LastEtoile.TEchan=TIndex;
+                        LastEtoile.TFin=Tfin;
+                        LastEtoile.InpCurrNoisRMS=Vrms/FPulse.GainElectrique();
+                        LastEtoile.X2Calc=X2Min;
+                        LastEtoile.NoDiode=NoPhDiod;
                         LastEtoile.FitForme=SuccesAFine;
-                        LastEtoile.X2Calc=X2Min;
-                        LastEtoile.InpCurrNoisRMS=Vrms/FPulse.GainElectrique();
-                        LastEtoile.TEchan=TIndex;
-                        LastEtoile.InpCurrent= AmplRes/FPulse.GainElectrique();
-                        LastEtoile.NoDiode=NoPhDiod;
                         LastEtoile.LargMiHauteur=widthHalfMax;
+/*                      cout<<setprecision(9);
+                        cout<<TIndex<<" "<<Tfin<<" "<<AmplRes<<" "<<AmplFin<<endl;
+                        cout<<TIndex-Tfin<<" "<<(AmplRes-AmplFin)/AmplFin<<endl;
+*/
+                }
+        }
 …
 // On calcule le Vrms sur le prepulse
         Vrms=0.;
+        for (int i=0; i<Prepulselong; i++) Vrms+=PhDiodPiste[i]*PhDiodPiste[i];
+        Vrms=sqrt(Vrms)/Prepulselong;
+        for (int i=0; i<Prepulselong; i++) {
+                Vrms+=PhDiodPiste[i]*PhDiodPiste[i];
+        }
+        Vrms=sqrt(Vrms/Prepulselong);
         if(Vrms<SSTLsbVal) Vrms=SSTLsbVal;              // Pour eviter les VRMS=0.
 …
 // Maintenant que l'on se doute que il y a une etoile
 // On veut connaitre son flux et son temps de passage
 // Nécessite d'avoir un ordre d'idee de la vitesse de rotation du ballon.
+// Nécessite de connaitre la vitesse de rotation du ballon.
                 // BUGG devra etre extrait des headers acquisition. Ca peut changer!!!!!!
         int IndexDebutEvt=MaxIndex-SSTLargEvtIndex();
 …
                 // Pointeur origine du tableau d'echantillon sur lequel sera calculé le fit
         double DeltaTemps=0.;
+// Pour chaque decalage en temps, on estime un amplitude et calcule le X2
+// Pour chaque decalage en temps, on estime une amplitude et calcule le X2
         for(int k=0; k<SSTnbPasFit; k++) {
                 DeltaTemps=TFExcursion-k*TPasTFit;
 …
         X2Min=1.e99;
         int X2MinIndex=-1;
         for(int k=0; k<SSTnbPasFit; k++)
+        for(int k=0; k<SSTnbPasFit; k++) {
                 if (X2[k]<X2Min) {
                         X2Min=X2[k];
                         X2MinIndex=k;
+                }
+        }
         if((X2MinIndex>0)&&(X2MinIndex<(SSTnbPasFit-1))) {
 // On fitte une parabole sur les trois X2 qui encadrent le min
         double Di;      // Index (non entier) du minimum du X2
+                double Di;      // Index (non entier) du minimum du X2
                 double X=(X2[X2MinIndex+1]-X2Min)/(X2Min-X2[X2MinIndex-1]);
                 Di=0.5*(X+1)/(X-1);
+// Le X2 est mimimum pour le n° de pas de fit X2MinIndex+Di
 // On calcule l'amplitude estimee a ce temps
                 double Som1=0.; pFitArray=pFitArray0;
                 DeltaTemps=TFExcursion-(X2MinIndex-Di)*TPasTFit;
                 for(int i=0; i<SSTLongIndexEvt(); i++, pFitArray++)
                         Som1+=(*pFitArray)*FPulse.PulseShape(i*PasTempsEch-DeltaTemps);
                 AmplRes=Som1/AmplNorm;
+//Le X2 est mimimum pour le n° de pas de fit X2MinIndex+Di
+                                Som1+=(*pFitArray)*FPulse.PulseShape(i*PasTempsEch-DeltaTemps);
+                AmplFin=Som1/AmplNorm;
 #ifdef PEtoileDebug
                 double* PulseTrak=new double [SSTLongIndexEvt()];
+                double *PulseTrak=new double [SSTLongIndexEvt()];
 #endif
+// On calcule le X2 correspondant
                 double X2Som1=0.; double dummy=0.;pFitArray=pFitArray0;
                 for(int i=0;i<SSTLongIndexEvt();i++,pFitArray++) {
 #ifdef PEtoileDebug
                         PulseTrak[i]=FPulse.PulseShape(i*PasTempsEch-DeltaTemps);
 #endif
                         dummy=(*pFitArray)-AmplRes*FPulse.PulseShape(i*PasTempsEch-DeltaTemps);         // Différence
+                        dummy=(*pFitArray)-AmplFin*FPulse.PulseShape(i*PasTempsEch-DeltaTemps);         // Différence
                         X2Som1+=dummy*dummy;                                                                                                            // Au carré
+                }
                 X2Min=X2Som1/(SSTLongIndexEvt()*Vrms*Vrms);
+                TIndex=indexDebutPiste+IndexDebutEvt+DeltaTemps/PasTempsEch;
+                                        //Ouff!!!!
+                Tfin= indexDebutPiste+IndexDebutEvt+DeltaTemps/PasTempsEch;
+                        //Ouff!!!!
+/*
+                cout<<setprecision(9)<<endl;
+                cout<<Tfin<<" "<<indexDebutPiste<<" "<<IndexDebutEvt<<" "<<DeltaTemps/PasTempsEch<<endl;
+*/
                 SuccesAFine=true;
+#ifdef PEtoileDebug
+                delete[] PulseTrak;
+#endif
+        }
         else {                          // Echec, c'est tres bizarre! On doit etre dans le bruit

trunk/Poubelle/archTOI.old/pisteetoile.h

-              r394
+              r534
 #include "formepulse.h"
 /*      // D. Yvon, CE Saclay, 08/99
+/*      // D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
         // Analyse the Photodiode track to find star pulse.
         // Caracterise the time and amplitude of the star Pulse
         We analyse a set of 32 samples. The 16 first are considered as prepulse sample.
+        We analyse a set of 48 samples. The 32 first are considered as prepulse sample.
         On the prepulse:
                 We compute baseline and substract it.
 …
                 Look for trigger
         If triggered :
+                Look for maximum. If max found 20< Max< 28
+                Look for maximum.
+                Compute Full width half maximum
+                Fit a paroabole on max of the pulse
+                return max fitte and time of max fitte
+        if analfine
                 Iterate time fit on pulse area and compute amplitude and X2
                 Find best X2.
 …
                 Build SSTStar.
         That's it. Test code tstpisteetoile.cc can be found in tstroutines folder.
+        That's it.
         When  done push the track eigth samples forward and iterate analysis.
 */
 …
 //      Parametrisation du Trigger
 #define SeuilGachette (6.)                              // On déclenche à SeuilGachette sigmas
 #define SSTLsbVal (10./4096)
+#define SSTLsbVal (10./4096.)
 //      Parametrisation Analyse Fine
 #define TFExcursion (15.e-3)                    // Secondes: Excursion du fit autour de MaxIndex
 #define TPasTFit (2.5e-3)                               // Secondes:
+#define TPasTFit (1.5e-3)                               // Secondes:
                         // Secondes:
 …
         inline int SSTLargEvtIndex();
         inline int SSTLongIndexEvt();
+        static TransFuncElec Trapani;
+        static FormePulse FPulse;
 protected:
         double PasTempsEch;     //Secondes
 // Idendificateur, repères et gestion
         short NoPhDiod;
         int indexDebutPiste;    // A actualiser dans push() et Fill()! BUGG!!!!
+        long indexDebutPiste;   // A actualiser dans push() et Fill()
         bool DiodVivante;               // Si false, On ne fait rien
 …
         double Vrms;
         short MaxIndex;                 // Index de l'amplitude Max en analyse Trig()
+        double TIndex;
         double AmplMax;
+        double AmplRes;
         int widthHalfMax;
+        int SSTnbPasFit;
 //      double integ;
 // Variables Resultat d'analyse fine
+        double AmplRes;
+        double TIndex;
+        double AmplFin;
+        double Tfin;
         double X2Min;
         bool SuccesAFine;
 …
         double* X2;             //Tableau des X2 calculés
+        static TransFuncElec Trapani;
+        static FormePulse FPulse;
 };

trunk/Poubelle/archTOI.old/sstetoile.cc

-              r394
+              r534
 #include "sstetoile.h"
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 void SSTEtoile::print(ofstream& str) {
         str<<setprecision(11);
+        str<<InpCurrent<<'\t'<<TEchan<<'\t'<<NoDiode<<'\t'<<LargMiHauteur;
+        str<<'\t'<<InpCurrNoisRMS<<endl;
+        //      <<'\t'<<X2Calc<<'\t'<<FitForme
+        str<<InpCurrent<<'\t'<<InpCurFin<<'\t'<<TEchan<<'\t'<<TFin<<'\t'<<NoDiode<<'\t'<<LargMiHauteur;
+        str<<'\t'<<InpCurrNoisRMS<<'\t'<<X2Calc<<'\t'<<FitForme<<endl;
         return;
+}
 string SSTEtoile::printHeader() {
         string a("InpCurrent\tTempsEchan\tNoDiode\tLarMiHauteur\tICurrRMS\tX2Calc\tFitForme");
+        string a("InpCurrent\tInpCurFin\tTempsEchan\tTFin\tNoDiode\tLarMiHauteur\tICurrRMS\tX2Calc\tFitForme");
         return a;
+}

trunk/Poubelle/archTOI.old/sstetoile.h

-              r394
+              r534
 #ifndef SSTEtoile_H
 #define SSTEtoile_H
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 #include <stdio.h>
 …
 public:
         double InpCurrent;              // En Ampère
+        double InpCurFin;               // En ampère,   Analyse fine
 //      double IInteg;                  // En Ampère
+        double TEchan;                  // Temps en indexEchantillon
+        double TEchan;                  // Temps en indexEchantillon analyse trig
+        double TFin;                    // Temps en indexEchantillon analyse fine
         double InpCurrNoisRMS;  // Bruit en Ampère
         double X2Calc;                  // Sans dimension
+        double X2Calc;                  // Sans dimension, analyse fine
         int NoDiode;
         bool FitForme;                  // Fit de forme?
 …
         void print(ofstream& str);
         string printHeader();
+        static string printHeader();
 };

trunk/Poubelle/archTOI.old/ssthandler.cc

-              r416
+              r534
 #include <math.h>
+#include <iostream>
 #include "ssthandler.h"
 #include "pisteetoile.h"
+#ifdef SST_DUMPLOG
 ofstream SSTHandler::sstchass("SSTChassLogFile");
 bool     SSTHandler::sstchassinit=false;
-#endif
 // diodpermut[i] = channel de la diode i
 …
+}
+// D. Yvon, CEA/DAPNIA/SPP Saclay, 08/99
+#ifdef SSTStatLog
+        static ofstream sststat("SSTStatLog");
+        static int compteurBlock=0;
+        static int nbStar=0;
+        #include "ssthardware.h"
+        #include "archparam.h"
+#endif
 int SSTHandler::FindStars(block_type_sst* blk) {
         NbStarInBlock=0;
 …
                                         //On empile sur la map.
                                 NbStarInBlock++;
 #ifdef SST_DEBUG
+#ifdef SSTDUMP
                         // On écrit les pistes ayant déclénchées
                                 int NoDiodeEvt=LastStar.NoDiode;
+                                if(pPisteDump[NoDiodeEvt]->is_open())
+                                        for(int noSamp=0; noSamp<NbSampleBlock;noSamp++)
+                                                (*pPisteDump[NoDiodeEvt])<<noStarDet<<'\t'<<NoDiode<<'\t'<<NoFirstSpInBlock+noSamp<<'\t'<<Diodedata[noSamp]<<endl;
+                                if(NoDiodeEvt<NbChanDump) {
+                                        if(pPisteDump[NoDiodeEvt]->is_open())
+                                                for(int noSamp=0; noSamp<NbSampleBlock;noSamp++)
+                                                        (*pPisteDump[NoDiodeEvt])<<noStarDet<<'\t'<<NoDiode<<'\t'<<NoFirstSpInBlock+noSamp<<'\t'<<Diodedata[noSamp]<<endl;
+                                }
 #endif
                                 noStarDet++;
 …
                 LastBlockStarVec.push_back((*StarHistReIter).second);
+#ifdef SST_DUMPLOG
+        if(NbStarInBlock>0) {
+#ifdef SSTChassLog
+        if(NbStarInBlock>0) {
         // On écrit les étoiles detectées
                 vector<SSTEtoile>::reverse_iterator StarVecRevIter;
 …
                         (*StarVecRevIter).print(sstchass);
+        }
+#endif
+#ifdef SSTStatLog
+        int nbBlockincircle=(int) (PeriodRotTeles/(archParam.acq.perEch*72.));
+//      if((compteurBlock%nbBlockincircle)==0)
+        {
+//      if((compteurBlock%5)==0) {
+                sststat<<numero_block(blk)<<'\t'<<nbStar<<'\t'<<noStarDet<<endl;
+                nbStar=0;
+        }
+                nbStar+=NbStarInBlock;
+                compteurBlock++;
 #endif
 …
         for(int i=0; i<NbPhotDiodBarette; i++) PisteBar[i]=new PisteEtoile(i);
         StarHistoryMap.clear();
-#ifdef SST_DUMPLOG
         if (!sstchassinit) {
           sstchass<<LastStar.printHeader()<<endl;
           sstchassinit = true;
+        }
+#endif
+#ifdef SST_DEBUG
+#ifdef SSTDUMP
         char s[32];
         string fileName;
 …
+                }
 */
                 pPisteDump[pistNumb]=new ofstream(fileName.c_str(),ios::out|ios::trunc);
+                pPisteDump[pistNumb]=new ofstream(fileName.c_str(),ios::out|ios::app);
                 if(!(*pPisteDump[pistNumb]).is_open()) {
 …
         delete[] PisteBar;
 #ifdef SST_DEBUG
+#ifdef SSTDUMP
         for (int pistNumb=0; pistNumb<NbChanDump; pistNumb++) {
                 pPisteDump[pistNumb]->close();
                 delete pPisteDump[pistNumb];
+        }
         delete  pPisteDump;
+        delete[] pPisteDump;
 #endif
 …
+{
   lastBlkNum = numero_block(blk);
+  if (prcTodo != 0) {
+    for (int i = 0; i<nb_per_block*2; i++) {
+      DecodeTMBlock(blk, i, diodeRaw[i]);
+    }
+  for (int i = 0; i<nb_per_block*2; i++) {
+    DecodeTMBlock(blk, i, diodeRaw[i]);
+  }
   if (prcTodo & permDiode) {

trunk/Poubelle/archTOI.old/ssthandler.h

-              r416
+              r534
 //#define SST_DEBUG
+//#define SST_DUMPLOG
+#define NbChanDump (10)
+//#define SSTDUMP
+#define SSTChassLog
+//#define SSTStatLog
+#define NbChanDump (3)
 #define NbPhotDiodBarette (46)
 #define NbEtInMapMax    (100)
-#define SSTChassLog (true)
 #define NbSampleBlock (2*nb_per_block)
 // Fin Zone findStar
 …
   void NeedProcess(int prcMask);
   void ProcessBlock(block_type_sst*);
-  bool NeedBlocks() {return prcTodo != 0;}
   static void Has2Bars(bool, int elecOffset=0); // elecOffset : which wired together
 …
   double getStarTime(int iSampl, int istar);
                         // return star time in this sample
-  static int getDiodPermut(int i) {return diodpermut[i];}
 protected:
 …
   PisteEtoile** PisteBar;
 #ifdef SST_DEBUG
+#ifdef SSTDUMP
         ofstream** pPisteDump;
 #endif
 …
   int noStarDet;
   int FirstSplNb(block_type_sst* blk); // Calcule l'index du premier echantillon du block analysé
   static ofstream sstchass;
   static bool     sstchassinit;

trunk/Poubelle/archTOI.old/ssthardware.h

-              r394
+              r534
 #ifndef SSTHARDWARE_SEEN
 #define SSTHARDWARE_SEEN
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 // Telescope SST
 …
 // Paramètres du Vol
 #define PeriodRotTeles (20.)                                            // Periode de rotation du telescope en seconde
+#define PeriodRotTeles (20.)                                    // Periode de rotation du telescope en seconde
+//#define PeriodRotTeles (30.)                                  // version Kiruna
 #endif

trunk/Poubelle/archTOI.old/templocator.cc

-              r436
+              r534
 #include "fitsio.h"
 #include "plgalcross.h"
+#include "archparam.h"
 #ifndef M_PI
 …
   // On trouve les croisements juste avant et juste apres notre sampleNum
-  int icross;
   for (icross=0; icross<nGalCross; icross++) {
     if (crossings[icross] > sampleNum) break;
 …
   if (fabs(azCr2-220) < fabs(azCr1-220)) azCross = azCr2;
   double rotSpeed = 360./(crossings[icross] - crossings[icross-1]); // °/sample
+  rotSpeed = 360./(crossings[icross] - crossings[icross-1]); // °/sample
   azimBolGC = azCross - (sampleNum - (crossings[icross-1]+12))*rotSpeed;
 …
   if (azimBolGC > 360) azimBolGC -= 360;
+}
+double TempLocator::getRotSpeed(int sampleNum) {
+  findGeomFromGC(sampleNum);
+  return rotSpeed / archParam.acq.perEch;
+}
+int  TempLocator::getCrossSamples(int sampleNum, int& SN1, int& SN2) {
+  findGeomFromGC(sampleNum);
+  if (icross == 0 || icross >= nGalCross) return -1;
+  SN1 = crossings[icross-1];
+  SN2 = crossings[icross];
+  return 0;
+}
 void TempLocator::getAltAzBolo(int sampleNum, int ibolo, double& elv, double& az) {
 …
   aa_hadec (lat * M_PI/180, elv * M_PI/180, az * M_PI/180, &ha, &dec);
   ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
+  while (ra>24) ra -= 24;
+  while (ra<0) ra += 24;
   dec = dec * 180. / M_PI;
   return ra;
 …
 double TempLocator::getDeltaCenter(int sampleNum) {
   return getDeltaBolo(sampleNum, 11);
+}
+double TempLocator::getAzimutCenter(int sampleNum) {
+  return getAzimutBolo(sampleNum, 11);
+}
+double TempLocator::getElvCenter(int sampleNum) {
+  return getElvBolo(sampleNum, 11);
+}

trunk/Poubelle/archTOI.old/templocator.h

-              r436
+              r534
   double getAlphaBolo(int sampleNum, int ibolo);
   double getDeltaBolo(int sampleNum, int ibolo);
+  double getAzimutCenter(int sampleNum);
+  double getElvCenter(int sampleNum);
   double getAlphaCenter(int sampleNum);
   double getDeltaCenter(int sampleNum);
+  double getRotSpeed(int sampleNum); // deg/sec
+  int    getCrossSamples(int sampleNum, int& SN1, int& SN2);
   void getAltAzBolo(int sampleNum, int ibolo, double& elv, double& az);
 …
   int xSampleNum; // pour le dernier calcul
   double azimBolGC;
+  double rotSpeed;
+  int icross;
   void findGeomFromGC(int sampleNum); // pour le bolo qui voit les xing

trunk/Poubelle/archTOI.old/toiiter.cc

-              r436
+              r534
 // Eric Aubourg         CEA/DAPNIA/SPP   juillet 1999
+#ifndef nb_max_bolo
+#define _archeops               //  Acquisition Archeops (avec transputer)
+#define programme
+extern "C" {
+#include "archeops.h"
+#include "arcunit.h"
+}
+#include "ark.h"
+#include "toiiter.h"
+#include "toimanager.h"
+#include "toiproducer.h"
+#include "archparam.h"
+#include "asigps.h"
+#include "toiauxgpsproducer.h"
+#include <iostream.h>
+#include <fstream.h>
+#define CHK_INITED if (initDone) \
+        throw ArchExc("Trying to modify TOIIter after init done.");
+TOIIter::TOIIter() {
+  mjdStart = -999999999;
+  mjdEnd   =  999999999;
+  utcStart = -999999999;
+  utcEnd   =  999999999;
+  sStart   = -999999999;
+  sEnd     =  999999999;
+  underSample = 1;
+  curSample = -1;
+  initDone = false;
+  incDone  = false;
+}
+#ifdef __MWERKS__
+#pragma mark -
 #endif
+#include "toiiter.h"
+#include "toiinterpolator.h"
+#include "archparam.h"
+#include <dirent.h>
+#include <iostream.h>
+#include "templocator.h"
+// Format bloc GPS
+// $GPGGA,hhmmss.ss,ddmm.mmmm,n,dddmm.mmmm,e,q,ss,y.y,a.a,z,
+TOIIter::TOIIter() {
+ // Toutes les valeurs initialisees par defaut. Le TOISvr les positionnera,
+ // puis lancera l'initialisation
+void TOIIter::addDirectory(string dir) {
+  CHK_INITED
+  fset.addDirectory(dir);
+}
+void TOIIter::addFile(string fn) {
+  CHK_INITED
+  fset.addFile(fn);
+}
+void TOIIter::setMJDInterval(double tStart, double tEnd) {
+  CHK_INITED
+  if (tStart>0) mjdStart = tStart;
+  if (tEnd>0)   mjdEnd = tEnd;
+}
+void TOIIter::setUTCInterval(double tStart, double tEnd) {
+  CHK_INITED
+  if (tStart>0) utcStart = tStart;
+  if (tEnd>0)   utcEnd = tEnd;
+}
+void TOIIter::setSNInterval(long tStart, long tEnd) {
+  CHK_INITED
+  if (tStart>0) sStart = tStart;
+  if (tEnd>0)   sEnd = tEnd;
+}
+void TOIIter::setUnderSample(int n) {
+  CHK_INITED
+  if (n<=1) n=1;
+  underSample = n;
+}
+int TOIIter::getUnderSample() {
+  return underSample;
+}
+void TOIIter::addTOI(TOI& toi, bool trg) {
+  CHK_INITED
+  TOIProducer* prod = TOIManager::findTOIProducer(toi);
+  if (!prod) throw ArchExc("Cannot produce " + toi.fullName());
+  prod->addTOI(toi, this);
+  request.push_back(TOIInfo(toi,prod,trg?1:0));
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
+void TOIIter::registerReqHandler(RequestHandler* h) {
+  CHK_INITED
+  handlers.push_back(h);
+}
+void TOIIter::readReq(istream& str) {
+  CHK_INITED
+   if (!incDone) defaultInclude();
+   string line;
+   while (str) {
+     getline(str,line);
+     if (!str) break;
+     if (line.substr(0,4)=="#END" && (line.length()==4 || line[5] == ' ')) break;
+     if (line[0] != '@' && line[0] != '#') continue;
+     if (!processRequest(line)) {
+       throw ArchExc("Unrecognized directive " + line);
+     }
+   }
+}
+bool TOIIter::processRequest(string line) {
+  if (line[0] == '#') {
+    int x = line.find(' ');
+    string keyw = line.substr(0, x);
+    string args = (x>0) ? line.substr(x) : string("");
+    bool handled = processOption(keyw,args);
+    for (list<RequestHandler*>::iterator i = handlers.begin();
+         i != handlers.end(); i++) {
+      handled |= (*i)->processOption(keyw,args);
+    }
+    return handled;
+  }
+  if (line[0] == '@') {
+    TOI toi(line.substr(1));
+    // if (kind  == sampleNum || kind == mjd || kind == mutc) notrig = true;
+    bool handled = processTOIReq(toi,line);
+    for (list<RequestHandler*>::iterator i = handlers.begin();
+         i != handlers.end(); i++) {
+      handled |= (*i)->processTOIReq(toi,line);
+    }
+    return handled;
+  }
+  return false;
+}
+bool TOIIter::processTOIReq(TOI& toi,string)
+{
+  TOI toi2 = toi;
+  bool trg=true;
+  if (toi.options.find("notrig") != toi.options.end()) {
+    toi2.options.erase("notrig");
+    trg=false;
+  }
+  addTOI(toi2, trg);
+  toi.unit = toi2.unit;
+  return true;
+}
+bool TOIIter::processOption(string key, string arg)
+{
+  if (arg.length()>0 && arg[0] == ' ') {
+    arg = arg.substr(arg.find_first_not_of(' '));
+  }
+  if (key == "#MJDRANGE") {
+    double tmin, tmax;
+    sscanf(arg.c_str(), "%lg %lg", &tmin, &tmax);
+    setMJDInterval(tmin, tmax);
+  } else if (key == "#UTCRANGE") {
+    double tmin, tmax;
+    sscanf(arg.c_str(), "%lg %lg", &tmin, &tmax);
+    setUTCInterval(tmin, tmax);
+  } else if (key == "#SNRANGE") {
+    long tmin, tmax;
+    sscanf(arg.c_str(), "%ld %ld", &tmin, &tmax);
+    setSNInterval(tmin, tmax);
+  } else if (key == "#PATH") {
+    addDirectory(arg);
+  } else if (key == "#FILE") {
+    addFile(arg);
+  } else if (key == "#UNDERSAMPLE") {
+    setUnderSample(atoi(arg.c_str()));
+  } else if (key == "#MJD0") {
+    double t0;
+    sscanf(arg.c_str(), "%lg", &t0);
+    archParam.acq.tBlock0 = t0;
+  } else if (key == "#UTCORIGIN") {
+    double t0;
+    sscanf(arg.c_str(), "%lg", &t0);
+    archParam.acq.utcOrigin = t0;
+  } else if (key == "#PERECH") {
+    double t0;
+    sscanf(arg.c_str(), "%lg", &t0);
+    archParam.acq.perEch = t0;
+  } else if (key == "#ASIGPS") {
+    ASIGPS* gps = new ASIGPS(arg);
+    TOIManager::registerProducer(new TOIAuxGPSProducer(gps));
+    //    gps->FitsDump("GPSDump.fits");
+  } else if (key == "#INCLUDE") {
+    ifstream f(arg.c_str());
+    readReq(f);
+  } else {
+    return false;
+  }
+  return true;
+}
+void TOIIter::defaultInclude() {
+  incDone = true;
+  processRequest("#REQVERSION V_270999");
+  processRequest("#MJD0 1376.8358818");
+  processRequest("#PERECH 0.005836818076");
+  processRequest("#UTCORIGIN 1376.5");
+  processRequest("#ASIGPS ASI_GPS_archeops1999.ascii");
+  processRequest("#COMMENT Archtoi V2 -- october 1999 -- Eric Aubourg CEA/DAPNIA");
+  processRequest("#COMMENT ***WARNING***");
+  processRequest("#COMMENT ***SOME TOI'S ARE PRELIMINARY***");
+  processRequest("#COMMENT gyroSpeed is not calibrated");
+  processRequest("#COMMENT azimut/alpha/delta use galaxy crossings and ASI GPS data");
+  processRequest("#COMMENT   and assume no pendulation");
+  processRequest("#COMMENT Focal plane center elevation found at 41.5 deg");
+  processRequest("#COMMENT   with Jupiter");
+  processRequest("#COMMENT boloMuV2 is not protected against glitches");
+}
+#ifdef __MWERKS__
+#pragma mark -
+#endif
+void TOIIter::init() {
+  if (initDone) return;
+  initDone = true;
+ file = NULL;
+ directory = "";
+  if (utcStart > 0) {
+    double t = (utcStart/24.) + archParam.acq.utcOrigin;
+    if (t > mjdStart) mjdStart=t;
+  }
+  if (utcEnd > 0) {
+    double t = (utcEnd/24.) + archParam.acq.utcOrigin;
+    if (t < mjdEnd) mjdEnd=t;
+  }
+  fset.setMJDRange(mjdStart-0.01, mjdEnd+0.01);
+  fset.setSNumRange(sStart-1000, sEnd+1000);
+  fset.init();
+  //curSample = fset.getSampleIndex();
+}
+bool TOIIter::isTrig(int column) {
+  return (request[column].third & triggering) != 0;
+}
+ files.clear();
+ isOnBoardRecorder = false;
+ imes=0;
+ tStart = -9.e99;
+ tEnd   = 9.e99;
+ utcStart = -9.e99;
+ utcEnd   = 9.e99;
+ sStart = -999999999;
+ sEnd   =  999999999;
+ //tBlock0 = -1;
+ // perEch = -1;
+ // Values for Trapani 99 = default values
+ //tBlock0 = 1376.8358818;
+ //perEch = 0.005836818076;
+ trigMask = 0;
+ rawIter = NULL;
+ interp  = NULL;
+ lastSample = -1;
+ maxLookAhead = 10000;
+ auxGPS = NULL;
+ initDone = false;
+ underSample = 1;
+}
+TOIIter::TOIIter(TOIIter const& x) {
+  directory = x.directory;
+  files = x.files;
+  if (x.initDone && x.curFile != x.files.end()) curFile = files.find(*(x.curFile));
+  isOnBoardRecorder = x.isOnBoardRecorder;
+  imes   = x.imes;
+TOI TOIIter::getKind(int column) {
+  return request[column].first;
+}
+long TOIIter::getSampleNum() {
+  return curSample;
+}
+bool TOIIter::next() {
+  if (!initDone) init();
+  for (int ii=0; ii<underSample; ii++)
+    if (!next1()) return false;
+  tStart   = x.tStart;
+  tEnd     = x.tEnd;
+  utcStart = x.utcStart;
+  utcEnd   = x.utcEnd;
+  sStart   = x.sStart;
+  sEnd     = x.sEnd;
+  return true;
+}
+bool TOIIter::next1() {
+// On tente de produire curSample+1 pour toutes les toi
+// Eventuellement en avancant sur le fichier...
+// Puis on regarde si une TOI triggering a eu une nouvelle
+// valeur.
+// Si on a epuise les fichiers de donnees, on s'arrete des qu'aucune
+// TOI n'a de valeurs apres curSample...
+  if (curSample <= 0) curSample = sStart-1;
+  if (curSample <= 0) curSample = fset.getSampleIndex()-1;
+  trigMask = x.trigMask;
+  infos = x.infos;
+  if (x.file)
+    file = new ArcheopsFile(*x.file);
+  else
+    file = NULL;
+  if (x.rawIter) {
+    rawIter = new TOIIter(*x.rawIter);
+    interp = new TOIInterpolator[infos.size()];
+    for (int i=0; i<infos.size(); i++)
+      interp[i] = x.interp[i];
+  } else {
+    rawIter = NULL;
+    interp  = NULL;
+  }
+  if (curSample >= sEnd || archParam.acq.SN2MJD(curSample) >= mjdEnd) return false;
+  lastSample = x.lastSample;
+  maxLookAhead = x.maxLookAhead;
+  auxGPS = x.auxGPS;
+  if (auxGPS) auxGPS = auxGPS->clone();
+  initDone = x.initDone;
+  underSample = x.underSample;
+}
+TOIIter::~TOIIter() {
+  delete file;
+  delete rawIter;
+  delete[] interp;
+  delete auxGPS;
+}
+#ifdef __MWERKS__
+#define filesep ':'
+#else
+#define filesep '/'
+#endif
+void TOIIter::Init() {
+ if (initDone) return;
+ initDone = true;
+ if (utcStart > 0) {
+   double t = (utcStart/24.) + archParam.acq.utcOrigin;
+   if (t > tStart) tStart=t;
+ }
+ if (utcEnd > 0) {
+   double t = (utcEnd/24.) + archParam.acq.utcOrigin;
+   if (t < tEnd) tEnd=t;
+ }
+ // On a soit un repertoire, soit une liste de fichiers....
+  if (directory == "") {
+    if (files.empty()) {   // Ni repertoire, ni fichiers
+      cerr << "toiiter : pas de repertoire, pas de fichiers" << endl;
+      exit(-1);
+    } else {
+      // On a deja une liste de fichiers
+  for (int i=0; i<request.size(); i++)
+    request[i].second->wontNeedEarlier(request[i].first, this, curSample);
+  bool endFound = false;
+  while(1) {
+    curSample++;
+    for (vector<TOIInfo>::iterator i = request.begin(); i != request.end(); i++) {
+      while ((*i).second->canGetValueLater(curSample, (*i).first)) {
+        if (! fset.next()) {endFound = true; break;} //return false;
+      }
+    }
+  } else { // On a un repertoire a explorer
+    // On cherche soit les fichiers dans le repertoire donne, soit des fichiers
+    // dans un sous-repertoire "arch-YY_MM_DD". Les fichiers ont un nom en
+    // "hYY_MM_DD-hh_mm_ss".
+    // Pour l'enregistreur de vol, les fichiers ont un nom en ARKxxxxxx.DAT
+    if (directory[directory.length()-1] != filesep)
+      directory += filesep;
+    DIR* dir = opendir(directory.c_str());
+    struct dirent* ent;
+    while ((ent = readdir(dir)) != NULL) {
+      // si c'est un repertoire, avec un nom de jour, il faut l'explorer...
+      if (!strncmp(ent->d_name, "arch-", 5)) {
+        double mjd = ArcheopsFile::decodeMJD(ent->d_name+5) - 2./24.; // ENTIER + .5 en temps local!
+        if (mjd >= tStart - 1. && mjd <= tEnd) {
+          string direc2 = directory + ent->d_name + filesep;
+          DIR* dir2 = opendir(direc2.c_str());
+          struct dirent* ent2;
+          while ((ent2 = readdir(dir2)) != NULL) {
+            if (*ent2->d_name == 'h') {
+              double mjd2 = ArcheopsFile::decodeMJD(ent->d_name+1) - 2./24.;
+              if (mjd2 >= tStart - 1./24. && mjd2 <= tEnd) {
+                files.insert(direc2 + ent2->d_name);
+              }
+            }
+          }
+        }
+      } else {
+        if (!isOnBoardRecorder && *ent->d_name == 'h') {
+          double mjd = ArcheopsFile::decodeMJD(ent->d_name+1) - 2./24.; // $CHECK$ UTCOffset
+          if (mjd >= tStart - 1./24. && mjd <= tEnd) {
+            files.insert(directory + ent->d_name);
+          }
+        } else if (isOnBoardRecorder) {
+          if (strncmp(ent->d_name, "ARK", 3) && strncmp(ent->d_name, "ark", 3)) continue;
+          char * sfx = ent->d_name + strlen(ent->d_name) - 4;
+          if (strcmp(sfx, ".DAT") && strcmp(sfx, ".dat")) continue;
+          files.insert(directory + ent->d_name);
+        }
+      }
+    bool found=false;
+    bool valuesAhead = false;
+    for (vector<TOIInfo>::iterator i = request.begin(); i != request.end(); i++) {
+      if (((*i).third & triggering) == 0) continue;
+      if ((*i).second->canGetValue(curSample, (*i).first)) {found=true;break;}
+      if ((*i).second->lastSampleNum((*i).first)>curSample) valuesAhead=true;
+    }
+    closedir(dir);
+  }
+  ScanFiles();
+  curFile = files.begin();
+  file = new ArcheopsFile((*curFile).c_str());
+  cout << "opening file " << (*curFile).c_str() << endl;
+  // On avance jusqu'a avoir au moins un bloc param et un bloc reglage,
+  // car on ne peut rien faire sans...
+  // Si on a des donnees de l'enregistreur de vol, pas de bloc param, et
+  // on en simule un
+  double oldTStart = tStart;
+  tStart = -9.e99; // pour init, on accepte des blocs avant tstart....
+  if (!file->lastParam()) {
+    if (isOnBoardRecorder) {
+      extern param_bolo  parametr;
+      block_type_param block;
+      block.param = parametr;
+      valide_block((block_type_modele*)&block, block_param, 0);
+      file->forceBlock((block_type_modele*)&block);
+    } else {
+      file->nextBlock(block_param_mask);
+    }
+  }
+  if (!file->lastReglage()) file->nextBlock(block_reglage_mask);
+  // On cherche un bloc GPS pour avoir la correspondance timestamp/UTC.
+  // Pour le moment, on se fonde sur le premier bloc GPS. On pourra faire
+  // mieux en prenant le min de tous les delta_T, a condition d'avoir un
+  // peu plus de details sur la facon dont le GPS est lu.
+  if (archParam.acq.tBlock0 < 0) {
+    archParam.acq.tBlock0 = file->getStartMJD();
+    file->pushMark();
+    if (file->lastGPS() || file->nextBlock(block_gps_mask)) {
+      // le temps du bloc courant, en secondes
+      double dt = file->blockNum() * file->perBlock();
+      archParam.acq.tBlock0 = file->getGPSMJD() - dt/86400.;
+    } else { // pas de bloc GPS...
+      archParam.acq.tBlock0 = file->getStartMJD();
+    }
+    file->popMark();
+  }
+  tStart = oldTStart;  // on restaure
+  if (archParam.acq.perEch < 0)
+    archParam.acq.perEch = file->perEchant();
+  bool hasInterp = false;
+  trigMask = 0;
+  for (vector<info>::iterator i = infos.begin(); i != infos.end(); i++) {
+    if ((*i).interpolated) hasInterp = true;
+    if ((*i).triggering) {
+      switch ((*i).kind) {
+      case boloMuV:
+      case boloMuV2:
+      case boloRawMuV:
+      case boloRawMuVCN:
+        trigMask |= block_bolo_mask;
+        break;
+      case boloMuV2T:
+      case boloRes:
+      case boloTemp:
+         trigMask |= block_bolo_mask | block_reglage_mask;
+        break;
+      case dilDAC:
+      case dilSwitch:
+        trigMask |= block_dilution_mask;
+      case boloGainAmpli:
+      case boloDACV:
+      case boloDACI:
+        trigMask |= block_reglage_mask;
+        break;
+      case gpsTime:
+      case longitude:
+      case latitude:
+      case altitude:
+        trigMask |= block_gps_mask;
+        break;
+      case azimuthBolo:
+        trigMask |= block_bolo_mask;
+        //file->needSSTProcessMask(SSTHandler::findPeriod);
+        //trigMask |= block_sst_mask;
+        break;
+      case sstStarCnt:
+      case sstStarZ:
+      case sstStarF:
+      case sstStarT:
+        file->needSSTProcessMask(SSTHandler::findStars);
+        trigMask |= block_sst_mask;
+        break;
+      case sstDiode:
+      case sstDiodeCN:
+        file->needSSTProcessMask(SSTHandler::permDiode);
+        trigMask |= block_sst_mask;
+        break;
+      case sstChannel:
+      case sstChannelCN:
+        trigMask |= block_sst_mask;
+        break;
+      case gyroRaw:
+      case gyroV:
+        trigMask |= block_gyro_mask;
+        break;
+      case gyroSpeed:
+        trigMask |= block_gyro_mask; // $CHECK$ + info to calibrate gyros
+        break;
+      case alphaRotAxis:
+      case deltaRotAxis:
+      case alphaSst:
+      case deltaSst:
+      case alphaFPAxis:
+      case deltaFPAxis:
+      case alphaBolo:
+      case deltaBolo:
+        trigMask |= block_bolo_mask;
+        //file->needSSTProcessMask(SSTHandler::findAxis);
+        //trigMask |= block_sst_mask;
+        break;
+      case alphaZenith:
+      case deltaZenith:
+        trigMask |= block_gps_mask;
+        break;
+      case voyantEVO:
+      case voyantEVF:
+      case commandeEV0:
+      case commandeEVF:
+      case commandeEVB:
+      case commandeEVV:
+      case pressEnt3He:
+      case debit3He:
+      case pressSor3He:
+      case pressEnt4He:
+      case debit4He:
+      case pressSor4He:
+      case pressAirVanne:
+      case pressPompChar:
+      case pressMembrane:
+      case pressExterne:
+      case tensPile10T:
+      case tensPileP18D:
+      case tensPileM18D:
+      case tensPile10B:
+      case tensPileP18B:
+      case tensPileM18B:
+      case tensPileCh:
+      case swPile5:
+      case swPile15:
+      case tempCaissH1:
+      case tempCaissH2:
+      case tempCaissB1:
+      case tempCaissB2:
+      case tempCaissTHe:
+      case tempCaissPiles:
+      case tempCaissDrv:
+      case pressHeBain:
+      case pressPirani:
+        trigMask |= block_dilution_mask;
+        break;
+      }
+    }
+  }
+  if (trigMask & (block_bolo_mask | block_sst_mask)) {
+    imes = 9999;
+  } else {
+    imes = 0;
+  }
+  if (hasInterp) {
+    rawIter = new TOIIter(*this);
+    interp = new TOIInterpolator[infos.size()];
+    for (int i=0; i<infos.size(); i++) {
+      rawIter->infos[i].interpolated = false;
+    }
+    delete file; file = NULL; // on ne travaille plus sur le fichier directement...
+  }
+}
+void TOIIter::ScanFiles() {
+  file1stSamp.clear();
+  cout << "Scanning all files" << endl;
+  // Petite astuce pour les STL non conformes comme celles de digital
+  // qui ne supportent pas files.erase(i) suivi de i++....
+  set<string> copy = files;
+  for (set<string>::iterator i = copy.begin(); i != copy.end(); i++) {
+    ArcheopsFile fich((*i).c_str());
+    if (fich.nextBlock()) {
+      file1stSamp[*i] = fich.blockNum()*72; // premier numsample
+      cout << "File " << *i << " 1st sample = " << fich.blockNum()*72 << endl;
+    } else {
+      cout << "File " << *i << " unrecoverable, skipping" << endl;
+      files.erase(*i);
+    }
+  }
+  cout << "Scan done" << endl;
+  // Et maintenant, on ne garde que ceux qui tombent dans l'intervalle...
+  copy = files;
+  string prev="";
+  for (set<string>::iterator i = copy.begin(); i != copy.end(); i++) {
+    double smp = file1stSamp[*i];
+    double t   = archParam.acq.tBlock0 + smp * archParam.acq.perEch/86400.;
+    if (t>tEnd || smp>sEnd) {  // premier echantillon apres tEnd
+     files.erase(*i);
+     prev = "";
+     continue;
+    }
+    if (t<tStart || smp<sStart) { // premier echantillon avant tStart -> on vire le precedent si existe
+      if (prev != "") {
+        files.erase(prev);
+      }
+    }
+    prev = *i;
+  }
+}
+bool TOIIter::NextFile() {
+  if (rawIter)
+    return rawIter->NextFile();
+  if (files.empty()) {
+    return false;
+  } else {
+    if (curFile == files.end()) return false;
+    curFile++;
+    if (curFile == files.end()) return false;
+    cout << "opening file " << (*curFile).c_str() << endl;
+    ArcheopsFile* newfile = new ArcheopsFile((*curFile).c_str());
+    newfile->grabLastBlocs(*file);
+    delete file;
+    file = newfile;
+    return true;
+  }
+}
+bool TOIIter::Next() {
+  for (int ii=0; ii<underSample; ii++)
+    if (!Next1()) return false;
+    if (found) break;
+    if (endFound && !valuesAhead) return false;
+  }
   return true;
+}
+bool TOIIter::Next1() {
+  if (!initDone) Init();
+  while (1) {
+    if (!NextSample()) return false; // end of files
+    double t = getMJD();
+    if (t < tStart) continue;
+    if (t > tEnd) return false;
+    return true;
+  }
+}
+bool TOIIter::NextSample() {
+  if (rawIter) {  // Delegation pour interpolation
+    // Trouve prochain sample disponible
+    for (int k=0; k<2; k++) {
+      long smp = 2147483647L;
+      for (int i=0; i<infos.size(); i++) {
+        long ss = interp[i].nextSample(lastSample+1);
+        if (ss > 0 && ss < smp) smp=ss;
+      }
+      if (smp != 2147483647L) {
+        lastSample = smp;
+        break;
+      }
+      if (!fetchAhead())  // tout le monde etait en bout de course,
+        return false;     // on lit un echantillon, ca suffit, d'ou le k<2
+    }
+    // Verifie que tous les interpolateurs ont assez de donnees pour
+    // trouver la valeur correspondante
+    for (int i=0; i<infos.size(); i++) {
+      //rif (infos[i].interpolated)
+        while (interp[i].needMoreDataFor(lastSample) &&
+               rawIter->getSampleIndex() - lastSample < maxLookAhead)
+          if (!fetchAhead()) return false;
+    }
+    // On est pret...
+    return true;
+  }
+  // trigger sur info indexee dans bloc bolo, bloc gyro ou bloc sst ?
+  if (trigMask & (block_bolo_mask | block_sst_mask | block_gyro_mask )) {
+    imes++;
+    if (imes < file->nEchBlock()) return true;
+    imes = 0;
+  }
+  // soit pas d'info indexee, soit fin bloc courant...
+  while (1) {
+    if (file->nextBlock(trigMask)) {
+      while (file->sameBlockNumAhead()) {  // tant que meme numero de bloc, on lit
+        if (!file->nextBlock()) {          // fin de fichier ?
+          if (NextFile()) file->nextBlock(); // fichier suivant
+          else return false;    // tout fini
+        }
+      }
+      return true;
+    }
+    if (!NextFile()) return false;
+  }
+}
+/* double TOIIter::getTime() {                // MJD
+   // le temps du bloc courant, en secondes
+   double dt = file->blockNum() * file->perBlock();
+   return tBlock0 + dt/86400. + imes*file->perEchant()/86400.;
+ }
+ */
+int TOIIter::getColTOI(TOI const& toi) {
+  for (int i=0; i<request.size(); i++)
+    if (request[i].first == toi) return i;
+  throw ArchExc("getColTOI : no such TOI " + toi.name);
+}
 bool TOIIter::canGetValue(int column) {
+  if (!initDone) Init();
+   if (column < 0 || column >= infos.size()) return false;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS &&
+       (kind == longitude || kind == latitude || kind == altitude || kind == tsid ||
+        kind == alphaZenith || kind == deltaZenith)) {
+         double dummy;
+         return auxGPS->getLocation(getMJD(), dummy, dummy, dummy) == 0;
+   }
+   if (rawIter) {
+       return interp[column].canGet(lastSample);
+   }
+   int index = infos[column].index;
+   switch (kind) {
+   case sampleNum:
+   case internalTime:
+   case mjd:
+   case mutc:
+     return true;
+   case boloMuV:
+   case boloMuV2:
+   case boloMuV2T:
+   case boloRes:
+   case boloRawMuVCN:
+   case boloTemp:
+     if (imes==0 && file->llastBolo()==NULL) return false;
+     return file->lastBolo() != NULL;
+   case boloRawMuV:
+     return file->lastBolo() != NULL;
+   case boloGainAmpli:
+   case boloDACV:
+   case boloDACI:
+     return file->lastReglage() != NULL;
+   case dilDAC:
+   case dilSwitch:
+     return file->lastDilution() != NULL;
+   case sstDiode:
+   case sstChannel:
+   case sstDiodeCN:
+   case sstChannelCN:
+     return file->lastSST() != NULL;
+   case sstStarCnt:
+   case sstStarZ:
+   case sstStarF:
+   case sstStarT:{
+     if (file->lastSST() == NULL) return false;
+     int n = file->getNumbStar(imes);
+     return (n > 0 && index < n);
+   }
+   case gyroRaw:
+   case gyroV:
+   case gyroSpeed:
+     return (file->lastGyro() != NULL);
+   case gpsTime:
+     return file->hasGPSTime();
+   case longitude:
+   case latitude:
+   case alphaZenith:
+   case deltaZenith:
+     return file->hasGPSPos();
+   case altitude:
+     return file->hasGPSAlt();
+   case tsid:
+     return file->hasGPSPos();
+   case azimuthBolo:
+   case alphaRotAxis:
+   case deltaRotAxis:
+   case alphaFPAxis:
+   case deltaFPAxis:
+   case alphaBolo:
+   case deltaBolo:
+     return true;
+     //return (file->lastGPS() != NULL && file->lastSST() != NULL);
+   case alphaSst:
+   case deltaSst:
+     return false;
+   case voyantEVO:
+   case voyantEVF:
+   case commandeEV0:
+   case commandeEVF:
+   case commandeEVB:
+   case commandeEVV:
+   case pressEnt3He:
+   case debit3He:
+   case pressSor3He:
+   case pressEnt4He:
+   case debit4He:
+   case pressSor4He:
+   case pressAirVanne:
+   case pressPompChar:
+   case pressMembrane:
+   case pressExterne:
+   case tensPile10T:
+   case tensPileP18D:
+   case tensPileM18D:
+   case tensPile10B:
+   case tensPileP18B:
+   case tensPileM18B:
+   case tensPileCh:
+   case swPile5:
+   case swPile15:
+   case tempCaissH1:
+   case tempCaissH2:
+   case tempCaissB1:
+   case tempCaissB2:
+   case tempCaissTHe:
+   case tempCaissPiles:
+   case tempCaissDrv:
+   case pressHeBain:
+   case pressPirani:
+     return file->lastDilution() != NULL;
+   }
+   return false;
+}
+   TOIInfo& info = request[column];
+   return(info.second->canGetValue(curSample, info.first));
+}
 double TOIIter::getValue(int column) {
+  if (!initDone) Init();
+   if (column < 0 || column >= infos.size()) return -1;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS &&
+       (kind == longitude || kind == latitude || kind == altitude || kind == tsid)) {
+         double lat,lon,alt;
+         if (auxGPS->getLocation(getMJD(), lat, lon, alt)) return -99999;
+         if (kind == longitude) return lon;
+         if (kind == latitude)  return lat;
+         if (kind == altitude)  return alt;
+         if (kind == tsid) {
+           tSid.setLongitude(lon);
+           return tSid.getLST(getMJD());
+         }
+   }
+   if (rawIter) {
+     if (infos[column].interpolated)
+       return interp[column].getIValue(lastSample);
+     else
+       return interp[column].getEValue(lastSample);
+   }
+   int index = infos[column].index;
+   switch (kind) {
+   case sampleNum:
+     return getSampleIndex();
+   case internalTime:
+     return getSampleIndex() * archParam.acq.perEch;
+   case mjd:
+     return getMJD();
+   case mutc:
+     return (getMJD()-archParam.acq.utcOrigin)*24.;
+   case boloMuV:
+     return file->getBoloMuV(index, imes);
+   case boloMuV2:
+     return file->getBoloMuV2(index, imes);
+   case boloRawMuV:
+     return file->getBoloRawMuV(index, imes);
+   case boloRawMuVCN:
+     return file->getBoloRawMuVCN(index, imes);
+   case boloGainAmpli:
+     return file->getGainAmpli(index);
+   case boloDACV:
+     return file->getDACV(index);
+   case boloDACI:
+     return file->getDACI(index);
+   case boloMuV2T:
+     return file->getBoloMuV2T(index, imes);
+   case boloRes:
+     return file->getBoloRes(index, imes);
+   case boloTemp:
+     return file->getBoloTemp(index, imes);
+   case dilDAC:
+     return file->getADCDil(index);
+   case dilSwitch:
+     return file->getSwitchDil();
+   case sstDiode:
+     return file->getSSTSignal(index, imes);
+   case sstChannel:
+     return file->getSSTRawSignal(index, imes);
+   case sstDiodeCN:
+     return file->getSSTSignalCN(index, imes);
+   case sstChannelCN:
+     return file->getSSTRawSignalCN(index, imes);
+   case sstStarCnt:
+     return file->getNumbStar(imes);
+   case sstStarZ:
+     return file->getSSTStarZ(index, imes);
+   case sstStarF:
+     return file->getSSTStarF(index, imes);
+   case sstStarT:
+     return file->getSSTStarT(index, imes);
+   case gyroRaw:
+     return file->getGyroRaw(index, imes);
+   case gyroV:
+     return file->getGyroTens(index, imes);
+   case gyroSpeed:
+     return file->getGyroSpeed(index, imes);
+   case gpsTime:
+     return file->getGPSUTC();
+   case longitude:
+     return file->getGPSLong();
+   case latitude:
+     return file->getGPSLat();
+   case altitude:
+     return file->getGPSAlt();
+   case tsid:
+     tSid.setLongitude(file->getGPSLong());
+     return tSid.getLST(getMJD());
+   case azimuthBolo:
+   case alphaRotAxis:
+   case deltaRotAxis:
+   case alphaFPAxis:
+   case deltaFPAxis:
+   case alphaZenith:
+   case deltaZenith:
+   case alphaBolo:
+   case deltaBolo: {
+        double lat, lon;
+        if (auxGPS) {
+          double alti;
+          auxGPS->getLocation(getMJD(), lat, lon, alti);
+        } else {
+          lon = file->getGPSLong();
+          lat = file->getGPSLat();
+        }
+        tSid.setLongitude(lon);
+        double ts = tSid.getLST(getMJD());
+        tempLocator.setEarthPos(lon, lat);
+        tempLocator.setTSid(ts);
+        if ((kind) == alphaRotAxis) return tempLocator.getAlphaZenith();
+        if ((kind) == deltaRotAxis) return tempLocator.getDeltaZenith();
+        if ((kind) == alphaZenith) return tempLocator.getAlphaZenith();
+        if ((kind) == deltaZenith) return tempLocator.getDeltaZenith();
+        if ((kind) == azimuthBolo) return tempLocator.getAzimutBolo(getSampleIndex(),index);
+        if ((kind) == alphaBolo) return tempLocator.getAlphaBolo(getSampleIndex(),index);
+        if ((kind) == deltaBolo) return tempLocator.getDeltaBolo(getSampleIndex(),index);
+        if ((kind) == alphaFPAxis) return tempLocator.getAlphaCenter(getSampleIndex());
+        if ((kind) == deltaFPAxis) return tempLocator.getDeltaCenter(getSampleIndex());
+        return 0;
+     }
+   case voyantEVO:
+        return voyant_EVO(file->lastDilution());
+   case voyantEVF:
+        return voyant_EVF(file->lastDilution());
+   case commandeEV0:
+        return commande_EVO(file->lastDilution());
+   case commandeEVF:
+        return commande_EVF(file->lastDilution());
+   case commandeEVB:
+        return commande_EVB(file->lastDilution());
+   case commandeEVV:
+        return commande_EVV(file->lastDilution());
+   case pressEnt3He:
+        return pression_entree_3He(file->lastDilution());
+   case debit3He:
+        return debit_3He(file->lastDilution());
+   case pressSor3He:
+        return pression_sortie_3He(file->lastDilution());
+   case pressEnt4He:
+        return pression_entree_4He(file->lastDilution());
+   case debit4He:
+        return debit_4He(file->lastDilution());
+   case pressSor4He:
+        return pression_sortie_4He(file->lastDilution());
+   case pressAirVanne:
+        return pression_air_vanne(file->lastDilution());
+   case pressPompChar:
+        return pression_pompe_charbon(file->lastDilution());
+   case pressMembrane:
+        return pression_membranne(file->lastDilution());
+   case pressExterne:
+        return pression_externe(file->lastDilution());
+   case tensPile10T:
+        return tension_pile_10T(file->lastDilution());
+   case tensPileP18D:
+        return tension_pile_p18D(file->lastDilution());
+   case tensPileM18D:
+        return tension_pile_m18D(file->lastDilution());
+   case tensPile10B:
+        return tension_pile_10B(file->lastDilution());
+   case tensPileP18B:
+        return tension_pile_p18B(file->lastDilution());
+   case tensPileM18B:
+        return tension_pile_m18B(file->lastDilution());
+   case tensPileCh:
+        return tension_pile_Ch(file->lastDilution());
+   case swPile5:
+        return switch_pile_5(file->lastDilution());
+   case swPile15:
+        return switch_pile_15(file->lastDilution());
+   case tempCaissH1:
+        return temperature_caisson_haut1(file->lastDilution());
+   case tempCaissH2:
+        return temperature_caisson_haut2(file->lastDilution());
+   case tempCaissB1:
+        return temperature_caisson_bas1(file->lastDilution());
+   case tempCaissB2:
+        return temperature_caisson_bas2(file->lastDilution());
+   case tempCaissTHe:
+        return temperature_caisson_tube_helium(file->lastDilution());
+   case tempCaissPiles:
+        return temperature_caisson_piles(file->lastDilution());
+   case tempCaissDrv:
+        return temperature_caisson_driver_moteur(file->lastDilution());
+   case pressHeBain:
+        return pression_helium_bain(file->lastDilution());
+   case pressPirani:
+        return pression_pirani(file->lastDilution());
+   }
+   return -1;
+}
+bool   TOIIter::newValue(int column) {
+  if (!initDone) Init();
+   if (column < 0 || column >= infos.size()) return false;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS &&
+       (kind == longitude || kind == latitude || kind == altitude)) {
+      return true;
+   }
+   if (rawIter) {
+       return interp[column].isNewValue(lastSample);
+   }
+   switch (kind) {
+     case sampleNum:
+     case internalTime:
+     case mjd:
+     case mutc:
+     case tsid:
+       return true;
+     case boloMuV:
+     case boloMuV2:
+     case boloMuV2T:
+     case boloRes:
+     case boloTemp:
+       return file->blockNum() == file->getBoloBlockNum();
+     case boloRawMuV:
+     case boloRawMuVCN:
+       return file->blockNum() == file->getBoloBlockNum();
+     case boloGainAmpli:
+       return file->blockNum() == file->getReglageBlockNum() && imes==0;
+     case boloDACI:
+     case boloDACV:
+       return file->blockNum() == file->getReglageBlockNum() && imes==0;
+     case dilDAC:
+     case dilSwitch:
+       return file->blockNum() == file->getDilutionBlockNum() && imes==0;
+     case sstChannel:
+     case sstDiode:
+     case sstStarCnt:
+     case sstStarZ:
+     case sstStarF:
+     case sstStarT:
+       return file->blockNum() == file->getSSTBlockNum();
+     case sstChannelCN:
+     case sstDiodeCN:
+       return file->blockNum() == file->getSSTCompBlockNum();
+     case gyroRaw:
+     case gyroV:
+     case gyroSpeed:
+       return file->blockNum() == file->getGyroBlockNum();
+     case gpsTime:
+       return file->blockNum() == file->getGPSBlockNum() && imes==0;
+     case longitude:
+       return file->blockNum() == file->getGPSBlockNum() && imes==0;
+     case latitude:
+       return file->blockNum() == file->getGPSBlockNum() && imes==0;
+     case altitude:
+       return file->blockNum() == file->getGPSBlockNum() && imes==0;
+     case azimuthBolo:
+     case alphaBolo:
+     case deltaBolo:
+     case alphaZenith:
+     case deltaZenith:
+     case alphaFPAxis:
+     case deltaFPAxis:
+       return true; // $CHECK$ with SSTHandler
+     case alphaRotAxis:
+       return true; // $CHECK$ with SSTHandler
+     case deltaRotAxis:
+       return true; // $CHECK$ with SSTHandler
+     case voyantEVO:
+     case voyantEVF:
+     case commandeEV0:
+     case commandeEVF:
+     case commandeEVB:
+     case commandeEVV:
+     case pressEnt3He:
+     case debit3He:
+     case pressSor3He:
+     case pressEnt4He:
+     case debit4He:
+     case pressSor4He:
+     case pressAirVanne:
+     case pressPompChar:
+     case pressMembrane:
+     case pressExterne:
+     case tensPile10T:
+     case tensPileP18D:
+     case tensPileM18D:
+     case tensPile10B:
+     case tensPileP18B:
+     case tensPileM18B:
+     case tensPileCh:
+     case swPile5:
+     case swPile15:
+     case tempCaissH1:
+     case tempCaissH2:
+     case tempCaissB1:
+     case tempCaissB2:
+     case tempCaissTHe:
+     case tempCaissPiles:
+     case tempCaissDrv:
+     case pressHeBain:
+     case pressPirani:
+       return file->blockNum() == file->getDilutionBlockNum() && imes==0;
+   }
+   return false;
+}
+bool   TOIIter::extendValue(int column) {
+   return (!infos[column].interpolated && !newValue(column));
+}
+bool   TOIIter::interpValue(int column) {
+   return (infos[column].interpolated && !newValue(column));
+}
+bool TOIIter::isTrig(int column) {
+   if (column < 0 || column >= infos.size()) return false;
+   return infos[column].triggering;
+}
+TOIKind TOIIter::getKind(int column) {
+   if (column < 0 || column >= infos.size()) return (TOIKind)-1;
+   return infos[column].kind;
+}
+int TOIIter::getIndex(int column) {
+   if (column < 0 || column >= infos.size()) return (TOIKind)-1;
+   return infos[column].index;
+}
+int TOIIter::getColTOI(TOIKind kind, int index) {
+   for (int i=0; i<infos.size(); i++)
+     if (infos[i].kind == kind && infos[i].index == index)
+       return i;
+   return -1;
+}
+bool TOIIter::canGetTOI(TOIKind kind, int index) {
+   int col = getColTOI(kind, index);
+   if (col<0) return false;
+   return canGetValue(col);
+}
+double TOIIter::getTOI(TOIKind kind, int index) {
+   int col = getColTOI(kind, index);
+   if (col<0) return -9.e99;
+   return getValue(col);
+}
+int TOIIter::getBlockSampleIndex() {
+   return imes;
+}
+int TOIIter::getUnderSampling() {
+  return underSample;
+}
+int TOIIter::getSampleIndex() {
+  if (!initDone) Init();
+   if (file) {
+     return file->blockNum() * file->nEchBlock() + imes;
+   } else {
+     return lastSample;
+   }
+}
+double TOIIter::getMJD() {
+  if (!initDone) Init();
+  int sample = getSampleIndex();
+  return archParam.acq.tBlock0 + sample*archParam.acq.perEch/86400.;
+}
+bool TOIIter::fetchAhead() { // Seulement si delegation
+  if (!rawIter) return false;
+  if (!rawIter->Next()) return false;
+  long sample = rawIter->getSampleIndex();
+  for (int i=0; i<infos.size(); i++) {
+    if (rawIter->canGetValue(i) && rawIter->newValue(i)) {
+      interp[i].enterValue(rawIter->getValue(i), sample);
+    }
+  }
+  return true;
+}
+block_type_param* TOIIter::lastParam() {
+  if (!initDone) Init();
+  if (file) return file->lastParam();
+  else return rawIter->lastParam();
+}
+   TOIInfo& info = request[column];
+   return(info.second->getValue(curSample, info.first));
+}
+bool TOIIter::canGetValue(TOI const& toi) {
+  return canGetValue(getColTOI(toi));
+}
+double TOIIter::getValue(TOI const& toi) {
+  return getValue(getColTOI(toi));
+}

trunk/Poubelle/archTOI.old/toiiter.h

-              r436
+              r534
 // toiiter.h
 // Eric Aubourg         CEA/DAPNIA/SPP   juillet 1999
+// Eric Aubourg         CEA/DAPNIA/SPP   septembre 1999
 #ifndef TOIITER_H
 #define TOIITER_H
 #include <vector>
+#include <list>
 #include <set>
 #include <map>
 #include <string>
-#include "archeopsfile.h"
-#include "auxgps.h"
-#include "tsid.h"
+enum TOIKind {
+  sampleNum,     // Numero d''echantillon (timer interne transputer)
+  internalTime,  // temps transputer depuis debut manip, en secondes.
+  mjd,           // UTC (en jours, MJD = JD - 2450000).
+  mutc,          // UTC en heure, avec soustraction d'origine
+  boloMuV,       // tension en microVolts, filtree avec filtre carre indexe par bolometre
+  boloMuV2,      // tension en microVolts, avec soustraction offset filtre
+  boloRawMuV,    // tension brute en microVolts, non filtree indexe par bolometre, sans chgt signe
+  boloTemp,      // Kelvins, need R(T) -> N/A
+  boloDACV,      // tension de reference du DAC (microVolts)
+  boloDACI,      // courant (microAmps)
+  boloRawMuVCN,  // compression noise on bolo (microVolts), no filter
+  boloMuV2T,     // total, microVolts, avec consigne DAC
+  boloRes,       // Resistance, ohms
+  boloGainAmpli, //
+  dilDAC,        // dilution DAC, index = 0-47
+  dilSwitch,     // dilution, mot de 32 bits
+  serviceTemp,   // service temperature, K
+  sstDiode,      // signal brut SST, indice = diode, 0-45
+  sstChannel,    // signal brut SST, indice = canal, 0-47
+  sstDiodeCN,    // compression noise
+  sstChannelCN,  // compression noise
+  sstStarCnt,    // nombre d''etoiles detectees durant ce sample
+  sstStarZ,      // numero de diode d''une eventuelle etoile,   i=etoile 0..n
+  sstStarF,      // flux d''une eventuelle etoile,             <0 si moins de i+1 etoiles
+  sstStarT,      // temps d''une eventuelle etoile,            <0 si moins de i+1 etoiles
+  gyroRaw,       // valeur brute du gyro, index = 0,1,2
+  gyroV,         // tension gyro, V, index = 0,1,2
+  gyroSpeed,     // vitesse gyro, deg/s, index = 0,1,2, calibree...
+  gpsTime,       // temps fourni par GPS
+  longitude,     // position ballon, degrees  +=EST
+  latitude,      // position ballon, degrees  +=NORD
+  altitude,      // position ballon, m
+  tsid,          // temps sideral en secondes
+  azimuthBolo,   // azimut
+  alphaRotAxis,  //
+  deltaRotAxis,  //
+  alphaZenith,   //
+  deltaZenith,   //
+  alphaSst,      //
+  deltaSst,      //
+  alphaFPAxis,   //
+  deltaFPAxis,   //
+  alphaBolo,     // sur le ciel, indexe par bolometre
+  deltaBolo,     // sur le ciel, indexe par bolometre
+  voyantEVO,
+  voyantEVF,
+  commandeEV0,
+  commandeEVF,
+  commandeEVB,
+  commandeEVV,
+  pressEnt3He,
+  debit3He,
+  pressSor3He,
+  pressEnt4He,
+  debit4He,
+  pressSor4He,
+  pressAirVanne,
+  pressPompChar,
+  pressMembrane,
+  pressExterne,
+  tensPile10T,
+  tensPileP18D,
+  tensPileM18D,
+  tensPile10B,
+  tensPileP18B,
+  tensPileM18B,
+  tensPileCh,
+  swPile5,
+  swPile15,
+  tempCaissH1,
+  tempCaissH2,
+  tempCaissB1,
+  tempCaissB2,
+  tempCaissTHe,
+  tempCaissPiles,
+  tempCaissDrv,
+  pressHeBain,
+  pressPirani
+};
+#include "toiabsorber.h"
+#include "requesthandler.h"
+#include "archfileset.h"
 class TOIInterpolator;
+class AuxGPS;
+class TOIIter {
+template <class T1, class T2, class T3>
+struct triple {
+    T1 first;
+    T2 second;
+    T3 third;
+        triple() : first(T1()), second(T2()), third(T3()) {}
+    triple (const T1& a, const T2& b, const T3& c) :
+           first (a), second (b), third (c) {}
+};
+class TOIIter : public TOIAbsorber /*, public RequestHandler*/ {
 public:
+   TOIIter(TOIIter const&);   // $CHECK$ prevoir operator = ?
+   ~TOIIter();
+   bool Next();
+  TOIIter();
+  // ~TOIIter();
+  void readReq(istream& s);
+  void defaultInclude();
+  bool processRequest(string line);
+  void registerReqHandler(RequestHandler*);
+  virtual bool processTOIReq(TOI& toi, string);
+  virtual bool processOption(string keyw, string args);
+  void addDirectory(string);
+  void addFile(string);
+  //void useAuxGPS(AuxGPS* gps);
+   bool    canGetValue(int column);  // data is available for that...
+   double  getValue(int column);
+   bool    newValue(int column);     // a juste change ?
+   bool    extendValue(int column);  // une valeur plus ancienne, etendue ?
+   bool    interpValue(int column);  // une valeur interpolee avec valeur future ?
+   bool    isTrig(int column);
+   TOIKind getKind(int column);
+   int     getIndex(int column);     // si plusieurs, exemple, tension bolo 2
+  void setMJDInterval(double tStart, double tEnd);
+  void setUTCInterval(double tStart, double tEnd);
+  void setSNInterval(long sStart, long sEnd);
+  void setUnderSample(int n);
+  int  getUnderSample();
+  void addTOI(TOI& toi, bool triggering=true);
+  void init();
+  bool next();
+   bool    canGetTOI(TOIKind kind, int index=0);
+   double  getTOI(TOIKind kind, int index=0);
+  long    getSampleNum();
+  bool    canGetValue(int column);  // data is available for that...
+  double  getValue(int column);
+  bool    canGetValue(TOI const& toi);  // data is available for that...
+  double  getValue(TOI const& toi);
+  bool    newValue(int column);     // a juste change ?
+  bool    extendValue(int column);  // une valeur plus ancienne, etendue ?
+  bool    interpValue(int column);  // une valeur interpolee avec valeur future ?
+  bool    isTrig(int column);
+  TOI     getKind(int column);
+  ArchFileSet& getFSet() {return fset;}
+protected:
+  ArchFileSet fset;
+  list<RequestHandler*> handlers;
+  enum {triggering = 1};
+  typedef ::triple<TOI,TOIProducer*,long> TOIInfo;
+  vector<TOIInfo> request;
+  double mjdStart, mjdEnd;  // MJD
+  double utcStart, utcEnd;  // UTC, will be converted towards tStart tEnd at init
+  long   sStart, sEnd;      // samplenum
+  int underSample;
+  long curSample;
+   int    getBlockSampleIndex();        // numero d''echantillon dans dernier bloc bolo
+   int    getSampleIndex();             // numero d''echantillon
+   int    getUnderSampling();
+   double getMJD();
+  int getColTOI(TOI const& toi);
+  bool next1();
+   block_type_param*         lastParam();
+   block_type_journal*       lastJournal();
+   block_type_reglage*       lastReglage();
+   block_type_dilution*      lastDilution();
+   block_type_gps*           lastGPS();
+   block_type_une_periode*   lastUnePeriode();
+   block_type_synchro_sol*   lastSynchroSol();
+   block_type_pointage_sol*  lastPointageSol();
+   block_type_bolo*          lastBolo();
+   block_type_gyro*          lastGyro();
+   block_type_sst*           lastSST();
+   block_type_bolo_comprime* lastBoloComp();
+   block_type_gyro_comprime* lastGyroComp();
+   block_type_sst_comprime*  lastSSTComp();
+   ArcheopsFile* currentFile() {if (!initDone) Init(); return file;}
+protected:
+   TOIIter();
+   void Init(); // After setting the options, opens the first file.
+   void ScanFiles(); // Parcourt tous les fichiers, cherche le premier echantillon...
+   bool Next1();
+   bool NextSample(); // no test of trange
+   bool NextFile();
+   ArcheopsFile* file;
+   string directory;   // If directory scanning
+   set<string> files;  // list of files to be read...
+   map<string, long> file1stSamp; // premier echantillon de chaque fichier
+   set<string>::iterator curFile;
+   bool  isOnBoardRecorder;
+   int imes;
+   double tStart, tEnd;      // MJD
+   double utcStart, utcEnd;  // UTC, will be converted towards tStart tEnd at init
+   long   sStart, sEnd;      // samplenum
+   //double tBlock0;  // MJD du bloc numero zero...
+   //double perEch;   // periode d''echantillonage en secondes
+   struct info {
+     TOIKind kind;
+     int index;
+     bool triggering;
+     bool interpolated;
+   };
+   int getColTOI(TOIKind kind, int index=0);
+   long trigMask;
+   vector<info>   infos;
+   friend class TOISvr;
+   // Gestion de l''interpolation, par delegation a un autre iterateur
+   TOIIter* rawIter; // iterateur sans interpolation, qui avance avec retard
+                     // en bufferisant les valeurs. S''il existe, on lui delegue
+                     // toutes les operations agissant sur le fichier.
+   TOIInterpolator* interp; // interpolateurs pour stocker les valeurs
+   long lastSample;  // last interpolated sample
+   long maxLookAhead;
+   bool fetchAhead(); // avance d''une TOI en nourrissant les interpolateurs
+   AuxGPS* auxGPS;
+   TSid tSid;
+   int underSample;
+   private:
+   bool initDone;
+private:
+  bool initDone;
+  bool incDone;
 };

trunk/Poubelle/archTOI.old/toisvr.cc

-              r465
+              r534
     else tsttoi(alphaBolo)
     else tsttoi(deltaBolo)
-    else tsttoi(voyantEVO)
-    else tsttoi(voyantEVF)
-    else tsttoi(commandeEV0)
-    else tsttoi(commandeEVF)
-    else tsttoi(commandeEVB)
-    else tsttoi(commandeEVV)
-    else tsttoi(pressEnt3He)
-    else tsttoi(debit3He)
-    else tsttoi(pressSor3He)
-    else tsttoi(pressEnt4He)
-    else tsttoi(debit4He)
-    else tsttoi(pressSor4He)
-    else tsttoi(pressAirVanne)
-    else tsttoi(pressPompChar)
-    else tsttoi(pressMembrane)
-    else tsttoi(pressExterne)
-    else tsttoi(tensPile10T)
-    else tsttoi(tensPileP18D)
-    else tsttoi(tensPileM18D)
-    else tsttoi(tensPile10B)
-    else tsttoi(tensPileP18B)
-    else tsttoi(tensPileM18B)
-    else tsttoi(tensPileCh)
-    else tsttoi(swPile5)
-    else tsttoi(swPile15)
-    else tsttoi(tempCaissH1)
-    else tsttoi(tempCaissH2)
-    else tsttoi(tempCaissB1)
-    else tsttoi(tempCaissB2)
-    else tsttoi(tempCaissTHe)
-    else tsttoi(tempCaissPiles)
-    else tsttoi(tempCaissDrv)
-    else tsttoi(pressHeBain)
-    else tsttoi(pressPirani)
     else {
       // cerr << "*Warning, unrecognized TOI " << line << endl;

trunk/Poubelle/archTOI.old/transfelec.cc

r394	r534
1	1	#include "transfelec.h"
2	2
	3	// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
3	4
4	5	// Parametres du circuit electrique de lecture SST

trunk/Poubelle/archTOI.old/transfelec.h

-              r394
+              r534
 //#include <math.h>
+// D. Yvon, CE Saclay, DAPNIA/SPP, 08/99
 #ifdef __MWERKS__

Context Navigation

Changeset 534 in Sophya for trunk

Legend:

trunk/Poubelle/archTOI.old/archeopsfile.cc

trunk/Poubelle/archTOI.old/archeopsfile.h

trunk/Poubelle/archTOI.old/archparam.cc

trunk/Poubelle/archTOI.old/archparam.h

trunk/Poubelle/archTOI.old/archtoi.cc

trunk/Poubelle/archTOI.old/archtoi.h

trunk/Poubelle/archTOI.old/archtoimain.cc

trunk/Poubelle/archTOI.old/archvers.h

trunk/Poubelle/archTOI.old/auxgps.h

trunk/Poubelle/archTOI.old/auxinterpgps.cc

trunk/Poubelle/archTOI.old/auxinterpgps.h

trunk/Poubelle/archTOI.old/dyffttools.cc

trunk/Poubelle/archTOI.old/dyffttools.h

trunk/Poubelle/archTOI.old/formepulse.cc

trunk/Poubelle/archTOI.old/formepulse.h

trunk/Poubelle/archTOI.old/pisteetoile.cc

trunk/Poubelle/archTOI.old/pisteetoile.h

trunk/Poubelle/archTOI.old/sstetoile.cc

trunk/Poubelle/archTOI.old/sstetoile.h

trunk/Poubelle/archTOI.old/ssthandler.cc

trunk/Poubelle/archTOI.old/ssthandler.h

trunk/Poubelle/archTOI.old/ssthardware.h

trunk/Poubelle/archTOI.old/templocator.cc

trunk/Poubelle/archTOI.old/templocator.h

trunk/Poubelle/archTOI.old/toiiter.cc

trunk/Poubelle/archTOI.old/toiiter.h

trunk/Poubelle/archTOI.old/toisvr.cc

trunk/Poubelle/archTOI.old/transfelec.cc

trunk/Poubelle/archTOI.old/transfelec.h

Download in other formats: