Changeset 534 in Sophya for trunk/Poubelle/archTOI.old/ssthandler.cc

Timestamp:

Nov 1, 1999, 10:59:29 AM (26 years ago)

Author:

ansari

Message:

V2

File:

• trunk/Poubelle/archTOI.old/ssthandler.cc (modified) (9 diffs)

trunk/Poubelle/archTOI.old/ssthandler.cc

@@ -12 +12 @@
 
 #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
@@ -180 +180 @@
 }
 
+
+// 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;
@@ -224 +237 @@
                                         //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++;
@@ -257 +272 @@
                 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;
@@ -266 +280 @@
                         (*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
 
@@ -276 +301 @@
         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;
@@ -301 +324 @@
                 }
 */                      
-                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()) {
@@ -317 +340 @@
         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
 
@@ -377 +400 @@
 {
   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) {