1 plotDiscreteWeight(Int_t len=5)
3 AliQuenchingWeights afq;
5 afq.PlotDiscreteWeights();
8 plotDiscreteWeightMacro(Int_t len=5)
10 AliQuenchingWeights afq;
13 TCanvas *c = new TCanvas("cdisc","Discrete Weight for Multiple Scattering",0,0,500,400);
14 TH2F *hframe = new TH2F("hdisc","",2,0,1.1,2,0,1);
16 hframe->SetXTitle("#hat{q} [GeV^{2}/fm]");
17 hframe->SetYTitle("Probability #Delta E = 0 , p_{0}");
20 TGraph *gq=new TGraph(20);
22 for(Double_t q=0.1;q<=1.05;q+=0.05){
24 afq.CalcMult(1,1.0, q, len, cont, disc);
25 //cout << " " << q << " " << disc << endl;
26 gq->SetPoint(i,q,disc);i++;
28 gq->SetMarkerStyle(20);
31 TGraph *gg=new TGraph(20);
33 for(Double_t q=0.05;q<=1.05;q+=0.05){
35 afq.CalcMult(2,1.0, q, 5., cont, disc);
36 //cout << " " << q << " " << disc << endl;
37 gg->SetPoint(i,q,disc);i++;
39 gg->SetMarkerStyle(24);
42 TLegend *l1a = new TLegend(0.5,0.6,.95,0.8);
44 l1a->SetBorderSize(0);
46 sprintf(label,"L = %d fm",len);
47 l1a->AddEntry(gq,label,"");
48 l1a->AddEntry(gq,"quark","pl");
49 l1a->AddEntry(gg,"gluon","pl");