PWG2/FLOW/papers/PRL105_252302/v2edepAlice.C

   1 // S. Voloshin   03-oct-2010\r
   2 // modified from plotV2edepBevalacRhicSqrtS.C - from Art   \r
   3 {\r
   4   float m_p=0.940;\r
   5   \r
   6   //  gROOT->Reset();\r
   7  \r
   8   //  gROOT->SetStyle("Bold");  \r
   9   //gROOT->SetStyle("Plain");  \r
  10   //  gStyle->SetTitleTextColor(kBlack);\r
  11   //  gStyle->SetLabelColor(kBlack,"X");\r
  12   //  gStyle->SetLabelColor(kBlack,"Y");\r
  13 \r
  14   int cCERES=kOrange+7;\r
  15   int cPHENIX=kYellow+3;\r
  16   int cFOPI=kGreen+2;\r
  17   int cE895=kMagenta+2;\r
  18   int cAlice=2;\r
  19   \r
  20   TCanvas *canvas = new TCanvas("v2edep","v2edep",10,10,750,600);\r
  21   canvas->cd();\r
  22 \r
  23   hist = new TH1F("v2 vs sqrt(s)","v2 vs sqrt(s)", 1, 1., 10000.);\r
  24   hist->SetLineColor(0);\r
  25   \r
  26   TAxis *axis = hist->GetXaxis();\r
  27   axis->SetTitle("#sqrt{s_{NN}} (GeV)");\r
  28   //axis->CenterTitle(kTRUE);\r
  29   axis->SetTitleSize(0.05);\r
  30   axis->SetLabelSize(0.04);\r
  31   axis->SetTitleOffset(1.5);\r
  32   axis = hist->GetYaxis();\r
  33   axis->SetTitle("v_{2}(y^{*}=0)");\r
  34   axis->SetTitle("v_{2}");\r
  35   axis->SetTitleOffset(1.6);\r
  36   //axis->CenterTitle(kTRUE);\r
  37   axis->SetTitleSize(0.05);\r
  38   axis->SetLabelSize(0.04);\r
  39  \r
  40   hist->SetStats(0); \r
  41   //  hist->SetTitle("Elliptic Flow");\r
  42   hist->SetTitle("");\r
  43   hist->SetMaximum(0.08); \r
  44   hist->SetMinimum(-0.085); \r
  45   //hist->SetLabelOffset(-0.01,"X");\r
  46   gPad->SetTopMargin(.1);\r
  47   gPad->SetBottomMargin(.2);\r
  48   gPad->SetLeftMargin(.2);\r
  49   gPad->SetRightMargin(.1);\r
  50   hist->Draw(); \r
  51   canvas->SetLogx(); \r
  52   hist->Draw(); \r
  53 \r
  54   TLine *tline=new TLine(1.,0.,10000.,0.);\r
  55   tline->SetLineWidth(1);\r
  56   tline->SetLineStyle(7); //wide dash\r
  57   tline->Draw("same");\r
  58   \r
  59   \r
  60   \r
  61   // CERES ----------------------------------------------------------------------------\r
  62   float xxr[3]={8.7,12.3,17.}; //sqrt(s)\r
  63   //float xxr[3]={9.57,18.7}; \r
  64   //float yyr[3]={0.028,0.031}; \r
  65   //float yyr[3]={0.028,0.031}; // taken from QM01 and INPC01\r
  66   //float yyr[3]={0.028,0.03,0.031}; // QM01 and INPC01, 80GeV interpolated\r
  67   //float yyr[3]={0.025,0.035,0.04}; // values integrated over phase space\r
  68   float yyr[3]={0.026,0.035,0.04}; // changed to those values after request by H. Appleshaeuser\r
  69   //CERES data from yugoslavian conference\r
  70   //float xxr[3]={9.57,13.5,18.7}; \r
  71   //float yyr[3]={0.028,0.033,0.036}; \r
  72    \r
  73   float exr[3]={0.,0.,0.}; \r
  74   float exr[3]={0.,0.,0.}; \r
  75   float eyr[3]={0.,0.,0.}; \r
  76   //float eyr[3]={0.0025,0.0025,0.0025}; \r
  77   float eyr[3]={0.005,0.005,0.005};  // changed to those values after request by H. Appleshaeuser\r
  78 \r
  79   //float eyr[3]={0.004,0.005,0.004}; \r
  80 \r
  81   ceres = new TGraphAsymmErrors(3,xxr,yyr,exr,exr,eyr,eyr); \r
  82   ceres->SetTitle(""); \r
  83   ceres->SetMarkerColor(cCERES); \r
  84   ceres->SetMarkerStyle(kOpenCircle); \r
  85   ceres->SetMarkerSize(1.2);\r
  86   ceres->SetLineWidth(2); \r
  87   ceres->SetLineColor(cCERES); \r
  88   ceres->Draw("p"); \r
  89   \r
  90   // NA49 --------------------------------------------------------------------------------\r
  91   float xxd[4]={8.7,17.,8.7,17.}; //sqrt(s)\r
  92   //  float yyd[3]={0.029867,0.0324}; //mod method with Full Acceptance (40)\r
  93   //float yyd[4]={0.023676,0.0324,0.023676,0.0324}; //mod method with Cos Cut (40)\r
  94   float yyd[4]={0.023676,0.0309,0.023676,0.0309}; //mod method with Cos Cut (40) '96 and '00 data (158)\r
  95   //float yyd[4]={0.0452,0.0379,0.0452,0.0379}; //mod method with Cos Cut (40) centrality 5;\r
  96   //float yyd[4]={0.0411,0.0406,0.0411,0.0406}; //mod method with Cos Cut (40) centrality 6;\r
  97   //float yyd[4]={0.0265,0.0341,0.0265,0.0341}; //mod method with Cos Cut (40) centrality 4;\r
  98   //float yyd[4]={0.0206,0.0276,0.0206,0.0276}; //mod method with Cos Cut (40) centrality 3;\r
  99   //float yyd[4]={-0.00495,0.0189,-0.00495,0.0189}; //mod method with Cos Cut (40) centrality 2;\r
 100   //hist->SetTitle("std method: centrality 2");\r
 101   float exd[4]={0.,0.,0.}; \r
 102   //  float eyd[3]={0.001734,0.000716}; //mod method with Full Acceptance (40)\r
 103   //float eyd[4]={0.00276,0.000716,0.005,0.005}; //mod method with Cos Cut (40)\r
 104   float eyd[4]={0.00276,0.00054,0.005,0.005}; //mod method with Cos Cut (40) '96 and '00 data (158)\r
 105   data = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); \r
 106   data->SetTitle(""); \r
 107   data->SetMarkerColor(kRed); \r
 108   data->SetMarkerStyle(kFullSquare); \r
 109   data->SetMarkerSize(1.2); \r
 110   data->SetLineWidth(2); \r
 111   data->SetLineColor(kRed); \r
 112   //data->Draw("p"); \r
 113   \r
 114   // NA49 \r
 115   float xxd[4]={8.7,17.}; //sqrt(s)\r
 116   float yyd[4]={0.023676,0.0324}; //mod method with Cos Cut (40)\r
 117   float exd[4]={0.,0.,0.}; \r
 118   float eyd[4]={0.01,0.005}; //mod method with Cos Cut (40)\r
 119   dataSys = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); \r
 120   dataSys->SetTitle(""); \r
 121   //dataSys->SetMarkerColor(kRed); \r
 122   //dataSys->SetMarkerStyle(kFullSquare); \r
 123   dataSys->SetMarkerSize(1.2); \r
 124   dataSys->SetLineWidth(2); \r
 125   dataSys->SetLineColor(kRed); \r
 126   //dataSys->Draw("[]"); \r
 127   \r
 128   \r
 129   // NA49 cumulant\r
 130   float xxd[2]={7.91,15.45};  //sqrt(s)\r
 131   float yyd[2]={0.02184,0.0291746}; //order 2\r
 132   //float yyd[2]={0.034883,0.027746}; //order 4\r
 133   float exd[3]={0.,0.,0.}; \r
 134   //float eyd[3]={0.0105334,0.0026345}; //order 4\r
 135   float eyd[3]={0.001173,0.000526268}; //order 2\r
 136   cumul = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); \r
 137   cumul->SetTitle(""); \r
 138   cumul->SetMarkerColor(kRed); \r
 139   cumul->SetMarkerStyle(kFullSquare); \r
 140   cumul->SetMarkerSize(1.2); \r
 141   cumul->SetLineWidth(2); \r
 142   cumul->SetLineColor(kBlue); \r
 143   cumul->Draw("p"); \r
 144   \r
 145   // STAR ------------------------------------------------------------------------------\r
 146   float xxd[2]={130.,200.}; //sqrt(s)\r
 147   float yyd[2]={0.0426,0.0478}; \r
 148   float exd[2]={0.}; \r
 149   float exd[2]={0.}; \r
 150   float eyd[2]={0.}; \r
 151   float eyd[2]={0.0026,0.0026}; \r
 152   star = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); \r
 153   star->SetTitle(""); \r
 154   star->SetMarkerColor(kRed); \r
 155   star->SetMarkerStyle(30); \r
 156   star->SetMarkerSize(1.6); \r
 157   star->SetLineWidth(2); \r
 158   star->SetLineColor(kRed); \r
 159   star->Draw("p"); \r
 160   \r
 161   \r
 162   // Phenix --------------------------------------------------------------------------\r
 163   float xxd[1]={220}; //sqrt(s) \r
 164   float yyd[1]={0.054}; \r
 165   float exd[1]={0.}; \r
 166   float exd[1]={0.}; \r
 167   float eyd[1]={0.}; \r
 168   float eyd[1]={0.0041}; \r
 169   phenix = new TGraphAsymmErrors(1,xxd,yyd,exd,exd,eyd,eyd); \r
 170   phenix->SetTitle(""); \r
 171   phenix->SetMarkerColor(cPHENIX); \r
 172   phenix->SetMarkerStyle(kOpenSquare); \r
 173   phenix->SetMarkerSize(1.2); \r
 174   phenix->SetLineWidth(2); \r
 175   phenix->SetLineColor(cPHENIX); \r
 176   phenix->Draw("p"); \r
 177   \r
 178 \r
 179   // Phobos -----------------------------------------------------------------------\r
 180   float xxd[2]={117.,180.}; //sqrt(s)\r
 181   float yyd[2]={0.048,0.051}; \r
 182   float exd[2]={0.}; \r
 183   float exd[2]={0.}; \r
 184   float eyd[2]={0.}; \r
 185   float eyd[2]={0.005,0.005}; \r
 186   phobos = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); \r
 187   phobos->SetTitle(""); \r
 188   phobos->SetMarkerColor(kBlue); \r
 189   phobos->SetMarkerStyle(28); \r
 190   phobos->SetMarkerSize(1.2); \r
 191   phobos->SetLineWidth(2); \r
 192   phobos->SetLineColor(kBlue); \r
 193   phobos->Draw("p"); \r
 194 \r
 195 \r
 196   //E877 -------------------------------------------------------------------------\r
 197   //float xxhh[4]={4.75,17.,130.,200.}; \r
 198   //float yyhh[4]={8.2,9.3,10.,9.4}; \r
 199   float xxhh[1]={4.75}; // sqrt(s)\r
 200   float yyhh[1]={0.019}; \r
 201   float exhh[1]={0.}; \r
 202   float exhh[1]={0.}; \r
 203   float eyhh[1]={0.}; \r
 204   float eyhh[1]={0.002}; \r
 205   e877 = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh,eyhh); \r
 206   e877->SetTitle(""); \r
 207   e877->SetMarkerColor(kBlue); \r
 208   e877->SetMarkerStyle(34); \r
 209   e877->SetMarkerSize(1.2); \r
 210   e877->SetLineWidth(2); \r
 211   e877->SetLineColor(kBlue); \r
 212   e877->Draw("p"); \r
 213   \r
 214   //E895 ------------------------------------------------------------------------\r
 215   //float xxhq[4]={4.75,17.,130.,200.}; \r
 216   //float yyhq[4]={9.1.,10.,9.8,8.4}; \r
 217   float xxhq[4]={2.68,3.32,3.83,4.24}; // sqrt(s)\r
 218   float yyhq[4]={-0.05,-0.005,0.01,0.015}; \r
 219   float exhq[4]={0.,0.,0.}; \r
 220   float eyhq[4]={0.004,0.003,0.003,0.004}; \r
 221   e895 = new TGraphAsymmErrors(4,xxhq,yyhq,exhq,exhq,eyhq,eyhq); \r
 222   e895->SetTitle(""); \r
 223   e895->SetMarkerColor(cE895); \r
 224   e895->SetMarkerStyle(kFullTriangleUp); \r
 225   e895->SetMarkerSize(1.2); \r
 226   e895->SetLineWidth(2); \r
 227   e895->SetLineColor(cE895); \r
 228   e895->Draw("p"); \r
 229   //Plastic Ball ----------------------------------------------------\r
 230   //float xxhq[4]={4.75,17.,130.,200.}; \r
 231   //float yyhq[4]={9.1.,10.,9.8,8.4}; \r
 232   float xxPB[6]={.15,.2,.4,.6,.8,1.2}; \r
 233   float yyPB[6]={-.035,-.08,-.09,-.06,-.04,-.03}; \r
 234   float exPB[6]={0.,0.,0.}; \r
 235   float eyPB[6]={0.,0.,0.}; \r
 236   plasticBall = new TGraphAsymmErrors(6,xxPB,yyPB,exPB,exPB,eyPB,eyPB); \r
 237   plasticBall->SetTitle(""); \r
 238   plasticBall->SetMarkerColor(kBlue); \r
 239   plasticBall->SetMarkerStyle(kFullTriangleDown); \r
 240   plasticBall->SetMarkerSize(1.2); \r
 241   plasticBall->SetLineWidth(1); \r
 242   plasticBall->SetLineColor(kBlue); \r
 243   //plasticBall->Draw("p");\r
 244 \r
 245 \r
 246   //Fopi------------------------------------------------------------\r
 247   float xxFopi[10]={0.09,0.12,0.15,0.25,0.4,0.6,0.8,1.0,1.2,1.49}; \r
 248   float yyFopi[10]={0.07456,0.02847,-0.00774,-0.05784,-0.08200,-0.07087,-0.06845,-0.06327,-0.05523,-0.04344}; \r
 249   float exFopi[10]={0.,0.,0.}; \r
 250   float eyFopi[10]={0.00746,0.00285,0.00077,0.00578,0.00656,0.00850,0.00684,0.00823,0.00828,0.00956}; \r
 251   for (int i=0;i<10;i++) xxFopi[i] = sqrt((xxFopi[i]+m_p)*m_p*2. +2.*m_p*m_p);\r
 252     \r
 253   fopi = new TGraphAsymmErrors(10,xxFopi,yyFopi,exFopi,exFopi,eyFopi,eyFopi); \r
 254   fopi->SetTitle(""); \r
 255   fopi->SetMarkerColor(cFOPI); \r
 256   fopi->SetMarkerStyle(kFullTriangleDown); \r
 257   fopi->SetMarkerSize(1.2); \r
 258   fopi->SetLineWidth(2); \r
 259   fopi->SetLineColor(cFOPI);\r
 260   fopi->Draw("p"); \r
 261 \r
 262 \r
 263   //EOS --------------------------------------------------------------------\r
 264   //float xxhq[4]={4.75,17.,130.,200.}; \r
 265   //float yyhq[4]={9.1.,10.,9.8,8.4}; \r
 266   float xxhq[4]={2.35}; \r
 267   float yyhq[4]={-0.065}; \r
 268   float exhq[4]={0.}; \r
 269   float eyhq[4]={0.007}; \r
 270   eos = new TGraphAsymmErrors(1,xxhq,yyhq,exhq,exhq,eyhq,eyhq); \r
 271   eos->SetTitle(""); \r
 272   eos->SetMarkerColor(kRed); \r
 273   eos->SetMarkerStyle(3); \r
 274   eos->SetMarkerSize(1.2); \r
 275   eos->SetLineWidth(2); \r
 276   eos->SetLineColor(kRed); \r
 277   eos->Draw("p"); \r
 278 \r
 279 \r
 280 \r
 281 // ALICE measured v2{4} + statistical error\r
 282   float v24ALICE = 0.073;\r
 283   float v24eALICE = 0.001;\r
 284   float corrpt = 0.88;\r
 285   float xxhh[1]={2760.}; // sqrt(s)\r
 286   float yyhh[1]={v24ALICE*corrpt}; \r
 287   float exhh[1]={0.};  \r
 288 \r
 289 \r
 290 // ALICE systematic error\r
 291 // hijing 15.3 %\r
 292 // therminator 6.1%\r
 293   float v24eALICE2 = 0.004;\r
 294   float eyhh2[1]={v24eALICE2};\r
 295 \r
 296   //draw with systematic error\r
 297   //alice = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh,eyhh); \r
 298   alice = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh2,eyhh2); \r
 299   alice->SetTitle(""); \r
 300   alice->SetMarkerColor(cAlice); \r
 301   alice->SetMarkerStyle(20); \r
 302   alice->SetMarkerSize(1.2); \r
 303   alice->SetLineWidth(2); \r
 304   alice->SetLineColor(cAlice); \r
 305   alice->Draw("p"); \r
 306   \r
 307   //------------------------------------------------------------------  \r
 308   TLegend *legend = new TLegend(.68,.24,.86,.66);\r
 309   // TLegend *legend = new TLegend(.45,.23,.62,.52);\r
 310   legend->Clear();\r
 311   legend->SetBorderSize(1);\r
 312   legend->SetFillColor(kWhite);\r
 313   legend->SetTextSize(0.028);\r
 314   //legend->AddEntry(e895,"E895","P");\r
 315 //   legend->AddEntry(e877,"E877","P");\r
 316 //   legend->AddEntry(data,"NA49 pions","P");\r
 317 //   legend->AddEntry(ceres,"CERES charged particles","P");\r
 318 //   legend->AddEntry(star,"STAR charged particles","P");\r
 319   //legend->SetHeader("A+A:");\r
 320   //legend->AddEntry(data,"NA49 std/mod","P");\r
 321   legend->AddEntry(alice,"ALICE","p");\r
 322   legend->AddEntry(star,"STAR","P");\r
 323   legend->AddEntry(phobos,"PHOBOS","P");\r
 324   legend->AddEntry(phenix,"PHENIX","P");\r
 325   legend->AddEntry(cumul,"NA49","P");\r
 326   legend->AddEntry(ceres,"CERES","P");\r
 327   legend->AddEntry(e877,"E877","P");\r
 328   legend->AddEntry(eos,"EOS","P");\r
 329   legend->AddEntry(e895,"E895","P");\r
 330   legend->AddEntry(fopi,"FOPI","P");\r
 331   legend->SetTextFont(22); // 22 = Times New Roman (bold)\r
 332 \r
 333   legend->Draw();\r
 334   //--------------------------- insert logo -----------\r
 335   TImage *ps = TImage::Open("macros/ALICElogo.png");\r
 336   float xlogo=0.75;\r
 337   float ylogo=0.76;\r
 338   float wlogo=0.08;\r
 339   TPad *aliceLogo = new TPad("aliceLogo", "aliceLogo",xlogo,ylogo,xlogo+wlogo,ylogo+wlogo);\r
 340   //  aliceLogo->Draw();\r
 341   aliceLogo->cd();\r
 342   //  ps->Draw("same"); \r
 343   //-------------------------- print -----------\r
 344   canvas->Print("v2edep.png");\r
 345   canvas->Print("v2edep.pdf");\r
 346   canvas->Print("v2edep.eps");\r
 347   \r
 348   \r
 349 } \r
 350  \r