void check1()
{
gStyle->SetOptStat(1002211);
gStyle->SetOptFit(1111);



TCanvas *c1 = new TCanvas("c1","Without constraint",800,600);
c1->Divide(4,1);

TH1D*h10 = new TH1D("h10","First Derivative (3 Bins)", 92,44,136);
TH1D*h11 = new TH1D("h11","First Derivative (5 Bins)",92,44,136);
TH1D*h12 = new TH1D("h12","First Derivative (7 Bins)", 92,44,136);
TH1D*h13 = new TH1D("h13","First Derivative (9 Bins)",92,44,136);

//Range

Double_t range1 = 44;
Double_t range2 = 136;

// Mass Limits
Double_t LowerLimit_Mass = 44;
Double_t UpperLimit_Mass = 136;


TF1 *f12_h2 = new TF1("f12_h2","[0]*exp(-0.5*((x-[1])/[2])**2)",range1,range2);
f12_h2->SetFillColor(4);

f12_h2->FixParameter(0,11000);
f12_h2->FixParameter(1, 88);
f12_h2->FixParameter(2, 10);

// First function
c1->cd(1);
gStyle->SetOptStat(1002211);
gStyle->SetOptFit(1111);
f12_h2->Draw("");
auto h2_I = (TH1D*)f12_h2->GetHistogram();
h2_I->Draw("HIST");

TF1 *f12_A1_i= new TF1("f12_A1_i","[0]*exp(-0.5*((x-(88-[1]))/[2])**2) + [0]*exp(-0.5*((x-(88+[1]))/[2])**2)",range1,range2);

f12_A1_i->FixParameter(0,2906);
f12_A1_i->FixParameter(1, 8);
f12_A1_i->FixParameter(2, 5);

// Second function
c1->cd(2);
f12_A1_i->Draw("");
auto h1_A1_i = (TH1D*)f12_A1_i->GetHistogram();
h1_A1_i->Draw("HIST");


TF1 *f12_A1_ii= new TF1("f12_A1_ii","[0]*exp(-0.5*((x-(88-[1]))/[2])**2) + [0]*exp(-0.5*((x-(88+[1]))/[2])**2)",range1,range2);

f12_A1_ii->FixParameter(0,200);
f12_A1_ii->FixParameter(1, 24);
f12_A1_ii->FixParameter(2, 5);

// Third function
c1->cd(3);
f12_A1_ii->Draw("");
auto h1_A1_ii = (TH1D*)f12_A1_ii->GetHistogram();
h1_A1_ii->Draw("HIST");

TH1D*h1 = new TH1D("h1","h1",92,44,136); 
h1->Add(h2_I);
h1->Add(h1_A1_i);
h1->Add(h1_A1_ii);

// Total function 
c1->cd(4);
h1->Draw("HIST");

}