3 TGeant3 *geant3 = (TGeant3*)gMC;
4 Float_t tkin[100], valPb[100], valW[100], val1[100], val2[100];
5 Float_t val3[100], val4[100], val5[100], pcut[5];
12 for (i=0; i< kdim; i++) {
13 tkin[i]=Float_t(i+1)*2.;
16 geant3->Gftmat( 4, ipart, chmeca, kdim, tkin, val1, pcut, ixst);
18 geant3->Gftmat(25, ipart, chmeca, kdim, tkin, val2, pcut, ixst);
20 geant3->Gftmat(28, ipart, chmeca, kdim, tkin, val3, pcut, ixst);
22 geant3->Gftmat(16, ipart, chmeca, kdim, tkin, val4, pcut, ixst);
24 geant3->Gftmat(13, ipart, chmeca, kdim, tkin, val5, pcut, ixst);
26 for (i=0; i< kdim; i++) {
27 valPb[i]=(225.*val1[i]+153.*val2[i]+15.*val3[i]+20*val4[i])/1000.;
30 for (i=0; i< kdim; i++) {
31 valW[i]=(225.*val1[i]+153.*val2[i]+35*val5[i])/1000.;
35 TGraph* dedx1 = new TGraph(kdim, tkin, valPb);
36 TGraph* dedx2 = new TGraph(kdim, tkin, valW);
38 TCanvas *c1=new TCanvas("c1","dedx",400,10,600,700);
39 dedx1->SetFillColor(42);
40 dedx1->SetMarkerColor(4);
41 dedx1->SetMarkerStyle(21);
43 dedx1->GetHistogram()->SetXTitle("Kinetic Energy (GeV)");
44 dedx1->GetHistogram()->SetYTitle("Mean Energy Loss (GeV)");
46 TCanvas *c2=new TCanvas("c2","dedx",400,10,600,700);
47 dedx2->SetFillColor(42);
48 dedx2->SetMarkerColor(4);
49 dedx2->SetMarkerStyle(21);
51 dedx2->GetHistogram()->SetXTitle("Kinetic Energy (GeV)");
52 dedx2->GetHistogram()->SetYTitle("Mean Energy Loss (GeV)");