RawAmpWatcher.h

Go to the documentation of this file.

#ifndef RAWAMPWATCHER_H
#define RAWAMPWATCHER_H

#include "TCanvas.h"
#include "TSpectrum.h"
#include "TTree.h"

#include "RecalEnergy.h"
#include "FitFunctions.h"

#include "SpectraViewerWatcher.h"

class BashColor;

class RawAmpWatcher : public SpectraViewerWatcher
{
public:

    double TabArea[38];
    double TabOffset[38];
    double TabGain[38];
    double TabMean[38];
    double TabE[38];
    double TabFWHM[38];
    double TabFWHM_10[38];
    double TabLTProp[38];
    double TabRTProp[38];

    TCanvas *fCanvas;

    double fWidth;
    double fThreshold;
    double fxmin;
    double fxmax;

    RawAmpWatcher(const char* name, const char *title);
    virtual ~RawAmpWatcher(){}

    virtual void  Exec(const char *option);

    virtual void ShowSpectra(int LibraryNumber=1, const char *DrawOpt=""); //*MENU*
    virtual void ShowHist(int HistNumber, int LibraryNumber = 1, bool same = false);  //*MENU*

    void FitAllRoot(const char *FitType = "Dino", const char *source = "60Co");

    void FitSelectedRoot(const char *FitType = "Dino", const char *source = "60Co", bool all=false);

    bool fUseOffset;
    void UseOffset(bool offset = false){fUseOffset = offset;}

    void SearchPeak(int NPeaks=2, double width=50, double threshold=0.1, double xmin=0., double xmax=30000.);
    bool fDefinePersoBackground;
    int fNumberIterations;
    int fDirection;
    int fFilterOrder;
    bool fSmoothing;
    int fsmoothingWindow;
    bool fCompton;

    double fFirstLibGain;

    void SetFirstLibGain(double gain){fFirstLibGain = gain;}

    void SetPersoBackground(int NumberIterations = 100, int Direction = TSpectrum::kBackDecreasingWindow,int FilterOrder = TSpectrum::kBackOrder8,bool Smoothing = true,int smoothingWindow = TSpectrum::kBackSmoothing3,bool Compton = false);


public:
    Int_t fLasteventType;
    Int_t fLasteventX;
    Int_t fLasteventY;

    void HandleMovement(Int_t eventType, Int_t eventX, Int_t eventY, TObject *select);


public:

    GWRecal *fGWRecal;
    ofstream fFileOutGWRecal;

    FitFunctions *fFitFunction;

    TH1 *fCurrentHist;

    bool fSourceDefine;

    int GetCurrentHistAndLib();

    void PrintHelp(); //*MENU={Hierarchy="GWRecal/PrintHelp"}*

    void FitAllRecal(); //*MENU={Hierarchy="GWRecal/1) Calibrate/A) FitAll"}*
    void FitSelectedRecal(bool PrintHeader=true); //*MENU={Hierarchy="GWRecal/1) Calibrate/B) FitSelected"}*
    void SpecificFit(int NbrOfSpectra=38, int FirstHist=0, int Step=1); //*MENU={Hierarchy="GWRecal/1) Calibrate/C) SpecificFit"}*

    void SetChannelOffset(int Off){fGWRecal->SetChannelOffset(Off);} //*MENU={Hierarchy="GWRecal/2) Configure/A) Channel/a) SetChannelOffset"}*
    void SetGlobalChannelLimits(int ChFrom=50, int ChTo=32768); //*MENU={Hierarchy="GWRecal/2) Configure/A) Channel/b) SetGlobalChannelLimits"}*
    void SetHistChannelLimits(int SpecNbr=0, int lib=1, int ChFrom=50, int ChTo=32768); //*MENU={Hierarchy="GWRecal/2) Configure/A) Channel/c) SetHistChannelLimits"}*
    void SetSelectedHistChannelLimits(int ChFrom=50, int ChTo=32768);



    void SetSource(const char * SourcesNames); //*MENU={Hierarchy="GWRecal/2) Configure/B) TheoreticalPeaks/a) SetSource"}*
    void AddPeak(double EPeak){fGWRecal->AddPeak(EPeak);} //*MENU={Hierarchy="GWRecal/2) Configure/B) TheoreticalPeaks/b) AddPeak"}*
    void RemoveClosestPeakTo(double EPeak){fGWRecal->RemoveClosestPeakTo(EPeak);} //*MENU={Hierarchy="GWRecal/2) Configure/B) TheoreticalPeaks/c) RemoveClosestPeakTo"}*
    void SetRefPeak(double EPeak){fGWRecal->SetRefPeak(EPeak);} //*MENU={Hierarchy="GWRecal/2) Configure/B) TheoreticalPeaks/d) SetRefPeak"}*


    void SetGlobalPeaksLimits(float DefFWHM=15, float DefAmpli=5); //*MENU={Hierarchy="GWRecal/2) Configure/C) PeakSearch/a) SetGlobalPeaksLimits"}*
    void SetHistPeaksLimits(int SpecNbr=0, int lib=1, float FWHM=15, float Ampli=5);//*MENU={Hierarchy="GWRecal/2) Configure/C) PeakSearch/b) SetHistPeaksLimits"}*
    void SetSelectedHistPeaksLimits(float FWHM=15, float Ampli=5);
    void UseFirstDerivativeSearch(){fGWRecal->UseFirstDerivativeSearch();} //*MENU={Hierarchy="GWRecal/2) Configure/C) PeakSearch/c) UseFirstDerivativeSearch"}*
    void UseSecondDerivativeSearch(){fGWRecal->UseSecondDerivativeSearch();} //*MENU={Hierarchy="GWRecal/2) Configure/C) PeakSearch/d) UseSecondDerivativeSearch"}*

    void UseFlatBackGround(){fGWRecal->UseFlatBackGround();} //*MENU={Hierarchy="GWRecal/2) Configure/D) BackGround/a) UseFlatBackGround"}*
    void UseAffineBackGround(){fGWRecal->UseAffineBackGround();} //*MENU={Hierarchy="GWRecal/2) Configure/D) BackGround/b) UseAffineBackGround"}*
    void NoBackGround(){fGWRecal->NoBackGround();} //*MENU={Hierarchy="GWRecal/2) Configure/D) BackGround/c) NoBackGround"}*

    void UnsetLeftTail(){fGWRecal->UnsetLeftTail();} //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/a) UnsetLeftTail"}*
    void UnsetRightTail(){fGWRecal->UnsetRightTail();} //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/b) UnsetRightTail"}*
    void UseLinearCalib(){fGWRecal->UseLinearCalib();} //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/c) UseLinearCalib"}*
    void UseAffineCalib(){fGWRecal->UseAffineCalib();} //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/d) UseAffineCalib"}*
    void UseParabolicCalib(){fGWRecal->UseParabolicCalib();} //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/e) UseParabolicCalib"}*
    void SetGain(float Gain); //*MENU={Hierarchy="GWRecal/2) Configure/E) FitOptions/f) SetGain"}*

    void PrintResForXTalk(int VerboseMod=0){fGWRecal->PrintResForXTalk(VerboseMod);} //*MENU={Hierarchy="GWRecal/3) VerboseOptions/A) PrintResForXTalk"}*
    void SetVerbosityLevel(int Verb=0){fGWRecal->SetVerbosityLevel(Verb);} //*MENU={Hierarchy="GWRecal/3) VerboseOptions/B) SetVerbosityLevel"}*
    void SetFullVerbosity(){fGWRecal->SetFullVerbosity();} //*MENU={Hierarchy="GWRecal/3) VerboseOptions/C) SetFullVerbosity"}*

    void GWRecalFit(){fGWRecal->SetFileName(fFileName); fGWRecal->StartCalib();}



public:
    TTree *fCalibResultsTree;
    int InitTree();

    TString PlotVarExp, PlotVarSel, PlotVarOpt;
    void GetPlotParams(const char *varexp = "FW05:spec", const char *selection="", const char* option="");
    void PlotCalibResultsAsTreeLaucher(); //*MENU={Hierarchy="GWRecal/4) PlotResults/A) PlotCalibResults"}*
    void PrintPlotHelp(); //*MENU={Hierarchy="GWRecal/4) PlotResults/B) PrintPlotHelp"}*
    void PlotCalibResultsAsTree(int LibNbr);

    virtual void WriteCalibFiles(const char* ConfDirectory = ""); //*MENU={Hierarchy="WriteCalibFiles"}*

    ClassDef(RawAmpWatcher,1)
};

#endif