{
//=========Macro generated from canvas: c1/Cos_theta dists for reconstructed top1
//=========  (Sat Jun 13 00:25:32 2009) by ROOT version5.21/02
   TCanvas *c1 = new TCanvas("c1", "Cos_theta dists for reconstructed top1",0,44,900,600);
   gStyle->SetOptStat(0);
   gStyle->SetOptTitle(0);

   c1->Range(-1.453333,-933.3904,1.213333,4252.112);
   c1->SetFillColor(10);
   c1->SetBorderMode(0);
   c1->SetBorderSize(2);
   c1->SetTickx(1);
   c1->SetTicky(1);
   c1->SetLeftMargin(0.17);
   c1->SetRightMargin(0.08);
   c1->SetTopMargin(0.08);
   c1->SetBottomMargin(0.18);
   c1->SetFrameLineWidth(2);
   c1->SetFrameBorderMode(0);
   c1->SetFrameLineWidth(2);
   c1->SetFrameBorderMode(0);
   
   TH1 *hCosT1Obs = new TH1D("hCosT1Obs","d#sigma_{obs} (#theta_{t1})  w/o charge-ID corr.",20,-1,1);
   hCosT1Obs->SetBinContent(1,1588.903);
   hCosT1Obs->SetBinContent(2,1474.753);
   hCosT1Obs->SetBinContent(3,1425.544);
   hCosT1Obs->SetBinContent(4,1345.152);
   hCosT1Obs->SetBinContent(5,1295.997);
   hCosT1Obs->SetBinContent(6,1324.99);
   hCosT1Obs->SetBinContent(7,1363.399);
   hCosT1Obs->SetBinContent(8,1381.063);
   hCosT1Obs->SetBinContent(9,1386.779);
   hCosT1Obs->SetBinContent(10,1444.67);
   hCosT1Obs->SetBinContent(11,1499.458);
   hCosT1Obs->SetBinContent(12,1609.798);
   hCosT1Obs->SetBinContent(13,1668.936);
   hCosT1Obs->SetBinContent(14,1848.95);
   hCosT1Obs->SetBinContent(15,2038.37);
   hCosT1Obs->SetBinContent(16,2154.986);
   hCosT1Obs->SetBinContent(17,2413.237);
   hCosT1Obs->SetBinContent(18,2654.825);
   hCosT1Obs->SetBinContent(19,2818.187);
   hCosT1Obs->SetBinContent(20,3141.676);
   hCosT1Obs->SetEntries(144737);
   hCosT1Obs->SetStats(0);
   hCosT1Obs->SetLineWidth(2);
   hCosT1Obs->SetMarkerStyle(20);
   hCosT1Obs->SetMarkerSize(0.7);
   hCosT1Obs->GetXaxis()->SetTitle("cos #theta_{t}");
   hCosT1Obs->GetXaxis()->SetNdivisions(506);
   hCosT1Obs->GetXaxis()->SetLabelFont(42);
   hCosT1Obs->GetXaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetXaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetXaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetXaxis()->SetTitleFont(42);
   hCosT1Obs->GetYaxis()->SetNdivisions(506);
   hCosT1Obs->GetYaxis()->SetLabelFont(42);
   hCosT1Obs->GetYaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetYaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetYaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetYaxis()->SetTitleOffset(1.1);
   hCosT1Obs->GetYaxis()->SetTitleFont(42);
   hCosT1Obs->GetYaxis()->SetTitle("dN/dcos#theta_{t}");
   hCosT1Obs->GetZaxis()->SetLabelFont(42);
   hCosT1Obs->GetZaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetZaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetZaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetZaxis()->SetTitleOffset(1.1);
   hCosT1Obs->GetZaxis()->SetTitleFont(42);
   hCosT1Obs->SetMaximum(4000);
   hCosT1Obs->SetMinimum(0);

   hCosT1Obs->Draw("E2");
   
   TH1 *hCosT1Obs = new TH1D("hCosT1Obs","d#sigma_{obs} (#theta_{t1})  w/o charge-ID corr.",20,-1,1);
   hCosT1Obs->SetBinContent(1,1588.903);
   hCosT1Obs->SetBinContent(2,1474.753);
   hCosT1Obs->SetBinContent(3,1425.544);
   hCosT1Obs->SetBinContent(4,1345.152);
   hCosT1Obs->SetBinContent(5,1295.997);
   hCosT1Obs->SetBinContent(6,1324.99);
   hCosT1Obs->SetBinContent(7,1363.399);
   hCosT1Obs->SetBinContent(8,1381.063);
   hCosT1Obs->SetBinContent(9,1386.779);
   hCosT1Obs->SetBinContent(10,1444.67);
   hCosT1Obs->SetBinContent(11,1499.458);
   hCosT1Obs->SetBinContent(12,1609.798);
   hCosT1Obs->SetBinContent(13,1668.936);
   hCosT1Obs->SetBinContent(14,1848.95);
   hCosT1Obs->SetBinContent(15,2038.37);
   hCosT1Obs->SetBinContent(16,2154.986);
   hCosT1Obs->SetBinContent(17,2413.237);
   hCosT1Obs->SetBinContent(18,2654.825);
   hCosT1Obs->SetBinContent(19,2818.187);
   hCosT1Obs->SetBinContent(20,3141.676);
   hCosT1Obs->SetEntries(144737);
   hCosT1Obs->SetStats(0);
   hCosT1Obs->SetLineWidth(2);
   hCosT1Obs->SetMarkerStyle(20);
   hCosT1Obs->SetMarkerSize(0.7);
   hCosT1Obs->GetXaxis()->SetTitle("cost #theta_{t}");
   hCosT1Obs->GetXaxis()->SetNdivisions(506);
   hCosT1Obs->GetXaxis()->SetLabelFont(42);
   hCosT1Obs->GetXaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetXaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetXaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetXaxis()->SetTitleFont(42);
   hCosT1Obs->GetYaxis()->SetNdivisions(506);
   hCosT1Obs->GetYaxis()->SetLabelFont(42);
   hCosT1Obs->GetYaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetYaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetYaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetYaxis()->SetTitleOffset(1.1);
   hCosT1Obs->GetYaxis()->SetTitleFont(42);
   hCosT1Obs->GetZaxis()->SetLabelFont(42);
   hCosT1Obs->GetZaxis()->SetLabelOffset(0.015);
   hCosT1Obs->GetZaxis()->SetLabelSize(0.06);
   hCosT1Obs->GetZaxis()->SetTitleSize(0.07);
   hCosT1Obs->GetZaxis()->SetTitleOffset(1.1);
   hCosT1Obs->GetZaxis()->SetTitleFont(42);
   hCosT1Obs->Draw("sameE");
   
   TH1 *hCosT1ObsMisChT = new TH1D("hCosT1ObsMisChT","d#sigma_{obs} (#theta_{t1})  w/o charge-ID corr.",20,-1,1);
   hCosT1ObsMisChT->SetBinContent(1,760.4787);
   hCosT1ObsMisChT->SetBinContent(2,652.2153);
   hCosT1ObsMisChT->SetBinContent(3,586.7313);
   hCosT1ObsMisChT->SetBinContent(4,540.5063);
   hCosT1ObsMisChT->SetBinContent(5,462.0874);
   hCosT1ObsMisChT->SetBinContent(6,446.7458);
   hCosT1ObsMisChT->SetBinContent(7,417.9211);
   hCosT1ObsMisChT->SetBinContent(8,390.2318);
   hCosT1ObsMisChT->SetBinContent(9,353.5446);
   hCosT1ObsMisChT->SetBinContent(10,321.8037);
   hCosT1ObsMisChT->SetBinContent(11,299.9149);
   hCosT1ObsMisChT->SetBinContent(12,287.3448);
   hCosT1ObsMisChT->SetBinContent(13,262.3811);
   hCosT1ObsMisChT->SetBinContent(14,269.1206);
   hCosT1ObsMisChT->SetBinContent(15,258.9037);
   hCosT1ObsMisChT->SetBinContent(16,237.2905);
   hCosT1ObsMisChT->SetBinContent(17,270.7557);
   hCosT1ObsMisChT->SetBinContent(18,263.9929);
   hCosT1ObsMisChT->SetBinContent(19,256.0966);
   hCosT1ObsMisChT->SetBinContent(20,271.5411);
   hCosT1ObsMisChT->SetEntries(30562);
   hCosT1ObsMisChT->SetStats(0);
   hCosT1ObsMisChT->SetFillColor(2);
   hCosT1ObsMisChT->SetFillStyle(3001);
   hCosT1ObsMisChT->SetLineColor(2);
   hCosT1ObsMisChT->SetLineWidth(2);
   hCosT1ObsMisChT->SetMarkerStyle(20);
   hCosT1ObsMisChT->SetMarkerSize(0.7);
   hCosT1ObsMisChT->GetXaxis()->SetNdivisions(506);
   hCosT1ObsMisChT->GetXaxis()->SetLabelFont(42);
   hCosT1ObsMisChT->GetXaxis()->SetLabelOffset(0.015);
   hCosT1ObsMisChT->GetXaxis()->SetLabelSize(0.06);
   hCosT1ObsMisChT->GetXaxis()->SetTitleSize(0.07);
   hCosT1ObsMisChT->GetXaxis()->SetTitleFont(42);
   hCosT1ObsMisChT->GetYaxis()->SetNdivisions(506);
   hCosT1ObsMisChT->GetYaxis()->SetLabelFont(42);
   hCosT1ObsMisChT->GetYaxis()->SetLabelOffset(0.015);
   hCosT1ObsMisChT->GetYaxis()->SetLabelSize(0.06);
   hCosT1ObsMisChT->GetYaxis()->SetTitleSize(0.07);
   hCosT1ObsMisChT->GetYaxis()->SetTitleOffset(1.1);
   hCosT1ObsMisChT->GetYaxis()->SetTitleFont(42);
   hCosT1ObsMisChT->GetZaxis()->SetLabelFont(42);
   hCosT1ObsMisChT->GetZaxis()->SetLabelOffset(0.015);
   hCosT1ObsMisChT->GetZaxis()->SetLabelSize(0.06);
   hCosT1ObsMisChT->GetZaxis()->SetTitleSize(0.07);
   hCosT1ObsMisChT->GetZaxis()->SetTitleOffset(1.1);
   hCosT1ObsMisChT->GetZaxis()->SetTitleFont(42);
   hCosT1ObsMisChT->Draw("][same");
   
   TH1 *hCosT1ObsMissBT = new TH1D("hCosT1ObsMissBT","d#sigma_{obs} (#theta_{t1})  w/o charge-ID corr.",20,-1,1);
   hCosT1ObsMissBT->SetBinContent(1,188.0752);
   hCosT1ObsMissBT->SetBinContent(2,170.3413);
   hCosT1ObsMissBT->SetBinContent(3,150.3884);
   hCosT1ObsMissBT->SetBinContent(4,138.4142);
   hCosT1ObsMissBT->SetBinContent(5,131.964);
   hCosT1ObsMissBT->SetBinContent(6,123.7492);
   hCosT1ObsMissBT->SetBinContent(7,124.2344);
   hCosT1ObsMissBT->SetBinContent(8,119.6621);
   hCosT1ObsMissBT->SetBinContent(9,109.9043);
   hCosT1ObsMissBT->SetBinContent(10,112.2594);
   hCosT1ObsMissBT->SetBinContent(11,108.9144);
   hCosT1ObsMissBT->SetBinContent(12,118.0097);
   hCosT1ObsMissBT->SetBinContent(13,109.4013);
   hCosT1ObsMissBT->SetBinContent(14,119.0831);
   hCosT1ObsMissBT->SetBinContent(15,134.7174);
   hCosT1ObsMissBT->SetBinContent(16,132.803);
   hCosT1ObsMissBT->SetBinContent(17,146.1582);
   hCosT1ObsMissBT->SetBinContent(18,157.7385);
   hCosT1ObsMissBT->SetBinContent(19,170.7357);
   hCosT1ObsMissBT->SetBinContent(20,194.2873);
   hCosT1ObsMissBT->SetEntries(10960);
   hCosT1ObsMissBT->SetStats(0);
   hCosT1ObsMissBT->SetFillColor(4);
   hCosT1ObsMissBT->SetFillStyle(3001);
   hCosT1ObsMissBT->SetLineColor(4);
   hCosT1ObsMissBT->SetLineWidth(2);
   hCosT1ObsMissBT->SetMarkerStyle(20);
   hCosT1ObsMissBT->SetMarkerSize(0.7);
   hCosT1ObsMissBT->GetXaxis()->SetNdivisions(506);
   hCosT1ObsMissBT->GetXaxis()->SetLabelFont(42);
   hCosT1ObsMissBT->GetXaxis()->SetLabelOffset(0.015);
   hCosT1ObsMissBT->GetXaxis()->SetLabelSize(0.06);
   hCosT1ObsMissBT->GetXaxis()->SetTitleSize(0.07);
   hCosT1ObsMissBT->GetXaxis()->SetTitleFont(42);
   hCosT1ObsMissBT->GetYaxis()->SetNdivisions(506);
   hCosT1ObsMissBT->GetYaxis()->SetLabelFont(42);
   hCosT1ObsMissBT->GetYaxis()->SetLabelOffset(0.015);
   hCosT1ObsMissBT->GetYaxis()->SetLabelSize(0.06);
   hCosT1ObsMissBT->GetYaxis()->SetTitleSize(0.07);
   hCosT1ObsMissBT->GetYaxis()->SetTitleOffset(1.1);
   hCosT1ObsMissBT->GetYaxis()->SetTitleFont(42);
   hCosT1ObsMissBT->GetZaxis()->SetLabelFont(42);
   hCosT1ObsMissBT->GetZaxis()->SetLabelOffset(0.015);
   hCosT1ObsMissBT->GetZaxis()->SetLabelSize(0.06);
   hCosT1ObsMissBT->GetZaxis()->SetTitleSize(0.07);
   hCosT1ObsMissBT->GetZaxis()->SetTitleOffset(1.1);
   hCosT1ObsMissBT->GetZaxis()->SetTitleFont(42);
   hCosT1ObsMissBT->Draw("][same");
   
   TLegend *leg = new TLegend(0.25,0.65,0.7,0.875,NULL,"brNDC");
   leg->SetBorderSize(0);
   leg->SetLineColor(1);
   leg->SetLineStyle(1);
   leg->SetLineWidth(2);
   leg->SetTextSize(0.05);
   leg->SetFillColor(10);
   leg->SetFillStyle(1001);
   TLegendEntry *entry=leg->AddEntry("hCosT1Obs","Total 500fb^{-1}","lp");
   entry->SetFillStyle(1001);
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(20);
   entry->SetMarkerSize(0.7);
   entry=leg->AddEntry("hCosT1ObsMisChT","Wrong charge","f");
   entry->SetFillColor(2);
   entry->SetFillStyle(3001);
   entry->SetLineColor(2);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(20);
   entry->SetMarkerSize(0.7);
   entry=leg->AddEntry("hCosT1ObsMissBT","Combinatorial Background","f");
   entry->SetFillColor(4);
   entry->SetFillStyle(3001);
   entry->SetLineColor(4);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(20);
   entry->SetMarkerSize(0.7);
   leg->Draw();
   c1->Modified();
   c1->cd();
   c1->SetSelected(c1);
   c1->SaveAs("cosThetaTop.eps");
   c1->SaveAs("cosThetaTop.gif");

}