elosstest.cxx

{
	gSystem->Load("/media/commondata/work/Eclipse/TELoss/Debug/libTELoss.so");

//	TELoss	target;
	TELoss	mSi;
	TELoss	mSi3;

//	target.SetEL(1, 0.0014625);
//	target.AddEL(1., 3., 1.);
//	target.SetZP(1., 1.);
//	target.SetEtab(1000, 0., 50.);

	mSi.SetEL(1, 2.321);
	mSi.AddEL(14., 28.086, 1);
//	mSi.SetZP(1., 1.);		//protony
	mSi.SetZP(2., 4.);		//alfy
	mSi.SetEtab(100000, 200.);	//deltaE
	mSi.SetDeltaEtab(300);

	mSi3.SetEL(1, 2.321);
	mSi3.AddEL(14., 28.086, 1);
//	mSi3.SetZP(1., 1.);		//protony
	mSi3.SetZP(2., 3.);		//3He
	mSi3.SetEtab(100000, 200.);	//deltaE
	mSi3.SetDeltaEtab(300);

//	TCanvas *c1 = new TCanvas("jmeno", "titulek", 617, 0, 1058, 972);

//	mSi.PrintdEEtoFile("out.txt");
//	TGraph *gr1 = new TGraph("out.txt");
//	gr1->Draw("A*");
//
//	mSi3.PrintdEEtoFile("out3.txt");
//	TGraph *gr2 = new TGraph("out3.txt");
//	gr2->SetMarkerColor(2);
//	gr2->Draw("*");


//	mSi.PrintREtoFile();
//	TGraph *gr2 = new TGraph("outputRE.txt");
//	gr2->Draw("A*");

//	char out[20], gr[10];
//	TObjArray glist(0);
//	TMultiGraph *mg = new TMultiGraph();
//	TGraph *g;

//	for (Int_t i = 0; i < 30; i++) {
//		sprintf(out, "output%d.txt", i);
//		cout << out << endl;
//		mSi.SetDeltaEtab(300+i);
//		mSi.PrintdEEtoFile(out);
//		g = new TGraph(out);
//		glist.Add(g);
		
//		cout << out << endl;
//	}
	


//	cout << mSi.GetE(7.686, 10) << endl;
//	cout << mSi.GetE(7.686, 1) << endl;

//	Double_t finalEnergy = -1.;

//	Double_t Energy = steelwin.GetE0(finalEnergy, 30.);
//	Energy = target.GetE0(Energy, 2000.);

//	printf("%f\n", Energy);

}