PWGCF/EBYE/Fluctuations/macro/drawEnergyDependence.C

   1 void drawEnergyDependence() {
   2   //Macro to draw the energy dependence of the nu_dyn
   3
   4   //================================ALICE================================//
   5   Double_t gNuDynALICEData[1] = {-0.000991451};
   6   Double_t gNuDynALICEDataStatError[1] = {3.06304e-05};
   7
   8   Double_t gALICEEnergy[1] = {2760.};
   9   Double_t gALICEEnergyError[1] = {10.};
  10
  11   TGraphErrors *grNuDynALICEData = new TGraphErrors(1,gALICEEnergy,
  12                                                     gNuDynALICEData,
  13                                                     gALICEEnergyError,
  14                                                     gNuDynALICEDataStatError);
  15   grNuDynALICEData->SetMarkerStyle(20);
  16   grNuDynALICEData->SetMarkerColor(2);
  17   grNuDynALICEData->SetMarkerSize(1.4);
  18
  19   //================================CERES================================//
  20   Double_t gNuDynCERESData[3] = {-0.00006,-0.00069,-0.00108};
  21   Double_t gNuDynCERESDataStatError[3] = {5.0e-05,5.0e-05,5.0e-05};
  22
  23   Double_t gCERESEnergy[3] = {7.8,12.3,17.3};
  24   Double_t gCERESEnergyError[3] = {0.5,0.5,0.5};
  25
  26   TGraphErrors *grNuDynCERESData = new TGraphErrors(3,gCERESEnergy,
  27                                                     gNuDynCERESData,
  28                                                     gCERESEnergyError,
  29                                                     gNuDynCERESDataStatError);
  30   grNuDynCERESData->SetMarkerStyle(23);
  31   grNuDynCERESData->SetMarkerColor(1);
  32   grNuDynCERESData->SetMarkerSize(1.4);
  33
  34   //================================STAR================================//
  35   Double_t gNuDynSTARData[4] = {-0.00113,-0.00147,-0.00122,-0.00163};
  36   Double_t gNuDynSTARDataStatError[4] = {0.00026,0.00018,0.00014,7e-05};
  37
  38   Double_t gSTAREnergy[4] = {19.6,62.4,130.,200.};
  39   Double_t gSTAREnergyError[4] = {0.5,0.5,0.5,0.5};
  40
  41   TGraphErrors *grNuDynSTARData = new TGraphErrors(4,gSTAREnergy,
  42                                                    gNuDynSTARData,
  43                                                    gSTAREnergyError,
  44                                                    gNuDynSTARDataStatError);
  45   grNuDynSTARData->SetMarkerStyle(29);
  46   grNuDynSTARData->SetMarkerColor(3);
  47   grNuDynSTARData->SetMarkerSize(1.6);
  48
  49   //================================PHENIX================================//
  50   Double_t gNuDynPHENIXData[1] = {-0.00086};
  51   Double_t gNuDynPHENIXDataStatError[1] = {0.00038};
  52
  53   Double_t gPHENIXEnergy[1] = {150.};
  54   Double_t gPHENIXEnergyError[4] = {0.5};
  55
  56   TGraphErrors *grNuDynPHENIXData = new TGraphErrors(1,gPHENIXEnergy,
  57                                                    gNuDynPHENIXData,
  58                                                    gPHENIXEnergyError,
  59                                                    gNuDynPHENIXDataStatError);
  60   grNuDynPHENIXData->SetMarkerStyle(21);
  61   grNuDynPHENIXData->SetMarkerColor(4);
  62   grNuDynPHENIXData->SetMarkerSize(1.2);
  63
  64   //_____________________________________________________//
  65   //Draw the results
  66   //_____________________________________________________//
  67   TLatex *latex = new TLatex();
  68   latex->SetTextSize(0.035);
  69   //====================================//
  70   //results vs centrality
  71   TH2F *gEmpty1 = new TH2F("gEmpty1",
  72                            ";#sqrt{s_{NN}} (GeV);",
  73                            100,1,10000,10000,-0.003,0.5e-03);
  74   gEmpty1->SetStats(kFALSE);
  75   gEmpty1->GetYaxis()->SetTitleOffset(1.5);
  76   gEmpty1->GetXaxis()->SetTitleOffset(1.5);
  77   gEmpty1->GetYaxis()->SetNdivisions(10);
  78   gEmpty1->GetXaxis()->SetNdivisions(10);
  79
  80   TF1 *f1 = new TF1("f1","0",0,10000);
  81   f1->SetLineColor(1); f1->SetLineStyle(3); f1->SetLineWidth(3);
  82
  83   //============================================================//
  84   //nu_{dyn.}
  85   TCanvas *c1 = new TCanvas("c1","Energy dependence: nu_dyn",
  86                             0,0,500,500);
  87   c1->SetFillColor(10); c1->SetHighLightColor(10);
  88   c1->SetLeftMargin(0.15); c1->SetBottomMargin(0.15);
  89   c1->SetGridx(); c1->SetGridy(); c1->SetLogx();
  90   gEmpty1->GetYaxis()->SetTitle("#nu_{dyn.}");
  91   gEmpty1->DrawCopy();
  92   f1->Draw("same");
  93   grNuDynALICEData->Draw("P");
  94   grNuDynCERESData->Draw("P");
  95   grNuDynSTARData->Draw("P");
  96   grNuDynPHENIXData->Draw("P");
  97
  98   DrawMarker(7., -0.0019, 20, 1.6, 2);
  99   DrawMarker(7., -0.0021, 29, 1.6, 3);
 100   DrawMarker(7., -0.0023, 21, 1.6, 4);
 101   DrawMarker(7., -0.0025, 23, 1.6, 1);
 102   latex->DrawLatex(10.,-0.00195,"ALICE");
 103   latex->DrawLatex(10.,-0.00215,"STAR");
 104   latex->DrawLatex(10.,-0.00235,"PHENIX");
 105   latex->DrawLatex(10.,-0.00255,"CERES");
 106
 107   c1->SaveAs("energyDependenceNuDyn.png");
 108 }
 109
 110 //_______________________________________________________________//
 111 void DrawMarker(Double_t x, Double_t y, Int_t style,
 112                 Double_t size, Int_t color) {
 113   TMarker *m = new TMarker(x,y,style);
 114   m->SetMarkerSize(size);
 115   m->SetMarkerColor(color);
 116   m->Draw();
 117 }