void ProcEvent(Int_t event, Int_t px, Int_t py, TObject *sel)
{
// print event type and current cursor position
TCanvas *c = (TCanvas *) gTQSender;
TPad *pad = (TPad *) c->GetSelectedPad();
if(!pad) return;
gROOT->SetEditHistograms(kFALSE);
//printf("event=%d, px=%d, py=%d\n", event, px, py);
//Float_t x = pad->AbsPixeltoX(px);
//Float_t y = pad->AbsPixeltoY(py);
//x = pad->PadtoX(x);
//y = pad->PadtoY(y);
//printf("x=%.3g, y=%.3g\n",x,y);
if(event==kButton1Double){
pad->Pop();
//printf("%s %d\n",pad->GetName(),pad->GetNumber());
pad->cd();
TCanvas *c_blow;
TIter next(pad->GetListOfPrimitives());
if((TCanvas*)gROOT->GetListOfCanvases()->FindObject("c_blow")){
c_blow=(TCanvas*)gROOT->GetListOfCanvases()->FindObject("c_blow");
c_blow->Clear();
//printf("show %d\n",c_blow->GetUniqueID());
//c_blow->GetCanvasImp()->Show();
//c_blow->RaiseWindow();
//c_blow->GetCanvasImp()->RaiseWindow();
//c_blow->Flush();
//gVirtualX->RaiseWindow(c_blow->GetUniqueID());
//c_blow->Show();
}
else{
// c_blow = new TCanvas("c_blow","blowup",750,0,743,525);
c_blow = new TCanvas("c_blow","blowup",100,50,600,600);
c_blow->SetFillColor(10); //white
//c_blow->ToggleEventStatus();
//c_blow->SetCrosshair();
c_blow->Draw();
}
c_blow->cd();
gROOT->SetSelectedPad((TPad*)c_blow);
TPad *clone =(TPad*)pad->Clone();
clone->SetPad(0,0,1,1);
clone->Draw();
clone->Modified();
clone->Update();
c_blow->Show();
}
}
// ______________________________________________________________________________________
void plotVsDeltaEta_nice(const Char_t* name = "ratioVsDeltaEta") {
gROOT->LoadMacro("include/toolsEtaNice.C++");
SetupStyle();
SetGlobals();
// -----------------------------------------------------
TFile *inFiles[nEtaSets][nDataSets][nMoments];
TGraphErrors *graphs[nEtaSets][nDataSets][nMoments];
for (int idxEta = 0 ; idxEta < nEtaSets; ++idxEta)
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet)
for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
inFiles[idxEta][idxDataSet][idxMoment] = TFile::Open(Form("output/%s/%s/Moments_%s.root",
aEtaSets[idxEta], aDataSets[idxDataSet], aMoments[idxMoment]));
graphs[idxEta][idxDataSet][idxMoment] = static_cast<TGraphErrors*>((inFiles[idxEta][idxDataSet][idxMoment]->Get(aMoments[idxMoment]))->Clone());
if (inFiles[idxEta][idxDataSet][idxMoment])
(inFiles[idxEta][idxDataSet][idxMoment])->Close();
}
// -----------------------------------------------------
TLegend *legRat[nDataSets];
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
legRat[idxDataSet] = new TLegend(0.12, 0.12, 0.70, 0.495);
legRat[idxDataSet]->SetTextAlign(12);
legRat[idxDataSet]->SetTextSize(0.06);
legRat[idxDataSet]->SetFillColor(1182);
legRat[idxDataSet]->SetLineColor(0);
legRat[idxDataSet]->SetBorderSize(0);
}
// -----------------------------------------------------
TGraphErrors *etaGraphs[9][nDataSets][nMoments];
for (int idxCent = 0; idxCent < nCent; idxCent++) {
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
for (int idxMoment = 4 ; idxMoment < nMoments; ++idxMoment) {
float y[nEtaSets];
float ey[nEtaSets];
for (int idxEta = 0 ; idxEta < nEtaSets; ++idxEta) {
y[idxEta] = graphs[idxEta][idxDataSet][idxMoment]->GetY()[idxCent];
ey[idxEta] = graphs[idxEta][idxDataSet][idxMoment]->GetEY()[idxCent];
}
etaGraphs[idxCent][idxDataSet][idxMoment] = new TGraphErrors(nEtaSets, aEtaSetBinWidth, y, 0, ey);
PrepareGraph(etaGraphs[idxCent][idxDataSet][idxMoment]);
} // for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
} // for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
} // for (int idxCent = 0; idxCent < 9; idxCent++) {
// -----------------------------------------------------
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
TPad *pad = SetupCanvas(Form("canDeltaEta_Ratio_%s", aDataSets[idxDataSet]), aDataSetsTitle[idxDataSet], "#Delta#eta", 0.48);
for (int idxMoment = 4 ; idxMoment < nMoments; ++idxMoment) {
pad->cd(idxMoment-3);
for (int idxCent = 0; idxCent < 9; idxCent++) {
if (idxCent != 0 && idxCent != 1 && idxCent != 8 && idxCent != 4 )
continue;
TGraphErrors *g = etaGraphs[idxCent][idxDataSet][idxMoment];
if (idxCent == 0)
ShiftGraphX(g, -0.015);
else if (idxCent == 1)
ShiftGraphX(g, -0.005);
else if (idxCent == 4)
ShiftGraphX(g, 0.005);
else if (idxCent == 8)
ShiftGraphX(g, 0.015);
ConfigGraph(g, idxMoment, idxCent);
if (idxCent == 0) {
g->Draw("AP");
if (idxMoment == 5) {
// TLine *line0 = new TLine(aMinX, 0, aMaxX, 0);
// line0->SetLineColor(kGray+1);
// line0->SetLineStyle(2);
// line0->SetLineWidth(2);
// line0->Draw();
}
else if (idxMoment == 6) {
TLine *line1 = new TLine(aMinX, 1, aMaxX, 1);
line1->SetLineColor(kGray+1);
line1->SetLineStyle(2);
line1->SetLineWidth(2);
line1->Draw();
legRat[idxDataSet]->AddEntry(line1, "Poisson", "l");
}
}
g->Draw("PSAME");
if (idxMoment == 4)
legRat[idxDataSet]->AddEntry(etaGraphs[idxCent][idxDataSet][4], Form("%s", cent1[idxCent]), "pl");
} // for (int idxCent = 0; idxCent < 9; idxCent++) {
} // for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
pad->cd(2);
TLatex *texb_3 = new TLatex(0.05, 0.55, "Au+Au collisions #sqrt{#it{s}_{NN}} = 14.5 GeV");
texb_3->SetTextSize(0.07);
texb_3->Draw("same");
TLatex *texb_3a = new TLatex(0.05, 0.49, "Net-Charge, 0.2 < #it{p}_{T} (GeV/#it{c}) < 2.0, 0-5%");
texb_3a->SetTextSize(0.07);
texb_3a->Draw("same");
TLatex *texb_3b = new TLatex(0.05,0.44, "statistical errors only");
texb_3b->SetTextSize(0.06);
texb_3b->Draw("same");
pad->cd(1);
TLatex *texb_4 = new TLatex(0.7, 19, "STAR Preliminary");
texb_4->SetTextSize(0.07);
texb_4->Draw("same");
pad->cd(3);
legRat[idxDataSet]->Draw();
pad->Modified();
} // for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
// -----------------------------------------------------
SaveCanvas(name);
// -----------------------------------------------------
TFile *fOut = TFile::Open("STAR_QM2015_Preliminary.root", "UPDATE");
fOut->cd();
TList* list = new TList;
for (int idxMoment = 4; idxMoment < nMoments; ++idxMoment) {
for (int idxCent = 0; idxCent < nCent; ++idxCent) {
if (idxCent > 1)
continue;
if (idxCent == 0)
ShiftGraphX(etaGraphs[idxCent][0][idxMoment], 0.015);
else if (idxCent == 1)
ShiftGraphX(etaGraphs[idxCent][0][idxMoment], 0.005);
else if (idxCent == 4)
ShiftGraphX(etaGraphs[idxCent][0][idxMoment], -0.005);
else if (idxCent == 8)
ShiftGraphX(etaGraphs[idxCent][0][idxMoment], -0.015);
etaGraphs[idxCent][0][idxMoment]->SetName(Form("Net-Charge_%s_DeltaEta_14.5GeV_%s_stat", aMoments[idxMoment], cent[idxCent]));
list->Add(etaGraphs[idxCent][0][idxMoment]);
}
}
list->Write("Net-Charge_VsDeltaEta", TObject::kSingleKey);
fOut->Close();
}
//################################################################################################################################
TCanvas *drawRatioPlotWithPurity(TH1D *prediction, TH1D *sr, TString xTitle, TString filename, double purity, bool LowECaloRegion){
TCanvas *c = new TCanvas("c"+filename,"c",0,0,500,500);
c->cd();
float y = 0.3;
TPad *pad1 = new TPad("pad1", "Control Plots 1", 0.01, y, 0.99, 0.99);
TPad *padRatio = new TPad("rp1", "Ratio1", 0.01, 0.01, 0.99, y-0.01);
TH1D *ratio = 0;
ratio = (TH1D*) sr->Clone();
ratio->Divide(prediction);
ratio->GetYaxis()->SetTitle((TString) sr->GetName() + "/" + (TString) prediction->GetName());
ratio->SetTitle("");
ratio->SetLabelSize(0.1,"X");
ratio->SetLabelSize(0.1,"Y");
ratio->SetTitleOffset(0.5,"Y");
ratio->SetTitleSize(0.11,"Y");
ratio->SetLineColor(kBlack);
ratio->SetMarkerColor(kBlack);
padRatio->cd();
ratio->GetYaxis()->SetRangeUser(0,2);
ratio->Draw();
// Draw line at one!
float xmin = ratio->GetXaxis()->GetXmin();
float xmax = ratio->GetXaxis()->GetXmax();
TLine *line = new TLine(xmin,1,xmax,1);
line->SetLineWidth(2);
line->Draw("same");
padRatio->Modified();
TLegend *leg = new TLegend(0.5,0.8,0.85,0.9);
leg->AddEntry(sr,sr->GetName(),"l");
leg->AddEntry(prediction,prediction->GetName(),"pel");
pad1->cd();
//pad1->SetLogy();
sr->SetLineColor(kRed);
sr->SetMarkerColor(kRed);
prediction->SetLineColor(kBlack);
prediction->SetMarkerColor(kBlack);
sr->SetTitle("");
prediction->SetTitle("");
prediction->GetXaxis()->SetTitle(xTitle);
sr->GetXaxis()->SetTitle(xTitle);
prediction->SetTitleSize(0.07,"X");
prediction->GetXaxis()->SetTitleOffset(0.7);
sr->SetTitleSize(0.07,"X");
sr->GetXaxis()->SetTitleOffset(0.7);
prediction->Draw("e");
sr->Draw("e same");
double maximum =prediction->GetMaximum()*1.2;
double minimum = 0.00001;
if(sr->GetMinimum()!=0) minimum=sr->GetMinimum()*0.9;
if(prediction->GetMinimum()!=0) minimum=prediction->GetMinimum()*0.9;
if(sr->GetMaximum()>prediction->GetMaximum()){
maximum=sr->GetMaximum()*1.2;
}
//leg->Draw("same");
//***********************
TLatex* info = new TLatex();
info->SetTextFont(132);
info-> SetNDC();
info->SetTextSize(0.06);
TString AuxString = Form("Purity of fakes in CR = %4.1f ",purity*100);
//info->DrawLatex(0.4,0.7,AuxString);
TLatex* info1 = new TLatex();
info1->SetTextFont(132);
info1-> SetNDC();
info1->SetTextSize(0.06);
if(LowECaloRegion) AuxString = "E_{calo}<10GeV";
else AuxString = "E_{calo}>10GeV";
//info1->DrawLatex(0.2,0.82,AuxString);
//***********************
// Draw both pads to canvas
c->cd();
pad1->Draw();
padRatio->SetGridy();
c -> SetBottomMargin(0.55);
//c->Modified();
padRatio->Draw();
c->SaveAs(filename);
return c;
}
//################################################################################################################################
//################################################################################################################################
//################################################################################################################################
TCanvas *drawRatioPlot(TH1D *prediction, TH1D *sr, TString xTitle, TString filename){
TCanvas *c = new TCanvas("c"+filename,"c",0,0,500,500);
c->cd();
float y = 0.3;
TPad *pad1 = new TPad("pad1", "Control Plots 1", 0.01, y, 0.99, 0.99);
TPad *padRatio = new TPad("rp1", "Ratio1", 0.01, 0.01, 0.99, y-0.01);
TH1D *ratio = 0;
ratio = (TH1D*) sr->Clone();
ratio->Divide(prediction);
ratio->GetYaxis()->SetTitle((TString) sr->GetName() + "/" + (TString) prediction->GetName());
ratio->SetTitle("");
ratio->SetLabelSize(0.1,"X");
ratio->SetLabelSize(0.1,"Y");
ratio->SetTitleOffset(0.5,"Y");
ratio->SetTitleSize(0.15,"Y");
padRatio->cd();
ratio->GetYaxis()->SetRangeUser(0,2);
ratio->Draw();
// Draw line at one!
float xmin = ratio->GetXaxis()->GetXmin();
float xmax = ratio->GetXaxis()->GetXmax();
TLine *line = new TLine(xmin,1,xmax,1);
line->SetLineWidth(2);
line->Draw("same");
padRatio->Modified();
TLegend *leg = new TLegend(0.6,0.8,0.9,0.9);
leg->AddEntry(sr,"bkg","l");
leg->AddEntry(prediction,prediction->GetName(),"pel");
pad1->cd();
//pad1->SetLogy();
sr->SetLineColor(kRed);
sr->SetMarkerColor(kRed);
sr->SetTitle("");
prediction->SetTitle("");
prediction->GetXaxis()->SetTitle(xTitle);
sr->GetXaxis()->SetTitle(xTitle);
prediction->SetTitleSize(0.07,"X");
prediction->GetXaxis()->SetTitleOffset(0.7);
sr->SetTitleSize(0.07,"X");
sr->GetXaxis()->SetTitleOffset(0.7);
prediction->Draw("e");
sr->Draw("e same");
if(sr->GetMaximum()>prediction->GetMaximum()){
sr->Draw("e");
prediction->Draw("e same");
}
leg->Draw("same");
// Draw both pads to canvas
c->cd();
pad1->Draw();
padRatio->SetGridy();
padRatio->Draw();
c->SaveAs(filename);
return c;
}
void limit()
{
//=========Macro generated from canvas: limit/limit
//========= (Thu Apr 27 14:38:33 2017) by ROOT version6.02/05
TCanvas *limit = new TCanvas("limit", "limit",0,0,600,600);
gStyle->SetOptFit(1);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
limit->SetHighLightColor(2);
limit->Range(0,0,1,1);
limit->SetFillColor(0);
limit->SetBorderMode(0);
limit->SetBorderSize(2);
limit->SetLeftMargin(0.16);
limit->SetRightMargin(0.04);
limit->SetTopMargin(0.06);
limit->SetBottomMargin(0.12);
limit->SetFrameFillStyle(0);
limit->SetFrameBorderMode(0);
// ------------>Primitives in pad: pad
TPad *pad = new TPad("pad", "pad",0,0,1,1);
pad->Draw();
pad->cd();
pad->Range(-59.99999,-2.046818,1565,15.01);
pad->SetFillColor(0);
pad->SetBorderMode(0);
pad->SetLogy();
pad->SetBorderSize(2);
pad->SetTickx(1);
pad->SetTicky(1);
pad->SetLeftMargin(0.16);
pad->SetRightMargin(0.04);
pad->SetTopMargin(0.06);
pad->SetBottomMargin(0.12);
pad->SetFrameFillStyle(0);
pad->SetFrameBorderMode(0);
pad->SetFrameFillStyle(0);
pad->SetFrameBorderMode(0);
TH1F *tmp01 = new TH1F("tmp01","Graph",100,200,1500);
tmp01->SetMinimum(0.);
tmp01->SetMaximum(13.98659);
tmp01->SetStats(0);
tmp01->SetLineStyle(0);
tmp01->SetMarkerStyle(20);
tmp01->GetXaxis()->SetTitle("m_{H} (GeV)");
tmp01->GetXaxis()->SetLabelFont(42);
tmp01->GetXaxis()->SetLabelOffset(0.01);
tmp01->GetXaxis()->SetTitleSize(0.05);
tmp01->GetXaxis()->SetTickLength(0.02);
tmp01->GetXaxis()->SetTitleOffset(1.08);
tmp01->GetXaxis()->SetTitleFont(42);
tmp01->GetYaxis()->SetTitle("95% CL limit on #sigma/#sigma_{SM}");
tmp01->GetYaxis()->SetLabelFont(42);
tmp01->GetYaxis()->SetLabelOffset(0.007);
tmp01->GetYaxis()->SetTitleSize(0.05);
tmp01->GetYaxis()->SetTickLength(0.02);
tmp01->GetYaxis()->SetTitleOffset(1.56);
tmp01->GetYaxis()->SetTitleFont(42);
tmp01->GetZaxis()->SetLabelFont(42);
tmp01->GetZaxis()->SetLabelOffset(0.007);
tmp01->GetZaxis()->SetTitleSize(0.05);
tmp01->GetZaxis()->SetTickLength(0.02);
tmp01->GetZaxis()->SetTitleFont(42);
tmp01->Draw("AXIS");
TH1F *tmp02 = new TH1F("tmp02","Graph",100,200,1500);
tmp02->SetMinimum(0);
tmp02->SetMaximum(13.98659);
tmp02->SetStats(0);
tmp02->SetLineStyle(0);
tmp02->SetMarkerStyle(20);
tmp02->GetXaxis()->SetTitle("m_{H} (GeV)");
tmp02->GetXaxis()->SetLabelFont(42);
tmp02->GetXaxis()->SetLabelOffset(0.01);
tmp02->GetXaxis()->SetTitleSize(0.05);
tmp02->GetXaxis()->SetTickLength(0.02);
tmp02->GetXaxis()->SetTitleOffset(1.08);
tmp02->GetXaxis()->SetTitleFont(42);
tmp02->GetYaxis()->SetTitle("95% CL limit on #sigma/#sigma_{SM}");
tmp02->GetYaxis()->SetLabelFont(42);
tmp02->GetYaxis()->SetLabelOffset(0.007);
tmp02->GetYaxis()->SetTitleSize(0.05);
tmp02->GetYaxis()->SetTickLength(0.02);
tmp02->GetYaxis()->SetTitleOffset(1.56);
tmp02->GetYaxis()->SetTitleFont(42);
tmp02->GetZaxis()->SetLabelFont(42);
tmp02->GetZaxis()->SetLabelOffset(0.007);
tmp02->GetZaxis()->SetTitleSize(0.05);
tmp02->GetZaxis()->SetTickLength(0.02);
tmp02->GetZaxis()->SetTitleFont(42);
tmp02->Draw("AXIGSAME");
Double_t Graph0_fx3001[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph0_fy3001[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
Double_t Graph0_felx3001[27] = {
4.032346e-316,
0,
0,
0,
6.906313e-310,
4.032343e-316,
2.121996e-314,
4.032582e-316,
0,
4.032588e-316,
4.032343e-316,
0,
4.032584e-316,
1.508018e-316,
1.508219e-316,
4.032589e-316,
2.124327e-314,
4.032593e-316,
5.582942e-322,
6.906313e-310,
0,
4.032344e-316,
4.032583e-316,
4.032588e-316,
1.58101e-322,
4.032293e-316,
5.582942e-322};
Double_t Graph0_fely3001[27] = {
2.370331,
1.057434,
0.405159,
0.1845837,
0.09618759,
0.05231667,
0.03192902,
0.02331161,
0.01754379,
0.01435089,
0.01183319,
0.009315491,
0.008041382,
0.006484985,
0.005966187,
0.004928589,
0.00440979,
0.004119873,
0.004119873,
0.003570557,
0.003570557,
0.00302124,
0.00302124,
0.002471924,
0.002471924,
0.002471924,
0.002471924};
Double_t Graph0_fehx3001[27] = {
4.031643e-316,
0,
0,
0,
6.906313e-310,
4.032106e-316,
2.121996e-314,
4.032349e-316,
0,
4.032354e-316,
4.032106e-316,
0,
4.032351e-316,
1.508018e-316,
1.508219e-316,
4.032356e-316,
2.124327e-314,
4.03236e-316,
5.582942e-322,
6.906313e-310,
0,
4.032107e-316,
4.03235e-316,
4.032354e-316,
1.58101e-322,
4.032222e-316,
5.582942e-322};
Double_t Graph0_fehy3001[27] = {
5.693069,
2.47175,
0.8989364,
0.4081147,
0.2108151,
0.1170805,
0.07362081,
0.05438244,
0.04146756,
0.03243252,
0.02825768,
0.02361592,
0.01990321,
0.01752287,
0.01712826,
0.0153559,
0.01456667,
0.01341494,
0.01192601,
0.0124739,
0.01149977,
0.01055484,
0.01051202,
0.008670811,
0.008670811,
0.008635777,
0.008635777};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(27,Graph0_fx3001,Graph0_fy3001,Graph0_felx3001,Graph0_fehx3001,Graph0_fely3001,Graph0_fehy3001);
grae->SetName("Graph0");
grae->SetTitle("Graph");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
grae->SetFillColor(ci);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3001 = new TH1F("Graph_Graph3001","Graph",100,70,1630);
Graph_Graph3001->SetMinimum(0);
Graph_Graph3001->SetMaximum(11.53875);
Graph_Graph3001->SetDirectory(0);
Graph_Graph3001->SetStats(0);
Graph_Graph3001->SetLineStyle(0);
Graph_Graph3001->SetMarkerStyle(20);
Graph_Graph3001->GetXaxis()->SetLabelFont(42);
Graph_Graph3001->GetXaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetXaxis()->SetTickLength(0.02);
Graph_Graph3001->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph3001->GetXaxis()->SetTitleFont(42);
Graph_Graph3001->GetYaxis()->SetLabelFont(42);
Graph_Graph3001->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph3001->GetYaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetYaxis()->SetTickLength(0.02);
Graph_Graph3001->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph3001->GetYaxis()->SetTitleFont(42);
Graph_Graph3001->GetZaxis()->SetLabelFont(42);
Graph_Graph3001->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph3001->GetZaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetZaxis()->SetTickLength(0.02);
Graph_Graph3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_Graph3001);
grae->Draw("3");
Double_t Graph1_fx3002[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph1_fy3002[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
Double_t Graph1_felx3002[27] = {
4.028644e-316,
2.608667e-321,
4.022558e-316,
0,
2.420922e-322,
4.028601e-316,
0,
2.11372e-314,
0,
4.940656e-323,
8.744962e-322,
0,
0,
4.940656e-323,
0,
6.906314e-310,
5.33769e+160,
0,
4.021573e-316,
6.906314e-310,
0,
1.188318e-312,
4.027998e-316,
0,
1.58101e-322,
4.027629e-316,
1.630417e-322};
Double_t Graph1_fely3002[27] = {
1.490716,
0.6650269,
0.2516418,
0.1146438,
0.05974151,
0.03228919,
0.01945674,
0.01420552,
0.01069075,
0.009417772,
0.007765532,
0.006113291,
0.005277157,
0.003647804,
0.00335598,
0.002772331,
0.002480507,
0.002317429,
0.002832413,
0.002231598,
0.002231598,
0.001888275,
0.001888275,
0.001544952,
0.001544952,
0.001544952,
0.002162933};
Double_t Graph1_fehx3002[27] = {
1.179415e-312,
2.608667e-321,
4.022558e-316,
0,
2.420922e-322,
4.028634e-316,
0,
2.113722e-314,
0,
4.940656e-323,
8.744962e-322,
0,
0,
4.940656e-323,
0,
6.906314e-310,
2.139459e+161,
0,
4.021573e-316,
6.906314e-310,
0,
1.188318e-312,
4.027998e-316,
0,
1.58101e-322,
4.027675e-316,
1.630417e-322};
Double_t Graph1_fehy3002[27] = {
2.351828,
1.053343,
0.3849801,
0.1768222,
0.09017253,
0.05078118,
0.03126747,
0.02303264,
0.01790216,
0.012869,
0.0113431,
0.009505779,
0.008205643,
0.007395977,
0.0068043,
0.006212621,
0.006088058,
0.005138258,
0.004204029,
0.005262822,
0.003643492,
0.004453157,
0.003768056,
0.004204029,
0.004204029,
0.003643492,
0.003643492};
grae = new TGraphAsymmErrors(27,Graph1_fx3002,Graph1_fy3002,Graph1_felx3002,Graph1_fehx3002,Graph1_fely3002,Graph1_fehy3002);
grae->SetName("Graph1");
grae->SetTitle("Graph");
ci = TColor::GetColor("#00ff00");
grae->SetFillColor(ci);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3002 = new TH1F("Graph_Graph3002","Graph",100,-8.010917e+160,2.406781e+161);
Graph_Graph3002->SetMinimum(0);
Graph_Graph3002->SetMaximum(7.86335);
Graph_Graph3002->SetDirectory(0);
Graph_Graph3002->SetStats(0);
Graph_Graph3002->SetLineStyle(0);
Graph_Graph3002->SetMarkerStyle(20);
Graph_Graph3002->GetXaxis()->SetLabelFont(42);
Graph_Graph3002->GetXaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetXaxis()->SetTickLength(0.02);
Graph_Graph3002->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph3002->GetXaxis()->SetTitleFont(42);
Graph_Graph3002->GetYaxis()->SetLabelFont(42);
Graph_Graph3002->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph3002->GetYaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetYaxis()->SetTickLength(0.02);
Graph_Graph3002->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph3002->GetYaxis()->SetTitleFont(42);
Graph_Graph3002->GetZaxis()->SetLabelFont(42);
Graph_Graph3002->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph3002->GetZaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetZaxis()->SetTickLength(0.02);
Graph_Graph3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_Graph3002);
grae->Draw("3");
Double_t Graph2_fx1[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph2_fy1[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
TGraph *graph = new TGraph(27,Graph2_fx1,Graph2_fy1);
graph->SetName("Graph2");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff0000");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
graph->SetMarkerStyle(20);
TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",100,70,1630);
Graph_Graph1->SetMinimum(0);
Graph_Graph1->SetMaximum(5.276123);
Graph_Graph1->SetDirectory(0);
Graph_Graph1->SetStats(0);
Graph_Graph1->SetLineStyle(0);
Graph_Graph1->SetMarkerStyle(20);
Graph_Graph1->GetXaxis()->SetLabelFont(42);
Graph_Graph1->GetXaxis()->SetTitleSize(0.05);
Graph_Graph1->GetXaxis()->SetTickLength(0.02);
Graph_Graph1->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph1->GetXaxis()->SetTitleFont(42);
Graph_Graph1->GetYaxis()->SetLabelFont(42);
Graph_Graph1->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph1->GetYaxis()->SetTitleSize(0.05);
Graph_Graph1->GetYaxis()->SetTickLength(0.02);
Graph_Graph1->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph1->GetYaxis()->SetTitleFont(42);
Graph_Graph1->GetZaxis()->SetLabelFont(42);
Graph_Graph1->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph1->GetZaxis()->SetTitleSize(0.05);
Graph_Graph1->GetZaxis()->SetTickLength(0.02);
Graph_Graph1->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph1);
graph->Draw("l");
Double_t Graph3_fx2[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph3_fy2[27] = {
4.79203,
2.152691,
0.8377939,
0.3861777,
0.1995838,
0.1090425,
0.0660598,
0.04759021,
0.03596779,
0.02751563,
0.02238263,
0.01815241,
0.01512282,
0.01228406,
0.01162986,
0.009836298,
0.008643669,
0.007866618,
0.007089044,
0.006702567,
0.005960677,
0.005744891,
0.005145786,
0.004956617,
0.004917011,
0.004476165,
0.004343825};
graph = new TGraph(27,Graph3_fx2,Graph3_fy2);
graph->SetName("Graph3");
graph->SetTitle("Graph");
graph->SetFillColor(1);
graph->SetLineWidth(2);
graph->SetMarkerStyle(20);
TH1F *Graph_Graph2 = new TH1F("Graph_Graph2","Graph",100,70,1630);
Graph_Graph2->SetMinimum(0);
Graph_Graph2->SetMaximum(5.270799);
Graph_Graph2->SetDirectory(0);
Graph_Graph2->SetStats(0);
Graph_Graph2->SetLineStyle(0);
Graph_Graph2->SetMarkerStyle(20);
Graph_Graph2->GetXaxis()->SetLabelFont(42);
Graph_Graph2->GetXaxis()->SetTitleSize(0.05);
Graph_Graph2->GetXaxis()->SetTickLength(0.02);
Graph_Graph2->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph2->GetXaxis()->SetTitleFont(42);
Graph_Graph2->GetYaxis()->SetLabelFont(42);
Graph_Graph2->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph2->GetYaxis()->SetTitleSize(0.05);
Graph_Graph2->GetYaxis()->SetTickLength(0.02);
Graph_Graph2->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph2->GetYaxis()->SetTitleFont(42);
Graph_Graph2->GetZaxis()->SetLabelFont(42);
Graph_Graph2->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph2->GetZaxis()->SetTitleSize(0.05);
Graph_Graph2->GetZaxis()->SetTickLength(0.02);
Graph_Graph2->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph2);
graph->Draw("pl");
TH1F *tmp0_copy3 = new TH1F("tmp0_copy3","Graph",100,200,1500);
tmp0_copy3->SetMinimum(0);
tmp0_copy3->SetMaximum(11.53875);
tmp0_copy3->SetDirectory(0);
tmp0_copy3->SetStats(0);
tmp0_copy3->SetLineStyle(0);
tmp0_copy3->SetMarkerStyle(20);
tmp0_copy3->GetXaxis()->SetLabelFont(42);
tmp0_copy3->GetXaxis()->SetTitleSize(0.05);
tmp0_copy3->GetXaxis()->SetTickLength(0.02);
tmp0_copy3->GetXaxis()->SetTitleOffset(1.08);
tmp0_copy3->GetXaxis()->SetTitleFont(42);
tmp0_copy3->GetYaxis()->SetLabelFont(42);
tmp0_copy3->GetYaxis()->SetLabelOffset(0.007);
tmp0_copy3->GetYaxis()->SetTitleSize(0.05);
tmp0_copy3->GetYaxis()->SetTickLength(0.02);
tmp0_copy3->GetYaxis()->SetTitleOffset(1.56);
tmp0_copy3->GetYaxis()->SetTitleFont(42);
tmp0_copy3->GetZaxis()->SetLabelFont(42);
tmp0_copy3->GetZaxis()->SetLabelOffset(0.007);
tmp0_copy3->GetZaxis()->SetTitleSize(0.05);
tmp0_copy3->GetZaxis()->SetTickLength(0.02);
tmp0_copy3->GetZaxis()->SetTitleFont(42);
tmp0_copy3->Draw("sameaxis");
TH1F *tmp0_copy4 = new TH1F("tmp0_copy4","Graph",100,200,1500);
tmp0_copy4->SetMinimum(0);
tmp0_copy4->SetMaximum(11.53875);
tmp0_copy4->SetDirectory(0);
tmp0_copy4->SetStats(0);
tmp0_copy4->SetLineStyle(0);
tmp0_copy4->SetMarkerStyle(20);
tmp0_copy4->GetXaxis()->SetLabelFont(42);
tmp0_copy4->GetXaxis()->SetTitleSize(0.05);
tmp0_copy4->GetXaxis()->SetTickLength(0.02);
tmp0_copy4->GetXaxis()->SetTitleOffset(1.08);
tmp0_copy4->GetXaxis()->SetTitleFont(42);
tmp0_copy4->GetYaxis()->SetLabelFont(42);
tmp0_copy4->GetYaxis()->SetLabelOffset(0.007);
tmp0_copy4->GetYaxis()->SetTitleSize(0.05);
tmp0_copy4->GetYaxis()->SetTickLength(0.02);
tmp0_copy4->GetYaxis()->SetTitleOffset(1.56);
tmp0_copy4->GetYaxis()->SetTitleFont(42);
tmp0_copy4->GetZaxis()->SetLabelFont(42);
tmp0_copy4->GetZaxis()->SetLabelOffset(0.007);
tmp0_copy4->GetZaxis()->SetTitleSize(0.05);
tmp0_copy4->GetZaxis()->SetTickLength(0.02);
tmp0_copy4->GetZaxis()->SetTitleFont(42);
tmp0_copy4->Draw("sameaxig");
TLegend *leg = new TLegend(0.495,0.825,0.945,0.925,NULL,"NBNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("Graph3","Observed","LP");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph2","Expected","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph1","#pm1#sigma Expected","F");
ci = TColor::GetColor("#00ff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph0","#pm2#sigma Expected","F");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TPave *pave = new TPave(0.16,0.81,0.96,0.94,1,"brNDC");
pave->SetFillColor(0);
pave->Draw();
leg = new TLegend(0.495,0.825,0.945,0.925,NULL,"NBNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
entry=leg->AddEntry("Graph3","Observed","LP");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph2","Expected","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph1","#pm1#sigma Expected","F");
ci = TColor::GetColor("#00ff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph0","#pm2#sigma Expected","F");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.196,0.9113,"CMS");
tex->SetNDC();
tex->SetTextAlign(13);
tex->SetTextSize(0.048);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.196,0.8537,"Internal");
tex->SetNDC();
tex->SetTextAlign(13);
tex->SetTextFont(52);
tex->SetTextSize(0.03648);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.96,0.952,"");
tex->SetNDC();
tex->SetTextAlign(31);
tex->SetTextFont(42);
tex->SetTextSize(0.036);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.16,0.952,"");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.036);
tex->SetLineWidth(2);
tex->Draw();
pad->Modified();
limit->cd();
limit->Modified();
limit->cd();
limit->SetSelected(limit);
}
void zj_ptj_8TeV()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Thu Apr 28 10:19:34 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,800,800);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: bottomPad
TPad *bottomPad = new TPad("bottomPad", "bottomPad",0.005,0.05,0.995,0.995);
bottomPad->Draw();
bottomPad->cd();
bottomPad->Range(0,0,1,1);
bottomPad->SetFillColor(0);
bottomPad->SetFillStyle(4000);
bottomPad->SetBorderMode(0);
bottomPad->SetBorderSize(2);
bottomPad->SetTickx(1);
bottomPad->SetTicky(1);
bottomPad->SetLeftMargin(0.14);
bottomPad->SetRightMargin(0.05);
bottomPad->SetTopMargin(0.05);
bottomPad->SetBottomMargin(0.16);
bottomPad->SetFrameLineWidth(2);
bottomPad->SetFrameBorderMode(0);
bottomPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(6.172839,-0.2132695,104.9383,0.4330016);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.33);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
Double_t xAxis19[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_toterr__19 = new TH1F("z_ptj_toterr__19","z_ptj",5, xAxis19);
z_ptj_toterr__19->SetBinContent(1,0.3608738);
z_ptj_toterr__19->SetBinContent(2,0.225932);
z_ptj_toterr__19->SetBinContent(3,0.1253807);
z_ptj_toterr__19->SetBinContent(4,0.06703913);
z_ptj_toterr__19->SetBinContent(5,0.01561868);
z_ptj_toterr__19->SetBinContent(6,0.780934);
z_ptj_toterr__19->SetBinError(1,0.02073388);
z_ptj_toterr__19->SetBinError(2,0.01400757);
z_ptj_toterr__19->SetBinError(3,0.009350811);
z_ptj_toterr__19->SetBinError(4,0.004623786);
z_ptj_toterr__19->SetBinError(5,0.001395185);
z_ptj_toterr__19->SetBinError(6,0.06975923);
z_ptj_toterr__19->SetEntries(12);
z_ptj_toterr__19->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetLineColor(ci);
z_ptj_toterr__19->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetMarkerColor(ci);
z_ptj_toterr__19->SetMarkerStyle(0);
z_ptj_toterr__19->GetXaxis()->SetNdivisions(4);
z_ptj_toterr__19->GetXaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetXaxis()->SetLabelOffset(999);
z_ptj_toterr__19->GetXaxis()->SetLabelSize(0);
z_ptj_toterr__19->GetXaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetXaxis()->SetTitleFont(132);
z_ptj_toterr__19->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_toterr__19->GetYaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetYaxis()->SetLabelSize(0.05);
z_ptj_toterr__19->GetYaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetYaxis()->SetTitleOffset(1.3);
z_ptj_toterr__19->GetYaxis()->SetTitleFont(132);
z_ptj_toterr__19->GetZaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetZaxis()->SetLabelSize(0.05);
z_ptj_toterr__19->GetZaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetZaxis()->SetTitleFont(132);
z_ptj_toterr__19->Draw("e2");
Double_t xAxis20[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_staterr__20 = new TH1F("z_ptj_staterr__20","z_ptj",5, xAxis20);
z_ptj_staterr__20->SetBinContent(1,0.3608738);
z_ptj_staterr__20->SetBinContent(2,0.225932);
z_ptj_staterr__20->SetBinContent(3,0.1253807);
z_ptj_staterr__20->SetBinContent(4,0.06703913);
z_ptj_staterr__20->SetBinContent(5,0.01561868);
z_ptj_staterr__20->SetBinContent(6,0.780934);
z_ptj_staterr__20->SetBinError(1,0.007287406);
z_ptj_staterr__20->SetBinError(2,0.005663236);
z_ptj_staterr__20->SetBinError(3,0.002951952);
z_ptj_staterr__20->SetBinError(4,0.002192817);
z_ptj_staterr__20->SetBinError(5,0.0004728113);
z_ptj_staterr__20->SetBinError(6,0.02364057);
z_ptj_staterr__20->SetEntries(12);
z_ptj_staterr__20->SetStats(0);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetLineColor(ci);
z_ptj_staterr__20->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetMarkerColor(ci);
z_ptj_staterr__20->SetMarkerStyle(0);
z_ptj_staterr__20->GetXaxis()->SetNdivisions(4);
z_ptj_staterr__20->GetXaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetXaxis()->SetLabelOffset(999);
z_ptj_staterr__20->GetXaxis()->SetLabelSize(0);
z_ptj_staterr__20->GetXaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetXaxis()->SetTitleFont(132);
z_ptj_staterr__20->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_staterr__20->GetYaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetYaxis()->SetLabelSize(0.05);
z_ptj_staterr__20->GetYaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetYaxis()->SetTitleFont(132);
z_ptj_staterr__20->GetZaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetZaxis()->SetLabelSize(0.05);
z_ptj_staterr__20->GetZaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetZaxis()->SetTitleFont(132);
z_ptj_staterr__20->Draw("e2same");
Double_t ptj_toterrs_fx3022[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t ptj_toterrs_fy3022[5] = {
0.3671394,
0.2404275,
0.1334798,
0.06537649,
0.01594908};
Double_t ptj_toterrs_felx3022[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fely3022[5] = {
0.02997783,
0.02293181,
0.009217184,
0.004503167,
0.001643067};
Double_t ptj_toterrs_fehx3022[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fehy3022[5] = {
0.02512764,
0.01751375,
0.008693834,
0.005222099,
0.001753864};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(5,ptj_toterrs_fx3022,ptj_toterrs_fy3022,ptj_toterrs_felx3022,ptj_toterrs_fehx3022,ptj_toterrs_fely3022,ptj_toterrs_fehy3022);
grae->SetName("ptj_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph_ptj_toterrs30153022 = new TH1F("Graph_Graph_ptj_toterrs30153022","",100,17.325,88.425);
Graph_Graph_ptj_toterrs30153022->SetMinimum(0.01287542);
Graph_Graph_ptj_toterrs30153022->SetMaximum(0.4300631);
Graph_Graph_ptj_toterrs30153022->SetDirectory(0);
Graph_Graph_ptj_toterrs30153022->SetStats(0);
Graph_Graph_ptj_toterrs30153022->SetLineWidth(2);
Graph_Graph_ptj_toterrs30153022->SetMarkerStyle(0);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetNdivisions(4);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelOffset(999);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelSize(0);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_ptj_toterrs30153022);
grae->Draw("p1p");
Double_t _fx3023[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t _fy3023[5] = {
0.3671394,
0.2404275,
0.1334798,
0.06537649,
0.01594908};
Double_t _felx3023[5] = {
0,
0,
0,
0,
0};
Double_t _fely3023[5] = {
0.01985385,
0.01348796,
0.006928003,
0.003412992,
0.001084033};
Double_t _fehx3023[5] = {
0,
0,
0,
0,
0};
Double_t _fehy3023[5] = {
0.01921597,
0.01338518,
0.006596172,
0.003531927,
0.001381032};
grae = new TGraphAsymmErrors(5,_fx3023,_fy3023,_felx3023,_fehx3023,_fely3023,_fehy3023);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph_Graph30163023 = new TH1F("Graph_Graph_Graph30163023","",100,17.325,88.425);
Graph_Graph_Graph30163023->SetMinimum(0.01337855);
Graph_Graph_Graph30163023->SetMaximum(0.4235044);
Graph_Graph_Graph30163023->SetDirectory(0);
Graph_Graph_Graph30163023->SetStats(0);
Graph_Graph_Graph30163023->SetLineWidth(2);
Graph_Graph_Graph30163023->SetMarkerStyle(0);
Graph_Graph_Graph30163023->GetXaxis()->SetNdivisions(4);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelOffset(999);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelSize(0);
Graph_Graph_Graph30163023->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph30163023->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_Graph30163023->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30163023->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph30163023->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30163023->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph30163023);
grae->Draw("p1p");
Double_t ptj_toterrs_fx3024[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t ptj_toterrs_fy3024[5] = {
0.389299,
0.2509263,
0.1443755,
0.07365244,
0.01874484};
Double_t ptj_toterrs_felx3024[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fely3024[5] = {
0.0283978,
0.01982285,
0.01219084,
0.00636368,
0.001684181};
Double_t ptj_toterrs_fehx3024[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fehy3024[5] = {
0.03041728,
0.02088176,
0.0127053,
0.007720362,
0.001934453};
grae = new TGraphAsymmErrors(5,ptj_toterrs_fx3024,ptj_toterrs_fy3024,ptj_toterrs_felx3024,ptj_toterrs_fehx3024,ptj_toterrs_fely3024,ptj_toterrs_fehy3024);
grae->SetName("ptj_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_ptj_toterrs30173024 = new TH1F("Graph_Graph_ptj_toterrs30173024","",100,17.175,72.075);
Graph_Graph_ptj_toterrs30173024->SetMinimum(0.01535459);
Graph_Graph_ptj_toterrs30173024->SetMaximum(0.4599818);
Graph_Graph_ptj_toterrs30173024->SetDirectory(0);
Graph_Graph_ptj_toterrs30173024->SetStats(0);
Graph_Graph_ptj_toterrs30173024->SetLineWidth(2);
Graph_Graph_ptj_toterrs30173024->SetMarkerStyle(0);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetNdivisions(4);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelOffset(999);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelSize(0);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_ptj_toterrs30173024);
grae->Draw("p1p");
Double_t _fx3025[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t _fy3025[5] = {
0.389299,
0.2509263,
0.1443755,
0.07365244,
0.01874484};
Double_t _felx3025[5] = {
0,
0,
0,
0,
0};
Double_t _fely3025[5] = {
0.02243218,
0.01539977,
0.009090877,
0.00492242,
0.001303763};
Double_t _fehx3025[5] = {
0,
0,
0,
0,
0};
Double_t _fehy3025[5] = {
0.02411537,
0.01689699,
0.0100922,
0.005544287,
0.001554499};
grae = new TGraphAsymmErrors(5,_fx3025,_fy3025,_felx3025,_fehx3025,_fely3025,_fehy3025);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph30183025 = new TH1F("Graph_Graph_Graph30183025","",100,17.175,72.075);
Graph_Graph_Graph30183025->SetMinimum(0.01569697);
Graph_Graph_Graph30183025->SetMaximum(0.4530117);
Graph_Graph_Graph30183025->SetDirectory(0);
Graph_Graph_Graph30183025->SetStats(0);
Graph_Graph_Graph30183025->SetLineWidth(2);
Graph_Graph_Graph30183025->SetMarkerStyle(0);
Graph_Graph_Graph30183025->GetXaxis()->SetNdivisions(4);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelOffset(999);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelSize(0);
Graph_Graph_Graph30183025->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph30183025->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_Graph30183025->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30183025->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph30183025->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30183025->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph30183025);
grae->Draw("p1p");
TLegend *leg = new TLegend(0.65,0.57,0.95,0.82,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("data_stat","Data_{stat}","f");
ci = TColor::GetColor("#ff9933");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("data_tot","Data_{tot}","f");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("powheg_temp","POWHEG","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("amcatnlo_temp","aMC@NLO","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(24);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.5,0.85,0.85,0.94,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb, #sqrt{s} = 8 TeV");
pt->Draw();
Double_t xAxis21[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_toterr_copy__21 = new TH1F("z_ptj_toterr_copy__21","z_ptj",5, xAxis21);
z_ptj_toterr_copy__21->SetBinContent(1,0.3608738);
z_ptj_toterr_copy__21->SetBinContent(2,0.225932);
z_ptj_toterr_copy__21->SetBinContent(3,0.1253807);
z_ptj_toterr_copy__21->SetBinContent(4,0.06703913);
z_ptj_toterr_copy__21->SetBinContent(5,0.01561868);
z_ptj_toterr_copy__21->SetBinContent(6,0.780934);
z_ptj_toterr_copy__21->SetBinError(1,0.02073388);
z_ptj_toterr_copy__21->SetBinError(2,0.01400757);
z_ptj_toterr_copy__21->SetBinError(3,0.009350811);
z_ptj_toterr_copy__21->SetBinError(4,0.004623786);
z_ptj_toterr_copy__21->SetBinError(5,0.001395185);
z_ptj_toterr_copy__21->SetBinError(6,0.06975923);
z_ptj_toterr_copy__21->SetEntries(12);
z_ptj_toterr_copy__21->SetDirectory(0);
z_ptj_toterr_copy__21->SetStats(0);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetLineColor(ci);
z_ptj_toterr_copy__21->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetMarkerColor(ci);
z_ptj_toterr_copy__21->SetMarkerStyle(0);
z_ptj_toterr_copy__21->GetXaxis()->SetNdivisions(4);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelOffset(999);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelSize(0);
z_ptj_toterr_copy__21->GetXaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetXaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_toterr_copy__21->GetYaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetYaxis()->SetLabelSize(0.05);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleOffset(1.3);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->GetZaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetZaxis()->SetLabelSize(0.05);
z_ptj_toterr_copy__21->GetZaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetZaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: lowerPad
TPad *lowerPad = new TPad("lowerPad", "lowerPad",0.005,0.05,0.995,0.995);
lowerPad->Draw();
lowerPad->cd();
lowerPad->Range(6.172839,-0.1935898,104.9383,4.832052);
lowerPad->SetFillColor(0);
lowerPad->SetFillStyle(4000);
lowerPad->SetBorderMode(0);
lowerPad->SetBorderSize(2);
lowerPad->SetTickx(1);
lowerPad->SetTicky(1);
lowerPad->SetLeftMargin(0.14);
lowerPad->SetRightMargin(0.05);
lowerPad->SetTopMargin(0.665);
lowerPad->SetBottomMargin(0.14);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
Double_t xAxis22[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp0__22 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp0__22","z_ptj",5, xAxis22);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(1,0.05745466);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(2,0.06199906);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(3,0.07457933);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(4,0.06897145);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(5,0.08932795);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(6,0.06975923);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleOffset(1.3);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->Draw("e2");
Double_t xAxis23[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp1__23 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp1__23","z_ptj",5, xAxis23);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(1,0.02019378);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(2,0.02506611);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(3,0.02354391);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(4,0.03270951);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(5,0.03027217);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(6,0.02364057);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetStats(0);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->Draw("e2same");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fx3026[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fy3026[5] = {
1.017362,
1.064159,
1.064596,
0.9751989,
1.021154};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_felx3026[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fely3026[5] = {
0.08307013,
0.1014988,
0.07351357,
0.06717222,
0.1051988};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fehx3026[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fehy3026[5] = {
0.06962999,
0.07751779,
0.06933948,
0.07789628,
0.1122927};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp2_fx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fy3026,zj_ptj_8TeV_400_1001_0ratiocomp2_felx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fehx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fely3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fehy3026);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp2");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026","",100,17.325,88.425);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026);
grae->Draw("p1");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fx3027[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fy3027[5] = {
1.017362,
1.064159,
1.064596,
0.9751989,
1.021154};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_felx3027[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fely3027[5] = {
0.05501606,
0.05969921,
0.05525572,
0.05091044,
0.06940621};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fehx3027[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fehy3027[5] = {
0.05324843,
0.05924428,
0.05260914,
0.05268456,
0.08842185};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp3_fx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fy3027,zj_ptj_8TeV_400_1001_0ratiocomp3_felx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fehx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fely3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fehy3027);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp3");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027","",100,17.325,88.425);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027);
grae->Draw("p1");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fx3028[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fy3028[5] = {
1.078768,
1.110628,
1.151497,
1.098648,
1.200155};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_felx3028[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fely3028[5] = {
0.07869176,
0.08773815,
0.0972306,
0.09492486,
0.1078312};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fehx3028[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fehy3028[5] = {
0.08428785,
0.09242501,
0.1013337,
0.115162,
0.1238551};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp4_fx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fy3028,zj_ptj_8TeV_400_1001_0ratiocomp4_felx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fehx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fely3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fehy3028);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp4");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028","",100,17.175,72.075);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028);
grae->Draw("p1");
Double_t xAxis24[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24","z_ptj",5, xAxis24);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(1,0.05745466);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(2,0.06199906);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(3,0.07457933);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(4,0.06897145);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(5,0.08932795);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(6,0.06975923);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetDirectory(0);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleOffset(1.3);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->Draw("sameaxis");
lowerPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
void probabilityDirection(){
int _vertex = 0;
int _pdg[MAXN];
float _probability[MAXN];
float _chi2[MAXN];
int _generation[MAXN];
gStyle->SetCanvasPreferGL(kTRUE);
TCanvas * c1 = new TCanvas("c1", "The 3d view",0,0,1000,500);
gStyle->SetCanvasColor(kWhite);
gStyle->SetPadColor(kWhite);
TPad *boxPad = new TPad("box", "box", 0.02, 0.02, 0.48, 0.92);
TPad *box1Pad = new TPad("box1", "box1", 0.52, 0.02, 0.98, 0.92);
boxPad->Draw();
box1Pad->Draw();
int nbin = 40;
TH1F * probhist = new TH1F("Phi","Phi",nbin,0,1.01 );
TH1F * probhist_non = new TH1F("Teta","Teta",nbin,0,1.01 );
probhist_non->Sumw2();
TH1F * chihist= new TH1F("PhiD","PhiD",40,0,30 );
TH1F * chihist_non = new TH1F("TetaD","TetaD",40,0,30 );
chihist_non->Sumw2();
TChain* T = new TChain("TaggedVertices");
T->Add("TrashRecoTest.root");
TChain* T2 = new TChain("Stats");
T2->Add("TrashMCTest.root");
int mTotalNumberOfEvents = T->GetEntries();
T->SetBranchAddress("numberOfTagged", &_vertex);
T->SetBranchAddress("PDG", _pdg);
T->SetBranchAddress("generation", _generation);
T->SetBranchAddress("chi2", _chi2);
T->SetBranchAddress("probability", _probability);
int _btotalnumber = 0;
int _bbartotalnumber = 0;
int _bnumber = 0;
int _bbarnumber = 0;
float _bdistance = 0.0;
float _bbardistance = 0.0;
float _bIPdistance = 0.0;
float _bbarIPdistance = 0.0;
T2->SetBranchAddress("btotalnumber", &_btotalnumber);
T2->SetBranchAddress("bbartotalnumber", &_bbartotalnumber);
T2->SetBranchAddress("bnumber", &_bnumber);
T2->SetBranchAddress("bbarnumber", &_bbarnumber);
T2->SetBranchAddress("bdistance", &_bdistance);
T2->SetBranchAddress("bbardistance", &_bbardistance);
T2->SetBranchAddress("bIPdistance", &_bIPdistance);
T2->SetBranchAddress("bbarIPdistance", &_bbarIPdistance);
cout << "mTotalNumberOfEvents: " << mTotalNumberOfEvents << '\n';
int counter = 0;
int totalcounter = 0;
int num = 2;
float distcut = .5;
for (unsigned int i = 0; i < mTotalNumberOfEvents; i += 1)
{
T->GetEntry(i);
T2->GetEntry(i);
vector<int> bs;
vector<int> bbars;
for (int j = 0; j < _vertex; j++)
{
if (_pdg[j] == 5)
{
bs.push_back(j);
}
else//if (_pdg[j] == -5)
{
bbars.push_back(j);
}
}
for (int k = 0; k < bbars.size(); k++)
{
/*if (
_bbarIPdistance < 0.50
|| _bbartotalnumber > 5
|| _bbartotalnumber < 3
|| _bbarnumber < 2
|| _bbartotalnumber - _bbarnumber < 2
)
{
continue;
}*/
if (bbars.size() > 1)
{
std::cout << "Nb: " << i << '\n';
probhist->Fill(_probability[bbars[k]]);
chihist->Fill(_chi2[bbars[k]]);
}
else
{
probhist_non->Fill(_probability[bbars[0]]);
chihist_non->Fill(_chi2[bbars[0]]);
}
}
for (int m = 0; m < bs.size(); m++)
{
/*if (
_bIPdistance < 0.50
|| _btotalnumber > 5
//|| _btotalnumber < 3
|| _bnumber < 2
|| _btotalnumber - _bnumber < 2
)
{
continue;
}*/
if (bs.size() > 1)
{
probhist->Fill(_probability[bs[m]]);
chihist->Fill(_chi2[bs[m]]);
}
else
{
probhist_non->Fill(_probability[bs[0]]);
chihist_non->Fill(_chi2[bs[0]]);
}
}
if ( bbars.size() > 0)
{
totalcounter++;
}
if ( bs.size() > 0)
{
totalcounter++;
}
}
boxPad->cd();
probhist_non->SetMarkerStyle(20);
chihist_non->SetMarkerStyle(20);
chihist_non->SetMarkerColor(kBlack);
chihist_non->SetMarkerSize(0.7);
chihist_non->GetXaxis()->SetTitle("#chi^2");
chihist->SetFillColor(kYellow);
chihist->GetXaxis()->SetTitle("#chi^{2}");
probhist_non->SetMarkerColor(kBlack);
probhist_non->SetMarkerSize(0.7);
probhist_non->GetXaxis()->SetTitle("P(V)");
//probhist->GetYaxis()->SetRange(0,0.1);
//probhist_non->GetYaxis()->SetRange(0,0.1);
probhist->SetMaximum(2000);
probhist_non->SetMaximum(2000);
probhist_non->SetMinimum(0);
probhist->SetMinimum(0);
probhist->SetFillColor(kYellow);
probhist->GetXaxis()->SetTitle("P(V)");
probhist->DrawNormalized("same");
probhist_non->DrawNormalized("ex0psame");
boxPad->SetGrid();
box1Pad->cd();
box1Pad->SetLogy();
chihist->DrawNormalized("");
chihist_non->DrawNormalized("ex0psame",1);
box1Pad->SetGrid();
box1Pad->Modified();
boxPad->Modified();
std::cout << "Statistics: " << probhist->GetEntries() << " with ternary; " << probhist_non->GetEntries() << " without ternary; " << totalcounter << " total; " << counter << " discarded\n";
std::cout << "Mean: " << probhist->GetMean() << " with ternary; " << probhist_non->GetMean() << " without ternary;\n";
}
// ______________________________________________________________________________________
void plotVsEta_nice(const Char_t* name = "ratioVsEta") {
gROOT->LoadMacro("include/toolsEtaNice.C++");
SetupStyle();
SetGlobals();
// -----------------------------------------------------
TFile *inFiles[nEtaSets][nDataSets][nMoments];
TGraphErrors *graphs[nEtaSets][nDataSets][nMoments];
for (int idxEta = 0 ; idxEta < nEtaSets; ++idxEta)
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet)
for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
inFiles[idxEta][idxDataSet][idxMoment] = TFile::Open(Form("output/%s/%s/Moments_%s.root",
aEtaSets[idxEta], aDataSets[idxDataSet], aMoments[idxMoment]));
graphs[idxEta][idxDataSet][idxMoment] = static_cast<TGraphErrors*>((inFiles[idxEta][idxDataSet][idxMoment]->Get(aMoments[idxMoment]))->Clone());
if (inFiles[idxEta][idxDataSet][idxMoment])
(inFiles[idxEta][idxDataSet][idxMoment])->Close();
}
// -----------------------------------------------------
TLegend *legRat[nDataSets];
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
legRat[idxDataSet] = new TLegend(0.14, 0.12, 0.78, 0.27);
legRat[idxDataSet]->SetTextAlign(12);
legRat[idxDataSet]->SetTextSize(0.06);
legRat[idxDataSet]->SetFillColor(1182);
legRat[idxDataSet]->SetLineColor(0);
legRat[idxDataSet]->SetBorderSize(0);
legRat[idxDataSet]->SetNColumns(3);
}
// -----------------------------------------------------
TGraphErrors *etaGraphs[9][nEtaSuperSets][nDataSets][nMoments];
for (int idxCent = 0; idxCent < nCent; idxCent++) {
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
for (int idxMoment = 4 ; idxMoment < nMoments; ++idxMoment) {
int idxEta = -1;
for (int idxEtaSuper = 0 ; idxEtaSuper < nEtaSuperSets; ++idxEtaSuper) {
float x[aEtaSuperSets[idxEtaSuper]];
float y[aEtaSuperSets[idxEtaSuper]];
float ex[aEtaSuperSets[idxEtaSuper]];
float ey[aEtaSuperSets[idxEtaSuper]];
for (int idxEtaSub = 0 ; idxEtaSub < aEtaSuperSets[idxEtaSuper]; ++idxEtaSub) {
++idxEta;
x[idxEtaSub] = aEtaSetBinCenter[idxEta];
ex[idxEtaSub] = aEtaSetBinWidth[idxEta]/2.;
y[idxEtaSub] = graphs[idxEta][idxDataSet][idxMoment]->GetY()[idxCent];
ey[idxEtaSub] = graphs[idxEta][idxDataSet][idxMoment]->GetEY()[idxCent];
} // for (int idxEtaSub = 0 ; idxEtaSub < aEtaSuperSets[idxEtaSuper]; ++idxEtaSub) {
etaGraphs[idxCent][idxEtaSuper][idxDataSet][idxMoment] = new TGraphErrors(aEtaSuperSets[idxEtaSuper], x, y, ex, ey);
PrepareGraph(etaGraphs[idxCent][idxEtaSuper][idxDataSet][idxMoment]);
} // for (int idxEtaSuper = 0 ; idxEtaSuper < nEtaSuperSets; ++idxEtaSuper) {
} // for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
} // for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
} // for (int idxCent = 0; idxCent < 9; idxCent++) {
// -----------------------------------------------------
for (int idxCent = 0; idxCent < 1; idxCent++) {
for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
TPad *pad = SetupCanvas(Form("canEta_Ratio_%s_%s", aDataSets[idxDataSet], cent[idxCent]), aDataSetsTitle[idxDataSet], "#eta", 0.45);
for (int idxEtaSuper = 0 ; idxEtaSuper < nEtaSuperSets; ++idxEtaSuper) {
if (idxEtaSuper == 1 || idxEtaSuper == 3 || idxEtaSuper == 5 || idxEtaSuper == 7 || idxEtaSuper == 8 )
continue;
for (int idxMoment = 4 ; idxMoment < nMoments; ++idxMoment) {
pad->cd(idxMoment-3);
TGraphErrors *g = etaGraphs[idxCent][idxEtaSuper][idxDataSet][idxMoment];
// if (idxCent == 0)
// ShiftGraphX(g, -0.015);
// else if (idxCent == 1)
// ShiftGraphX(g, -0.005);
// else if (idxCent == 4)
// ShiftGraphX(g, 0.005);
// else if (idxCent == 8)
// ShiftGraphX(g, 0.015);
ConfigGraph(g, idxMoment, idxEtaSuper);
if (idxEtaSuper == 0) {
g->Draw("AP");
TLine *line05 = new TLine( 0.5, aMinY[idxMoment], 0.5, aMaxY[idxMoment]);
line05->SetLineColor(kGray+1);
line05->SetLineStyle(3);
line05->Draw();
TLine *line50 = new TLine(-0.5, aMinY[idxMoment], -0.5, aMaxY[idxMoment]);
line50->SetLineColor(kGray+1);
line50->SetLineStyle(3);
line50->Draw();
TLine *line00 = new TLine(0., aMinY[idxMoment],0, aMaxY[idxMoment]);
line00->SetLineColor(kGray+1);
line00->SetLineStyle(3);
line00->Draw();
if (idxMoment == 5) {
// TLine *line0 = new TLine(aMinX, 0, aMaxX, 0);
// line0->SetLineColor(kGray+1);
// line0->SetLineStyle(2);
// line0->SetLineWidth(2);
// line0->Draw();
}
else if (idxMoment == 6) {
TLine *line1 = new TLine(aMinX, 1, aMaxX, 1);
line1->SetLineColor(kGray+1);
line1->SetLineStyle(2);
line1->SetLineWidth(2);
line1->Draw();
}
}
g->Draw("PSAME");
} // for (int idxMoment = 0 ; idxMoment < nMoments; ++idxMoment) {
if (idxCent == 0) {
legRat[idxDataSet]->AddEntry(etaGraphs[0][idxEtaSuper][idxDataSet][4],
Form("#Delta#eta = %.1f", aEtaSuperSetsBinWidth[idxEtaSuper]), "pl");
}
} // for (int idxEta = 0 ; idxEta < nEta; ++idxEta) {
TLine *linex = new TLine(-0.55, 1, 0.55, 1);
linex->SetLineColor(kGray+1);
linex->SetLineStyle(2);
linex->SetLineWidth(2);
legRat[idxDataSet]->AddEntry(linex, "Poisson", "l");
pad->cd(1);
TLatex *texb_3 = new TLatex(-0.49, 9.3, "Au+Au collisions #sqrt{#it{s}_{NN}} = 14.5 GeV");
texb_3->SetTextSize(0.07);
texb_3->Draw("same");
TLatex *texb_3a = new TLatex(-0.49,8.9, "Net-Charge, 0.2 < #it{p}_{T} (GeV/#it{c}) < 2.0, 0-5%");
texb_3a->SetTextSize(0.07);
texb_3a->Draw("same");
pad->cd(2);
TLatex *texb_3b = new TLatex(0.49,0.42, "statistical errors only");
texb_3b->SetTextSize(0.06);
texb_3b->SetTextAlign(31);
texb_3b->Draw("same");
TLatex *texb_3c = new TLatex(0.49,0.38, "horizontal error: #Delta#eta width");
texb_3c->SetTextSize(0.06);
texb_3c->SetTextAlign(31);
texb_3c->Draw("same");
TLatex *texb_4 = new TLatex(-0.49, 0.42, "STAR Preliminary");
texb_4->SetTextSize(0.07);
texb_4->Draw("same");
pad->cd(3);
legRat[idxDataSet]->Draw("lt");
pad->Modified();
} // for (int idxDataSet = 0 ; idxDataSet < nDataSets; ++idxDataSet) {
} // for (int idxCent = 0; idxCent < 9; idxCent++) {
// -----------------------------------------------------
SaveCanvas(name);
}
// draw the same thing but after reco
void genPlots02(std::string fullPath, int nOverlay = 500, bool custBinning = false)
{
const int fVerbose(1);
setTDRStyle();
gStyle->SetOptStat(112211);
gStyle->SetPalette(1);
// Canvas
c = new TCanvas("c2","c2",1000,600);
const unsigned int nPadX = 1;
const unsigned int nPadY = 1;
c->Divide(nPadX,nPadY);
const unsigned int nPads=nPadX*nPadY;
for(unsigned int i=1; i<=nPads; i++)
{
TPad* pad= (TPad*)c->cd(i);
pad->SetTopMargin(0.10);
pad->SetRightMargin(0.20);
pad->SetLeftMargin(0.15);
}
// Open file
TFile *f = TFile::Open(fullPath.c_str());
if (f==0)
{
cout << "File " << fullPath << " not found -- exiting" << endl;
return;
}
if(fVerbose>0)
cout << "Succesfully opened file " << fullPath << endl;
// Get TTree
TTree* t = (TTree*) f->Get("events");
if(fVerbose>0) cout << "Got TTree with " << t->GetEntries() << " entries" << endl;
// Do a cut, if needed
//t->Draw(">>lst","chi2lb>.1&&mlb>5.61&&mlb<5.63");
//t->Draw(">>lst","chi2lb>.1&&isSig==1");
t->Draw(">>lst","(rid1m&4)==4&&(rid2m&4)==4&&mjp>2.895&&mjp<3.295&&prob1m>0.1&&prob2m>0.1&&ptjp>2&&probjp>0.005&&ml0>1.101&&ml0<1.129&&probpr>0.02&&probpi>0.02&&rptpr>rptpi&&ptl0>3&&rptpr>1&&rptpi>0.5&&probl0>0.02&&alphal0<0.3&&d3l0>1&&d3l0/d3El0>10&&problb>0.001&&alphalb<0.3");
TEventList *lst;
lst = (TEventList*)gDirectory->Get("lst");
t->SetEventList(lst);
if(fVerbose>0) cout << "Got TTree with " << t->GetEntries() << " entries" << endl;
// Do plots
c->cd(1);
//doPlot2d(t,"hrzL0vtx", "vrl0:TMath::Abs(vzl0)",30,0,300,30,0,120,"Tit","|z|","r","cm","cm");
if (custBinning)
{
double newbinsX[]={0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50};
const int newbinsX_size = sizeof(newbinsX)/sizeof(double);
std::vector<double> binvecX(newbinsX,newbinsX+newbinsX_size);
//double newbinsY[]={0,1,2,3,4,5,6,7,8,9,10};
double newbinsY[]={0,0.5,1,2,4,8,16,32};
const int newbinsY_size = sizeof(newbinsY)/sizeof(double);
std::vector<double> binvecY(newbinsY,newbinsY+newbinsY_size);
doPlot2d(t,"hrzL0vtxreco", "vrl0:TMath::Abs(vzl0)",binvecX, binvecY,"#Lambda vertices","|z|","r","cm","cm");
}
else
{
doPlot2d(t,"hrzL0vtxreco", "vrl0:TMath::Abs(vzl0)",30,0,50,30,0,30,"#Lambda vertices","|z|","r","cm","cm");
}
// add tracker
TPad* pad;
pad = (TPad*)c->cd(1);
pad->Modified();
pad->Update();
repositionPalette("hrzL0vtxreco");
pad->Update();
pad->SetLogz();
drawTracker(pad);
if (nOverlay<=0) return;
int maxN = nOverlay;
if (maxN > t->GetEntries()) maxN = t->GetEntries();
double vrl0,vzl0,ppr,ppi,etapr,etapi;
t->SetBranchAddress("vrl0",&vrl0);
t->SetBranchAddress("vzl0",&vzl0);
t->SetBranchAddress("ppr",&ppr);
t->SetBranchAddress("etapr",&etapr);
t->SetBranchAddress("ppi",&ppi);
t->SetBranchAddress("etapi",&etapi);
double scalepr = 4;
double scalepi = 8;
{ // reference indicator
const double x1pr = 0; const double y1pr = -3;
const double x2pr = scalepr; const double y2pr = y1pr;
const double versatz = 14;
const double x1pi = x1pr+versatz; const double y1pi = -3;
const double x2pi = x2pr+versatz+scalepi; const double y2pi = y1pi;
TArrow *a;
a = new TArrow(x1pr,y1pr,x2pr,y2pr,.01,">");
a->SetLineColor(24);
a->Draw();
TLatex tl;
tl.SetTextSize(20);
tl.SetTextFont(4);
tl.DrawLatex(x1pr,y2pr-1.2,"p(p) / 1 GeV");
a = new TArrow(x1pi,y1pi,x2pi,y2pi,.01,">");
a->SetLineColor(20);
a->Draw();
tl.SetTextSize(20);
tl.SetTextFont(4);
tl.DrawLatex(x1pi,y2pi-1.2,"p(#pi) / 1 GeV");
}
for (int i = 0; i!=maxN; i++)
{
t->GetEntry(i);
const double thetapr = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapr)));
const double thetapi = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapi)));
const double x1=TMath::Abs(vzl0);
const double y1=vrl0;
const double x2pr=x1+scalepr*ppr*TMath::Cos(thetapr);
const double y2pr=y1+scalepr*ppr*TMath::Sin(thetapr);
const double x2pi=x1+scalepi*ppi*TMath::Cos(thetapi);
const double y2pi=y1+scalepi*ppi*TMath::Sin(thetapi);
TArrow *a;
a = new TArrow(x1,y1,x2pr,y2pr,.01,">");
a->SetLineColor(24);
a->Draw();
a = new TArrow(x1,y1,x2pi,y2pi,.01,">");
a->SetLineColor(20);
a->Draw();
TMarker *m = new TMarker(x1,y1,7);
m->SetMarkerColor(28);
m->Draw();
}
}
void SingleAnalysis::DrawGraphs()
{
std::cout << "Drawing graphs" << std::endl;
// special case for erec and deviation
m_pERecCanvas->cd();
TPad *pERecPad = new TPad("ERecPad", "ERecPad", 0.02, 0.37, 0.95, 0.99);
pERecPad->Draw();
pERecPad->cd();
pERecPad->Range(-19.23077,0.01,96.15385,1428.571);
pERecPad->SetFillColor(0);
pERecPad->SetBorderMode(0);
pERecPad->SetBorderSize(2);
pERecPad->SetTickx(1);
pERecPad->SetTicky(1);
pERecPad->SetLeftMargin(0.1666667);
pERecPad->SetRightMargin(0.01);
pERecPad->SetTopMargin(0.02);
pERecPad->SetBottomMargin(0.0);
pERecPad->SetFrameLineWidth(1);
pERecPad->SetFrameBorderMode(0);
m_canvasMultiGraphMap[E_REC].second->Draw("ap");
this->PostDrawMultiGraph(E_REC, m_canvasMultiGraphMap[E_REC].second);
TLegend *pLegend = pERecPad->BuildLegend();
pLegend->SetBorderSize(0);
pLegend->SetLineColor(0);
pLegend->SetLineStyle(0);
pLegend->SetLineWidth(1);
pLegend->SetFillColor(0);
pLegend->SetFillStyle(0);
pLegend->SetTextSize(0.04);
pERecPad->Modified();
m_pERecCanvas->cd();
TPad *pERecDeviationPad = new TPad("ERecDeviationPad", "ERecDeviationPad", 0.02, 0, 0.95, 0.37);
pERecDeviationPad->Draw();
pERecDeviationPad->cd();
pERecDeviationPad->Range(-19.23077,-0.3993507,96.15385,0.25);
pERecDeviationPad->SetFillColor(0);
pERecDeviationPad->SetBorderMode(0);
pERecDeviationPad->SetBorderSize(2);
pERecDeviationPad->SetTickx(1);
pERecDeviationPad->SetTicky(1);
pERecDeviationPad->SetLeftMargin(0.1666667);
pERecDeviationPad->SetRightMargin(0.01);
pERecDeviationPad->SetTopMargin(0);
pERecDeviationPad->SetBottomMargin(0.23);
pERecDeviationPad->SetFrameLineWidth(1);
pERecDeviationPad->SetFrameBorderMode(0);
m_canvasMultiGraphMap[E_REC_DEVIATION].second->Draw("ap");
TF1 *pLin = new TF1("lin", "0", 0, 90);
pLin->SetLineColor(kBlack);
pLin->SetLineWidth(1);
pLin->SetLineStyle(2);
pLin->Draw("same");
this->PostDrawMultiGraph(E_REC_DEVIATION, m_canvasMultiGraphMap[E_REC_DEVIATION].second);
pERecDeviationPad->Modified();
for(CanvasMultiGraphMap::iterator graphIter = m_canvasMultiGraphMap.begin(), graphEndIter = m_canvasMultiGraphMap.end() ;
graphEndIter != graphIter ; ++graphIter)
{
if(graphIter->first == E_REC || graphIter->first == E_REC_DEVIATION)
continue;
graphIter->second.first->cd();
graphIter->second.second->Draw("ap");
this->PostDrawMultiGraph(graphIter->first, graphIter->second.second);
TLegend *pLegend = graphIter->second.first->BuildLegend();
pLegend->SetBorderSize(0);
pLegend->SetLineColor(0);
pLegend->SetLineStyle(1);
pLegend->SetLineWidth(1);
pLegend->SetFillColor(0);
pLegend->SetFillStyle(0);
pLegend->SetTextSize(0.04);
}
}
void plot(int mass) {
double myQCDRelUncert = 0.038;
double myEWKRelUncert = 0.131;
double myFakesRelUncert = 0.238;
double delta = 1.4;
double br = 0.05;
bool debug = false;
bool log = false;
double ymin = 0.001;
double ymax = 48;
static bool bMessage = false;
if (!bMessage) {
cout << "Values used as relative uncertainty (please check):" << endl;
cout << " QCD: " << myQCDRelUncert << endl;
cout << " EWK genuine tau: " << myEWKRelUncert << endl;
cout << " EWK fake tau: " << myFakesRelUncert << endl << endl;
bMessage = true;
}
cout << "Processing mass point: " << mass << " GeV/c2" << endl;
gStyle->SetOptFit(1);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
gStyle->SetTitleFont(43, "xyz");
gStyle->SetTitleSize(33, "xyz");
gStyle->SetLabelFont(43, "xyz");
gStyle->SetLabelSize(27, "xyz");
//std::string infile = "EPS_data_nodeltaphi/hplus_100.root";
//std::string infile = "EPS_data_deltaphi160/hplus_100.root";
std::stringstream s;
s << "lands_histograms_hplushadronic_m" << mass << ".root";
std::string infile = s.str();
// Canvas
TCanvas *myCanvas = new TCanvas("myCanvas", "",0,0,600,600);
myCanvas->SetHighLightColor(2);
myCanvas->Range(0,0,1,1);
myCanvas->SetFillColor(0);
myCanvas->SetBorderMode(0);
myCanvas->SetBorderSize(2);
if (log)
myCanvas->SetLogy();
myCanvas->SetTickx(1);
myCanvas->SetTicky(1);
myCanvas->SetLeftMargin(0.16);
myCanvas->SetRightMargin(0.05);
myCanvas->SetTopMargin(0.05);
myCanvas->SetBottomMargin(0.08);
myCanvas->SetFrameFillStyle(0);
myCanvas->SetFrameBorderMode(0);
myCanvas->SetFrameFillStyle(0);
myCanvas->SetFrameBorderMode(0);
myCanvas->cd();
Int_t ci;
TFile* f = TFile::Open(infile.c_str());
s.str("");
s << "HW" << mass << "_1";
TH1* hw = (TH1*)f->Get(s.str().c_str());
s.str("");
s << "HH" << mass << "_1";
TH1* hh = (TH1*)f->Get(s.str().c_str());
TH1* data = (TH1*)f->Get("data_obs");
data->SetLineWidth(2);
data->SetMarkerStyle(20);
data->SetMarkerSize(1.2);
TH1* ewktau = (TH1*)f->Get("EWK_Tau");
ci = TColor::GetColor("#993399");
ewktau->SetFillColor(ci);
ewktau->SetLineWidth(0);
TH1* ewkDY = (TH1*)f->Get("EWK_DYx");
TH1* ewkVV = (TH1*)f->Get("EWK_VVx");
ewktau->Add(ewkDY);
ewktau->Add(ewkVV);
//TH1* qcd = (TH1*)f->Get("QCDInv");
TH1* qcd = (TH1*)f->Get("QCD");
ci = TColor::GetColor("#ffcc33");
qcd->SetFillColor(ci);
qcd->SetLineWidth(0);
TH1* fakett = (TH1*)f->Get("fake_tt");
ci = TColor::GetColor("#669900");
fakett->SetFillColor(ci);
fakett->SetLineWidth(0);
TH1* fakeW = (TH1*)f->Get("fake_W");
ci = TColor::GetColor("#cc3300");
fakeW->SetFillColor(ci);
fakeW->SetLineWidth(0);
TH1* faket = (TH1*)f->Get("fake_t");
TH1F *hFrame = new TH1F("hFrame","",20,0,400);
hFrame->SetMinimum(ymin);
if (log)
hFrame->SetMaximum(ymax*1.5);
else
hFrame->SetMaximum(ymax);
hFrame->SetDirectory(0);
hFrame->SetStats(0);
hFrame->SetLineStyle(0);
hFrame->SetMarkerStyle(20);
hFrame->SetXTitle("Transverse mass (#tau jet, E_{T}^{miss}), (GeV/c^{2})");
if (paperStatus)
hFrame->SetXTitle("Transverse mass (#tau_{h}, E_{T}^{miss}), (GeV/c^{2})");
hFrame->SetYTitle("Events / 20 GeV/c^{2}");
hFrame->GetXaxis()->SetTitleSize(0);
hFrame->GetXaxis()->SetLabelSize(0);
hFrame->GetYaxis()->SetTitleFont(43);
hFrame->GetYaxis()->SetTitleSize(27);
hFrame->GetYaxis()->SetTitleOffset(1.3);
// signal
hh->Scale(br*br);
hw->Scale(2*br*(1.0-br));
TH1* signal = (TH1*)hh->Clone();
signal->Add(hw);
ci = TColor::GetColor("#ff3399");
signal->SetLineColor(ci);
signal->SetLineStyle(2);
signal->SetLineWidth(2);
// Fakes
TH1* fakes = (TH1*)(fakett->Clone());
fakes->Add(fakeW);
fakes->Add(faket);
// stacked backgrounds
THStack *exp = new THStack();
exp->SetName("exp");
exp->SetTitle("exp");
exp->Add(fakes);
exp->Add(ewktau);
exp->Add(qcd);
exp->Add(signal);
TH1* hExpBkg = (TH1*)fakes->Clone();
hExpBkg->Add(ewktau);
hExpBkg->Add(qcd);
// uncertainty
TH1* uncert = (TH1*)fakeW->Clone();
uncert->Add(fakett);
uncert->Add(ewktau);
uncert->Add(qcd);
uncert->SetFillColor(1);
uncert->SetFillStyle(3344);
uncert->SetLineColor(0);
uncert->SetLineStyle(0);
uncert->SetLineWidth(0);
TH1* hExpBkgTotalUncert = (TH1*)uncert->Clone();
hExpBkgTotalUncert->SetFillStyle(3354);
TH1* hAgreement = (TH1*)data->Clone();
hAgreement->Divide(hExpBkg);
TGraphErrors* hAgreementRelUncert = new TGraphErrors(hAgreement->GetNbinsX());
hAgreementRelUncert->SetLineWidth(2);
hAgreementRelUncert->SetLineColor(kBlack);
for (int i = 1; i <= hFrame->GetNbinsX(); ++i) {
double myQCDTotalUncert = TMath::Power(qcd->GetBinError(i), 2)
+ TMath::Power(qcd->GetBinContent(i)*myQCDRelUncert, 2);
double myEWKTotalUncert = TMath::Power(ewktau->GetBinError(i), 2)
+ TMath::Power(ewktau->GetBinContent(i)*myEWKRelUncert, 2);
double myFakesTotalUncert = TMath::Power(fakes->GetBinError(i), 2)
+ TMath::Power(fakes->GetBinContent(i)*myFakesRelUncert, 2);
hExpBkgTotalUncert->SetBinError(i, TMath::Sqrt(myQCDTotalUncert + myEWKTotalUncert + myFakesTotalUncert));
if (hExpBkg->GetBinContent(i) > 0) {
hAgreementRelUncert->SetPoint(i-1, hExpBkg->GetBinCenter(i), data->GetBinContent(i) / hExpBkg->GetBinContent(i));
double myUncertData = 0;
if (data->GetBinContent(i) > 0)
myUncertData = TMath::Power(data->GetBinError(i) / data->GetBinContent(i), 2);
double myUncertBkg = (myQCDTotalUncert + myEWKTotalUncert + myFakesTotalUncert) / TMath::Power(hExpBkg->GetBinContent(i), 2);
hAgreementRelUncert->SetPointError(i-1, 0, data->GetBinContent(i) / hExpBkg->GetBinContent(i) * TMath::Sqrt(myUncertData + myUncertBkg));
} else {
hAgreementRelUncert->SetPoint(i-1, hExpBkg->GetBinCenter(i), 0);
hAgreementRelUncert->SetPointError(i-1, 0, 0);
}
if (debug) {
cout << "Point: " << hAgreementRelUncert->GetX()[i-1]-10 << "-" << hAgreementRelUncert->GetX()[i-1]+10
<< " GeV/c2, agreement: " << hAgreementRelUncert->GetY()[i-1] << ", uncert: " << hAgreement->GetBinError(i) << ", " << hAgreementRelUncert->GetErrorY(i-1) << endl;
cout << " bkg. stat. uncert. " << hExpBkg->GetBinError(i) << " (i.e. " << hExpBkg->GetBinError(i) / hExpBkg->GetBinContent(i) * 100.0 << " %)"
<< ", stat+syst uncert. " << TMath::Sqrt(myQCDTotalUncert + myEWKTotalUncert + myFakesTotalUncert)
<< " (i.e. " << TMath::Sqrt(myQCDTotalUncert + myEWKTotalUncert + myFakesTotalUncert) / hExpBkg->GetBinContent(i) * 100.0 << " %)" << endl;
}
}
// Agreement pad
TPad* pad = new TPad("ratiopad","ratiopad",0.,0.,1.,.3);
pad->Draw();
pad->cd();
pad->Range(0,0,1,1);
pad->SetFillColor(0);
pad->SetFillStyle(4000);
pad->SetBorderMode(0);
pad->SetBorderSize(2);
pad->SetTickx(1);
pad->SetTicky(1);
pad->SetLeftMargin(0.16);
pad->SetRightMargin(0.05);
pad->SetTopMargin(0);
pad->SetBottomMargin(0.34);
pad->SetFrameFillStyle(0);
pad->SetFrameBorderMode(0);
// Plot here ratio
if (1.0-delta > 0)
hAgreement->SetMinimum(1.0-delta);
else
hAgreement->SetMinimum(0.);
hAgreement->SetMaximum(1.0+delta);
hAgreement->GetXaxis()->SetLabelOffset(0.007);
hAgreement->GetXaxis()->SetLabelFont(43);
hAgreement->GetXaxis()->SetLabelSize(27);
hAgreement->GetYaxis()->SetLabelFont(43);
hAgreement->GetYaxis()->SetLabelSize(27);
hAgreement->GetYaxis()->SetLabelOffset(0.007);
hAgreement->GetYaxis()->SetNdivisions(505);
hAgreement->GetXaxis()->SetTitleFont(43);
hAgreement->GetYaxis()->SetTitleFont(43);
hAgreement->GetXaxis()->SetTitleSize(33);
hAgreement->GetYaxis()->SetTitleSize(33);
hAgreement->SetTitleSize(27, "xyz");
hAgreement->GetXaxis()->SetTitleOffset(3.2);
hAgreement->GetYaxis()->SetTitleOffset(1.3);
hAgreement->SetXTitle(hFrame->GetXaxis()->GetTitle());
hAgreement->SetYTitle("Data/#Sigmabkg");
hAgreement->Draw("e2");
// Plot line at zero
TH1* hAgreementLine = dynamic_cast<TH1*>(hAgreement->Clone());
for (int i = 1; i <= hAgreementLine->GetNbinsX(); ++i) {
hAgreementLine->SetBinContent(i,1.0);
hAgreementLine->SetBinError(i,0.0);
}
hAgreementLine->SetLineColor(kRed);
hAgreementLine->SetLineWidth(2);
hAgreementLine->SetLineStyle(3);
hAgreementLine->Draw("hist same");
hAgreement->Draw("same");
hAgreementRelUncert->Draw("[]");
pad->RedrawAxis();
myCanvas->cd();
TPad* plotpad = new TPad("plotpad", "plotpad",0,0.3,1.,1.);
plotpad->Draw();
plotpad->cd();
plotpad->Range(0,0,1,1);
plotpad->SetFillColor(0);
plotpad->SetFillStyle(4000);
plotpad->SetBorderMode(0);
plotpad->SetBorderSize(2);
//if (logy)
// plotpad->SetLogy();
plotpad->SetTickx(1);
plotpad->SetTicky(1);
plotpad->SetLeftMargin(0.16);
plotpad->SetRightMargin(0.05);
plotpad->SetTopMargin(0.065);
plotpad->SetBottomMargin(0.0);
plotpad->SetFrameFillStyle(0);
plotpad->SetFrameBorderMode(0);
hFrame->GetXaxis()->SetTitleSize(0);
hFrame->GetXaxis()->SetLabelSize(0);
hFrame->GetYaxis()->SetTitleFont(43);
hFrame->GetYaxis()->SetTitleSize(33);
hFrame->GetYaxis()->SetTitleOffset(1.3);
// Draw objects
hFrame->Draw();
exp->Draw("hist same");
uncert->Draw("E2 same");
hExpBkgTotalUncert->Draw("E2 same");
// Data
data->Draw("same");
//signal->Draw("same");
TLegend *leg = new TLegend(0.53,0.6,0.87,0.91,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(63);
leg->SetTextSize(18);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(kWhite);
//leg->SetFillStyle(4000); // enabling this will cause the plot to be erased from the pad
TLegendEntry* entry = leg->AddEntry(data, "Data", "P");
s.str("");
s << "with H^{#pm}#rightarrow#tau^{#pm}#nu";
entry = leg->AddEntry(signal, s.str().c_str(), "L");
entry = leg->AddEntry(qcd, "QCD (meas.)", "F");
entry = leg->AddEntry(ewktau, "EWK genuine #tau (meas.)", "F");
entry = leg->AddEntry(fakes, "EWK fake #tau (MC)", "F");
entry = leg->AddEntry(uncert, "stat. uncert.", "F");
entry = leg->AddEntry(hExpBkgTotalUncert, "stat. #oplus syst. uncert.", "F");
leg->Draw();
string myTitle = "CMS Preliminary";
if (paperStatus)
myTitle = "CMS";
TLatex *tex = new TLatex(0.62,0.945,myTitle.c_str());
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(27);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.2,0.945,"#sqrt{s} = 7 TeV");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(27);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.43,0.945,"2.2 fb^{-1}");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(27);
tex->SetLineWidth(2);
tex->Draw();
s.str("");
s << "m_{H^{#pm}} = " << mass << " GeV/c^{2}";
tex = new TLatex(0.28,0.865,s.str().c_str());
tex->SetNDC();
tex->SetTextFont(63);
tex->SetTextSize(20);
tex->SetLineWidth(2);
tex->Draw();
s.str("");
s << "BR(t#rightarrowbH^{#pm})=" << setprecision(2) << br;
tex = new TLatex(0.28,0.805,s.str().c_str());
tex->SetNDC();
tex->SetTextFont(63);
tex->SetTextSize(20);
tex->SetLineWidth(2);
tex->Draw();
plotpad->RedrawAxis();
plotpad->Modified();
s.str("");
s << "mT_datadriven_m" << mass << ".png";
myCanvas->Print(s.str().c_str());
s.str("");
s << "mT_datadriven_m" << mass << ".C";
myCanvas->Print(s.str().c_str());
s.str("");
s << "mT_datadriven_m" << mass << ".eps";
myCanvas->Print(s.str().c_str());
}
TCanvas *drawRatioPlot(TH1D *prediction, TH1D *sr, TString xTitle, TString filename){
TCanvas *c = new TCanvas("c"+filename,"c",0,0,500,500);
float y = 0.3;
TPad *pad1 = new TPad("pad1", "Control Plots 1", 0.01, y, 0.99, 0.99);
TPad *padRatio = new TPad("rp1", "Ratio1", 0.01, 0.01, 0.99, y-0.01);
TH1D *ratio = 0;
//ratio = (TH1D*) sr->Clone();
//ratio->Divide(prediction);
ratio = (TH1D*) prediction->Clone();
ratio->Divide(sr);
for(int i=1; i<=ratio->GetNbinsX();i++){
cout<<"Scaling Factor in "<<i<<". bin ="<<ratio->GetBinContent(i)<<" +/- "<<ratio->GetBinError(i)<<endl;
}
ratio->GetYaxis()->SetTitle("Pred/Bkg");
ratio->SetTitle("");
ratio->SetLabelSize(0.1,"X");
ratio->SetLabelSize(0.1,"Y");
ratio->SetTitleOffset(0.5,"Y");
ratio->SetTitleSize(0.15,"Y");
padRatio->cd();
//ratio->GetYaxis()->SetRangeUser(0,2);
ratio->Draw();
// Draw line at one!
float xmin = ratio->GetXaxis()->GetXmin();
float xmax = ratio->GetXaxis()->GetXmax();
TLine *line = new TLine(xmin,1,xmax,1);
line->SetLineWidth(2);
line->Draw("same");
padRatio->Modified();
TLegend *leg = new TLegend(0.6,0.8,0.9,0.9);
leg->AddEntry(sr,"bkg","l");
leg->AddEntry(prediction,"bkg prediction","pel");
pad1->cd();
pad1->SetLogy();
sr->SetLineColor(kRed);
sr->SetMarkerColor(kRed);
sr->SetMarkerStyle(20);
sr->SetTitle("");
prediction->SetMarkerStyle(20);
prediction->SetTitle("");
prediction->GetXaxis()->SetTitle(xTitle);
sr->GetXaxis()->SetTitle(xTitle);
prediction->SetTitleSize(0.07,"X");
prediction->GetXaxis()->SetTitleOffset(0.7);
sr->SetTitleSize(0.07,"X");
sr->GetXaxis()->SetTitleOffset(0.7);
prediction->Draw("e");
sr->Draw("e same");
double maximum =prediction->GetMaximum()*2.5;
double minimum = 0.00001;
if(sr->GetMinimum()!=0) minimum=sr->GetMinimum()*0.5;
if(prediction->GetMinimum()!=0) minimum=prediction->GetMinimum()*0.5;
if(sr->GetMaximum()>prediction->GetMaximum()){
maximum=sr->GetMaximum()*2.5;
}
prediction->GetYaxis()->SetRangeUser(minimum,maximum);
leg->Draw("same");
// Draw both pads to canvas
c->cd();
pad1->Draw();
padRatio->SetGridy();
padRatio->Draw();
c->SaveAs(filename);
return c;
}
void wmj_fit_stack()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Tue May 3 07:35:30 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,700,500);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(-0.1730123,-5358.989,1.06279,32919.5);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.14);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3005[20] = {
0.025,
0.075,
0.125,
0.175,
0.225,
0.275,
0.325,
0.375,
0.425,
0.475,
0.525,
0.575,
0.625,
0.675,
0.725,
0.775,
0.825,
0.875,
0.925,
0.975};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3005[20] = {
0,
0,
9,
26,
118,
315,
881,
1812,
3086,
4575,
5856,
6537,
6824,
6348,
5662,
4769,
4941,
6726,
13179,
28019};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3005[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3005[20] = {
0,
0,
2.943461,
5.066015,
10.37428,
17.25528,
29.18586,
42.07053,
55.05403,
67.14059,
76.02614,
80.35326,
82.10902,
79.17591,
74.74792,
68.55976,
69.79403,
81.51372,
114.3009,
166.8895};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3005[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3005[20] = {
1.147874,
1.147874,
4.110204,
6.164324,
11.37428,
18.25528,
30.18586,
43.07053,
56.05403,
68.14059,
77.02614,
81.35326,
83.10902,
80.17591,
75.74792,
69.55976,
70.79403,
82.51372,
115.3009,
167.8895};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(20,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3005);
grae->SetName("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]");
grae->SetTitle("Histogram of combData_plot__muminus_ISO_");
grae->SetFillColor(1);
grae->SetFillStyle(0);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005 = new TH1F("Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMinimum(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetDirectory(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetStats(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetLineWidth(2);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005);
grae->Draw("e1pa");
Double_t h5_fx7[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t h5_fy7[47] = {
0,
0,
0,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
3.789762,
3.789762,
71.12063,
71.12063,
211.4394,
211.4394,
857.7007,
857.7007,
1798.534,
1798.534,
3099.609,
3099.609,
4425.901,
4425.901,
5601.584,
5601.584,
6533.687,
6533.687,
6910.208,
6910.208,
6392.171,
6392.171,
5777.193,
5777.193,
5020.919,
5020.919,
5207.834,
5207.834,
6798.159,
6798.159,
13033.98,
13033.98,
27939.27,
27939.27,
0,
0,
0,
0};
TGraph *graph = new TGraph(47,h5_fx7,h5_fy7);
graph->SetName("h5");
graph->SetTitle("Projection of ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#3333ff");
graph->SetFillColor(ci);
ci = TColor::GetColor("#3333ff");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#3333ff");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_h57 = new TH1F("Graph_h57","Projection of ",100,-0.1292683,1.129268);
Graph_h57->SetMinimum(0);
Graph_h57->SetMaximum(30733.2);
Graph_h57->SetDirectory(0);
Graph_h57->SetStats(0);
Graph_h57->SetLineWidth(2);
Graph_h57->SetMarkerStyle(0);
Graph_h57->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_h57->GetXaxis()->SetRange(11,90);
Graph_h57->GetXaxis()->SetNdivisions(1005);
Graph_h57->GetXaxis()->SetLabelFont(132);
Graph_h57->GetXaxis()->SetLabelOffset(0.02);
Graph_h57->GetXaxis()->SetLabelSize(0.05);
Graph_h57->GetXaxis()->SetTitleSize(0.06);
Graph_h57->GetXaxis()->SetTitleFont(132);
Graph_h57->GetYaxis()->SetTitle("Events / (0.05)");
Graph_h57->GetYaxis()->SetLabelFont(132);
Graph_h57->GetYaxis()->SetLabelSize(0.05);
Graph_h57->GetYaxis()->SetTitleSize(0.06);
Graph_h57->GetYaxis()->SetTitleFont(132);
Graph_h57->GetZaxis()->SetLabelFont(132);
Graph_h57->GetZaxis()->SetLabelSize(0.05);
Graph_h57->GetZaxis()->SetTitleSize(0.06);
Graph_h57->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_h57);
graph->Draw("f");
Double_t sumBkgPdf_totbkg_fitW_fx8[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t sumBkgPdf_totbkg_fitW_fy8[47] = {
0,
0,
0,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
3.789757,
3.789757,
71.12063,
71.12063,
211.4394,
211.4394,
855.0132,
855.0132,
1794.513,
1794.513,
3088.88,
3088.88,
4399.772,
4399.772,
5568.272,
5568.272,
6460.972,
6460.972,
6790.725,
6790.725,
6142.168,
6142.168,
5333.208,
5333.208,
4062.953,
4062.953,
3103.771,
3103.771,
2394.858,
2394.858,
1734.507,
1734.507,
1650.883,
1650.883,
0,
0,
0,
0};
graph = new TGraph(47,sumBkgPdf_totbkg_fitW_fx8,sumBkgPdf_totbkg_fitW_fy8);
graph->SetName("sumBkgPdf_totbkg_fitW");
graph->SetTitle("Projection of ");
ci = TColor::GetColor("#ffff66");
graph->SetFillColor(ci);
ci = TColor::GetColor("#ffff66");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#ffff66");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618 = new TH1F("Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618","Projection of ",100,-0.1292683,1.129268);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMinimum(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMaximum(7469.797);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetDirectory(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetStats(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetLineWidth(2);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetRange(11,90);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618);
graph->Draw("f");
Double_t sumPdf_Norm[muminus_ISO_]_fx9[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t sumPdf_Norm[muminus_ISO_]_fy9[47] = {
0,
0,
0,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
3.789756,
3.789756,
71.12054,
71.12054,
211.4391,
211.4391,
857.6998,
857.6998,
1798.532,
1798.532,
3099.606,
3099.606,
4425.897,
4425.897,
5601.578,
5601.578,
6533.68,
6533.68,
6910.201,
6910.201,
6392.164,
6392.164,
5777.187,
5777.187,
5020.914,
5020.914,
5207.829,
5207.829,
6798.152,
6798.152,
13033.97,
13033.97,
27939.24,
27939.24,
0,
0,
0,
0};
graph = new TGraph(47,sumPdf_Norm[muminus_ISO_]_fx9,sumPdf_Norm[muminus_ISO_]_fy9);
graph->SetName("sumPdf_Norm[muminus_ISO_]");
graph->SetTitle("Projection of ");
ci = TColor::GetColor("#ff0000");
graph->SetFillColor(ci);
graph->SetFillStyle(0);
ci = TColor::GetColor("#ff0000");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#ff0000");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839 = new TH1F("Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839","Projection of ",100,-0.1292683,1.129268);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMinimum(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMaximum(30733.17);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetDirectory(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetStats(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetLineWidth(2);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetRange(11,90);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839);
graph->Draw("");
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3006[20] = {
0.025,
0.075,
0.125,
0.175,
0.225,
0.275,
0.325,
0.375,
0.425,
0.475,
0.525,
0.575,
0.625,
0.675,
0.725,
0.775,
0.825,
0.875,
0.925,
0.975};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3006[20] = {
0,
0,
9,
26,
118,
315,
881,
1812,
3086,
4575,
5856,
6537,
6824,
6348,
5662,
4769,
4941,
6726,
13179,
28019};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3006[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3006[20] = {
0,
0,
2.943461,
5.066015,
10.37428,
17.25528,
29.18586,
42.07053,
55.05403,
67.14059,
76.02614,
80.35326,
82.10902,
79.17591,
74.74792,
68.55976,
69.79403,
81.51372,
114.3009,
166.8895};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3006[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3006[20] = {
1.147874,
1.147874,
4.110204,
6.164324,
11.37428,
18.25528,
30.18586,
43.07053,
56.05403,
68.14059,
77.02614,
81.35326,
83.10902,
80.17591,
75.74792,
69.55976,
70.79403,
82.51372,
115.3009,
167.8895};
grae = new TGraphAsymmErrors(20,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3006);
grae->SetName("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]");
grae->SetTitle("Histogram of combData_plot__muminus_ISO_");
grae->SetFillColor(1);
grae->SetFillStyle(0);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006 = new TH1F("Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMinimum(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetDirectory(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetStats(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetLineWidth(2);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006);
grae->Draw("e1p");
TLegend *leg = new TLegend(0.2,0.53,0.6,0.78,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]","Data, #sqrt{s}=8 TeV","lep");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("h5","Electroweak","f");
ci = TColor::GetColor("#3333ff");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
ci = TColor::GetColor("#3333ff");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("sumBkgPdf_totbkg_fitW","QCD","f");
ci = TColor::GetColor("#ffff66");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
ci = TColor::GetColor("#ffff66");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("sumPdf_Norm[muminus_ISO_]","Fit","l");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.2,0.8,0.4,0.87,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb");
pt->Draw();
TH1F *Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy = new TH1F("Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMinimum(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetDirectory(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetStats(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetLineWidth(2);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
void test3()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Fri May 6 10:50:13 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,700,500);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(1.567901,0.02811051,4.654321,0.1367104);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.14);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
TH1F *test3_red_0_0_norm__1 = new TH1F("test3_red_0_0_norm__1","ETA",10,2,4.5);
test3_red_0_0_norm__1->SetBinContent(1,0.1040398);
test3_red_0_0_norm__1->SetBinContent(2,0.1186439);
test3_red_0_0_norm__1->SetBinContent(3,0.1256071);
test3_red_0_0_norm__1->SetBinContent(4,0.1229555);
test3_red_0_0_norm__1->SetBinContent(5,0.118442);
test3_red_0_0_norm__1->SetBinContent(6,0.1127067);
test3_red_0_0_norm__1->SetBinContent(7,0.09951538);
test3_red_0_0_norm__1->SetBinContent(8,0.08332527);
test3_red_0_0_norm__1->SetBinContent(9,0.06651872);
test3_red_0_0_norm__1->SetBinContent(10,0.0482457);
test3_red_0_0_norm__1->SetBinError(1,0.001708584);
test3_red_0_0_norm__1->SetBinError(2,0.001736969);
test3_red_0_0_norm__1->SetBinError(3,0.001484444);
test3_red_0_0_norm__1->SetBinError(4,0.001426268);
test3_red_0_0_norm__1->SetBinError(5,0.001336676);
test3_red_0_0_norm__1->SetBinError(6,0.00122978);
test3_red_0_0_norm__1->SetBinError(7,0.00125216);
test3_red_0_0_norm__1->SetBinError(8,0.00113555);
test3_red_0_0_norm__1->SetBinError(9,0.0009819959);
test3_red_0_0_norm__1->SetBinError(10,0.0009418311);
test3_red_0_0_norm__1->SetEntries(58129);
test3_red_0_0_norm__1->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm__1->SetFillColor(ci);
test3_red_0_0_norm__1->SetFillStyle(0);
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm__1->SetLineColor(ci);
test3_red_0_0_norm__1->SetLineWidth(2);
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm__1->SetMarkerColor(ci);
test3_red_0_0_norm__1->SetMarkerStyle(20);
test3_red_0_0_norm__1->GetXaxis()->SetTitle("#eta");
test3_red_0_0_norm__1->GetXaxis()->SetNdivisions(1005);
test3_red_0_0_norm__1->GetXaxis()->SetLabelFont(132);
test3_red_0_0_norm__1->GetXaxis()->SetLabelOffset(0.02);
test3_red_0_0_norm__1->GetXaxis()->SetLabelSize(0.05);
test3_red_0_0_norm__1->GetXaxis()->SetTitleSize(0.06);
test3_red_0_0_norm__1->GetXaxis()->SetTitleFont(132);
test3_red_0_0_norm__1->GetYaxis()->SetNdivisions(505);
test3_red_0_0_norm__1->GetYaxis()->SetLabelFont(132);
test3_red_0_0_norm__1->GetYaxis()->SetLabelSize(0.05);
test3_red_0_0_norm__1->GetYaxis()->SetTitleSize(0.06);
test3_red_0_0_norm__1->GetYaxis()->SetTitleFont(132);
test3_red_0_0_norm__1->GetZaxis()->SetLabelFont(132);
test3_red_0_0_norm__1->GetZaxis()->SetLabelSize(0.05);
test3_red_0_0_norm__1->GetZaxis()->SetTitleSize(0.06);
test3_red_0_0_norm__1->GetZaxis()->SetTitleFont(132);
test3_red_0_0_norm__1->Draw("");
Double_t xAxis1[11] = {2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4, 4.25, 4.5};
TH1F *test3_green_0_1_norm__2 = new TH1F("test3_green_0_1_norm__2","ETA",10, xAxis1);
test3_green_0_1_norm__2->SetBinContent(1,0.1214835);
test3_green_0_1_norm__2->SetBinContent(2,0.1235428);
test3_green_0_1_norm__2->SetBinContent(3,0.1252758);
test3_green_0_1_norm__2->SetBinContent(4,0.1207572);
test3_green_0_1_norm__2->SetBinContent(5,0.1150848);
test3_green_0_1_norm__2->SetBinContent(6,0.1056844);
test3_green_0_1_norm__2->SetBinContent(7,0.09450962);
test3_green_0_1_norm__2->SetBinContent(8,0.07910584);
test3_green_0_1_norm__2->SetBinContent(9,0.06467657);
test3_green_0_1_norm__2->SetBinContent(10,0.04987957);
test3_green_0_1_norm__2->SetEntries(10);
test3_green_0_1_norm__2->SetStats(0);
ci = TColor::GetColor("#00ff00");
test3_green_0_1_norm__2->SetFillColor(ci);
test3_green_0_1_norm__2->SetFillStyle(0);
ci = TColor::GetColor("#00ff00");
test3_green_0_1_norm__2->SetLineColor(ci);
test3_green_0_1_norm__2->SetLineWidth(2);
ci = TColor::GetColor("#00ff00");
test3_green_0_1_norm__2->SetMarkerColor(ci);
test3_green_0_1_norm__2->SetMarkerStyle(20);
test3_green_0_1_norm__2->GetXaxis()->SetTitle("#eta");
test3_green_0_1_norm__2->GetXaxis()->SetNdivisions(1005);
test3_green_0_1_norm__2->GetXaxis()->SetLabelFont(132);
test3_green_0_1_norm__2->GetXaxis()->SetLabelOffset(0.02);
test3_green_0_1_norm__2->GetXaxis()->SetLabelSize(0.05);
test3_green_0_1_norm__2->GetXaxis()->SetTitleSize(0.06);
test3_green_0_1_norm__2->GetXaxis()->SetTitleFont(132);
test3_green_0_1_norm__2->GetYaxis()->SetNdivisions(505);
test3_green_0_1_norm__2->GetYaxis()->SetLabelFont(132);
test3_green_0_1_norm__2->GetYaxis()->SetLabelSize(0.05);
test3_green_0_1_norm__2->GetYaxis()->SetTitleSize(0.06);
test3_green_0_1_norm__2->GetYaxis()->SetTitleFont(132);
test3_green_0_1_norm__2->GetZaxis()->SetLabelFont(132);
test3_green_0_1_norm__2->GetZaxis()->SetLabelSize(0.05);
test3_green_0_1_norm__2->GetZaxis()->SetTitleSize(0.06);
test3_green_0_1_norm__2->GetZaxis()->SetTitleFont(132);
test3_green_0_1_norm__2->Draw("same");
TLegend *leg = new TLegend(0.2,0.7,0.4,0.85,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("MuonTrackingMC2015RW_Sim09/ETA_tot","T&P","lep");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
ci = TColor::GetColor("#ff0000");
entry->SetMarkerColor(ci);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("MuonUbsTrackingMC2015/ETA_tot","Truth","lep");
ci = TColor::GetColor("#00ff00");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
ci = TColor::GetColor("#00ff00");
entry->SetMarkerColor(ci);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TH1F *test3_red_0_0_norm_copy__3 = new TH1F("test3_red_0_0_norm_copy__3","ETA",10,2,4.5);
test3_red_0_0_norm_copy__3->SetBinContent(1,0.1040398);
test3_red_0_0_norm_copy__3->SetBinContent(2,0.1186439);
test3_red_0_0_norm_copy__3->SetBinContent(3,0.1256071);
test3_red_0_0_norm_copy__3->SetBinContent(4,0.1229555);
test3_red_0_0_norm_copy__3->SetBinContent(5,0.118442);
test3_red_0_0_norm_copy__3->SetBinContent(6,0.1127067);
test3_red_0_0_norm_copy__3->SetBinContent(7,0.09951538);
test3_red_0_0_norm_copy__3->SetBinContent(8,0.08332527);
test3_red_0_0_norm_copy__3->SetBinContent(9,0.06651872);
test3_red_0_0_norm_copy__3->SetBinContent(10,0.0482457);
test3_red_0_0_norm_copy__3->SetBinError(1,0.001708584);
test3_red_0_0_norm_copy__3->SetBinError(2,0.001736969);
test3_red_0_0_norm_copy__3->SetBinError(3,0.001484444);
test3_red_0_0_norm_copy__3->SetBinError(4,0.001426268);
test3_red_0_0_norm_copy__3->SetBinError(5,0.001336676);
test3_red_0_0_norm_copy__3->SetBinError(6,0.00122978);
test3_red_0_0_norm_copy__3->SetBinError(7,0.00125216);
test3_red_0_0_norm_copy__3->SetBinError(8,0.00113555);
test3_red_0_0_norm_copy__3->SetBinError(9,0.0009819959);
test3_red_0_0_norm_copy__3->SetBinError(10,0.0009418311);
test3_red_0_0_norm_copy__3->SetEntries(58129);
test3_red_0_0_norm_copy__3->SetDirectory(0);
test3_red_0_0_norm_copy__3->SetStats(0);
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm_copy__3->SetFillColor(ci);
test3_red_0_0_norm_copy__3->SetFillStyle(0);
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm_copy__3->SetLineColor(ci);
test3_red_0_0_norm_copy__3->SetLineWidth(2);
ci = TColor::GetColor("#ff0000");
test3_red_0_0_norm_copy__3->SetMarkerColor(ci);
test3_red_0_0_norm_copy__3->SetMarkerStyle(20);
test3_red_0_0_norm_copy__3->GetXaxis()->SetTitle("#eta");
test3_red_0_0_norm_copy__3->GetXaxis()->SetNdivisions(1005);
test3_red_0_0_norm_copy__3->GetXaxis()->SetLabelFont(132);
test3_red_0_0_norm_copy__3->GetXaxis()->SetLabelOffset(0.02);
test3_red_0_0_norm_copy__3->GetXaxis()->SetLabelSize(0.05);
test3_red_0_0_norm_copy__3->GetXaxis()->SetTitleSize(0.06);
test3_red_0_0_norm_copy__3->GetXaxis()->SetTitleFont(132);
test3_red_0_0_norm_copy__3->GetYaxis()->SetNdivisions(505);
test3_red_0_0_norm_copy__3->GetYaxis()->SetLabelFont(132);
test3_red_0_0_norm_copy__3->GetYaxis()->SetLabelSize(0.05);
test3_red_0_0_norm_copy__3->GetYaxis()->SetTitleSize(0.06);
test3_red_0_0_norm_copy__3->GetYaxis()->SetTitleFont(132);
test3_red_0_0_norm_copy__3->GetZaxis()->SetLabelFont(132);
test3_red_0_0_norm_copy__3->GetZaxis()->SetLabelSize(0.05);
test3_red_0_0_norm_copy__3->GetZaxis()->SetTitleSize(0.06);
test3_red_0_0_norm_copy__3->GetZaxis()->SetTitleFont(132);
test3_red_0_0_norm_copy__3->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
// draws the decay vertices of the Lambda_s
void genDisplay(std::string fullPath, int iGen, int iReco = -1, bool savePdf = false)
{
const int fVerbose(0);
setTDRStyle();
//gStyle->SetOptStat(112211);
gStyle->SetOptStat(0);
gStyle->SetPalette(1);
// Canvases
crz = new TCanvas("crz","crz",1000,600);
TPad* padrz = (TPad*)crz->cd(1);
setPadMargins(padrz, 0.02, 0.05, 0.02, 0.05);
crphi = new TCanvas("crphi","crphi",600,600);
TPad* padrphi = (TPad*)crphi->cd(1);
setPadMargins(padrphi, 0.02, 0.05, 0.02, 0.05);
// Open file
TFile *f = TFile::Open(fullPath.c_str());
if (f==0)
{
cout << "File " << fullPath << " not found -- exiting" << endl;
return;
}
if(fVerbose>0)
cout << "Succesfully opened file " << fullPath << endl;
// Get TTree with GenEvents
TTree* tgen = (TTree*) f->Get("genevents");
if(fVerbose>0) cout << "Got TTree with " << tgen->GetEntries() << " entries" << endl;
// Do a cut, if needed
tgen->Draw(">>lst","ptmu1>3&&ptmu2>3&&TMath::Abs(etamu1)<2.5&&TMath::Abs(etamu2)<2.5");
TEventList *lst;
lst = (TEventList*)gDirectory->Get("lst");
tgen->SetEventList(lst);
// Get TTree with iRecoEvents
TTree* treco = (TTree*) f->Get("events");
if(fVerbose>0) cout << "Got TTree with " << treco->GetEntries() << " entries" << endl;
// Do plots
const int nbins(30);
const double rangeR(30), rangeZ(100); // standard
//const double rangeR(60), rangeZ(200); // zoomed out
//const double rangeR(100), rangeZ(200); // zoomed in
//const double rangeR(12), rangeZ(30); // zoomed further in
padrz->cd();
TH2F *hrz = new TH2F("hdisp","", nbins,-rangeZ,rangeZ,nbins,-rangeR,rangeR);
hrz->Draw();
padrphi->cd();
TH2F *hrphi = new TH2F("hdisp","", nbins,-rangeR,rangeR,nbins,-rangeR,rangeR);
hrphi->Draw();
// add tracker
padrz->Modified();
padrz->Update();
drawTrackerRZ(padrz);
padrphi->Modified();
padrphi->Update();
drawTrackerRPhi(padrphi);
// set branch addresses for gentree
double vrl0,vxl0,vyl0,vzl0,vrlb,vxlb,vylb,vzlb;
double pmu1,phimu1,etamu1,pmu2,phimu2,etamu2;
double ppr,ppi,phipr,phipi,etapr,etapi;
tgen->SetBranchAddress("vrl0",&vrl0);
tgen->SetBranchAddress("vxl0",&vxl0);
tgen->SetBranchAddress("vyl0",&vyl0);
tgen->SetBranchAddress("vzl0",&vzl0);
tgen->SetBranchAddress("vrlb",&vrlb);
tgen->SetBranchAddress("vxlb",&vxlb);
tgen->SetBranchAddress("vylb",&vylb);
tgen->SetBranchAddress("vzlb",&vzlb);
tgen->SetBranchAddress("pmu1",&pmu1);
tgen->SetBranchAddress("etamu1",&etamu1);
tgen->SetBranchAddress("phimu1",&phimu1);
tgen->SetBranchAddress("pmu2",&pmu2);
tgen->SetBranchAddress("etamu2",&etamu2);
tgen->SetBranchAddress("phimu2",&phimu2);
tgen->SetBranchAddress("ppr",&ppr);
tgen->SetBranchAddress("etapr",&etapr);
tgen->SetBranchAddress("phipr",&phipr);
tgen->SetBranchAddress("ppi",&ppi);
tgen->SetBranchAddress("etapi",&etapi);
tgen->SetBranchAddress("phipi",&phipi);
int genrun, genls, genevt;
tgen->SetBranchAddress("run",&genrun);
tgen->SetBranchAddress("LS",&genls);
tgen->SetBranchAddress("event",&genevt);
// set branch addresses for recotree
double rvrl0,rvxl0,rvyl0,rvzl0,rvrlb,rvxlb,rvylb,rvzlb;
double rpmu1,retamu1,rphimu1,rpmu2,retamu2,rphimu2;
double rppr,rppi,retapr,retapi,rphipr,rphipi;
treco->SetBranchAddress("vrl0",&rvrl0);
treco->SetBranchAddress("vxl0",&rvxl0);
treco->SetBranchAddress("vyl0",&rvyl0);
treco->SetBranchAddress("vzl0",&rvzl0);
treco->SetBranchAddress("vrlb",&rvrlb);
treco->SetBranchAddress("vxlb",&rvxlb);
treco->SetBranchAddress("vylb",&rvylb);
treco->SetBranchAddress("vzlb",&rvzlb);
treco->SetBranchAddress("rpt1m",&rpmu1);
treco->SetBranchAddress("reta1m",&retamu1);
treco->SetBranchAddress("rphi1m",&rphimu1);
treco->SetBranchAddress("rpt2m",&rpmu2);
treco->SetBranchAddress("reta2m",&retamu2);
treco->SetBranchAddress("rphi2m",&rphimu2);
treco->SetBranchAddress("rptpr",&rppr);
treco->SetBranchAddress("retapr",&retapr);
treco->SetBranchAddress("rphipr",&rphipr);
treco->SetBranchAddress("rptpi",&rppi);
treco->SetBranchAddress("retapi",&retapi);
treco->SetBranchAddress("rphipi",&rphipi);
drawArrowLegend(padrz);
drawArrowLegend(padrphi,.8);
// write the event indices
writeLatex(padrz,.1,.9,("iGen: " + toString(iGen)).c_str());
writeLatex(padrphi,.1,.9,("iGen: " + toString(iGen)).c_str());
if (iReco >= 0)
{
writeLatex(padrz,.1,.85,("iReco: " + toString(iReco)).c_str());
writeLatex(padrphi,.1,.85,("iReco: " + toString(iReco)).c_str());
}
if (iGen<0) return;
if (iGen>tgen->GetEntries()) return;
if (iReco>treco->GetEntries()) return;
tgen->GetEntry(iGen);
cout << "iGen: " << iGen << " - " << genrun << " " << genls << " " << genevt << endl;
padrz->cd();
drawEventZR(vrlb, vzlb, vrl0, vzl0, pmu1, etamu1, pmu2, etamu2, ppr, etapr, ppi, etapi, 0);
padrphi->cd();
drawEventRPhi(vxlb, vylb, vxl0, vyl0, pmu1, phimu1, pmu2, phimu2, ppr, phipr, ppi, phipi, 0);
if(iReco>=0)
{
treco->GetEntry(iReco);
cout << "iReco: " << iReco << endl;
padrz->cd();
drawEventZR(rvrlb, rvzlb, rvrl0, rvzl0, rpmu1, retamu1, rpmu2, retamu2, rppr, retapr, rppi, retapi, 1);
padrphi->cd();
drawEventRPhi(rvxlb, rvylb, rvxl0, rvyl0, rpmu1, rphimu1, rpmu2, rphimu2, rppr, rphipr, rppi, rphipi, 1);
cout << phimu1 << " " << rphimu1 << endl;
}
if(savePdf)
{
padrz->Update();
padrphi->Update();
std::string filename = "genDisplay_" + toString(iGen);
if(iReco >= 0) filename += "_" + toString(iReco);
crz->SaveAs((filename+"_rz.pdf").c_str());
crphi->SaveAs((filename+"_rphi.png").c_str()); // weiss der Geier warum auf PDF die Kreise und die Skala fehlen...
}
return;
}
void genPlotsRatio0102(std::string fullPath, bool custBinning = false)
{
genPlots01(fullPath, 0, custBinning);
TH2F *h1 = (TH2F*)gDirectory->GetList()->FindObject("hrzL0vtx");
genPlots02(fullPath, 0, custBinning);
TH2F *h2 = (TH2F*)gDirectory->GetList()->FindObject("hrzL0vtxreco");
if (0==h1)
{
cout << "Histo hrzL0vtx not present" << endl;
return;
}
if (0==h2)
{
cout << "Histo hrzL0vtxreco not present" << endl;
return;
}
TH2F *hratio;
if (custBinning)
{
double newbinsX[]={0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50};
const int newbinsX_size = sizeof(newbinsX)/sizeof(double);
double newbinsY[]={0,0.5,1,2,4,8,16,32};
const int newbinsY_size = sizeof(newbinsY)/sizeof(double);
hratio = new TH2F("hratio","hratio",newbinsX_size-1,newbinsX,newbinsY_size-1,newbinsY);
}
else
{
hratio = new TH2F("hratio","hratio",30,0,50,30,0,30);
}
c = new TCanvas("c3","c3",1000,600);
const unsigned int nPadX = 1;
const unsigned int nPadY = 1;
c->Divide(nPadX,nPadY);
const unsigned int nPads=nPadX*nPadY;
for(unsigned int i=1; i<=nPads; i++)
{
TPad* pad= (TPad*)c->cd(i);
pad->SetTopMargin(0.10);
pad->SetRightMargin(0.20);
pad->SetLeftMargin(0.15);
//pad->SetLogy();
}
for(int i=1; i<=h1->GetNbinsX(); i++)
for(int j=1; j<=h1->GetNbinsY(); j++)
{
if (h1->GetBinContent(i,j) != 0.)
hratio->SetBinContent(i,j,h2->GetBinContent(i,j) / h1->GetBinContent(i,j));
else
hratio->SetBinContent(i,j,0);
}
hratio->Draw("COLZ");
// add tracker
TPad* pad;
pad = (TPad*)c->cd(1);
pad->Modified();
pad->Update();
repositionPalette("hratio");
pad->Update();
drawTracker(pad);
}
TCanvas *drawRatioPlot(TH1D *prediction, TH1D *sr, TH1D *data, TString xTitle, TString filename, double ecaloCut){
gStyle -> SetPadLeftMargin(0.20);
data->SetMarkerStyle(20);
data->SetMarkerColor(kGreen);
data->SetLineColor(kGreen);
TCanvas *c = new TCanvas("c"+filename,"c",0,0,500,500);
float y = 0.3;
TPad *pad1 = new TPad("pad1", "Control Plots 1", 0.01, y, 0.99, 0.99);
TPad *padRatio = new TPad("rp1", "Ratio1", 0.01, 0.01, 0.99, y-0.01);
pad1->SetNumber(100);
pad1->SetTicks(0,1);
cout<<"number pad1 = "<<pad1->GetNumber()<<endl;
cout<<"number padRatio = "<<padRatio->GetNumber()<<endl;
TH1D *ratio = 0;
//ratio = (TH1D*) sr->Clone();
//ratio->Divide(prediction);
ratio = (TH1D*) prediction->Clone();
ratio->Divide(data);
for(int i=1; i<=ratio->GetNbinsX();i++){
if(ratio->GetBinContent(i) != 0){
cout<<"N in CR in "<<i<<". bin ="<<prediction->GetBinContent(i)<<endl;
cout<<"N in SR in "<<i<<". bin ="<<sr->GetBinContent(i)<<endl;
cout<<"Rel. difference in "<<i<<". bin ="<<(1./ratio->GetBinContent(i)-1.)*100<<"%"<<endl;
}
else if(sr->GetBinContent(i) == 0 && prediction->GetBinContent(i) !=0) cout<<"Scaling Factor in "<<i<<". bin <"<<prediction->GetBinContent(i)/1.15<<" +/- "<<ratio->GetBinError(i)<<endl;
else if(sr->GetBinContent(i) != 0 && prediction->GetBinContent(i) ==0) cout<<"Scaling Factor in "<<i<<". bin <"<<(sr->GetEntries()/prediction->GetEntries())/sr->GetBinContent(i)<<" +/- "<<ratio->GetBinError(i)<<endl;
}
ratio->GetYaxis()->SetTitle("#frac{CR (MC)}{CR (data)}");
ratio->SetTitle("");
ratio->SetLabelSize(0.1,"X");
ratio->SetLabelSize(0.1,"Y");
ratio->SetTitleOffset(0.5,"Y");
ratio->SetTitleSize(0.15,"Y");
padRatio->cd();
//ratio->GetYaxis()->SetRangeUser(0,2);
ratio->Draw("e");
// Draw line at one!
float xmin = ratio->GetXaxis()->GetXmin();
float xmax = ratio->GetXaxis()->GetXmax();
TLine *line = new TLine(xmin,1,xmax,1);
line->SetLineWidth(2);
line->Draw("same");
padRatio->Modified();
TLegend *leg = new TLegend(0.5,0.7,0.9,0.9);
leg->AddEntry(sr,"SR (MC)","l");
leg->AddEntry(prediction,"lepton CR (MC)","pel");
pad1->cd();
pad1->SetLogy();
// pad1->SetLogx();
sr->SetLineColor(kRed);
sr->SetMarkerColor(kRed);
sr->SetMarkerStyle(20);
sr->SetTitle("");
prediction->SetMarkerStyle(20);
prediction->SetTitle("");
prediction->GetXaxis()->SetTitle(xTitle);
sr->GetXaxis()->SetTitle(xTitle);
prediction->SetTitleSize(0.07,"X");
prediction->GetXaxis()->SetTitleOffset(0.7);
sr->SetTitleSize(0.07,"X");
sr->GetXaxis()->SetTitleOffset(0.7);
double maximum = 0;
double minimum = 1000000;
if(sr->GetMinimum()!=0 && sr->GetMinimum()<minimum){
minimum=sr->GetMinimum()*0.5;
}
if(prediction->GetMinimum()!=0 && prediction->GetMinimum()<minimum){
minimum=prediction->GetMinimum()*0.5;
}
if(data->GetMinimum()!=0 && data->GetMinimum()<minimum){
minimum=data->GetMinimum()*0.5;
}
if(sr->GetMaximum()>maximum){
maximum=sr->GetMaximum()*2.5;
}
if(prediction->GetMaximum()>maximum){
maximum=prediction->GetMaximum()*2.5;
}
if(data->GetMaximum()>maximum){
maximum=data->GetMaximum()*2.5;
}
prediction->GetYaxis()->SetRangeUser(minimum,maximum);
ratio->GetYaxis()->SetRangeUser(0,2);
prediction->Draw("e");
sr->Draw("e same");
leg->AddEntry(data,"lepton CR (data)","pel");
data->Draw("e same");
leg->Draw("same");
TLatex* info1 = new TLatex();
info1->SetTextFont(132);
info1-> SetNDC();
info1->SetTextSize(0.06);
info1->DrawLatex(0.55,0.62,Form("E_{calo}<%.0fGeV",ecaloCut));
// Draw both pads to canvas
c->cd();
pad1->Draw();
padRatio->SetGridy();
padRatio->Draw();
c->SaveAs(filename);
return c;
}
void zj_dphi_8TeV_log()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Thu Apr 28 10:19:34 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,800,800);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: bottomPad
TPad *bottomPad = new TPad("bottomPad", "bottomPad",0.005,0.05,0.995,0.995);
bottomPad->Draw();
bottomPad->cd();
bottomPad->Range(0,0,1,1);
bottomPad->SetFillColor(0);
bottomPad->SetFillStyle(4000);
bottomPad->SetBorderMode(0);
bottomPad->SetBorderSize(2);
bottomPad->SetTickx(1);
bottomPad->SetTicky(1);
bottomPad->SetLeftMargin(0.14);
bottomPad->SetRightMargin(0.05);
bottomPad->SetTopMargin(0.05);
bottomPad->SetBottomMargin(0.16);
bottomPad->SetFrameLineWidth(2);
bottomPad->SetFrameBorderMode(0);
bottomPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(-0.5429913,-2.062865,3.335518,1.327192);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetLogy();
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.33);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
TH1F *z_dphi_toterr__37 = new TH1F("z_dphi_toterr__37","z_dphi",6,0,3.141593);
z_dphi_toterr__37->SetBinContent(1,0.261202);
z_dphi_toterr__37->SetBinContent(2,0.3162337);
z_dphi_toterr__37->SetBinContent(3,0.3852335);
z_dphi_toterr__37->SetBinContent(4,0.7150939);
z_dphi_toterr__37->SetBinContent(5,1.796506);
z_dphi_toterr__37->SetBinContent(6,7.258657);
z_dphi_toterr__37->SetBinContent(7,3.800624);
z_dphi_toterr__37->SetBinError(1,0.03375302);
z_dphi_toterr__37->SetBinError(2,0.03634611);
z_dphi_toterr__37->SetBinError(3,0.0427511);
z_dphi_toterr__37->SetBinError(4,0.05747346);
z_dphi_toterr__37->SetBinError(5,0.1278198);
z_dphi_toterr__37->SetBinError(6,0.3295685);
z_dphi_toterr__37->SetBinError(7,0.1725617);
z_dphi_toterr__37->SetEntries(14);
z_dphi_toterr__37->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetLineColor(ci);
z_dphi_toterr__37->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetMarkerColor(ci);
z_dphi_toterr__37->SetMarkerStyle(0);
z_dphi_toterr__37->GetXaxis()->SetNdivisions(4);
z_dphi_toterr__37->GetXaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetXaxis()->SetLabelOffset(999);
z_dphi_toterr__37->GetXaxis()->SetLabelSize(0);
z_dphi_toterr__37->GetXaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetXaxis()->SetTitleFont(132);
z_dphi_toterr__37->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_toterr__37->GetYaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetYaxis()->SetLabelSize(0.05);
z_dphi_toterr__37->GetYaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetYaxis()->SetTitleOffset(1.3);
z_dphi_toterr__37->GetYaxis()->SetTitleFont(132);
z_dphi_toterr__37->GetZaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetZaxis()->SetLabelSize(0.05);
z_dphi_toterr__37->GetZaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetZaxis()->SetTitleFont(132);
z_dphi_toterr__37->Draw("e2");
TH1F *z_dphi_staterr__38 = new TH1F("z_dphi_staterr__38","z_dphi",6,0,3.141593);
z_dphi_staterr__38->SetBinContent(1,0.261202);
z_dphi_staterr__38->SetBinContent(2,0.3162337);
z_dphi_staterr__38->SetBinContent(3,0.3852335);
z_dphi_staterr__38->SetBinContent(4,0.7150939);
z_dphi_staterr__38->SetBinContent(5,1.796506);
z_dphi_staterr__38->SetBinContent(6,7.258657);
z_dphi_staterr__38->SetBinContent(7,3.800624);
z_dphi_staterr__38->SetBinError(1,0.01764622);
z_dphi_staterr__38->SetBinError(2,0.01951211);
z_dphi_staterr__38->SetBinError(3,0.02108379);
z_dphi_staterr__38->SetBinError(4,0.03021233);
z_dphi_staterr__38->SetBinError(5,0.04899558);
z_dphi_staterr__38->SetBinError(6,0.1002775);
z_dphi_staterr__38->SetBinError(7,0.05250516);
z_dphi_staterr__38->SetEntries(14);
z_dphi_staterr__38->SetStats(0);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetLineColor(ci);
z_dphi_staterr__38->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetMarkerColor(ci);
z_dphi_staterr__38->SetMarkerStyle(0);
z_dphi_staterr__38->GetXaxis()->SetNdivisions(4);
z_dphi_staterr__38->GetXaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetXaxis()->SetLabelOffset(999);
z_dphi_staterr__38->GetXaxis()->SetLabelSize(0);
z_dphi_staterr__38->GetXaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetXaxis()->SetTitleFont(132);
z_dphi_staterr__38->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_staterr__38->GetYaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetYaxis()->SetLabelSize(0.05);
z_dphi_staterr__38->GetYaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetYaxis()->SetTitleFont(132);
z_dphi_staterr__38->GetZaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetZaxis()->SetLabelSize(0.05);
z_dphi_staterr__38->GetZaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetZaxis()->SetTitleFont(132);
z_dphi_staterr__38->Draw("e2same");
Double_t dphi_toterrs_fx3043[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t dphi_toterrs_fy3043[6] = {
0.1926519,
0.2207154,
0.3546338,
0.6864677,
2.01515,
7.779723};
Double_t dphi_toterrs_felx3043[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fely3043[6] = {
0.03415733,
0.03577312,
0.05682396,
0.09598937,
0.1897309,
0.5023379};
Double_t dphi_toterrs_fehx3043[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fehy3043[6] = {
0.04673612,
0.05211068,
0.07287587,
0.1207667,
0.143423,
0.3949306};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(6,dphi_toterrs_fx3043,dphi_toterrs_fy3043,dphi_toterrs_felx3043,dphi_toterrs_fehx3043,dphi_toterrs_fely3043,dphi_toterrs_fehy3043);
grae->SetName("dphi_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_dphi_toterrs3043 = new TH1F("Graph_dphi_toterrs3043","",100,0.07853982,3.220132);
Graph_dphi_toterrs3043->SetMinimum(0.1426451);
Graph_dphi_toterrs3043->SetMaximum(8.976269);
Graph_dphi_toterrs3043->SetDirectory(0);
Graph_dphi_toterrs3043->SetStats(0);
Graph_dphi_toterrs3043->SetLineWidth(2);
Graph_dphi_toterrs3043->SetMarkerStyle(0);
Graph_dphi_toterrs3043->GetXaxis()->SetNdivisions(4);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelOffset(999);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelSize(0);
Graph_dphi_toterrs3043->GetXaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetXaxis()->SetTitleFont(132);
Graph_dphi_toterrs3043->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_dphi_toterrs3043->GetYaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetYaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3043->GetYaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetYaxis()->SetTitleFont(132);
Graph_dphi_toterrs3043->GetZaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetZaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3043->GetZaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_dphi_toterrs3043);
grae->Draw("p1p");
Double_t _fx3044[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t _fy3044[6] = {
0.1926519,
0.2207154,
0.3546338,
0.6864677,
2.01515,
7.779723};
Double_t _felx3044[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fely3044[6] = {
0.02873042,
0.03223617,
0.04948444,
0.07740101,
0.1117696,
0.3251484};
Double_t _fehx3044[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fehy3044[6] = {
0.03948655,
0.04420709,
0.06698785,
0.1028379,
0.1104173,
0.2636435};
grae = new TGraphAsymmErrors(6,_fx3044,_fy3044,_felx3044,_fehx3044,_fely3044,_fehy3044);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph3044 = new TH1F("Graph_Graph3044","",100,0.07853982,3.220132);
Graph_Graph3044->SetMinimum(0.1475293);
Graph_Graph3044->SetMaximum(8.831311);
Graph_Graph3044->SetDirectory(0);
Graph_Graph3044->SetStats(0);
Graph_Graph3044->SetLineWidth(2);
Graph_Graph3044->SetMarkerStyle(0);
Graph_Graph3044->GetXaxis()->SetNdivisions(4);
Graph_Graph3044->GetXaxis()->SetLabelFont(132);
Graph_Graph3044->GetXaxis()->SetLabelOffset(999);
Graph_Graph3044->GetXaxis()->SetLabelSize(0);
Graph_Graph3044->GetXaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetXaxis()->SetTitleFont(132);
Graph_Graph3044->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_Graph3044->GetYaxis()->SetLabelFont(132);
Graph_Graph3044->GetYaxis()->SetLabelSize(0.05);
Graph_Graph3044->GetYaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetYaxis()->SetTitleFont(132);
Graph_Graph3044->GetZaxis()->SetLabelFont(132);
Graph_Graph3044->GetZaxis()->SetLabelSize(0.05);
Graph_Graph3044->GetZaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph3044);
grae->Draw("p1p");
Double_t dphi_toterrs_fx3045[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t dphi_toterrs_fy3045[6] = {
0.3071206,
0.3788296,
0.5235754,
0.9026085,
2.09773,
8.042891};
Double_t dphi_toterrs_felx3045[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fely3045[6] = {
0.03223837,
0.05049376,
0.0513308,
0.07859791,
0.1671821,
0.6048482};
Double_t dphi_toterrs_fehx3045[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fehy3045[6] = {
0.05069997,
0.04670219,
0.05741781,
0.0972528,
0.2140117,
0.6377428};
grae = new TGraphAsymmErrors(6,dphi_toterrs_fx3045,dphi_toterrs_fy3045,dphi_toterrs_felx3045,dphi_toterrs_fehx3045,dphi_toterrs_fely3045,dphi_toterrs_fehy3045);
grae->SetName("dphi_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_dphi_toterrs3045 = new TH1F("Graph_dphi_toterrs3045","",100,0,3.063053);
Graph_dphi_toterrs3045->SetMinimum(0.247394);
Graph_dphi_toterrs3045->SetMaximum(9.521209);
Graph_dphi_toterrs3045->SetDirectory(0);
Graph_dphi_toterrs3045->SetStats(0);
Graph_dphi_toterrs3045->SetLineWidth(2);
Graph_dphi_toterrs3045->SetMarkerStyle(0);
Graph_dphi_toterrs3045->GetXaxis()->SetNdivisions(4);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelOffset(999);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelSize(0);
Graph_dphi_toterrs3045->GetXaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetXaxis()->SetTitleFont(132);
Graph_dphi_toterrs3045->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_dphi_toterrs3045->GetYaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetYaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3045->GetYaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetYaxis()->SetTitleFont(132);
Graph_dphi_toterrs3045->GetZaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetZaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3045->GetZaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_dphi_toterrs3045);
grae->Draw("p1p");
Double_t _fx3046[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t _fy3046[6] = {
0.3071206,
0.3788296,
0.5235754,
0.9026085,
2.09773,
8.042891};
Double_t _felx3046[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fely3046[6] = {
0.02710968,
0.03295216,
0.04124188,
0.06402796,
0.1338949,
0.4694553};
Double_t _fehx3046[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fehy3046[6] = {
0.03469016,
0.04020747,
0.048986,
0.07362466,
0.1488417,
0.5027875};
grae = new TGraphAsymmErrors(6,_fx3046,_fy3046,_felx3046,_fehx3046,_fely3046,_fehy3046);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3046 = new TH1F("Graph_Graph3046","",100,0,3.063053);
Graph_Graph3046->SetMinimum(0.2520099);
Graph_Graph3046->SetMaximum(9.372245);
Graph_Graph3046->SetDirectory(0);
Graph_Graph3046->SetStats(0);
Graph_Graph3046->SetLineWidth(2);
Graph_Graph3046->SetMarkerStyle(0);
Graph_Graph3046->GetXaxis()->SetNdivisions(4);
Graph_Graph3046->GetXaxis()->SetLabelFont(132);
Graph_Graph3046->GetXaxis()->SetLabelOffset(999);
Graph_Graph3046->GetXaxis()->SetLabelSize(0);
Graph_Graph3046->GetXaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetXaxis()->SetTitleFont(132);
Graph_Graph3046->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_Graph3046->GetYaxis()->SetLabelFont(132);
Graph_Graph3046->GetYaxis()->SetLabelSize(0.05);
Graph_Graph3046->GetYaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetYaxis()->SetTitleFont(132);
Graph_Graph3046->GetZaxis()->SetLabelFont(132);
Graph_Graph3046->GetZaxis()->SetLabelSize(0.05);
Graph_Graph3046->GetZaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph3046);
grae->Draw("p1p");
TLegend *leg = new TLegend(0.2,0.57,0.5,0.82,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("data_stat","Data_{stat}","f");
ci = TColor::GetColor("#ff9933");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("data_tot","Data_{tot}","f");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("powheg_temp","POWHEG","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("amcatnlo_temp","aMC@NLO","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(24);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.2,0.8,0.55,0.95,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb, #sqrt{s} = 8 TeV");
pt->Draw();
TH1F *z_dphi_toterr_copy__39 = new TH1F("z_dphi_toterr_copy__39","z_dphi",6,0,3.141593);
z_dphi_toterr_copy__39->SetBinContent(1,0.261202);
z_dphi_toterr_copy__39->SetBinContent(2,0.3162337);
z_dphi_toterr_copy__39->SetBinContent(3,0.3852335);
z_dphi_toterr_copy__39->SetBinContent(4,0.7150939);
z_dphi_toterr_copy__39->SetBinContent(5,1.796506);
z_dphi_toterr_copy__39->SetBinContent(6,7.258657);
z_dphi_toterr_copy__39->SetBinContent(7,3.800624);
z_dphi_toterr_copy__39->SetBinError(1,0.03375302);
z_dphi_toterr_copy__39->SetBinError(2,0.03634611);
z_dphi_toterr_copy__39->SetBinError(3,0.0427511);
z_dphi_toterr_copy__39->SetBinError(4,0.05747346);
z_dphi_toterr_copy__39->SetBinError(5,0.1278198);
z_dphi_toterr_copy__39->SetBinError(6,0.3295685);
z_dphi_toterr_copy__39->SetBinError(7,0.1725617);
z_dphi_toterr_copy__39->SetEntries(14);
z_dphi_toterr_copy__39->SetDirectory(0);
z_dphi_toterr_copy__39->SetStats(0);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetLineColor(ci);
z_dphi_toterr_copy__39->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetMarkerColor(ci);
z_dphi_toterr_copy__39->SetMarkerStyle(0);
z_dphi_toterr_copy__39->GetXaxis()->SetNdivisions(4);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelOffset(999);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelSize(0);
z_dphi_toterr_copy__39->GetXaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetXaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_toterr_copy__39->GetYaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetYaxis()->SetLabelSize(0.05);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleOffset(1.3);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->GetZaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetZaxis()->SetLabelSize(0.05);
z_dphi_toterr_copy__39->GetZaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetZaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: lowerPad
TPad *lowerPad = new TPad("lowerPad", "lowerPad",0.005,0.05,0.995,0.995);
lowerPad->Draw();
lowerPad->cd();
lowerPad->Range(-0.5429913,-0.4371795,3.335518,5.614103);
lowerPad->SetFillColor(0);
lowerPad->SetFillStyle(4000);
lowerPad->SetBorderMode(0);
lowerPad->SetBorderSize(2);
lowerPad->SetTickx(1);
lowerPad->SetTicky(1);
lowerPad->SetLeftMargin(0.14);
lowerPad->SetRightMargin(0.05);
lowerPad->SetTopMargin(0.665);
lowerPad->SetBottomMargin(0.14);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp0__40 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp0__40","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(1,0.1292219);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(2,0.1149343);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(3,0.1109745);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(4,0.08037191);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(5,0.07114914);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(6,0.04540351);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(7,0.1725617);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleOffset(1.3);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->Draw("e2");
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp1__41 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp1__41","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(1,0.06755777);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(2,0.06170157);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(3,0.0547299);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(4,0.04224946);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(5,0.02727271);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(6,0.01381488);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(7,0.05250516);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetStats(0);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->Draw("e2same");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fx3047[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fy3047[6] = {
0.7375589,
0.6979502,
0.9205685,
0.9599686,
1.121705,
1.071786};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_felx3047[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fely3047[6] = {
0.1307698,
0.1131224,
0.1475052,
0.1342332,
0.1056111,
0.06920535};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehx3047[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehy3047[6] = {
0.1789271,
0.1647854,
0.1891732,
0.1688822,
0.07983441,
0.05440822};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fy3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_felx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fely3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehy3047);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp2");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047","",100,0.07853982,3.220132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047);
grae->Draw("p1");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fx3048[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fy3048[6] = {
0.7375589,
0.6979502,
0.9205685,
0.9599686,
1.121705,
1.071786};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_felx3048[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fely3048[6] = {
0.1099931,
0.1019378,
0.1284531,
0.1082389,
0.06221498,
0.04479457};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehx3048[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehy3048[6] = {
0.1511724,
0.1397925,
0.1738889,
0.1438103,
0.06146226,
0.03632125};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fy3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_felx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fely3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehy3048);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp3");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048","",100,0.07853982,3.220132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048);
grae->Draw("p1");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fx3049[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fy3049[6] = {
1.175797,
1.197942,
1.359112,
1.262224,
1.167673,
1.108041};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_felx3049[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fely3049[6] = {
0.1234231,
0.1596723,
0.1332459,
0.1099127,
0.09305962,
0.08332785};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehx3049[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehy3049[6] = {
0.1941025,
0.1476825,
0.1490468,
0.136,
0.1191266,
0.08785962};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fy3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_felx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fely3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehy3049);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp4");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049","",100,0,3.063053);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049);
grae->Draw("p1");
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(1,0.1292219);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(2,0.1149343);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(3,0.1109745);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(4,0.08037191);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(5,0.07114914);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(6,0.04540351);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(7,0.1725617);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetDirectory(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleOffset(1.3);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->Draw("sameaxis");
lowerPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}