BOOL CDlgHisto::OnInitDialog()
{
CDialog::OnInitDialog();
DrawHisto();
return TRUE; // return TRUE unless you set the focus to a control
// EXCEPTION: OCX Property Pages should return FALSE
}
void PlotJetScheme()
{
gROOT->LoadMacro("$LOCAL/include/TUntilsOpHisto.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsAliFigs.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsPlotStd.h"); SetStyle();
//=============================================================================
TFile *file = TFile::Open("data/incl/AnalysisOutputs_jet_py8_JetR05_SjeR01.root", "READ");
TH1D *h5py8I = (TH1D*)file->Get("hJet"); h5py8I->SetName("h5py8I"); h5py8I->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_py6_JetR05_SjeR01.root", "READ");
TH1D *h5py6I = (TH1D*)file->Get("hJet"); h5py6I->SetName("h5py6I"); h5py6I->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_qpy_JetR05_SjeR01.root", "READ");
TH1D *h5qpyI = (TH1D*)file->Get("hJet"); h5qpyI->SetName("h5qpyI"); h5qpyI->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_vac_JetR05_SjeR01.root", "READ");
TH1D *h5vacI = (TH1D*)file->Get("hJet"); h5vacI->SetName("h5vacI"); h5vacI->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_med_JetR05_SjeR01.root", "READ");
TH1D *h5medI = (TH1D*)file->Get("hJet"); h5medI->SetName("h5medI"); h5medI->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_py8_JetR04_SjeR01.root", "READ");
TH1D *h4py8I = (TH1D*)file->Get("hJet"); h4py8I->SetName("h4py8I"); h4py8I->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_py6_JetR04_SjeR01.root", "READ");
TH1D *h4py6I = (TH1D*)file->Get("hJet"); h4py6I->SetName("h4py6I"); h4py6I->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_qpy_JetR04_SjeR01.root", "READ");
TH1D *h4qpyI = (TH1D*)file->Get("hJet"); h4qpyI->SetName("h4qpyI"); h4qpyI->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_vac_JetR04_SjeR01.root", "READ");
TH1D *h4vacI = (TH1D*)file->Get("hJet"); h4vacI->SetName("h4vacI"); h4vacI->SetDirectory(0);
file->Close();
file = TFile::Open("data/incl/AnalysisOutputs_jet_med_JetR04_SjeR01.root", "READ");
TH1D *h4medI = (TH1D*)file->Get("hJet"); h4medI->SetName("h4medI"); h4medI->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_py8_JetR05_SjeR01.root", "READ");
TH1D *h5py8M = (TH1D*)file->Get("hJet"); h5py8M->SetName("h5py8M"); h5py8M->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_py6_JetR05_SjeR01.root", "READ");
TH1D *h5py6M = (TH1D*)file->Get("hJet"); h5py6M->SetName("h5py6M"); h5py6M->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_qpy_JetR05_SjeR01.root", "READ");
TH1D *h5qpyM = (TH1D*)file->Get("hJet"); h5qpyM->SetName("h5qpyM"); h5qpyM->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_vac_JetR05_SjeR01.root", "READ");
TH1D *h5vacM = (TH1D*)file->Get("hJet"); h5vacM->SetName("h5vacM"); h5vacM->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_med_JetR05_SjeR01.root", "READ");
TH1D *h5medM = (TH1D*)file->Get("hJet"); h5medM->SetName("h5medM"); h5medM->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_py8_JetR04_SjeR01.root", "READ");
TH1D *h4py8M = (TH1D*)file->Get("hJet"); h4py8M->SetName("h4py8M"); h4py8M->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_py6_JetR04_SjeR01.root", "READ");
TH1D *h4py6M = (TH1D*)file->Get("hJet"); h4py6M->SetName("h4py6M"); h4py6M->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_qpy_JetR04_SjeR01.root", "READ");
TH1D *h4qpyM = (TH1D*)file->Get("hJet"); h4qpyM->SetName("h4qpyM"); h4qpyM->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_vac_JetR04_SjeR01.root", "READ");
TH1D *h4vacM = (TH1D*)file->Get("hJet"); h4vacM->SetName("h4vacM"); h4vacM->SetDirectory(0);
file->Close();
file = TFile::Open("data/mass/AnalysisOutputs_jet_med_JetR04_SjeR01.root", "READ");
TH1D *h4medM = (TH1D*)file->Get("hJet"); h4medM->SetName("h4medM"); h4medM->SetDirectory(0);
file->Close();
//=============================================================================
TH1F *hfm = 0;
TCanvas *can = 0;
TLegend *leg = 0;
TLatex *tex = 0;
TLine *line = 0;
const Float_t dflx = -40., dfux = 200.;
const Float_t dfly = 8e-10, dfuy = 8e3;
const Float_t dlsx = 0.05, dlsy = 0.05;
const Float_t dtsx = 0.06, dtsy = 0.06;
const Float_t dtox = 1.10, dtoy = 1.00;
const TString stnx = "#it{p}_{T,jet}^{ch} (GeV/#it{c})";
const TString stny = "d#sigma/d#it{p}_{T} (mb/GetV/#it{c})";
TH1D *hSE = new TH1D("hSE", "", 10., 0., 10.);
hSE->SetLineWidth(2);
hSE->SetLineColor(wcl[0]);
hSE->SetMarkerStyle(wmk[0]);
hSE->SetMarkerColor(wcl[0]);
TH1D *hME = new TH1D("hME", "", 10., 0., 10.);
hME->SetLineWidth(2);
hME->SetLineColor(wcl[0]);
hME->SetMarkerStyle(wmk[2]);
hME->SetMarkerColor(wcl[0]);
//=============================================================================
const Double_t dBin[] = { 0., 10., 20., 30., 40., 50., 60., 80., 100., 120., 150., 200., 250., 300., 350., 400., 500. };
const Int_t nBin = sizeof(dBin) / sizeof(Double_t) -1;
TH1D *h5py8Irb = h5py8I->Rebin(nBin, "h5py8Irb", dBin); TH1D *h4py8Irb = h4py8I->Rebin(nBin, "h4py8Irb", dBin);
TH1D *h5py6Irb = h5py6I->Rebin(nBin, "h5py6Irb", dBin); TH1D *h4py6Irb = h4py6I->Rebin(nBin, "h4py6Irb", dBin);
TH1D *h5qpyIrb = h5qpyI->Rebin(nBin, "h5qpyIrb", dBin); TH1D *h4qpyIrb = h4qpyI->Rebin(nBin, "h4qpyIrb", dBin);
TH1D *h5vacIrb = h5vacI->Rebin(nBin, "h5vacIrb", dBin); TH1D *h4vacIrb = h4vacI->Rebin(nBin, "h4vacIrb", dBin);
TH1D *h5medIrb = h5medI->Rebin(nBin, "h5medIrb", dBin); TH1D *h4medIrb = h4medI->Rebin(nBin, "h4medIrb", dBin);
TH1D *h5py8Mrb = h5py8M->Rebin(nBin, "h5py8Mrb", dBin); TH1D *h4py8Mrb = h4py8M->Rebin(nBin, "h4py8Mrb", dBin);
TH1D *h5py6Mrb = h5py6M->Rebin(nBin, "h5py6Mrb", dBin); TH1D *h4py6Mrb = h4py6M->Rebin(nBin, "h4py6Mrb", dBin);
TH1D *h5qpyMrb = h5qpyM->Rebin(nBin, "h5qpyMrb", dBin); TH1D *h4qpyMrb = h4qpyM->Rebin(nBin, "h4qpyMrb", dBin);
TH1D *h5vacMrb = h5vacM->Rebin(nBin, "h5vacMrb", dBin); TH1D *h4vacMrb = h4vacM->Rebin(nBin, "h4vacMrb", dBin);
TH1D *h5medMrb = h5medM->Rebin(nBin, "h5medMrb", dBin); TH1D *h4medMrb = h4medM->Rebin(nBin, "h4medMrb", dBin);
//=============================================================================
TH1D *h5py8R = (TH1D*)h5py8Irb->Clone("h5py8R"); h5py8R->Divide(h5py8Mrb);
TH1D *h5qpyR = (TH1D*)h5qpyIrb->Clone("h5qpyR"); h5qpyR->Divide(h5qpyMrb);
TH1D *h5medR = (TH1D*)h5medIrb->Clone("h5medR"); h5medR->Divide(h5medMrb);
TH1D *h4py8R = (TH1D*)h4py8Irb->Clone("h4py8R"); h4py8R->Divide(h4py8Mrb);
TH1D *h4qpyR = (TH1D*)h4qpyIrb->Clone("h4qpyR"); h4qpyR->Divide(h4qpyMrb);
TH1D *h4medR = (TH1D*)h4medIrb->Clone("h4medR"); h4medR->Divide(h4medMrb);
//=============================================================================
can = MakeCanvas("Jet5R"); can->SetGridx(); can->SetGridy();
hfm = can->DrawFrame(0., 0.95, 500., 1.2); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy,
"#it{p}_{T,jet} (GeV/#it{c})",
"Ratio: BIpt scheme / E scheme");
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(h5py8R, wcl[3], wmk[0], "SAME");
DrawHisto(h5qpyR, wcl[1], wmk[0], "SAME");
DrawHisto(h5medR, wcl[5], wmk[0], "SAME");
leg = new TLegend(0.62, 0.64, 0.98, 0.88); SetupLegend(leg);
leg->AddEntry(h5py8R, "PYTHIA 8", "LP")->SetTextSizePixels(24);
leg->AddEntry(h5qpyR, "Q-PYTHIA", "LP")->SetTextSizePixels(24);
leg->AddEntry(h5medR, "JEWEL+PYTHIA", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "2.76 TeV, Jets anti-#it{k}_{T}, #it{R}=0.5, |#eta_{jet}|<1");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("Jet4R"); can->SetGridx(); can->SetGridy();
hfm = can->DrawFrame(0., 0.95, 500., 1.2); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy,
"#it{p}_{T,jet} (GeV/#it{c})",
"Ratio: BIpt scheme / E scheme");
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(h4py8R, wcl[3], wmk[0], "SAME");
DrawHisto(h4qpyR, wcl[1], wmk[0], "SAME");
DrawHisto(h4medR, wcl[5], wmk[0], "SAME");
leg = new TLegend(0.62, 0.64, 0.98, 0.88); SetupLegend(leg);
leg->AddEntry(h4py8R, "PYTHIA 8", "LP")->SetTextSizePixels(24);
leg->AddEntry(h4qpyR, "Q-PYTHIA", "LP")->SetTextSizePixels(24);
leg->AddEntry(h4medR, "JEWEL+PYTHIA", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "2.76 TeV, Jets anti-#it{k}_{T}, #it{R}=0.4, |#eta_{jet}|<1");
CanvasEnd(can);
//=============================================================================
TH1D *h5qpyIAA = (TH1D*)h5qpyIrb->Clone("h5qpyIAA"); h5qpyIAA->Divide(h5py6Irb);
TH1D *h5medIAA = (TH1D*)h5medIrb->Clone("h5medIAA"); h5medIAA->Divide(h5vacIrb);
TH1D *h4qpyIAA = (TH1D*)h4qpyIrb->Clone("h4qpyIAA"); h4qpyIAA->Divide(h4py6Irb);
TH1D *h4medIAA = (TH1D*)h4medIrb->Clone("h4medIAA"); h4medIAA->Divide(h4vacIrb);
TH1D *h5qpyMAA = (TH1D*)h5qpyMrb->Clone("h5qpyMAA"); h5qpyMAA->Divide(h5py6Mrb);
TH1D *h5medMAA = (TH1D*)h5medMrb->Clone("h5medMAA"); h5medMAA->Divide(h5vacMrb);
TH1D *h4qpyMAA = (TH1D*)h4qpyMrb->Clone("h4qpyMAA"); h4qpyMAA->Divide(h4py6Mrb);
TH1D *h4medMAA = (TH1D*)h4medMrb->Clone("h4medMAA"); h4medMAA->Divide(h4vacMrb);
//=============================================================================
can = MakeCanvas("Jet5RAA"); can->SetGridx(); can->SetGridy();
hfm = can->DrawFrame(0., 0., 500., 1.2); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy,
"#it{p}_{T,jet} (GeV/#it{c})",
"#it{R}_{AA}");
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(h5qpyIAA, wcl[1], wmk[0], "SAME");
DrawHisto(h5medIAA, wcl[5], wmk[0], "SAME");
DrawHisto(h5qpyMAA, wcl[1], wmk[2], "SAME");
DrawHisto(h5medMAA, wcl[5], wmk[2], "SAME");
leg = new TLegend(0.42, 0.56, 0.98, 0.88); SetupLegend(leg);
leg->AddEntry(h5qpyIAA, "Q-PYTHIA, BIpt scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h5qpyMAA, "Q-PYTHIA, BIpt scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h5medIAA, "JEWEL+PYTHIA, E scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h5medMAA, "JEWEL+PYTHIA, E scheme", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "2.76 TeV, Jet anti-#it{k}_{T}, #it{R}=0.5, |#eta_{jet}|<1");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("Jet4RAA"); can->SetGridx(); can->SetGridy();
hfm = can->DrawFrame(0., 0., 500., 1.2); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy,
"#it{p}_{T,jet} (GeV/#it{c})",
"#it{R}_{AA}");
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(h4qpyIAA, wcl[1], wmk[0], "SAME");
DrawHisto(h4medIAA, wcl[5], wmk[0], "SAME");
DrawHisto(h4qpyMAA, wcl[1], wmk[2], "SAME");
DrawHisto(h4medMAA, wcl[5], wmk[2], "SAME");
leg = new TLegend(0.42, 0.56, 0.98, 0.88); SetupLegend(leg);
leg->AddEntry(h4qpyIAA, "Q-PYTHIA, BIpt scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h4qpyMAA, "Q-PYTHIA, BIpt scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h4medIAA, "JEWEL+PYTHIA, E scheme", "LP")->SetTextSizePixels(24);
leg->AddEntry(h4medMAA, "JEWEL+PYTHIA, E scheme", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "2.76 TeV, Jet anti-#it{k}_{T}, #it{R}=0.4, |#eta_{jet}|<1");
CanvasEnd(can);
//=============================================================================
return;
}
void PlotIsjMean()
{
gROOT->LoadMacro("$LOCAL/include/TUntilsOpHisto.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsAliFigs.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsPlotStd.h"); SetStyle();
//=============================================================================
TFile *file = TFile::Open("data/AnalysisOutputs_Isj.root", "READ");
TList *listPy8R4S1 = (TList*)file->Get("list_py8_JetR04_SjeR01");
TList *listPy8R4S2 = (TList*)file->Get("list_py8_JetR04_SjeR02");
TList *listPy8R4S3 = (TList*)file->Get("list_py8_JetR04_SjeR03");
TList *listPy8R5S1 = (TList*)file->Get("list_py8_JetR05_SjeR01");
TList *listPy8R5S2 = (TList*)file->Get("list_py8_JetR05_SjeR02");
TList *listPy8R5S3 = (TList*)file->Get("list_py8_JetR05_SjeR03");
TList *listPy8R5S4 = (TList*)file->Get("list_py8_JetR05_SjeR04");
TList *listVacR4S1 = (TList*)file->Get("list_vac_JetR04_SjeR01");
TList *listVacR4S2 = (TList*)file->Get("list_vac_JetR04_SjeR02");
TList *listVacR4S3 = (TList*)file->Get("list_vac_JetR04_SjeR03");
TList *listVacR5S1 = (TList*)file->Get("list_vac_JetR05_SjeR01");
TList *listVacR5S2 = (TList*)file->Get("list_vac_JetR05_SjeR02");
TList *listVacR5S3 = (TList*)file->Get("list_vac_JetR05_SjeR03");
TList *listVacR5S4 = (TList*)file->Get("list_vac_JetR05_SjeR04");
TList *listMedR4S1 = (TList*)file->Get("list_med_JetR04_SjeR01");
TList *listMedR4S2 = (TList*)file->Get("list_med_JetR04_SjeR02");
TList *listMedR4S3 = (TList*)file->Get("list_med_JetR04_SjeR03");
TList *listMedR5S1 = (TList*)file->Get("list_med_JetR05_SjeR01");
TList *listMedR5S2 = (TList*)file->Get("list_med_JetR05_SjeR02");
TList *listMedR5S3 = (TList*)file->Get("list_med_JetR05_SjeR03");
TList *listMedR5S4 = (TList*)file->Get("list_med_JetR05_SjeR04");
file->Close();
TProfile *hPy8R4S1 = (TProfile*)listPy8R4S1->FindObject("hJetIsj_mean"); hPy8R4S1->SetName("hPy8R4S1");
TProfile *hPy8R4S2 = (TProfile*)listPy8R4S2->FindObject("hJetIsj_mean"); hPy8R4S2->SetName("hPy8R4S2");
TProfile *hPy8R4S3 = (TProfile*)listPy8R4S3->FindObject("hJetIsj_mean"); hPy8R4S3->SetName("hPy8R4S3");
TProfile *hPy8R5S1 = (TProfile*)listPy8R5S1->FindObject("hJetIsj_mean"); hPy8R5S1->SetName("hPy8R5S1");
TProfile *hPy8R5S2 = (TProfile*)listPy8R5S2->FindObject("hJetIsj_mean"); hPy8R5S2->SetName("hPy8R5S2");
TProfile *hPy8R5S3 = (TProfile*)listPy8R5S3->FindObject("hJetIsj_mean"); hPy8R5S3->SetName("hPy8R5S3");
TProfile *hPy8R5S4 = (TProfile*)listPy8R5S4->FindObject("hJetIsj_mean"); hPy8R5S4->SetName("hPy8R5S4");
TProfile *hVacR4S1 = (TProfile*)listVacR4S1->FindObject("hJetIsj_mean"); hVacR4S1->SetName("hVacR4S1");
TProfile *hVacR4S2 = (TProfile*)listVacR4S2->FindObject("hJetIsj_mean"); hVacR4S2->SetName("hVacR4S2");
TProfile *hVacR4S3 = (TProfile*)listVacR4S3->FindObject("hJetIsj_mean"); hVacR4S3->SetName("hVacR4S3");
TProfile *hVacR5S1 = (TProfile*)listVacR5S1->FindObject("hJetIsj_mean"); hVacR5S1->SetName("hVacR5S1");
TProfile *hVacR5S2 = (TProfile*)listVacR5S2->FindObject("hJetIsj_mean"); hVacR5S2->SetName("hVacR5S2");
TProfile *hVacR5S3 = (TProfile*)listVacR5S3->FindObject("hJetIsj_mean"); hVacR5S3->SetName("hVacR5S3");
TProfile *hVacR5S4 = (TProfile*)listVacR5S4->FindObject("hJetIsj_mean"); hVacR5S4->SetName("hVacR5S4");
TProfile *hMedR4S1 = (TProfile*)listMedR4S1->FindObject("hJetIsj_mean"); hMedR4S1->SetName("hMedR4S1");
TProfile *hMedR4S2 = (TProfile*)listMedR4S2->FindObject("hJetIsj_mean"); hMedR4S2->SetName("hMedR4S2");
TProfile *hMedR4S3 = (TProfile*)listMedR4S3->FindObject("hJetIsj_mean"); hMedR4S3->SetName("hMedR4S3");
TProfile *hMedR5S1 = (TProfile*)listMedR5S1->FindObject("hJetIsj_mean"); hMedR5S1->SetName("hMedR5S1");
TProfile *hMedR5S2 = (TProfile*)listMedR5S2->FindObject("hJetIsj_mean"); hMedR5S2->SetName("hMedR5S2");
TProfile *hMedR5S3 = (TProfile*)listMedR5S3->FindObject("hJetIsj_mean"); hMedR5S3->SetName("hMedR5S3");
TProfile *hMedR5S4 = (TProfile*)listMedR5S4->FindObject("hJetIsj_mean"); hMedR5S4->SetName("hMedR5S4");
//=============================================================================
TH1F *hfm = 0;
TCanvas *can = 0;
TLegend *leg = 0;
TLatex *tex = 0;
TLine *line = 0;
const Float_t dflx = 0., dfux = 500.;
const Float_t dfly = 0., dfuy = 500.;
const Float_t dlsx = 0.05, dlsy = 0.05;
const Float_t dtsx = 0.06, dtsy = 0.06;
const Float_t dtox = 1.10, dtoy = 1.00;
const TString stnx = "#it{p}_{T,jet} (GeV/#it{c})";
const TString stny = "#it{p}_{T,subjet}^{mean} (GeV/#it{c})";
TH1D *hSE = new TH1D("hSE", "", 10., 0., 10.);
hSE->SetLineWidth(2);
hSE->SetLineColor(wcl[0]);
hSE->SetMarkerStyle(wmk[0]);
hSE->SetMarkerColor(wcl[0]);
TH1D *hME = new TH1D("hME", "", 10., 0., 10.);
hME->SetLineWidth(2);
hME->SetLineColor(wcl[0]);
hME->SetMarkerStyle(wmk[2]);
hME->SetMarkerColor(wcl[0]);
//=============================================================================
can = MakeCanvas("IsjPy8R04_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hPy8R4S1, wcl[1], wmk[0], "SAME");
DrawHisto(hPy8R4S2, wcl[3], wmk[0], "SAME");
DrawHisto(hPy8R4S3, wcl[4], wmk[0], "SAME");
leg = new TLegend(0.68, 0.66, 0.98, 0.90); SetupLegend(leg);
leg->AddEntry(hPy8R4S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hPy8R4S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hPy8R4S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA 8, tune 4C, pp #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.4, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("IsjPy8R05_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hPy8R5S1, wcl[1], wmk[0], "SAME");
DrawHisto(hPy8R5S2, wcl[3], wmk[0], "SAME");
DrawHisto(hPy8R5S3, wcl[4], wmk[0], "SAME");
DrawHisto(hPy8R5S4, wcl[2], wmk[0], "SAME");
leg = new TLegend(0.58, 0.74, 0.98, 0.90); SetupLegend(leg);
leg->SetNColumns(2);
leg->AddEntry(hPy8R5S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hPy8R5S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hPy8R5S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->AddEntry(hPy8R5S4, "#it{R}_{sub}=0.4", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA 8, tune 4C, pp #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.5, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("IsjVacR04_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hVacR4S1, wcl[1], wmk[0], "SAME");
DrawHisto(hVacR4S2, wcl[3], wmk[0], "SAME");
DrawHisto(hVacR4S3, wcl[4], wmk[0], "SAME");
leg = new TLegend(0.68, 0.66, 0.98, 0.90); SetupLegend(leg);
leg->AddEntry(hVacR4S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hVacR4S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hVacR4S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, vacuum, #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.4, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("IsjVacR05_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hVacR5S1, wcl[1], wmk[0], "SAME");
DrawHisto(hVacR5S2, wcl[3], wmk[0], "SAME");
DrawHisto(hVacR5S3, wcl[4], wmk[0], "SAME");
DrawHisto(hVacR5S4, wcl[2], wmk[0], "SAME");
leg = new TLegend(0.58, 0.74, 0.98, 0.90); SetupLegend(leg);
leg->SetNColumns(2);
leg->AddEntry(hVacR5S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hVacR5S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hVacR5S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->AddEntry(hVacR5S4, "#it{R}_{sub}=0.4", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, vacuum, #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.5, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("IsjMedR04_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hMedR4S1, wcl[1], wmk[0], "SAME");
DrawHisto(hMedR4S2, wcl[3], wmk[0], "SAME");
DrawHisto(hMedR4S3, wcl[4], wmk[0], "SAME");
leg = new TLegend(0.68, 0.66, 0.98, 0.90); SetupLegend(leg);
leg->AddEntry(hMedR4S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hMedR4S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hMedR4S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, medium, #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.4, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
can = MakeCanvas("IsjMedR05_mean"); can->SetGridx(); can->SetGridy(); //can->SetLogy();
hfm = can->DrawFrame(dflx,dfly,dfux,dfuy); SetupFrame(hfm,dlsx,dlsy,dtsx,dtsy,dtox,dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hMedR5S1, wcl[1], wmk[0], "SAME");
DrawHisto(hMedR5S2, wcl[3], wmk[0], "SAME");
DrawHisto(hMedR5S3, wcl[4], wmk[0], "SAME");
DrawHisto(hMedR5S4, wcl[2], wmk[0], "SAME");
leg = new TLegend(0.58, 0.74, 0.98, 0.90); SetupLegend(leg);
leg->SetNColumns(2);
leg->AddEntry(hMedR5S1, "#it{R}_{sub}=0.1", "LP")->SetTextSizePixels(24);
leg->AddEntry(hMedR5S2, "#it{R}_{sub}=0.2", "LP")->SetTextSizePixels(24);
leg->AddEntry(hMedR5S3, "#it{R}_{sub}=0.3", "LP")->SetTextSizePixels(24);
leg->AddEntry(hMedR5S4, "#it{R}_{sub}=0.4", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, medium, #sqrt{#it{s}_{NN}}=2.76 TeV");
tex->DrawLatex(0.16, 0.84, "Jet: anti-#it{k}_{T}, #it{R}=0.5, |#eta_{jet}|<1");
tex->DrawLatex(0.16, 0.76, "Sub-jet: #it{k}_{T}");
CanvasEnd(can);
//=============================================================================
return;
}
void CDlgHisto::OnKillfocusEdtNBins()
{
DrawHisto();
}
void PlotTrkPt()
{
gROOT->LoadMacro("$LOCAL/include/TUntilsOpHisto.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsAliFigs.h");
gROOT->LoadMacro("$LOCAL/include/TUntilsPlotStd.h"); SetStyle();
//=============================================================================
TFile *file = TFile::Open("data/AnalysisResults_vac_JetR05_skTR01.root", "READ");
TH2D *hTrkPtEtaVac = (TH2D*)file->Get("hTrkPtEta"); hTrkPtEtaVac->SetName("hTrkPtEtaVac"); hTrkPtEtaVac->SetDirectory(0);
file->Close();
file = TFile::Open("data/AnalysisResults_med_JetR05_skTR01.root", "READ");
TH2D *hTrkPtEtaMed = (TH2D*)file->Get("hTrkPtEta"); hTrkPtEtaMed->SetName("hTrkPtEtaMed"); hTrkPtEtaMed->SetDirectory(0);
file->Close();
//=============================================================================
TH1D *hTrkPtVac = hTrkPtEtaVac->ProjectionX("hTrkPtVac");
TH1D *hTrkPtMed = hTrkPtEtaMed->ProjectionX("hTrkPtMed");
hTrkPtVac->Rebin(10); hTrkPtVac->Scale(1./5.);
hTrkPtMed->Rebin(10); hTrkPtMed->Scale(1./5.);
//=============================================================================
TH1F *hfm = 0;
TCanvas *can = 0;
TLegend *leg = 0;
TLatex *tex = 0;
TLine *line = 0;
const Float_t dflx = 0., dfux = 300.;
const Float_t dfly = 8e-13, dfuy = 8e2;
const Float_t dlsx = 0.05, dlsy = 0.05;
const Float_t dtsx = 0.06, dtsy = 0.06;
const Float_t dtox = 1.10, dtoy = 1.00;
const TString stnx = "#it{p}_{T} (GeV/#it{c})";
const TString stny = "d#sigma/d#it{p}_{T} (mb/GeV/#it{c})";
TH1D *hSE = new TH1D("hSE", "", 10., 0., 10.);
hSE->SetLineWidth(2);
hSE->SetLineColor(wcl[0]);
hSE->SetMarkerStyle(wmk[0]);
hSE->SetMarkerColor(wcl[0]);
TH1D *hME = new TH1D("hME", "", 10., 0., 10.);
hME->SetLineWidth(2);
hME->SetLineColor(wcl[0]);
hME->SetMarkerStyle(wmk[2]);
hME->SetMarkerColor(wcl[0]);
//=============================================================================
can = MakeCanvas("TrkPt"); can->SetGridx(); can->SetGridy(); can->SetLogy();
hfm = can->DrawFrame(dflx, dfly, dfux, dfuy); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy, stnx, stny);
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hTrkPtVac, wcl[3], wmk[0], "SAME");
DrawHisto(hTrkPtMed, wcl[1], wmk[0], "SAME");
leg = new TLegend(0.62, 0.72, 0.98, 0.88); SetupLegend(leg);
leg->AddEntry(hTrkPtVac, "Vacuum", "LP")->SetTextSizePixels(24);
leg->AddEntry(hTrkPtMed, "Medium", "LP")->SetTextSizePixels(24);
leg->Draw();
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, 2.76 TeV, |#eta|<2.6");
CanvasEnd(can);
//=============================================================================
const Double_t dBins[] = { 0., 5., 10., 20., 30., 50., 80., 120., 160., 200., 250., 300., 400., 500. };
const Int_t nBins = sizeof(dBins) / sizeof(Double_t) - 1;
TH1D *hTrkPtVacRB = hTrkPtVac->Rebin(nBins, "hTrkPtVacRB", dBins);
TH1D *hTrkPtMedRB = hTrkPtMed->Rebin(nBins, "hTrkPtMedRB", dBins); hTrkPtMedRB->Divide(hTrkPtVacRB);
can = MakeCanvas("TrkPtRAA"); can->SetGridx(); can->SetGridy();
hfm = can->DrawFrame(dflx, 0., dfux, 2.); SetupFrame(hfm, dlsx, dlsy, dtsx, dtsy, dtox, dtoy, stnx, "#it{R}_{AA}");
hfm->GetXaxis()->SetNdivisions(510);
hfm->GetYaxis()->SetNdivisions(510);
DrawHisto(hTrkPtMedRB, wcl[0], wmk[0], "SAME");
tex = new TLatex();
tex->SetNDC();
tex->SetTextSizePixels(24);
tex->DrawLatex(0.16, 0.92, "PYTHIA+JEWEL, 2.76 TeV, |#eta|<2.6");
CanvasEnd(can);
//=============================================================================
return;
}
void DrawSpectraAndRatios() {
const Int_t nCent=7;
const Char_t *title[nCent]={
"0-5 %",
"5-10 %",
"10-20 %",
"20-40 %",
"40-60 %",
"60-80 %",
"80-90 %"
};
const Int_t colour[nCent]={2, 635, 797, 419, 4 , 6, 1 };
const Int_t marker[nCent]={22, 29, 34, 21, 23, 33, 20 };
const Float_t masize[nCent]={1.3, 1.6, 1.3, 1.2, 1.4, 1.8, 1.3};
const Float_t factor[nCent]={1.0, 1.0, 1.1, 1.5, 3.0,7.5,15.0}; //scale for drawing
const Float_t factor2[nCent]={1.0, 0.9, 0.9, 1.1, 2.0,5.5,15.0}; //scale for lin lambda drawing
const Float_t factor1[nCent]={1.0, 1/2., 1/4., 1/8., 1/16., 1/32., 1/64.}; //scale for log drawing
const Char_t *rNameL[2*nCent]={ // file name, histo name
"raw.root", "YieldLambda_0005",
"raw.root", "YieldLambda_0510",
"raw.root", "YieldLambda_1020",
"raw.root", "YieldLambda_2040",
"raw.root", "YieldLambda_4060",
"raw.root", "YieldLambda_6080",
"raw.root", "YieldLambda_8090"
};
const Char_t *eNameL[2*nCent]={ // file name, histo name
"eff.root", "eff_Lambda_comb_0005",
"eff.root", "eff_Lambda_comb_0510",
"eff.root", "eff_Lambda_comb_1020",
"eff.root", "eff_Lambda_comb_2040",
"eff.root", "eff_Lambda_comb_4060",
"eff.root", "eff_Lambda_comb_6080",
"eff.root", "eff_Lambda_comb_8090"
};
const Char_t *rNameK[2*nCent]={ // file name, histo name
"raw.root", "YieldK0Short_0005",
"raw.root", "YieldK0Short_0510",
"raw.root", "YieldK0Short_1020",
"raw.root", "YieldK0Short_2040",
"raw.root", "YieldK0Short_4060",
"raw.root", "YieldK0Short_6080",
"raw.root", "YieldK0Short_8090"
};
const Char_t *eNameK[2*nCent]={ // file name, histo name
"eff.root", "eff_K0s_comb_0005",
"eff.root", "eff_K0s_comb_0510",
"eff.root", "eff_K0s_comb_1020",
"eff.root", "eff_K0s_comb_2040",
"eff.root", "eff_K0s_comb_4060",
"eff.root", "eff_K0s_comb_6080",
"eff.root", "eff_K0s_comb_8090"
};
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
gStyle->SetLegendFillColor(0);
TH1 *raw=0;
TH1 *eff=0;
TString option(""), ratio("ratio");
TCanvas *c1=new TCanvas; c1->SetLogy();
c1->SetLeftMargin(0.13); c1->SetBottomMargin(0.13);
TCanvas *c1lin=new TCanvas;
c1lin->SetLeftMargin(0.13); c1lin->SetBottomMargin(0.13);
TCanvas *c2=new TCanvas; c2->SetLogy();
c2->SetLeftMargin(0.13); c2->SetBottomMargin(0.13);
TCanvas *c2lin=new TCanvas;
c2lin->SetLeftMargin(0.13); c2lin->SetBottomMargin(0.13);
TCanvas *c3=new TCanvas;
c3->SetLeftMargin(0.13); c3->SetBottomMargin(0.13);
TH1 *lkRatio[nCent]={0};
for (Int_t cent=0; cent<nCent; cent++) {
const Char_t *tit=title[cent];
Int_t col=colour[cent];
Int_t mar=marker[cent];
Float_t siz=masize[cent];
// Lambda
if (!GetHistos(rNameL+2*cent, eNameL+2*cent, raw, eff)) return;
TH1 *rawHl=MapHisto(raw);
TH1 *effHl=MapHisto(eff);
//Feed down
FeedDown(rawHl);
rawHl->Divide(effHl);
SetAttributes(rawHl,tit,col,mar,siz);
c1->cd();
DrawHisto(rawHl, option.Data(), sysEffLam, sysSigLam, factor1[cent]);
TH1 *linHl=(TH1*)rawHl->Clone();
SetAttributes(linHl,tit,col,mar,siz,0.,20.,32);
c1lin->cd();
DrawHisto(linHl, option.Data(), sysEffLam, sysSigLam, factor2[cent]);
// K0s
if (!GetHistos(rNameK+2*cent, eNameK+2*cent, raw, eff)) return;
TH1 *rawHk=MapHisto(raw);
TH1 *effHk=MapHisto(eff);
rawHk->Divide(effHk);
SetAttributes(rawHk,tit,col,mar,siz,1e-7);
c2->cd();
DrawHisto(rawHk, option.Data(), sysEffK0s, sysSigK0s, factor1[cent]);
TH1 *linHk=(TH1*)rawHk->Clone();
SetAttributes(linHk,tit,col,mar,siz,0.,120.,32);
c2lin->cd();
DrawHisto(linHk, option.Data(), sysEffK0s, sysSigK0s, factor[cent]);
// Lambda/K0s
TH1 *rawHlk=(TH1*)rawHl->Clone();
lkRatio[cent]=rawHlk;
TString name=ratio+rawHlk->GetName();
rawHlk->SetName(name.Data());
rawHlk->SetMaximum(1.7);
rawHlk->Divide(rawHk);
rawHlk->GetYaxis()->SetTitle("#Lambda/K^{0}_{S}");
c3->cd();
//if (cent!=1)
DrawRatio(rawHlk,option.Data());
option+="same";
}
for (Int_t cent=0; cent<nCent; cent++) {
//if (cent != 1)
lkRatio[cent]->Draw("same");
}
TLegend *leg=c1->BuildLegend(0.68,0.46,0.88,0.82,"Centrality:");
leg->SetBorderSize(0);
leg->SetFillColor(0);
TLegendEntry *entry=leg->AddEntry("NULL","systematic uncertainty","lpf");
Int_t ci = TColor::GetColor("#cccccc");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(10);
entry->SetMarkerColor(ci);
c1->cd();
TLatex * tex = new TLatex(0.5,0.65,"#Lambda");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.11);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.27,0.83,"Pb-Pb at #sqrt{s_{NN}}=2.76 TeV, |y|<0.5");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetLineWidth(2);
tex->Draw();
//Float_t offx=0.15, offy=0.16, sizx=0.22, sizy=0.22;
//DrawALICELogo(offx,offy,offx+sizx,offy+sizy);
leg=c1lin->BuildLegend(0.69,0.43,0.88,0.80,"Centrality:");
leg->SetBorderSize(0);
leg->SetFillColor(0);
entry=leg->AddEntry("NULL","systematic uncertainty","lpf");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(10);
entry->SetMarkerColor(ci);
entry=leg->AddEntry("NULL","BGBW fit","l");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
c1lin->cd();
tex=new TLatex(0.27,0.83,"Pb-Pb at #sqrt{s_{NN}}=2.76 TeV, |y|<0.5");
tex->SetNDC();
tex->Draw();
tex = new TLatex(0.5,0.65,"#Lambda");
tex->SetNDC();
tex->SetTextSize(0.11);
tex->Draw();
{
TFile::Open(" BWFitResults_Lambda0_stat.root");
const Char_t *name[nCent]={
"BWFit_0005",
"BWFit_0510","BWFit_1020","BWFit_2040",
"BWFit_4060","BWFit_6080","BWFit_8090"
};
DrawFit(name, factor2, nCent);
}
//Float_t offx1=0.70, offy1=0.18;
//DrawALICELogo(offx1,offy1,offx1+sizx,offy1+sizy);
leg=c2->BuildLegend(0.68,0.46,0.88,0.82,"Centrality:");
leg->SetBorderSize(0);
leg->SetFillColor(0);
entry=leg->AddEntry("NULL","systematic uncertainty","lpf");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(10);
entry->SetMarkerColor(ci);
c2->cd();
tex = new TLatex(0.5,0.65,"K^{0}_{S}");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.089);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.27,0.83,"Pb-Pb at #sqrt{s_{NN}}=2.76 TeV, |y|<0.5");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetLineWidth(2);
tex->Draw();
//DrawALICELogo(offx,offy,offx+sizx,offy+sizy);
leg=c2lin->BuildLegend(0.69,0.43,0.88,0.80,"Centrality:");
leg->SetBorderSize(0);
leg->SetFillColor(0);
entry=leg->AddEntry("NULL","systematic uncertainty","lpf");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(10);
entry->SetMarkerColor(ci);
entry=leg->AddEntry("NULL","BGBW fit","l");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
c2lin->cd();
tex = new TLatex(0.27,0.83,"Pb-Pb at #sqrt{s_{NN}}=2.76 TeV, |y|<0.5");
tex->SetNDC();
tex->Draw();
tex = new TLatex(0.5,0.65,"K^{0}_{S}");
tex->SetNDC();
tex->SetTextSize(0.089);
tex->Draw();
{
TFile::Open("BWFitResults_K0_stat.root");
const Char_t *name[nCent]={
"BWFit_0005",
"BWFit_0510","BWFit_1020","BWFit_2040",
"BWFit_4060","BWFit_6080","BWFit_8090"
};
DrawFit(name, factor, nCent);
}
//DrawALICELogo(offx1,offy1,offx1+sizx,offy1+sizy);
//leg=c3->BuildLegend(0.74,0.62,0.88,0.88,"Centrality:");
leg=c3->BuildLegend(0.55,0.55,0.88,0.88,"Centrality:");
leg->SetBorderSize(0);
leg->SetFillColor(0);
entry=leg->AddEntry("NULL","systematic uncertainty","lpf");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(10);
entry->SetMarkerColor(ci);
return;
}
