ROOT logo
//------------------------------------------------------------------------------
// makeCompYield.C
//
// create figure with ALICE vs. UA1 data
//------------------------------------------------------------------------------


void makeCompYield()
{
using namespace std;
 
TGraphErrors *graphYieldUa1 = new TGraphErrors(binsUa1,ptUa1,yieldUa1,errPtUa1,err2YieldUa1);
graphYieldUa1->SetMarkerColor(colorUa1);
graphYieldUa1->SetLineColor(colorUa1);
graphYieldUa1->SetMarkerStyle(26);

TGraphErrors *graphYieldAlice = new TGraphErrors(binsYieldAlice,ptYieldAlice,yieldAlice,0,err2YieldAlice);
graphYieldAlice->SetMarkerColor(colorAliceFit);
graphYieldAlice->SetLineColor(colorAliceFit);
graphYieldAlice->SetLineWidth(1);
graphYieldAlice->SetMarkerStyle(20);


TMultiGraph *all = new TMultiGraph("all","");


all->Add(graphYieldUa1);
//all->Add(graphCms);
//all->Add(graphAlice);
all->Add(graphYieldAlice);
//all->Draw("AP");

Double_t one[binsYieldAlice];  
for (Int_t i=0; i < binsYieldAlice; i++) { 
    one[i] = 1.0;
}

TGraphErrors *ratioUa1Alice = divide(graphYieldUa1,graphYieldAlice);
ratioUa1Alice->SetFillColor(colorUa1);
ratioUa1Alice->SetLineColor(colorUa1);
ratioUa1Alice->SetFillStyle(3345);


TGraphErrors *ratioAliceAlice = new TGraphErrors(binsYieldAlice,ptYieldAlice,one,0,relErr2YieldAlice);
ratioAliceAlice->SetFillColor(colorAliceErrors);
ratioAliceAlice->SetLineColor(20);

TCanvas *can4 = new TCanvas("can4","CompYield",520,700);

TPad *pad4_1 = new TPad("pad4_1","pad4_1",0.0,0.35,1.0,1.0);
setAttrib(pad4_1);

TPad *pad4_2 = new TPad("pad4_2","pad4_2",0.0,0.0,1.0,0.35);
setAttrib(pad4_2);

//
can4->cd();

pad4_1->Draw();
pad4_1->cd();
pad4_1->SetLogx();
pad4_1->SetLogy();
//all->Draw("AE4");
all->Draw("APZ");
all->GetXaxis()->SetTitle("p_{T} (GeV/c)");
all->GetXaxis()->SetTitleOffset(1.6);
all->GetYaxis()->SetTitleOffset(1.6);
all->GetYaxis()->SetTitle("1/N_{evt} 1/(2#pi p_{T}) (d^{2}N_{ch})/(dy dp_{T}) (GeV/c)^{-2}");
all->GetXaxis()->SetLimits(minPt,maxPt);
setAttrib(all);


graphYieldAlice->SetFillColor(2);

TLegend *l4_1 = new TLegend(0.200,0.03,0.65,0.203);
//TLegend *l4_1 = new TLegend(0.200,0.031,0.302,0.137);
l4_1->SetHeader("#sqrt{s} = 900 GeV");
l4_1->AddEntry(graphYieldAlice,"ALICE pp, NSD, | #eta | < 0.8","LP");
l4_1->AddEntry(graphYieldUa1,"UA1 p#bar{p}, NSD, | #eta | < 2.5","LP");
//l1->AddEntry(graphCms,"CMS |#eta|<2.4","LP");
//la->AddEntry(graphUA1,"UA1","p");
//la->AddEntry(graphCMS,"CMS","p");
//la->AddEntry(graphCMSred,"CMS, #eta<0.8","p");
//la->AddEntry(graphATLAS,"ATLAS","p");
l4_1->SetFillColor(0);
l4_1->SetLineColor(0);
l4_1->SetTextSize(legendTextSize);
l4_1->Draw();


//histo1->GetXaxis()->SetNdivisions(405);

can4->cd();
pad4_2->Draw();
pad4_2->cd();
pad4_2->SetLogx();

TMultiGraph *ratios = new TMultiGraph("ratios","");


ratios->Add(ratioAliceAlice);
//ratios->Add(ratioAtlasAlice);
ratios->Add(ratioUa1Alice);


//ratio->SetTitle("ratio ALICE/ATLAS (different #eta)");
ratios->Draw("AE3");
ratios->GetXaxis()->SetLimits(minPt,maxPt);
ratios->GetXaxis()->SetTitle("p_{T} (GeV/c)");
ratios->GetXaxis()->SetTitleOffset(1.4);
ratios->GetYaxis()->SetTitleOffset(0.9);
ratios->GetYaxis()->SetTitle("ratio");
ratios->GetYaxis()->SetLabelSize(0.07);
ratios->GetXaxis()->SetLabelSize(0.07);
ratios->GetXaxis()->SetTitleSize(0.07);
ratios->GetYaxis()->SetTitleSize(0.07);
ratios->GetYaxis()->CenterTitle();
setAttrib(ratios);
//ratios->GetYaxis()->SetRangeUser(0.2,1.5);
//ratios->GetYaxis()->SetLimits(0.2,1.5);

ratios->GetYaxis()->SetLabelOffset(0.023);

graphYieldAlice->SetFillColor(2);

TLegend *l4_2 = new TLegend(0.200,0.323,0.65,0.537);
l4_2->AddEntry(ratioAliceAlice,"ALICE uncertainties","F");
l4_2->AddEntry(ratioUa1Alice,"UA1 / ALICE","F");
//l2->AddEntry(ratioCmsAlice,"CMS / ALICE","F");
l4_2->SetFillColor(0);
l4_2->SetLineColor(0);
l4_2->SetTextSize(legendTextSize);
l4_2->Draw();


//aliceData->Draw("AE4");
TF1 *fOne = new TF1("one","1",0.1,20);
fOne->SetRange(0.1,20);
fOne->SetLineWidth(1);
fOne->Draw("SAME");
//ratio->GetYaxis()->SetTitle("ratio ALICE/ATLAS");
/*
histo1c->Draw();
histo1c->GetXaxis()->SetLabelSize(0.08);
histo1c->GetYaxis()->SetLabelSize(0.08);
histo1c->GetXaxis()->SetNdivisions(405);
histo1c->GetYaxis()->SetNdivisions(405);
*/

logoPrelim(can4);
}
 makeCompYield.C:1
 makeCompYield.C:2
 makeCompYield.C:3
 makeCompYield.C:4
 makeCompYield.C:5
 makeCompYield.C:6
 makeCompYield.C:7
 makeCompYield.C:8
 makeCompYield.C:9
 makeCompYield.C:10
 makeCompYield.C:11
 makeCompYield.C:12
 makeCompYield.C:13
 makeCompYield.C:14
 makeCompYield.C:15
 makeCompYield.C:16
 makeCompYield.C:17
 makeCompYield.C:18
 makeCompYield.C:19
 makeCompYield.C:20
 makeCompYield.C:21
 makeCompYield.C:22
 makeCompYield.C:23
 makeCompYield.C:24
 makeCompYield.C:25
 makeCompYield.C:26
 makeCompYield.C:27
 makeCompYield.C:28
 makeCompYield.C:29
 makeCompYield.C:30
 makeCompYield.C:31
 makeCompYield.C:32
 makeCompYield.C:33
 makeCompYield.C:34
 makeCompYield.C:35
 makeCompYield.C:36
 makeCompYield.C:37
 makeCompYield.C:38
 makeCompYield.C:39
 makeCompYield.C:40
 makeCompYield.C:41
 makeCompYield.C:42
 makeCompYield.C:43
 makeCompYield.C:44
 makeCompYield.C:45
 makeCompYield.C:46
 makeCompYield.C:47
 makeCompYield.C:48
 makeCompYield.C:49
 makeCompYield.C:50
 makeCompYield.C:51
 makeCompYield.C:52
 makeCompYield.C:53
 makeCompYield.C:54
 makeCompYield.C:55
 makeCompYield.C:56
 makeCompYield.C:57
 makeCompYield.C:58
 makeCompYield.C:59
 makeCompYield.C:60
 makeCompYield.C:61
 makeCompYield.C:62
 makeCompYield.C:63
 makeCompYield.C:64
 makeCompYield.C:65
 makeCompYield.C:66
 makeCompYield.C:67
 makeCompYield.C:68
 makeCompYield.C:69
 makeCompYield.C:70
 makeCompYield.C:71
 makeCompYield.C:72
 makeCompYield.C:73
 makeCompYield.C:74
 makeCompYield.C:75
 makeCompYield.C:76
 makeCompYield.C:77
 makeCompYield.C:78
 makeCompYield.C:79
 makeCompYield.C:80
 makeCompYield.C:81
 makeCompYield.C:82
 makeCompYield.C:83
 makeCompYield.C:84
 makeCompYield.C:85
 makeCompYield.C:86
 makeCompYield.C:87
 makeCompYield.C:88
 makeCompYield.C:89
 makeCompYield.C:90
 makeCompYield.C:91
 makeCompYield.C:92
 makeCompYield.C:93
 makeCompYield.C:94
 makeCompYield.C:95
 makeCompYield.C:96
 makeCompYield.C:97
 makeCompYield.C:98
 makeCompYield.C:99
 makeCompYield.C:100
 makeCompYield.C:101
 makeCompYield.C:102
 makeCompYield.C:103
 makeCompYield.C:104
 makeCompYield.C:105
 makeCompYield.C:106
 makeCompYield.C:107
 makeCompYield.C:108
 makeCompYield.C:109
 makeCompYield.C:110
 makeCompYield.C:111
 makeCompYield.C:112
 makeCompYield.C:113
 makeCompYield.C:114
 makeCompYield.C:115
 makeCompYield.C:116
 makeCompYield.C:117
 makeCompYield.C:118
 makeCompYield.C:119
 makeCompYield.C:120
 makeCompYield.C:121
 makeCompYield.C:122
 makeCompYield.C:123
 makeCompYield.C:124
 makeCompYield.C:125
 makeCompYield.C:126
 makeCompYield.C:127
 makeCompYield.C:128
 makeCompYield.C:129
 makeCompYield.C:130
 makeCompYield.C:131
 makeCompYield.C:132
 makeCompYield.C:133
 makeCompYield.C:134
 makeCompYield.C:135
 makeCompYield.C:136
 makeCompYield.C:137
 makeCompYield.C:138
 makeCompYield.C:139
 makeCompYield.C:140
 makeCompYield.C:141
 makeCompYield.C:142
 makeCompYield.C:143
 makeCompYield.C:144
 makeCompYield.C:145
 makeCompYield.C:146
 makeCompYield.C:147
 makeCompYield.C:148
 makeCompYield.C:149
 makeCompYield.C:150
 makeCompYield.C:151
 makeCompYield.C:152