00001 void getRdaPrel()
00002 {
00003
00004 ifstream frank("./datapoints/frank_pp05.dat");
00005 Float_t frx[100];
00006 Float_t fry[100];
00007 Float_t frex[100];
00008 Float_t frey[100];
00009 Int_t ifr=0;
00010 while(ifr<12){
00011 if(!frank.good()) break;
00012 frank>>frx[ifr]>>fry[ifr]>>frey[ifr];
00013 fry[ifr]*=0.85/26.1;
00014 frey[ifr]*=0.85/26.1;
00015 frex[ifr]=0.;
00016 ifr++;
00017 }
00018 TGraphErrors *frank_pp=new TGraphErrors(ifr,frx,fry,frex,frey);
00019 frank_pp->SetName("frank");
00020 frank_pp->SetMarkerStyle(25);
00021 frank_pp->SetMarkerColor(4);
00022
00023 TGraphErrors *andre_dau = new TGraphErrors(20);
00024 andre_dau->SetName("andre");
00025 andre_dau->SetFillColor(1);
00026 andre_dau->SetMarkerStyle(20);
00027 andre_dau->SetMarkerSize(1.4);
00028 andre_dau->SetPoint(1,1.75,0.00500784);
00029 andre_dau->SetPointError(1,0,0.000146705);
00030 andre_dau->SetPoint(2,2.25,0.00178272);
00031 andre_dau->SetPointError(2,0,5.55874e-05);
00032 andre_dau->SetPoint(3,2.75,0.000527415);
00033 andre_dau->SetPointError(3,0,2.38637e-05);
00034 andre_dau->SetPoint(4,3.25,0.000173036);
00035 andre_dau->SetPointError(4,0,1.17154e-05);
00036 andre_dau->SetPoint(5,3.75,7.47168e-05);
00037 andre_dau->SetPointError(5,0,1.44958e-05);
00038 andre_dau->SetPoint(6,4.25,2.65451e-05);
00039 andre_dau->SetPointError(6,0,4.05473e-06);
00040 andre_dau->SetPoint(7,4.75,1.0569e-05);
00041 andre_dau->SetPointError(7,0,2.9797e-06);
00042 andre_dau->SetPoint(8,5.25,4.08643e-06);
00043 andre_dau->SetPointError(8,0,1.23376e-06);
00044 andre_dau->SetPoint(9,6,1.29798e-06);
00045 andre_dau->SetPointError(9,0,4.30548e-07);
00046 andre_dau->SetPoint(10,6.5,7.02737e-07);
00047 andre_dau->SetPointError(10,0,1.00577e-07);
00048 andre_dau->SetPoint(11,7.5,2.8363e-07);
00049 andre_dau->SetPointError(11,0,5.20925e-08);
00050 andre_dau->SetPoint(12,8.5,1.22991e-07);
00051 andre_dau->SetPointError(12,0,2.57462e-08);
00052 andre_dau->SetPoint(13,9.5,5.68615e-08);
00053 andre_dau->SetPointError(13,0,6.45742e-09);
00054 andre_dau->SetPoint(14,10.5,2.77471e-08);
00055 andre_dau->SetPointError(14,0,3.97804e-09);
00056 andre_dau->SetPoint(15,11.5,1.61125e-08);
00057 andre_dau->SetPointError(15,0,4.34288e-09);
00058 andre_dau->SetPoint(16,12.5,9.93835e-09);
00059 andre_dau->SetPointError(16,0,2.42936e-09);
00060 andre_dau->SetPoint(17,13.5,5.44627e-09);
00061 andre_dau->SetPointError(17,0,3.26254e-09);
00062 andre_dau->SetPoint(18,14.5,3.15326e-09);
00063 andre_dau->SetPointError(18,0,2.1746e-09);
00064 andre_dau->SetPoint(19,15.5,1.56401e-09);
00065 andre_dau->SetPointError(19,0,1.13746e-09);
00066
00067 TMultiGraph *mgfit=new TMultiGraph();
00068 mgfit->Add(frank_pp);
00069 TF1 *fitQCD=new TF1("fitQCD","[0]*pow(1.+x,[1])*([2]+[3]*x+[4]*x*x+[5]*x*x*x)",1.5,16.);
00070 fitQCD->SetParameters(600.,-9.5,1.,0.,0.,0.);
00071 mgfit->Fit(fitQCD,"R0");
00072
00073 TMultiGraph *m=new TMultiGraph("m","none;p_{T} (GeV/c);1/(2#piNp_{T}) dN/dydp_{T}");
00074 gStyle->SetOptStat(0);
00075 gStyle->SetOptTitle(0);
00076 m->Add(andre_dau);
00077 m->Add(frank_pp);
00078 m->SetMaximum(.1);
00079 m->SetMinimum(1.e-11);
00080
00081 TCanvas *c=new TCanvas("c","c",400,400);
00082 gPad->SetLogy();
00083 m->Draw("ap");
00084 fitQCD->SetLineWidth(1);
00085 fitQCD->Draw("same");
00086
00087 TLegend *legg=new TLegend(.4,.6,.8,.8);
00088 legg->AddEntry(andre_dau,"d+Au preliminary","p");
00089 legg->AddEntry(frank_pp,"p+p preliminary","p");
00090 legg->Draw("same");
00091
00092 c->SaveAs("rda_prel.pdf");
00093
00094 TCanvas *c2=new TCanvas("c2","c2",400,150);
00095 TGraphErrors *andre_dau2=new TGraphErrors();
00096 andre_dau2->SetMarkerStyle(20);
00097 andre_dau2->SetMarkerSize(1.4);
00098 TGraphErrors *frank_pp2=new TGraphErrors();
00099
00100 for(int i=0;i<19;i++){
00101 double xx=0.;
00102 double yy=0;
00103 double exx=0.;
00104 double eyy=0.;
00105 if(i<12){
00106 frank_pp->GetPoint(i,xx,yy);
00107 eyy=frank_pp->GetErrorY(i);
00108 frank_pp2->SetPoint(i,xx,yy/fitQCD->Eval(xx));
00109 frank_pp2->SetPointError(i,exx,eyy/fitQCD->Eval(xx));
00110 }
00111 andre_dau->GetPoint(i,xx,yy);
00112 eyy=andre_dau->GetErrorY(i);
00113 andre_dau2->SetPoint(i,xx,yy/(7.5*fitQCD->Eval(xx)));
00114 andre_dau2->SetPointError(i,exx,eyy/(7.5*fitQCD->Eval(xx)));
00115 }
00116 frank_pp2->SetMarkerStyle(25);
00117 frank_pp2->SetMarkerColor(4);
00118
00119
00120
00121 TMultiGraph *m_rda=new TMultiGraph();
00122 m_rda->Add(frank_pp2);
00123 m_rda->Add(andre_dau2);
00124 m_rda->SetMinimum(0.01);
00125 m_rda->SetMaximum(3.);
00126 m_rda->Draw("ap");
00127
00128 c2->SaveAs("rda_prel2.pdf");
00129 }
