plotDiscreteWeight(Int_t len=5)
{
  AliQuenchingWeights afq;
  afq.InitMult();
  afq.PlotDiscreteWeights();
}

plotDiscreteWeightMacro(Int_t len=5)
{
  AliQuenchingWeights afq;
  afq.InitMult();

  TCanvas *c = new TCanvas("cdisc","Discrete Weight for Multiple Scattering",0,0,500,400);
  TH2F *hframe = new TH2F("hdisc","",2,0,1.1,2,0,1);
  hframe->SetStats(0);
  hframe->SetXTitle("#hat{q} [GeV^{2}/fm]");
  hframe->SetYTitle("Probability #Delta E = 0 , p_{0}");
  hframe->Draw();

  TGraph *gq=new TGraph(20);
  Int_t i=0;
  for(Double_t q=0.1;q<=1.05;q+=0.05){
    Double_t disc,cont;
    afq.CalcMult(1,1.0, q, len, cont, disc);
    //cout << " " << q << " " << disc << endl;
    gq->SetPoint(i,q,disc);i++;
  }
  gq->SetMarkerStyle(20);
  gq->Draw("pl");

  TGraph *gg=new TGraph(20);
  Int_t i=0;
  for(Double_t q=0.05;q<=1.05;q+=0.05){
    Double_t disc,cont;
    afq.CalcMult(2,1.0, q, 5., cont, disc);
    //cout << " " << q << " " << disc << endl;
    gg->SetPoint(i,q,disc);i++;
  }
  gg->SetMarkerStyle(24);
  gg->Draw("pl");

  TLegend *l1a = new TLegend(0.5,0.6,.95,0.8);
  l1a->SetFillStyle(0);
  l1a->SetBorderSize(0);
  Char_t label[100];
  sprintf(label,"L = %d fm",len);
  l1a->AddEntry(gq,label,"");
  l1a->AddEntry(gq,"quark","pl");
  l1a->AddEntry(gg,"gluon","pl");
  l1a->Draw();

  c->Update();

}