void getAiZ()
{
TFile *data1 = new TFile("results/CSFrame_ZJets_madgraph_0_inf_plus.root");
TFile *data2 = new TFile("results/CSFrame_ZJets_madgraph_0_inf_minus.root");
TString HistPA0 = "MC_CSFramePlots_0toinf/MC_CSWPlusA0Numerator";
TString HistPA1 = "MC_CSFramePlots_0toinf/MC_CSWPlusA1Numerator";
TString HistPA2 = "MC_CSFramePlots_0toinf/MC_CSWPlusA2Numerator";
TString HistPA3 = "MC_CSFramePlots_0toinf/MC_CSWPlusA3Numerator";
TString HistPA4 = "MC_CSFramePlots_0toinf/MC_CSWPlusA4Numerator";
TString HistPD = "MC_CSFramePlots_0toinf/MC_CSWPlusAiDenominator";
TString HistMA0 = "MC_CSFramePlots_0toinf/MC_CSWMinusA0Numerator";
TString HistMA1 = "MC_CSFramePlots_0toinf/MC_CSWMinusA1Numerator";
TString HistMA2 = "MC_CSFramePlots_0toinf/MC_CSWMinusA2Numerator";
TString HistMA3 = "MC_CSFramePlots_0toinf/MC_CSWMinusA3Numerator";
TString HistMA4 = "MC_CSFramePlots_0toinf/MC_CSWMinusA4Numerator";
TString HistMD = "MC_CSFramePlots_0toinf/MC_CSWMinusAiDenominator";
TH1D *hPA0 = (TH1D*)data1->Get(HistPA0);
TH1D *hPA1 = (TH1D*)data1->Get(HistPA1);
TH1D *hPA2 = (TH1D*)data1->Get(HistPA2);
TH1D *hPA3 = (TH1D*)data1->Get(HistPA3);
TH1D *hPA4 = (TH1D*)data1->Get(HistPA4);
TH1D *hPD = (TH1D*)data1->Get(HistPD);
TH1D *hMA0 = (TH1D*)data2->Get(HistMA0);
TH1D *hMA1 = (TH1D*)data2->Get(HistMA1);
TH1D *hMA2 = (TH1D*)data2->Get(HistMA2);
TH1D *hMA3 = (TH1D*)data2->Get(HistMA3);
TH1D *hMA4 = (TH1D*)data2->Get(HistMA4);
TH1D *hMD = (TH1D*)data2->Get(HistMD);
hPA0->Rebin(20);
hPA1->Rebin(20);
hPA2->Rebin(20);
hPA3->Rebin(20);
hPA4->Rebin(20);
hPD->Rebin(20);
hMA0->Rebin(20);
hMA1->Rebin(20);
hMA2->Rebin(20);
hMA3->Rebin(20);
hMA4->Rebin(20);
hMD->Rebin(20);
TCanvas * aCanvas = new TCanvas("CSAiZ", "CSAiZ", 1100, 500);
aCanvas->SetFillColor(0);
aCanvas->SetFrameFillColor(0);
aCanvas->SetFrameBorderMode(0);
aCanvas->Divide(2,1);
setPadMargins(aCanvas->cd(1));
TLegend * aLegend = new TLegend(0.23,0.92,0.4,0.62);
aLegend->SetFillColor(0);
aLegend->SetLineColor(0);
aLegend->SetTextSize(0.0625);
aLegend->AddEntry(hPA0,"A0", "L");
//aLegend->AddEntry(hPA1,"A1", "L");
aLegend->AddEntry(hPA2,"A2", "L");
aLegend->AddEntry(hPA3,"A3", "L");
//aLegend->AddEntry(hPA4,"A4", "L");
hPA0->Scale(20.0/3.0);
hPA0->Divide(hPD);
for(unsigned int i=0; i<hPA0->GetNbinsX(); i++) {
hPA0->SetBinContent(i, (2.0/3.0) + hPA0->GetBinContent(i));
}
hPA0->SetLineWidth(3);
hPA0->GetYaxis()->SetRangeUser(-1.,2.);
hPA0->GetXaxis()->SetTitle("P_{T}(Z) [GeV]");
hPA0->GetXaxis()->SetRangeUser(0,220);
hPA0->GetYaxis()->SetTitle("A_{i}");
hPA0->GetXaxis()->SetTitleSize(0.0625);
hPA0->GetYaxis()->SetTitleSize(0.0625);
hPA0->GetXaxis()->SetTitleOffset(1.00);
hPA0->GetYaxis()->SetTitleOffset(1.20);
hPA0->GetXaxis()->SetLabelSize(0.052);
hPA0->GetYaxis()->SetLabelSize(0.052);
hPA0->GetXaxis()->SetLabelOffset(0.01);
hPA0->GetYaxis()->SetLabelOffset(0.01);
hPA0->DrawCopy("hist");
hPA1->Scale(5.0);
hPA1->Divide(hPD);
hPA1->SetLineColor(kRed);
hPA1->SetLineWidth(3);
//hPA1->DrawCopy("histsame");//we don't have sensitivity to A1 - sign of cos(theta_CS)
hPA2->Scale(10.0);
hPA2->Divide(hPD);
hPA2->SetLineColor(kBlue);
hPA2->SetLineWidth(3);
hPA2->DrawCopy("histsame");
hPA3->Scale(4.0);
hPA3->Divide(hPD);
hPA3->SetLineColor(kGreen);
hPA3->SetLineWidth(3);
hPA3->DrawCopy("histsame");
hPA4->Scale(4.0);
hPA4->Divide(hPD);
hPA4->SetLineColor(kCyan);
hPA4->SetLineWidth(3);
//hPA4->DrawCopy("histsame");//we don't have sensitivity to A4 - sign of cos(theta_CS)
TLatex *tplus = labelLatex(0.25,0.25,"Lepton Plus Charge");
tplus->DrawClone("same");
aLegend->DrawClone("same");
setPadMargins(aCanvas->cd(2));
hMA0->Scale(20.0/3.0);
hMA0->Divide(hMD);
for(unsigned int i=0; i<hMA0->GetNbinsX(); i++) {
hMA0->SetBinContent(i, (2.0/3.0) + hMA0->GetBinContent(i));
}
hMA0->SetLineWidth(3);
hMA0->GetYaxis()->SetRangeUser(-1.,2.);
hMA0->GetXaxis()->SetTitle("P_{T}(Z) [GeV]");
hMA0->GetXaxis()->SetRangeUser(0,220);
hMA0->GetYaxis()->SetTitle("A_{i}");
hMA0->GetXaxis()->SetTitleSize(0.0625);
hMA0->GetYaxis()->SetTitleSize(0.0625);
hMA0->GetXaxis()->SetTitleOffset(1.00);
hMA0->GetYaxis()->SetTitleOffset(1.20);
hMA0->GetXaxis()->SetLabelSize(0.052);
hMA0->GetYaxis()->SetLabelSize(0.052);
hMA0->GetXaxis()->SetLabelOffset(0.01);
hMA0->GetYaxis()->SetLabelOffset(0.01);
hMA0->DrawCopy("hist");
hMA1->Scale(5.0);
hMA1->Divide(hMD);
hMA1->SetLineColor(kRed);
hMA1->SetLineWidth(3);
//hMA1->DrawCopy("histsame");
hMA2->Scale(10.0);
hMA2->Divide(hMD);
hMA2->SetLineColor(kBlue);
hMA2->SetLineWidth(3);
hMA2->DrawCopy("histsame");
hMA3->Scale(4.0);
hMA3->Divide(hMD);
hMA3->SetLineColor(kGreen);
hMA3->SetLineWidth(3);
hMA3->DrawCopy("histsame");
hMA4->Scale(4.0);
hMA4->Divide(hMD);
hMA4->SetLineColor(kCyan);
hMA4->SetLineWidth(3);
//hMA4->DrawCopy("histsame");
TLatex *tminus = labelLatex(0.25,0.25,"Lepton Minus Charge");
tminus->DrawClone("same");
aLegend->DrawClone("same");
aCanvas->Print(".png");
return;
}
void MCtoSTATerrors() {
TLatex *tplus = labelLatex(0.1364943,0.8114407,"MC sample");
TFile *file0 = TFile::Open("results/StandardCutflow/allmc.root");
TCanvas *c1 = new TCanvas();
//TCanvas *c2 = new TCanvas();
TH1D * LPplus = (TH1D*)file0->Get("RECO_PolPlots_50toinf/RECO_ICVarPFPlus");
TH1D * LPminus = (TH1D*)file0->Get("RECO_PolPlots_50toinf/RECO_ICVarPFMinus");
LPplus->Rebin(20);
LPminus->Rebin(20);
LPplus->Scale(0.3);
LPminus->Scale(0.3);
TH1D * newLPplus = (TH1D*)LPplus->Clone();
TH1D * newLPminus = (TH1D*)LPminus->Clone();
for(unsigned int i=1; i < LPplus->GetXaxis()->GetNbins(); i++) {
newLPplus->SetBinError(i, TMath::Sqrt(LPplus->GetBinContent(i)));
newLPminus->SetBinError(i, TMath::Sqrt(LPminus->GetBinContent(i)));
}
c1->cd();
newLPplus->GetYaxis()->SetTitle("Events / 300 nb^{-1}");
//newLPplus->GetXaxis()->SetTitle("LP(#mu^{+})");
newLPplus->GetXaxis()->SetTitle("LP(#mu)");
newLPplus->GetXaxis()->SetRangeUser(-1.0,2.0);
newLPplus->SetLineWidth(3);
newLPplus->SetLineStyle(2);
newLPplus->SetLineColor(kRed);
newLPplus->DrawCopy();
tplus->DrawClone("same");
// c2->cd();
// newLPminus->GetYaxis()->SetTitle("Events / 100 nb^{-1}");
// newLPminus->GetXaxis()->SetTitle("LP(#mu^{-})");
newLPminus->SetLineWidth(3);
newLPminus->SetLineStyle(2);
newLPminus->SetLineColor(kBlack);
newLPminus->DrawCopy("same");
// tplus->DrawClone("same");
leg = new TLegend(0.762931,0.720339,0.8649425,0.8622881,NULL,"brNDC");
leg->AddEntry(newLPplus,"#mu^{+}","l");
leg->AddEntry(newLPminus,"#mu^{-}","l");
leg->SetFillColor(kWhite);
leg->SetBorderSize(0);
leg->SetTextFont(62);
leg->DrawClone();
file0->Close();
return;
}
