Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/hanson/hanson.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 {
  3    gROOT->Reset();
  4 
  5    // Draw histograms fill by Geant4 TestEm11 simulation
  6    TFile f1("./19um-local.root");
  7    TH1D* h1 = (TH1D*) f1.Get("h12");
  8    h1->SetTitle("Angular distribution of 15.7 MeV e- after 19um Au foil");
  9    //h1->SetTitle("Angular distribution of 15.7 MeV e- after 9um Au foil");   
 10    h1->GetXaxis()->SetTitle("theta (deg)");
 11    h1->GetYaxis()->SetTitle("dN/dOmega");
 12    h1->SetStats(kFALSE);  // Eliminate statistics box
 13    h1->SetLineColor(kOrange);
 14    h1->Draw("HIST");
 15    
 16    TFile f2("./19um-opt3.root");
 17    TH1D* h2 = (TH1D*) f2.Get("h12");
 18    h2->SetStats(kFALSE);  // Eliminate statistics box
 19    h2->SetLineColor(kBlue);
 20    h2->Draw("SAME HIST");
 21    
 22    TFile f3("./19um-opt4.root");
 23    TH1D* h3 = (TH1D*) f3.Get("h12");
 24    h3->SetStats(kFALSE);  // Eliminate statistics box
 25    h3->SetLineColor(kGreen);
 26    h3->Draw("SAME HIST"); 
 27      
 28 /* data
 29 * angle distribution of  15.7 MeV electrons
 30 * transmitted through thin gold foils.
 31 * A.O.Hanson et al. Phys.Rev.84 (1951) page 634.
 32 */
 33 
 34    ifstream in;
 35    in.open("./19um-wwols.ascii");
 36    //in.open("./9um-wwols.ascii");
 37    
 38    TMarker *pt;
 39    Double_t x, y;
 40    // First indicate number of data
 41    int nbdata = 0;
 42    in >> nbdata;
 43    for ( int i = 0 ; i < nbdata ; i++ ) {
 44       in >> x >> y ;
 45       if (!in.good()) break;
 46       pt = new TMarker(x,y,23); // 23 for triangle-down TMarker
 47       pt->SetMarkerColor(kMagenta);
 48       pt->Draw();
 49    }
 50    in.close();
 51 
 52    // Print the histograms legend
 53    TLegend* legend = new TLegend(0.6,0.5,0.8,0.68);
 54    legend->AddEntry(h1,"local ","l");
 55    legend->AddEntry(h2,"opt3 ","l");
 56    legend->AddEntry(h3,"opt4","l");      
 57    legend->AddEntry(pt,"Hanson data","P");
 58    legend->Draw();
 59 }
 60