Changeset 645 in Sophya for trunk/Poubelle/archediab.old/archediab.sources/c

Timestamp:

Nov 25, 1999, 2:35:30 PM (26 years ago)

Author:

ansari

Message:

archediab 27

Location:

trunk/Poubelle/archediab.old/archediab.sources/c

Files:

• bit_block.c (modified) (6 diffs)
• dilution.c (modified) (3 diffs)
• tache_archeops_bit2.c (modified) (4 diffs)
• telecommande.c (modified) (1 diff)

trunk/Poubelle/archediab.old/archediab.sources/c/bit_block.c

@@ -9 +9 @@
 #include "archeops.h"
 #include "bit_block.h"
+#include "choix_param.h"
+#include "structure.h"
+#include "tm.h"
+#include "tache.h"
+#include "bolo.h"
+#include "carte_acqui.h"
+#include "carte_pci.h"
+ 
 
-
-
+#define cadence_taux    16
+#define mode_acquisition        0x00    //0xff en mode inverse, 0x00 en mode normal
 unsigned long   buffer;
 int     llgg,somme,type_du_block,numero;
 int     deb;    /*deb vaut zero si aucun debut n'a encore ete trouve*/
 int     compteur;
-int cpt_cherche_debut;
+int     cpt_cherche_debut;
+unsigned long   taux;
+unsigned long   Nvide, Nplein;
 def_long_block
 
@@ -25 +35 @@
 {
         int     i,m,bit;
+        
 //      unsigned char   bit;    /*c'est le bit lu*/
         
-        bit=un_bit();
+        bit=un_bit()^mode_acquisition;
         switch(bit)
                 {
@@ -38 +49 @@
                 }
         
+        
         compteur++;
         i=compteur%32;
 
-                
+        
+        
         
         switch(deb)
@@ -108 +121 @@
                                                 if(pbk->mot[m] == fin_block_mesure)
                                                 {
+                                                        /*
                                                         if(calcul_somme_block(pbk,llgg)!=somme)
                                                         {
                                                                 ecrit_erreur_bit(erreur_somme);
                                                                 
-                                                                /*on continue quand meme*/
+                                                                
                                                         }
+                                                        */
+                                                        
 //                                                      else    ecrit_erreur_bit(bloc_correct); 
                                                                 
                                                         deb=0;
-
+                                                        Nplein=compteur+1;
+                                                        compteur=0;
+                                                        //Calcul du taux de remplissage
+                                                        taux=taux+(((0xffff-taux)*Nplein)>>cadence_taux);
+                                                        tt->nb_lec_fofo_ext=((taux*1000)>>cadence_taux);
                                                         return(lit_bit_un_block);
                                                                                 
@@ -129 +149 @@
                                                 }
                                                 
+                                        
+                                                
+                                                
                                                 
                                         }
-                                }               
+                                        
+                                }
+                                
+                                                
                         }
                         break;
@@ -138 +164 @@
                 default :
                 
-                                
+                        if(compteur>1000)
+                        {
+                                Nvide=compteur;
+                                compteur=0;
+                                taux=taux-((taux*Nvide)>>cadence_taux);
+                                tt->nb_lec_fofo_ext=((taux*1000)>>cadence_taux);
+                        }       
                         if(buffer==debut_block_mesure)
                         {
                                 deb=1;
+                                Nvide=compteur-32;
                                 compteur=31;
                                 pbk->debut=buffer;

trunk/Poubelle/archediab.old/archediab.sources/c/dilution.c

@@ -208 +208 @@
 
 
+#define c(j,i)          (1e-4*(double)parametr.nom_coef[parametr.bolo[j].numero_nom_coef].coef[i])      
 
 // les temperatures sur les cartes modifiées
 {
-double I,V,R[4];
+double I,V,R[8],llR;
 int j,k;
 def_gains;
 k=0;
+
 for(j=0;(j<nb_max_bolo) && (k<4);j++)
         {
@@ -222 +224 @@
                 V=0.001*bol_micro_volt(blk->temperature[k],(double)parametr.bolo[j].bolo_gain*gain_ampli(gg->reglage.bolo[j]));
                 if(I>0.0000001)         R[k]=V/I;       else    R[k]=0;                         //      R  en ‡
-                k++;
+ //------  calcul de l'etalonnage en temperature des cartes temperature simplifiées  
+ // j=numero de bebo    T[k] = R[k+4] en Kelvin
+                if ((R[k]-c(j,6))>1.) 
+                        {if((log(R[k]-c(j,6))-c(j,0))>0.001) llR= log(log(R[k]-c(j,6))-c(j,0)) ; else llR=0;} 
+                else llR=0;
+                R[k+4]  = exp(  c(j,1) + c(j,2)* llR + c(j,3)* llR* llR + c(j,4)* llR* llR* llR + c(j,5)* llR* llR* llR* llR) ; 
+                if(R[k+4]>9999) R[k+4]=9999;
+                k++;            
                 }
         }
 
-if(R[0]<90.)    ecritD(fenetre_dilution,dil_temp,"germ= %6.2f‡  4K=%5.1f‡   1.6K=%5.1f‡   10K=%5.1f‡",R[0],R[1],R[2],R[3]);
-else    ecritD(fenetre_dilution,dil_temp,"germ= %6.3fk‡  4K=%5.1f‡   1.6K=%5.1f‡   10K=%5.1f‡",R[0]*0.001,R[1],R[2],R[3]);
+
+if(R[0]<90.)    ecritD(fenetre_dilution,dil_temp
+        , "ge=%6.2f‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fK  4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
+                ,R[0],R[1],R[2],R[3],R[4],R[5],R[6],R[7]);
+else    ecritD(fenetre_dilution,dil_temp
+        ,"ge=%6.3fk‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fmK 4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
+        ,R[0]*0.001,R[1],R[2],R[3],R[4]*1000.,R[5],R[6],R[7]);
+
+
 if(litD(fenetre_dilution,dil_trace_tempe,0))// trace les sondes de niveau helium
         {
@@ -240 +256 @@
                         strcpy(graph->xtitre,"minutes");
                         }
-        tracen(fenetre_temperature_dil,4,minutes,R);
+        tracen(fenetre_temperature_dil,8,minutes,R);
         }
 }

trunk/Poubelle/archediab.old/archediab.sources/c/tache_archeops_bit2.c

@@ -23 +23 @@
 //
 //
-//   lit directement les blocks recut sur la fifo et les envoie au programme principal
+//   lit directement les blocks recus sur la fifo et les envoie au programme principal
 //
 
@@ -34 +34 @@
 
 //*****************************************************************************************
+unsigned long   N_total_bits,M;
+unsigned short  Max,Min;
 
 void init_tache(void)
@@ -69 +71 @@
 
 
-char    un_bit(void)
+char    un_bit()
 {
 static int i;
 static unsigned long RR;
+unsigned short  Val;
+
 lit_carte;      if(rien_a_lire)         return(bit_vide);       //  fifo vide  ->  retour  ->  sort du while
+
+N_total_bits++;
+Val=((RR>>1) & 0x3f);   //permet de lire les bits de temps dans RR
+
+//Calcul de la nouvelle moyenne glissee du temps
+M=M-(M>>16)+(Val<<8);
+
+        
+//Calcul du Max
+if(N_total_bits==0) Max=Val;    //initialisation au premier passage
+if(Val>Max) Max=Val;    //changement de maximum
+if((N_total_bits%60000)==0)
+{
+        Max=Max-1;      //Diminution tous les 60000 bits
+        N_total_bits=0;
+}
+        
+//Calcul du Min
+if(N_total_bits==0) Min=Val;    //initialisation au premier passage
+if(Val<Min) Min=Val;    //changement de minimum
+if((N_total_bits%60000)==0) Min=Min+1;  //Augmentation tous les 60000 bits
+//tt->nb_lec_fofo_ext=M;
+
+
+        
 //if( RR & 1 )          return(bit_un);
 //else          return(bit_zero);
 if( RR & 1 )    return(bit_zero);
 else            return(bit_un);
+
 }
 
@@ -135 +165 @@
         }
 if(lec_fifo > tt->nb_lec_fofo_int)      tt->nb_lec_fofo_int=lec_fifo;
-if(!tt->nb_lec_fofo_ext)                {tt->nb_lec_fofo_ext=tt->nb_lec_fofo_int;tt->nb_lec_fofo_int=0;}
+//if(!tt->nb_lec_fofo_ext)              {tt->nb_lec_fofo_ext=tt->nb_lec_fofo_int;tt->nb_lec_fofo_int=0;}
 
 }

trunk/Poubelle/archediab.old/archediab.sources/c/telecommande.c

@@ -217 +217 @@
         case    tc_dir_reset_epld       :       //printf("tc_reset_epld:\n");
                 
-                        for(j=0;j<14;j++) tc_char(j);
+//                      for(j=0;j<14;j++) tc_char(j);
+                        for(j=0;j<15;j++) tc_char(j);   // reset epld+transputer
                         tc_char(1);tc_char(1);
                         tc_char(0x00);tc_char(0x00);    //  synchro