plot_jpsi_multiplicity_Haiwang(bool draw_alice = false)
{
	//=========Macro generated from canvas: c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc/c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc
	//=========  (Fri May  2 11:31:57 2014) by ROOT version5.30/03

	TCanvas *c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc = new TCanvas("c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc", "c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc",-292,49,1000,674);
	gStyle->SetOptStat(0);
	gStyle->SetOptTitle(0);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->Range(-0.6856555,-1.293858,4.800745,6.5);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetFillColor(0);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetBorderMode(0);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetBorderSize(2);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetLeftMargin(0.1375502);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetRightMargin(0.02911647);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetTopMargin(0.04383117);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetBottomMargin(0.1720779);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetFrameBorderMode(0);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetFrameBorderMode(0);

	TGraphErrors *gre = new TGraphErrors(5);
	gre->SetName("g_ALICE");
	gre->SetTitle("J/#psi Yield vs. FVTX multiplicity (Normalized to MB);#frac{N_{ch}}{<N_{ch}^{MB}>};#frac{N_{J/#psi}/N_{MB}}{<N_{J/#psi}/N_{MB}>}");
	gre->SetFillColor(1);

	gre->SetLineWidth(3);

	gre->SetMarkerStyle(21);
	gre->SetMarkerColor(kGray+1);
	gre->SetLineColor(kGray+1);
	gre->SetMarkerSize(2);
	gre->SetPoint(0,0.48,0.481);
	gre->SetPointError(0,0.0,0.037);
	gre->SetPoint(1,1.23,1.07);
	gre->SetPointError(1,0.0,0.1);
	gre->SetPoint(2,1.815173,1.61122);
	gre->SetPointError(2,0.0,0.15);
	gre->SetPoint(3,2.66,3.26);
	gre->SetPointError(3,0.0,0.25);
	gre->SetPoint(4,4.02,4.83);
	gre->SetPointError(4,0.0,0.57);

	TGraphErrors *gre_sys = (TGraphErrors*)gre->Clone("gre_sys");
	gre_sys->SetFillColor(kBlack);
	gre_sys->SetFillStyle(0);
	gre_sys->SetPointError(1,0.05, 0.007);
	gre_sys->SetPointError(2,0.05, 0.052);
	gre_sys->SetPointError(3,0.05, 0.275);
	gre_sys->SetPointError(4,0.05, 0.312);
	gre_sys->SetPointError(5,0.05, 0.198);
	gre_sys->SetMarkerStyle(21);
	gre_sys->SetMarkerColor(kBlack);
	gre_sys->SetLineColor(kGray+2);
	gre_sys->SetLineWidth(1);
	gre_sys->SetMarkerSize(1);

	TH1F *Graph_g_ALICE1 = new TH1F("Graph_g_ALICE1","J/#psi Yield vs. FVTX multiplicity (Normalized to MB)",100,0.06899999,4.641);
	Graph_g_ALICE1->SetMinimum(0);
	Graph_g_ALICE1->SetMaximum(6.5);
	Graph_g_ALICE1->SetDirectory(0);
	Graph_g_ALICE1->SetStats(0);

	Graph_g_ALICE1->SetLineColor(kGray);
	Graph_g_ALICE1->GetXaxis()->SetTitle("#frac{N_{ch}}{<N_{ch}^{MB}>}");
	Graph_g_ALICE1->GetXaxis()->SetNdivisions(309);
	Graph_g_ALICE1->GetXaxis()->SetLabelFont(42);
	Graph_g_ALICE1->GetXaxis()->SetLabelSize(0.05);
	Graph_g_ALICE1->GetXaxis()->SetTitleSize(0.05);
	Graph_g_ALICE1->GetXaxis()->SetTitleOffset(1.23);
	Graph_g_ALICE1->GetXaxis()->SetTitleFont(42);
	Graph_g_ALICE1->GetYaxis()->SetTitle("#frac{N_{J/#psi}/N_{MB}}{<N_{J/#psi}/N_{MB}>}");
	Graph_g_ALICE1->GetYaxis()->CenterTitle(true);
	Graph_g_ALICE1->GetYaxis()->SetNdivisions(309);
	Graph_g_ALICE1->GetYaxis()->SetLabelFont(42);
	Graph_g_ALICE1->GetYaxis()->SetLabelSize(0.05);
	Graph_g_ALICE1->GetYaxis()->SetTitleSize(0.05);
	Graph_g_ALICE1->GetYaxis()->SetTitleFont(42);
	Graph_g_ALICE1->GetZaxis()->SetLabelFont(42);
	Graph_g_ALICE1->GetZaxis()->SetLabelSize(0.035);
	Graph_g_ALICE1->GetZaxis()->SetTitleSize(0.035);
	Graph_g_ALICE1->GetZaxis()->SetTitleFont(42);
	gre->SetHistogram(Graph_g_ALICE1);

	TH1F *h_jpsi_mb_ratio_sys_tot__1 = new TH1F("h_jpsi_mb_ratio_sys_tot__1","ALumi_00_10",40,-0.1393602,11.00945);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(1,-0.2491511);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(2,0.2367082);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(3,0.4167507);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(4,0.6989338);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(5,1.053275);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(6,1.391411);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(7,1.739571);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(8,2.131778);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(9,2.602424);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(10,2.951751);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(11,3.397758);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(12,3.658449);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(13,4.231115);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(14,4.318662);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(15,4.960003);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(16,4.267373);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(17,5.751378);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(18,5.63537);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(19,4.86447);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(20,14.28462);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(21,4.451932);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(22,8.569675);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(23,9.796638);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(24,4.086076);
	h_jpsi_mb_ratio_sys_tot__1->SetBinContent(25,2.363666);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(1,0.2956494);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(2,0.04366202);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(3,0.04969899);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(4,0.07121363);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(5,0.08756496);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(6,0.09810783);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(7,0.1218351);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(8,0.1347862);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(9,0.1544787);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(10,0.2046136);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(11,0.1976428);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(12,0.2231408);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(13,0.2905741);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(14,0.3801666);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(15,0.5210414);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(16,0.5019844);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(17,0.4724341);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(18,0.5538025);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(19,1.281816);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(20,7.932367);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(21,6.499417);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(22,6.358303);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(23,5.741969);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(24,1.543781);
	h_jpsi_mb_ratio_sys_tot__1->SetBinError(25,1.933283);
	h_jpsi_mb_ratio_sys_tot__1->SetEntries(40);
	h_jpsi_mb_ratio_sys_tot__1->SetDirectory(0);
	h_jpsi_mb_ratio_sys_tot__1->SetFillStyle(0);
	h_jpsi_mb_ratio_sys_tot__1->SetMarkerColor(kBlack);

	TGraphErrors* gsys = new TGraphErrors();
	gsys->SetMarkerStyle(20);
	gsys->SetMarkerColor(kRed+1);
	gsys->SetLineColor(kRed+1);
	gsys->SetLineWidth(3);
	for (int i=1; i<=h_jpsi_mb_ratio_sys_tot__1->GetNbinsX(); i++)
	{
		float nch = h_jpsi_mb_ratio_sys_tot__1->GetXaxis()->GetBinCenter(i);
		float yield = h_jpsi_mb_ratio_sys_tot__1->GetBinContent(i);
		float yield_error = h_jpsi_mb_ratio_sys_tot__1->GetBinError(i);
		gsys->SetPoint(i-1, nch, yield);
		gsys->SetPointError(i-1, 0.05, yield_error);
	}
	gsys->SetFillColor(kRed-10);

	gsys->GetXaxis()->SetRangeUser(0.5,4.8);
	gsys->GetYaxis()->SetRangeUser(0.0,6.5);
	gsys->SetTitle("J/#psi Yield vs. FVTX multiplicity (Normalized to MB);#frac{N_{ch}}{<N_{ch}^{MB}>};#frac{N_{J/#psi}/N_{MB}}{<N_{J/#psi}/N_{MB}>}");
	gsys->SetHistogram(Graph_g_ALICE1);

	gsys->Draw("ae2");

	if(draw_alice)
	{
		//gre->Draw("p");
		gre_sys->Draw("e2same");
		gre->Draw("p");
	}

	TGraphErrors* gsys_cont = (TGraphErrors*)gsys->Clone("gsys_cont");
	gsys_cont->SetFillStyle(0);
	gsys_cont->SetFillColor(kBlack);
	gsys_cont->SetLineWidth(1);
	gsys_cont->Draw("e2");
	if(draw_alice)
	{
		gre_sys->Draw("e2");
		gre->Draw("p");
	}

	TH1F *h_jpsi_mb_ratio_stat__2 = new TH1F("h_jpsi_mb_ratio_stat__2","ALumi_00_10",40,-0.1393602,11.00945);
	h_jpsi_mb_ratio_stat__2->SetBinContent(1,-0.2491511);
	h_jpsi_mb_ratio_stat__2->SetBinContent(2,0.2367082);
	h_jpsi_mb_ratio_stat__2->SetBinContent(3,0.4167507);
	h_jpsi_mb_ratio_stat__2->SetBinContent(4,0.6989338);
	h_jpsi_mb_ratio_stat__2->SetBinContent(5,1.053275);
	h_jpsi_mb_ratio_stat__2->SetBinContent(6,1.391411);
	h_jpsi_mb_ratio_stat__2->SetBinContent(7,1.739571);
	h_jpsi_mb_ratio_stat__2->SetBinContent(8,2.131778);
	h_jpsi_mb_ratio_stat__2->SetBinContent(9,2.602424);
	h_jpsi_mb_ratio_stat__2->SetBinContent(10,2.951751);
	h_jpsi_mb_ratio_stat__2->SetBinContent(11,3.397758);
	h_jpsi_mb_ratio_stat__2->SetBinContent(12,3.658449);
	h_jpsi_mb_ratio_stat__2->SetBinContent(13,4.231115);
	h_jpsi_mb_ratio_stat__2->SetBinContent(14,4.318662);
	h_jpsi_mb_ratio_stat__2->SetBinContent(15,4.960003);
	h_jpsi_mb_ratio_stat__2->SetBinContent(16,4.267373);
	h_jpsi_mb_ratio_stat__2->SetBinContent(17,5.751378);
	h_jpsi_mb_ratio_stat__2->SetBinContent(18,5.63537);
	h_jpsi_mb_ratio_stat__2->SetBinContent(19,4.86447);
	h_jpsi_mb_ratio_stat__2->SetBinContent(20,14.28462);
	h_jpsi_mb_ratio_stat__2->SetBinContent(21,4.451932);
	h_jpsi_mb_ratio_stat__2->SetBinContent(22,8.569675);
	h_jpsi_mb_ratio_stat__2->SetBinContent(23,9.796638);
	h_jpsi_mb_ratio_stat__2->SetBinContent(24,4.086076);
	h_jpsi_mb_ratio_stat__2->SetBinContent(25,2.363666);
	h_jpsi_mb_ratio_stat__2->SetBinError(1,0.2956494);
	h_jpsi_mb_ratio_stat__2->SetBinError(2,0.01402017);
	h_jpsi_mb_ratio_stat__2->SetBinError(3,0.004105629);
	h_jpsi_mb_ratio_stat__2->SetBinError(4,0.006996293);
	h_jpsi_mb_ratio_stat__2->SetBinError(5,0.01058105);
	h_jpsi_mb_ratio_stat__2->SetBinError(6,0.01504929);
	h_jpsi_mb_ratio_stat__2->SetBinError(7,0.02069979);
	h_jpsi_mb_ratio_stat__2->SetBinError(8,0.02851407);
	h_jpsi_mb_ratio_stat__2->SetBinError(9,0.04013851);
	h_jpsi_mb_ratio_stat__2->SetBinError(10,0.05404271);
	h_jpsi_mb_ratio_stat__2->SetBinError(11,0.07540871);
	h_jpsi_mb_ratio_stat__2->SetBinError(12,0.1009367);
	h_jpsi_mb_ratio_stat__2->SetBinError(13,0.1422414);
	h_jpsi_mb_ratio_stat__2->SetBinError(14,0.1859517);
	h_jpsi_mb_ratio_stat__2->SetBinError(15,0.2643097);
	h_jpsi_mb_ratio_stat__2->SetBinError(16,0.3137726);
	h_jpsi_mb_ratio_stat__2->SetBinError(17,0.5042191);
	h_jpsi_mb_ratio_stat__2->SetBinError(18,0.6728657);
	h_jpsi_mb_ratio_stat__2->SetBinError(19,0.7879694);
	h_jpsi_mb_ratio_stat__2->SetBinError(20,2.377066);
	h_jpsi_mb_ratio_stat__2->SetBinError(21,1.680578);
	h_jpsi_mb_ratio_stat__2->SetBinError(22,2.85569);
	h_jpsi_mb_ratio_stat__2->SetBinError(23,6.917796);
	h_jpsi_mb_ratio_stat__2->SetBinError(24,4.078428);
	h_jpsi_mb_ratio_stat__2->SetBinError(25,2.363666);
	h_jpsi_mb_ratio_stat__2->SetEntries(40);
	h_jpsi_mb_ratio_stat__2->SetDirectory(0);
	h_jpsi_mb_ratio_stat__2->SetMarkerColor(kRed+1);
	h_jpsi_mb_ratio_stat__2->SetLineColor(kRed+1);

	h_jpsi_mb_ratio_stat__2->SetLineWidth(3);

	h_jpsi_mb_ratio_stat__2->SetMarkerStyle(20);
	h_jpsi_mb_ratio_stat__2->SetMarkerSize(2);
	h_jpsi_mb_ratio_stat__2->GetXaxis()->SetLabelFont(42);
	h_jpsi_mb_ratio_stat__2->GetXaxis()->SetLabelSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetXaxis()->SetTitleSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetXaxis()->SetTitleFont(42);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->CenterTitle(true);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->SetLabelFont(42);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->SetLabelSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->SetTitleSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->SetTitleOffset(1.3);
	h_jpsi_mb_ratio_stat__2->GetYaxis()->SetTitleFont(42);
	h_jpsi_mb_ratio_stat__2->GetZaxis()->SetLabelFont(42);
	h_jpsi_mb_ratio_stat__2->GetZaxis()->SetLabelSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetZaxis()->SetTitleSize(0.035);
	h_jpsi_mb_ratio_stat__2->GetZaxis()->SetTitleFont(42);
	h_jpsi_mb_ratio_stat__2->Draw("eX0same");

	if(draw_alice)
	{
	TLegend *leg = new TLegend(0.15,0.55,0.50,0.95,NULL,"brNDC");
		leg->SetBorderSize(0);
		leg->SetTextFont(62);
		leg->SetTextSize(0.030);
		leg->SetLineColor(1);
		leg->SetLineStyle(1);
		leg->SetLineWidth(1);
		leg->SetFillColor(0);
		leg->SetFillStyle(1001);
		//leg->AddEntry(h_jpsi_mb_ratio_stat__2,"PHENIX p+p #sqrt{s}=510 GeV J/#psi: 1.2<|y|<2.2 Chg. Mult. 1.2<|#eta|<2.2","P");
		//leg->AddEntry("g_ALICE","ALICE p+p #sqrt{s}=7 TeV  J/#psi: 2.5<y<4 Chg. Mult. |#eta|<1","p");
		leg->AddEntry(h_jpsi_mb_ratio_stat__2,"PHENIX p+p #sqrt{s}=510 GeV","P");
		leg->AddEntry((TObject*)0,"J/#psi: 1.2<|y|<2.2 Chg. Mult. 1.2<|#eta|<2.2","");
		leg->AddEntry((TObject*)0,"Chg. Mult. included J/#psi decay muons if reconstructed","");
		leg->AddEntry((TObject*)0,"<N_{ch}^{MB}> = 3.597 #pm 0.004(stat.)","");
		leg->AddEntry((TObject*)0,"","");
		leg->AddEntry("g_ALICE","ALICE p+p #sqrt{s}=7 TeV","p");
		leg->AddEntry((TObject*)0,"J/#psi: 2.5<|y|<4.0 Chg. Mult. |#eta|<1","");
		leg->AddEntry((TObject*)0,"<N_{ch}^{MB}> = 6.01 #pm 0.01(stat.) ^{+0.20}_{-0.12}(syst.)","");
		leg->Draw();
	}
	else
	{
		TPaveText *pav = new TPaveText(0.15,0.55,0.50,0.95,"NDC");
		pav->SetTextSize(0.028);
		pav->SetTextAlign(12);
		pav->SetFillColor(0);
		pav->SetBorderSize(0);
		pav->AddText("PHENIX p+p #sqrt{s}=510 GeV");
		pav->AddText("J/#psi: 1.2<|y|<2.2 Chg. Mult. 1.2<|#eta|<2.2");
		pav->AddText("Chg. Mult. included J/#psi decay muons if reconstructed");
		pav->AddText("<N_{ch}^{MB}> = 3.597 #pm 0.004(stat.)");
		pav->Draw();
	}

	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->Modified();
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->cd();
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SetSelected(c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->ToggleToolBar();

	std::string name = Form("result_one_vertex_sys_%d",draw_alice);
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SaveAs((name+".root").data());
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SaveAs((name+".eps").data());
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SaveAs((name+".png").data());
	c_compare_ALICE_sys_stath_jpsi_mb_ratio_sys_totc->SaveAs((name+".pdf").data());
}

//Kopeliovitch data points

/*
	 0.515800822047032, 0.32679738562091415
	 0.6500460436179051, 0.5446623093681915
	 0.8565685151494721, 0.871459694989107
	 1.0114322934213762, 1.0893246187363843
	 1.2798778159602893, 1.4814814814814805
	 1.5379466792444336, 1.8082788671023966
	 1.8579385935359258, 2.200435729847491
	 2.3327718257978987, 2.7886710239651418
	 2.7250185297486693, 3.2679738562091463
	 2.952092177077012, 3.5294117647058814
	 3.210116119758324, 3.812636165577339
	 3.4372122273880916, 4.0958605664487955
	 3.757159221076742, 4.444444444444443
	 4.149361004424674, 4.880174291938998
	 4.407362486804567, 5.1416122004357305

	 0.5051546391752576, 0
	 0.5982076679469035, 0.26143790849673276
	 0.6808840374637831, 0.45751633986928136
	 0.8255059182894222, 0.7407407407407415
	 0.9907912764189291, 1.0675381263616548
	 1.197358668553331, 1.4379084967320255
	 1.38324012308245, 1.7429193899782118
	 1.5588122992610558, 2.047930283224398
	 1.7447162140915908, 2.374727668845315
	 2.0234934752824385, 2.7886710239651373
	 2.2506794241178687, 3.1590413943355116
	 2.5397659636592307, 3.5729847494553324
	 2.78759292949711, 3.9651416122004344
	 3.014711497428291, 4.270152505446622
	 3.252184264312825, 4.618736383442261
	 3.4690160141949096, 4.945533769063175
	 3.7787211104372997, 5.359477124182998
	 4.057475911326728, 5.751633986928097
	 4.398153763223497, 6.209150326797382

*/