{

Double_t pHeight = 1.0;
Double_t bMargin = 0.25;
Double_t tMargin = 0.01;
Double_t Height;
Height = 5*pHeight + tMargin + bMargin;
Double_t pWidth = 3.0;
Double_t rMargin = 0.01;
Double_t lMargin = 0.055;
Double_t yAxisWidth = 0.15;
Double_t Width;
Width = pWidth + rMargin + yAxisWidth;

TFile fresA("/twist/data15/analysis/andr_analysis_fits/standard/nominal_A/nominal_A_standard.root");
fresA->cd();
TH1D* hpresA;
hpresA = (TH1D*)full_residuals_ptot->Clone();
TFile fresB("/twist/data15/analysis/andr_analysis_fits/standard/nominal_B/nominal_B_standard.root");
fresB->cd();
TH1D* hpresB;
hpresB = (TH1D*)full_residuals_ptot->Clone();
TFile fresLF("/twist/data15/analysis/andr_analysis_fits/standard/low_field/low_field_standard.root");
fresLF->cd();
TH1D* hpresLF;
hpresLF = (TH1D*)full_residuals_ptot->Clone();
TFile fresHF("/twist/data15/analysis/andr_analysis_fits/standard/high_field/high_field_standard.root");
fresHF->cd();
TH1D* hpresHF;
hpresHF = (TH1D*)full_residuals_ptot->Clone();
TFile fresCl("/twist/data15/analysis/andr_analysis_fits/standard/cloud/cloud_standard.root");
fresCl->cd();
TH1D* hpresCl;
hpresCl = (TH1D*)full_residuals_ptot->Clone();

Double_t Lbl = 0.15;

hpresA->SetTitle("");
hpresA->GetXaxis()->SetLabelSize(Lbl);
hpresA->GetYaxis()->SetLabelSize(Lbl*pHeight/(pHeight+tMargin));
hpresA->GetYaxis()->SetNdivisions(205,kTRUE);
hpresA->SetMaximum(3.9);
hpresA->SetMinimum(-3.9);

hpresB->SetTitle("");
hpresB->GetXaxis()->SetLabelSize(Lbl);
hpresB->GetYaxis()->SetLabelSize(Lbl);
hpresB->GetYaxis()->SetNdivisions(205,kTRUE);
hpresB->SetMaximum(3.9);
hpresB->SetMinimum(-3.9);

hpresLF->SetTitle("");
hpresLF->GetXaxis()->SetLabelSize(Lbl);
hpresLF->GetYaxis()->SetLabelSize(Lbl);
hpresLF->GetYaxis()->SetNdivisions(205,kTRUE);
hpresLF->SetMaximum(3.9);
hpresLF->SetMinimum(-3.9);

hpresHF->SetTitle("");
hpresHF->GetXaxis()->SetLabelSize(Lbl);
hpresHF->GetYaxis()->SetLabelSize(Lbl);
hpresHF->GetYaxis()->SetNdivisions(205,kTRUE);
hpresHF->SetMaximum(3.9);
hpresHF->SetMinimum(-3.9);

hpresCl->SetTitle("");
hpresCl->GetXaxis()->SetTitle("Momentum (MeV/c)");
hpresCl->GetXaxis()->SetTitleSize(Lbl*pHeight/(pHeight+bMargin));
hpresCl->GetXaxis()->CenterTitle();
hpresCl->GetXaxis()->SetLabelSize(Lbl*pHeight/(pHeight+bMargin));
hpresCl->GetYaxis()->SetLabelSize(Lbl*pHeight/(pHeight+bMargin));
hpresCl->GetYaxis()->SetNdivisions(205,kTRUE);
hpresCl->SetMaximum(3.9);
hpresCl->SetMinimum(-3.9);

TCanvas *cpres;
cpres = new TCanvas("cpres","Fit residuals vs. momentum",150*Width,150*Height);
cpres->SetFillColor(0);
cpres->SetBorderSize(0.0);
cpres->SetBorderMode(0);
cpres->SetFrameBorderMode(0);
cpres->Divide(1,6);
//----------------------------------------------------------------
cpres->cd(1);
TPad *p1;
p1 = gPad;
p1->SetPad("p1","Nominal A",yAxisWidth/Width,(4*pHeight+bMargin)/Height,1.0,1.0);
p1->SetFillColor(0);
p1->SetBorderMode(0);
p1->SetBorderSize(0);
p1->SetFrameBorderMode(0);
p1->SetTopMargin(tMargin);
p1->SetRightMargin(rMargin);
p1->SetLeftMargin(lMargin);
p1->SetBottomMargin(0.0);
TPaveText *plA;
plA = new TPaveText(0,2,16,4);
plA->SetFillStyle(0);
plA->SetBorderSize(0);
plA->SetTextSize(Lbl);
TText *lA = plA->AddText("Nominal A");
//----------------------------------------------------------------
cpres->cd(2);
TPad *p2;
p2 = gPad;
p2->SetPad("p2","Nominal B",yAxisWidth/Width,(3*pHeight+bMargin)/Height,1.0,(4*pHeight+bMargin)/Height);
p2->SetFillColor(0);
p2->SetBorderMode(0);
p2->SetBorderSize(0);
p2->SetFrameBorderMode(0);
p2->SetTopMargin(0.0);
p2->SetRightMargin(rMargin);
p2->SetLeftMargin(lMargin);
p2->SetBottomMargin(0.0);
TPaveText *plB;
plB = new TPaveText(0,2,16,4);
plB->SetFillStyle(0);
plB->SetBorderSize(0);
plB->SetTextSize(Lbl);
TText *lB = plB->AddText("Nominal B");
//----------------------------------------------------------------
cpres->cd(3);
TPad *p3;
p3 = gPad;
p3->SetPad("p3","Low Field",yAxisWidth/Width,(2*pHeight+bMargin)/Height,1.0,(3*pHeight+bMargin)/Height);
p3->SetFillColor(0);
p3->SetBorderMode(0);
p3->SetBorderSize(0);
p3->SetFrameBorderMode(0);
p3->SetTopMargin(0.0);
p3->SetRightMargin(rMargin);
p3->SetLeftMargin(lMargin);
p3->SetBottomMargin(0.0);
TPaveText *plLF;
plLF = new TPaveText(0,2,16,4);
plLF->SetFillStyle(0);
plLF->SetBorderSize(0);
plLF->SetTextSize(Lbl);
TText *lLF = plLF->AddText("Low Field");
//----------------------------------------------------------------
cpres->cd(4);
TPad *p4;
p4 = gPad;
p4->SetPad("p4","High Field",yAxisWidth/Width,(pHeight+bMargin)/Height,1.0,(2*pHeight+bMargin)/Height);
p4->SetFillColor(0);
p4->SetBorderMode(0);
p4->SetBorderSize(0);
p4->SetFrameBorderMode(0);
p4->SetTopMargin(0.0);
p4->SetRightMargin(rMargin);
p4->SetLeftMargin(lMargin);
p4->SetBottomMargin(0.0);
TPaveText *plHF;
plHF = new TPaveText(0,2,16,4);
plHF->SetFillStyle(0);
plHF->SetBorderSize(0);
plHF->SetTextSize(Lbl);
TText *lHF = plHF->AddText("High Field");
//----------------------------------------------------------------
cpres->cd(5);
TPad *p5;
p5 = gPad;
p5->SetPad("p5","Cloud",yAxisWidth/Width,0.0,1.0,(pHeight+bMargin)/Height);
p5->SetFillColor(0);
p5->SetBorderMode(0);
p5->SetBorderSize(0);
p5->SetFrameBorderMode(0);
p5->SetTopMargin(0.0);
p5->SetRightMargin(rMargin);
p5->SetLeftMargin(lMargin);
p5->SetBottomMargin(bMargin);
TPaveText *plCl;
plCl = new TPaveText(0,2,16,4);
plCl->SetFillStyle(0);
plCl->SetBorderSize(0);
plCl->SetTextSize(Lbl*pHeight/(pHeight+bMargin));
TText *lCl = plCl->AddText("Cloud");

//----------------------------------------------------------------
cpres->cd(6);
TPad *p6;
p6 = gPad;
p6->SetPad("p6","y axis",0.0,0.0,yAxisWidth/Width,1.0);
p6->SetFillColor(0);
p6->SetBorderMode(0);
p6->SetBorderSize(0);
p6->SetFrameBorderMode(0);
p6->Range(0.0,0.0,1.0,1.0);
TPaveText *plyAxis;
plyAxis = new TPaveText(0.0,0.524,1.0,0.525);
plyAxis->SetFillStyle(0);
plyAxis->SetBorderSize(0);
plyAxis->SetTextSize(7*Lbl);
TText *lyAxis = plyAxis->AddText("(Data-Fit)/#sigma");
lyAxis->SetTextAngle(90);

cpres->cd(1);
hpresA->Draw();
plA->Draw();
cpres->cd(2);
hpresB->Draw();
plB->Draw();
cpres->cd(3);
hpresLF->Draw();
plLF->Draw();
cpres->cd(4);
hpresHF->Draw();
plHF->Draw();
cpres->cd(5);
hpresCl->Draw();
plCl->Draw();
cpres->cd(6);
plyAxis->Draw();


}