MUON/dedx.C

   1 void dedx(Int_t imate)
   2 {
   3     TGeant3 *geant3 = (TGeant3*)gMC;
   4     Float_t   tkin[100], value[100], pcut[5];
   5     Int_t ipart=5;
   6     char chmeca[4];
   7     strcpy(chmeca,"LOSS");
   8     Int_t ixst, i;
   9     Int_t kdim=100;
  10
  11     for (i=0; i< kdim; i++) {
  12         tkin[i]=Float_t(i)*1.+1;
  13     }
  14     geant3->Gftmat(imate, ipart, chmeca, kdim, tkin, value, pcut, ixst);
  15     for (i=0; i< kdim; i++) {
  16         printf("\n Energy %f dE/dx %f", tkin[i], value[i]);
  17     }
  18     TGraph*  dedx = new TGraph(kdim, tkin, value);
  19     TCanvas *c1=new TCanvas("c1","dedx",400,10,600,700);
  20     dedx->SetFillColor(42);
  21     dedx->SetMarkerColor(4);
  22     dedx->SetMarkerStyle(21);
  23     dedx->Draw("AC");
  24     dedx->GetHistogram()->SetXTitle("Kinetic Energy (GeV)");
  25     dedx->GetHistogram()->SetYTitle("dE/dx   ");
  26  }
  27
  28
  29
  30
  31
  32
  33
  34
  35
  36