#include <TH1.h>
#include <TF1.h>

#include <Riostream.h>
#include <stdlib.h>
#include <stdio.h>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <string>
#include "TChain.h"
#include "TClassTable.h"
#include "TH1F.h"
#include "TF1.h"
#include "TCanvas.h"
#include "TPostScript.h"
#include "TSystem.h"
#include "TROOT.h"
#include "TLorentzVector.h"

#define MAX_HEIGHT  9000000.
#define MAX_AREA    9000000.

#define DEFAULT_SIGMA     10
#define DEFAULT_RANGE_MAG 8

Double_t MyQuadroPeaks(Double_t* x, Double_t* par);
Double_t MyQuadroTails(Double_t* x, Double_t* par);
Double_t FitMyQuadroPeaks(TH1F* spec, Double_t* par, Double_t* par_err, 
			  Double_t range1, Double_t range2, Double_t range3, Double_t range4,
			  Double_t min_area1, Double_t min_area2, Double_t min_area3, Double_t min_area4,
			  Double_t fit_low, Double_t fit_high, Int_t noBG,
			  Int_t color, Int_t integrate, TCanvas* canvas);