source: Sophya/trunk/Poubelle/archTOI.old/pisteetoile.h@ 614

#ifndef PisteEtoileSeen
#define PisteEtoileSeen
#include "sstetoile.h"
#include "transfelec.h"
#include "formepulse.h"

/*      // 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 48 samples. The 32 first are considered as prepulse sample.
        On the prepulse: 
                We compute baseline and substract it. 
                Then compute slope due to pile up and substract it
                Then compute Vrms.
                Look for trigger
        If triggered :
                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. 
                Fit parabole arond the minimum of X2 to get the true minimum.
                Build SSTStar.
                
        That's it. 
        When  done push the track eigth samples forward and iterate analysis. 
*/

class FormePulse;

#define PEtoileDebug


#define PhDiodTabLong (48)
#define Prepulselong (32)
#define Pousslong (8)


//      Parametrisation du Trigger 
#define SeuilGachette (6.)                              // On déclenche à SeuilGachette sigmas
#define SSTLsbVal (10./4096.)

//      Parametrisation Analyse Fine
#define TFExcursion (15.e-3)                    // Secondes: Excursion du fit autour de MaxIndex
#define TPasTFit (1.5e-3)                               // Secondes:
                        // Secondes:

class PisteEtoile {
public:
        PisteEtoile(short NoPhotDiod=0);
        PisteEtoile (PisteEtoile const & x);
        ~PisteEtoile();
        PisteEtoile& operator = (PisteEtoile const & x);
        void fill(int* pData, int inDebutPiste=0, int nData=PhDiodTabLong);
        void push(int* pData, int nbdata=Pousslong); 
                //"Pousse" les données dans PhDiodArray de nbdata
        bool traque() ; 
                // Traque l'etoile.     
        bool trig(); 
                //Trig sur les donnees? Si Oui declenche une analyse complète
        void demasque();
                // Met en Oeuvre un fit de forme pour determiner le temps et l'amplitude.

        bool isAlive()  {       return  DiodVivante;}
        void kill()             {       DiodVivante=false;}
        void restore()  {       DiodVivante=true;}
        
        SSTEtoile DonneEtoile() {return LastEtoile;}
                //Retourne la dernière Etoile détectee
        inline int SSTLargEvtIndex(); 
        inline int SSTLongIndexEvt();
        static TransFuncElec Trapani;
        static FormePulse FPulse;
        
protected:
        double PasTempsEch;     //Secondes 
// Idendificateur, repères et gestion 
        short NoPhDiod;
        long indexDebutPiste;   // A actualiser dans push() et Fill()
        bool DiodVivante;               // Si false, On ne fait rien 
        
// Stocke donnees brutes et traitees
        int* PhDiodArray;
        double* PhDiodPiste;
        
// Derniere Etoile Vue
        SSTEtoile LastEtoile;
        
// Variables d'analyse
        bool gachette;
        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 AmplFin;
        double Tfin;
        double X2Min;
        bool SuccesAFine;
        
// Utilitaires de calcul analyse fine
        double AmplNorm;
                // Normalisation d'amplitude sur FPulse dans l'ajustement iteratif
        double *pFitArray0;
                // Pointeur origine du tableau d'echantillon sur lequel sera calculé le fit
        double *pFitArray;
                // Pointeur courant du tableau d'echantillon 
        double* Ampl;   //Tableau des amplitudes calculées

        double* X2;             //Tableau des X2 calculés
        

        
};

#endif