1 void testAliGlauberQuenching() {
5 g.SetCentralityClass(0.00,0.10); // centrality (fraction of geometrical cross section)
7 AliQuenchingWeights afq;
12 afq.SampleEnergyLoss();
14 TCanvas *c = new TCanvas("cELD","Energy Loss Distribution",0,0,800,500);
17 for(Int_t itype=1;itype<=2;itype++){
23 sprintf(name,"Energy Loss Distribution - Quarks;E_{loss} (GeV);#");
24 sprintf(hname,"hQuarks");
26 sprintf(name,"Energy Loss Distribution - Gluons;E_{loss} (GeV);#");
27 sprintf(hname,"hGluons");
30 TH1F *h = new TH1F(hname,name,100,0,200);
32 for(Int_t i=0;i<10000;i++){
33 if(i % 100 == 0) cout << "." << flush;
36 Double_t loss=afq.GetELossRandom(itype,ell);
55 // second example using ell distribution
56 void testFastAliGlauberQuenching(Char_t *fname) {
59 TH1F *hEll=(TH1F*)f.Get("hEll");
62 AliQuenchingWeights afq;
64 afq.SetQTransport(1.);
67 afq.SampleEnergyLoss();
69 TCanvas *c = new TCanvas("cELD","Energy Loss Distribution",0,0,800,500);
72 for(Int_t itype=1;itype<=2;itype++){
77 sprintf(name,"Energy Loss Distribution - Quarks;E_{loss} (GeV);#");
79 sprintf(name,"Energy Loss Distribution - Gluons;E_{loss} (GeV);#");
82 sprintf(hname,"hQuarks");
84 sprintf(hname,"hGluons");
86 TH1F *h = afq.ComputeELossHisto(itype,1.,hEll);