Do not include from subdirectories
[u/mrichter/AliRoot.git] / FASTSIM / plotDiscreteWeight.C
1 plotDiscreteWeight(Int_t len=5)
2 {
3   AliQuenchingWeights afq;
4   afq.InitMult();
5   afq.PlotDiscreteWeights();
6 }
7
8 plotDiscreteWeightMacro(Int_t len=5)
9 {
10   AliQuenchingWeights afq;
11   afq.InitMult();
12
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);
15   hframe->SetStats(0);
16   hframe->SetXTitle("#hat{q} [GeV^{2}/fm]");
17   hframe->SetYTitle("Probability #Delta E = 0 , p_{0}");
18   hframe->Draw();
19
20   TGraph *gq=new TGraph(20);
21   Int_t i=0;
22   for(Double_t q=0.1;q<=1.05;q+=0.05){
23     Double_t disc,cont;
24     afq.CalcMult(1,1.0, q, len, cont, disc);
25     //cout << " " << q << " " << disc << endl;
26     gq->SetPoint(i,q,disc);i++;
27   }
28   gq->SetMarkerStyle(20);
29   gq->Draw("pl");
30
31   TGraph *gg=new TGraph(20);
32   Int_t i=0;
33   for(Double_t q=0.05;q<=1.05;q+=0.05){
34     Double_t disc,cont;
35     afq.CalcMult(2,1.0, q, 5., cont, disc);
36     //cout << " " << q << " " << disc << endl;
37     gg->SetPoint(i,q,disc);i++;
38   }
39   gg->SetMarkerStyle(24);
40   gg->Draw("pl");
41
42   TLegend *l1a = new TLegend(0.5,0.6,.95,0.8);
43   l1a->SetFillStyle(0);
44   l1a->SetBorderSize(0);
45   Char_t label[100];
46   sprintf(label,"L = %d fm",len);
47   l1a->AddEntry(gq,label,"");
48   l1a->AddEntry(gq,"quark","pl");
49   l1a->AddEntry(gg,"gluon","pl");
50   l1a->Draw();
51
52   c->Update();
53
54 }