Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/berger/berger.C

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 {
  2    gROOT->Reset();
  3 
  4    // Draw histograms fill by Geant4 TestEm5 simulation
  5    TFile f1("./local.root");
  6    TH1D* h1 = (TH1D*) f1.Get("h1");
  7    h1->SetTitle("Energy deposited in 530 um of Si by 1 MeV e-");
  8    h1->GetXaxis()->SetTitle("Ekine (MeV)");
  9    h1->GetYaxis()->SetTitle("nb / MeV");
 10    h1->SetStats(kFALSE);  // Eliminate statistics box
 11    h1->SetLineColor(kBlue);
 12    h1->Draw("HIST");
 13    
 14    TFile f2("./opt4.root");
 15    TH1D* h2 = (TH1D*) f2.Get("h1");
 16    h2->SetStats(kFALSE);  // Eliminate statistics box
 17    h2->SetLineColor(kRed);
 18    h2->Draw("SAME HIST");
 19 /*
 20 * e- 1 MeV in Silicon 530 um 
 21 * M.J. Berger et al. NIM 69 (p.181) 1969
 22 * distribution of energy deposition
 23 * (from 110 keV to 1.03 MeV by bin of 10 keV --> 93 bins) 
 24 */
 25 
 26   ifstream in;
 27   in.open("530um.ascii");
 28    
 29   TMarker *pt;
 30   Double_t x, y;   
 31   // First indicate number of data
 32   int nbdata = 0;
 33   in >> nbdata;
 34   for ( int i = 0 ; i < nbdata ; i++ ) {
 35       in >> x >> y ;
 36       if (!in.good()) break;
 37       pt = new TMarker(x,y,32); // 32 for open triangle-down
 38       pt->SetMarkerColor(kGreen);
 39       pt->Draw();
 40   }   
 41   in.close();
 42 
 43   // Print the histograms legend
 44   TLegend *legend = new TLegend(0.6,0.6,0.8,0.8);
 45   legend->AddEntry(h1,"local","l");
 46   legend->AddEntry(h2,"opt4 ","l");   
 47   legend->AddEntry(pt,"Berger data","P");
 48   legend->Draw();
 49 }
 50