Test routines and example for use of AliFastGlauber and AliQuenchingWeights

[u/mrichter/AliRoot.git] / FASTSIM / testAliQuenching.C

diff --git a/FASTSIM/testAliQuenching.C b/FASTSIM/testAliQuenching.C
new file mode 100644 (file)
index 0000000..daee853
--- /dev/null
+++ b/FASTSIM/testAliQuenching.C
@@ -0,0 +1,24 @@
+testAliQuenching()
+{
+  AliQuenchingWeights afq;
+  afq.InitMult();
+  afq.SetQTransport(1.);
+  //afq.SetECMethod(0);
+  afq.SampleEnergyLoss();
+
+  Int_t lengths[6]={2,4,6,8,10,12};
+  TCanvas *c = new TCanvas("cELW","Energy Loss Weights",0,0,800,500);
+  c->Divide(3,2);
+  for(Int_t len=0;len<6;len++){
+    c->cd(len+1);
+    TH1F *h = new TH1F(*afq.GetHisto(2,lengths[len]));
+    if(h) h->DrawCopy();
+  }
+  c->Update();
+
+  afq.PlotDiscreteWeights(4);
+  afq.PlotAvgELoss(4);
+  afq.PlotAvgELossVsPt(1.,4);
+}
+
+