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

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);
}
Commit	Line	Data
ffb3506e	1	testAliQuenching()
	2	{
	3	AliQuenchingWeights afq;
	4	afq.InitMult();
	5	afq.SetQTransport(1.);
	6	//afq.SetECMethod(0);
	7	afq.SampleEnergyLoss();
	8
	9	Int_t lengths[6]={2,4,6,8,10,12};
	10	TCanvas *c = new TCanvas("cELW","Energy Loss Weights",0,0,800,500);
	11	c->Divide(3,2);
	12	for(Int_t len=0;len<6;len++){
	13	c->cd(len+1);
	14	TH1F h = new TH1F(afq.GetHisto(2,lengths[len]));
	15	if(h) h->DrawCopy();
	16	}
	17	c->Update();
	18
	19	afq.PlotDiscreteWeights(4);
	20	afq.PlotAvgELoss(4);
	21	afq.PlotAvgELossVsPt(1.,4);
	22	}
	23
	24