QCalPeak.hh

Go to the documentation of this file.

#ifndef __QCALPEAK_HH_
#define __QCALPEAK_HH_

#include <utility>
#include <iostream>
#include <string>
#include <vector>
#include "TF1.h"
#include "TTree.h"
#include "TH1D.h"

class QCalPeak {
        public:
      QCalPeak();
                
                QCalPeak(const QCalPeak& orig);
                
                QCalPeak(std::string Name, double actualEnergy);
                
                QCalPeak(std::string Name, double actualEnergy, TF1* fitFunction);
                
                virtual ~QCalPeak();
                
                bool Found() const {return fFound;};
                
                double GetEnergy() const {return fActualEnergy;};
                
                std::pair<double, double> GetPeakMean();
                
                std::pair<double, double> GetPeakWidth();
                
                std::string GetName() {return fName;};
                
                void SetName(std::string name) {fName = name;};         
                
                double GetMinimum() {return fFitFunction->GetMinimum();};
                
                double GetMaximum() {return fFitFunction->GetMaximum();};
                
                void SetFitFunction(TF1 *newFitFunction);
                
                void SetRange(double newMin, double newMax);
                
                void SetMinimum(double newMin);
                
                void SetMaximum(double newMax);
                
                bool SetInitialValues(std::vector<double> init);
                
//              TF1* Fit(std::vector<double> data, TH1D *&temphist);
                TF1* Fit(std::vector<double> data);
                
                bool AutoInitialize(double ApproximateMean, std::vector<double> data, double NearestNeighbor, double ApproximateWidth=-1);
                
                std::string Print(bool toScreen=true);


        private:
                double fActualEnergy;
                
                std::string fName;
                
                TF1 *fFitFunction;
                
                bool fFound;

};

Double_t DefaultFit(Double_t *x, Double_t *par);
double chisquare(TF1* function, TH1D* histogram);

#endif

---