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

Timestamp:

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

Author:

ansari

Message:

V2

File:

• trunk/Poubelle/archTOI.old/pisteetoile.cc (modified) (8 diffs)

trunk/Poubelle/archTOI.old/pisteetoile.cc

@@ -8 +8 @@
 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.;
@@ -39 +38 @@
 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;
@@ -115 +113 @@
 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);
@@ -152 +150 @@
         
                 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;
+*/
                 }
         }
@@ -199 +204 @@
 // 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.
         
@@ -281 +288 @@
 // 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();
@@ -292 +297 @@
                 // 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;
@@ -314 +320 @@
         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