TF1NormSum.h

Go to the documentation of this file.

// @(#)root/hist:$Id$
// Authors: Lorenzo Moneta, Aurélie Flandi  27/08/14
 
/**********************************************************************
 *                                                                    *
 * Copyright (c) 2015  ROOT  Team, CERN/PH-SFT                        *
 *                                                                    *
 *                                                                    *
 **********************************************************************/
 
#ifndef ROOT_TF1NormSum__
#define ROOT_TF1NormSum__
 
#include <iostream>
#include "TF1.h"
#include <memory>
#include "TF1AbsComposition.h"
 
class TF1NormSum : public TF1AbsComposition {
 
protected:
   unsigned int fNOfFunctions;                         ///< Number of functions to add
   Double_t     fScale;                                ///< Fixed Scale parameter to normalize function (e.g. bin width)
   Double_t fXmin;                                     /// Minimal bound of range of NormSum
   Double_t fXmax;                                     /// Maximal bound of range of NormSum
   std::vector<std::unique_ptr<TF1>> fFunctions;       ///< Vector of size fNOfFunctions containing TF1 functions
   std::vector < Double_t  > fCoeffs;                  ///< Vector of size afNOfFunctions containing coefficients in front of each function
   std::vector < Int_t     > fCstIndexes;              ///< Vector with size of fNOfFunctions containing the index of the constant parameter/ function (the removed ones)
   std::vector< TString >    fParNames;                ///< Parameter names
 
   void InitializeDataMembers(const std::vector<TF1 *> &functions, const std::vector<Double_t> &coeffs,
                              Double_t scale); // acts as a constructor
 
public:
 
   TF1NormSum();
   TF1NormSum(const std::vector <TF1*>&functions, const std::vector <Double_t> &coeffs, Double_t scale = 1.);
   TF1NormSum(TF1* function1, TF1* function2, Double_t coeff1 = 1., Double_t coeff2 = 1., Double_t scale = 1.);
   TF1NormSum(TF1* function1, TF1* function2, TF1*function3, Double_t coeff1 = 1., Double_t coeff2 = 1., Double_t coeff3 = 1., Double_t scale = 1.);
   TF1NormSum(const TString &formula, Double_t xmin, Double_t xmax);
 
   // Copy constructor
   TF1NormSum(const TF1NormSum &nsum);
 
   TF1NormSum &operator=(const TF1NormSum &rhs);
 
   virtual ~TF1NormSum() {}
 
   double operator()(const Double_t *x, const Double_t *p);
 
   std::vector<double> GetParameters() const;
 
   void        SetScale(Double_t scale) { fScale = scale; };
 
   void SetParameters(const Double_t *params);
 
   void        SetParameters(Double_t p0, Double_t p1, Double_t p2=0., Double_t p3=0., Double_t p4=0.,
                                   Double_t p5=0., Double_t p6=0., Double_t p7=0., Double_t p8=0., Double_t p9=0., Double_t p10=0.);
 
   void SetRange(Double_t a, Double_t b);
 
   Int_t       GetNpar() const;
 
   Double_t    GetScale() const { return fScale; }
 
   const char* GetParName(Int_t ipar) const { return fParNames.at(ipar).Data(); }
 
   Double_t GetXmin() const { return fXmin; }
 
   Double_t GetXmax() const { return fXmax; }
 
   void GetRange(Double_t &a, Double_t &b) const;
 
   void Update();
 
   void Copy(TObject &obj) const;
 
   ClassDef(TF1NormSum, 1);
};
#endif /* defined(ROOT_TF1NormSum__) */