TruthGeneration/Vitev/plot_bjet_raa.C

0001 /*

0002 Here are preliminary results for the RAA of b-jets

0003 I have used the parameterization of Cesar (had convert it to d^2pT)

0004 I have used 10% g and 90 % b quarks

0005  3 radii  0.2, 0.4 and 0.6  - there are names 2, 4, 6 in the name of files

0006 two different couplings g= 2  (labeled g41.0)  and g=2.2 labeled g41.5

0007 I suggest using g=2.2 as default

0008 since ichi is exponential parameterizations and things will start falling steeply,

0009 I sugest not using results beyond  60 GeV

0010 This should give you a good sense of the expected suppression and p-T dependence

0011 as a function of radius and coupling between the jet and the media,

0012 Ivan Vitev, May 12, 2017

0013  */
0014 void plot_bjet_raa()
0015 {
0016   TGraph* gdef = new TGraph("R_SigR4.W0.0g41.5.SPEHIXdat");
0017   TGraph* ggeq2 = new TGraph("R_SigR4.W0.0g41.0.SPEHIXdat");
0018   gdef->SetLineColor(kBlack);
0019   gdef->SetLineWidth(4);
0020   ggeq2->SetLineWidth(4);
0021   ggeq2->SetLineColor(kBlue);
0022   TCanvas *c1 = new TCanvas("c1", "c1",10,45,700,502);
0023   gStyle->SetOptTitle(1);
0024   c1->SetHighLightColor(0);
0025   c1->Range(-9.413793,-0.2360656,62.7931,1.321311);
0026   c1->SetFillColor(0);
0027   c1->SetBorderMode(0);
0028   c1->SetBorderSize(0);
0029   c1->SetLeftMargin(0.1303725);
0030   c1->SetRightMargin(0.03868195);
0031   c1->SetTopMargin(0.07789474);
0032   c1->SetBottomMargin(0.1515789);
0033   c1->SetFrameBorderMode(0);
0034
0035   TH2F* hframe = new TH2F("hframe","Jet radii = 0.4;p_{T} [GeV/c];R_{AA}",
0036               100,0,60,100,0,1.2);
0037   hframe->GetXaxis()->SetLabelFont(132);
0038   hframe->GetXaxis()->SetLabelSize(0.06);
0039   hframe->GetXaxis()->SetTitleSize(0.07);
0040   hframe->GetXaxis()->SetTitleOffset(0.9);
0041   hframe->GetXaxis()->SetTitleFont(132);
0042   hframe->GetYaxis()->SetTitle("R_{AA}");
0043   hframe->GetYaxis()->SetNdivisions(310);
0044   hframe->GetYaxis()->SetLabelFont(132);
0045   hframe->GetYaxis()->SetLabelSize(0.06);
0046   hframe->GetYaxis()->SetTitleSize(0.08);
0047   hframe->GetYaxis()->SetTitleOffset(0.65);
0048   hframe->GetYaxis()->SetTitleFont(132);
0049   hframe->SetTitleFont(132);
0050   hframe->SetTitleSize(0.06315789);
0051   hframe->Draw();
0052
0053   gdef->Draw("L");
0054   ggeq2->Draw("L");
0055
0056   TLegend* l = new TLegend(0.3323782,0.6221053,0.6332378,0.9231579,
0057                "Au+Au#rightarrow b-jet @ #sqrt{s_{NN}}=200 GeV");
0058   l->SetBorderSize(0);
0059   l->SetFillStyle(0);
0060   l->SetTextFont(132);
0061   l->SetTextSize(0.06315789);
0062   l->AddEntry(ggeq2,"coupling g=2.0","L");
0063   l->AddEntry(gdef,"coupling g=2.2","L");
0064   l->Draw();
0065
0066 }