30template<
class Iterator>
35 std::multimap<TString, int > parMap;
36 for (Iterator it = begin; it != end; ++it) {
37 parMap.insert( std::make_pair( *it, std::distance(begin,it) ) );
39 for (
auto & elem : parMap) {
41 int n = parMap.count(
name);
43 std::pair <std::multimap<TString,int>::iterator, std::multimap<TString,int>::iterator> ret;
44 ret = parMap.equal_range(
name);
46 for (std::multimap<TString,int>::iterator it=ret.first; it!=ret.second; ++it) {
68 fFunctions = std::vector<std::unique_ptr<TF1>>(functions.size());
72 TF1 *
f = (
TF1*) functions[
n]->IsA()->New();
73 functions[
n]->Copy(*
f);
78 Fatal(
"InitializeDataMembers",
"Invalid input function -- abort");
92 for (
int i=0; i<npar; i++)
100 for (
int i=0; i < npar; i++)
136 fFunctions = std::vector<std::unique_ptr<TF1>>(0);
137 fCoeffs = std::vector < Double_t >(0) ;
155 std::vector<TF1 *> functions(2);
156 std::vector < Double_t > coeffs(2);
158 functions = {function1, function2};
159 coeffs = {coeff1, coeff2};
169 std::vector<TF1 *> functions(3);
170 std::vector < Double_t > coeffs(3);
172 functions = {function1, function2, function3};
173 coeffs = {coeff1, coeff2, coeff3};
190 Int_t noffunctions = (formula.
Tokenize(
"+")) -> GetEntries();
191 Int_t nofobj = arrayall -> GetEntries();
192 Int_t nofcoeffs = nofobj - noffunctions;
194 std::vector<TF1 *> functions(noffunctions);
195 std::vector < Double_t > coeffs(noffunctions);
196 std::vector < TString > funcstringall(nofobj);
197 std::vector < Int_t > indexsizetimes(nofcoeffs+1);
198 std::vector < Bool_t > isacoeff(nofobj);
200 for (
int i=0; i<nofobj; i++)
202 funcstringall[i] = ((
TObjString*)((*arrayall)[i])) -> GetString();
203 funcstringall[i].ReplaceAll(
" ",
"");
209 for (
int i=0; i<nofcoeffs+1; i++)
211 indexsizetimes[i] = ( ( ( (
TObjString*)(*arraytimes)[i] ) -> GetString() ).Tokenize(
"+") ) -> GetEntries();
212 while (k < indexsizetimes[i])
217 j = j+indexsizetimes[i];
218 if (j==nofobj) isacoeff[j-1] = 0;
219 else isacoeff[j-1] = 1;
223 Double_t old_xmin = 0.0, old_xmax = 0.0;
225 for (
int i=0; i<noffunctions; i++)
229 coeffs[i] = funcstringall[k].Atof();
236 functions[i] = (
TF1 *)(
gROOT->GetListOfFunctions()->FindObject(funcstringall[k]));
238 Error(
"TF1NormSum",
"Function %s does not exist", funcstringall[k].Data());
241 functions[i]->GetRange(old_xmin, old_xmax);
250 if (noffunctions > 0 && functions[0])
251 functions[0]->SetRange(old_xmin, old_xmax);
293 std::vector<double> params(
GetNpar() );
295 int nOfNonCstParams = 0;
299 offset += nOfNonCstParams;
301 for (
int j = 0; j <
fFunctions[
n]->GetNpar(); ++j) {
326 bool equalParams =
true;
331 for (
int i = 0; i < npar; ++i) {
336 equalParams &= (funcParams[i] == params[k+
fNOfFunctions+offset] );
356 const double params[] = {p0,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10};
380 Warning(
"SetRange",
"Invalid range: %f >= %f",
a,
b);
427 ((
TF1NormSum &)obj).fFunctions[
n] = std::unique_ptr<TF1>(
f);
void FixDuplicateNames(Iterator begin, Iterator end)
Function to find and rename duplicate parameters with the same name.
Class adding two functions: c1*f1+c2*f2.
Double_t fScale
Fixed Scale parameter to normalize function (e.g. bin width)
const char * GetParName(Int_t ipar) const
void GetRange(Double_t &a, Double_t &b) const
TF1NormSum & operator=(const TF1NormSum &rhs)
Operator =.
Double_t fXmax
Minimal bound of range of NormSum.
std::vector< TString > fParNames
Parameter names.
void SetParameters(const Double_t *params)
Initialize array of all parameters.
std::vector< double > GetParameters() const
Return array of parameters.
std::vector< std::unique_ptr< TF1 > > fFunctions
Maximal bound of range of NormSum.
void InitializeDataMembers(const std::vector< TF1 * > &functions, const std::vector< Double_t > &coeffs, Double_t scale)
void Copy(TObject &obj) const
Copy this to obj.
std::vector< Double_t > fCoeffs
Vector of size afNOfFunctions containing coefficients in front of each function.
Int_t GetNpar() const
Return the number of (non constant) parameters including the coefficients: for 2 functions: c1,...
void Update()
Update the component functions of the normalized sum.
unsigned int fNOfFunctions
Number of functions to add.
double operator()(const Double_t *x, const Double_t *p)
Overload the parenthesis to add the functions.
void SetRange(Double_t a, Double_t b)
std::vector< Int_t > fCstIndexes
Vector with size of fNOfFunctions containing the index of the constant parameter/ function (the remov...
static void InitStandardFunctions()
Create the basic function objects.
Collectable string class.
Mother of all ROOT objects.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
virtual void Fatal(const char *method, const char *msgfmt,...) const
Issue fatal error message.
TObjArray * Tokenize(const TString &delim) const
This function is used to isolate sequential tokens in a TString.
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString.
static long int sum(long int i)