TestEm11/sandia/egs.C

0001
0002 {
0003    gROOT->Reset();
0004
0005    // Draw histograms fill by Geant4 TestEm11 simulation
0006    //
0007    ///TFile f1("./Al_1033keV_opt3.root");
0008    TFile f1("./Ta_1000keV_opt3.root");
0009    TH1D* h1 = (TH1D*) f1.Get("8");
0010    h1->SetTitle("Depth dose distribution of 1000 keV e- in Ta");
0011    h1->GetXaxis()->SetTitle("Edep (Mev.cm2/g) along x/r0                    x/r0");
0012    h1->GetYaxis()->SetTitle("MeV*cm2/g");
0013    h1->SetStats(kFALSE);  // Eliminate statistics box
0014    h1->SetLineColor(kBlack);
0015    h1->Draw("HIST");
0016
0017 /* EGSnrc
0018 * Yann Perrot
0019 */
0020
0021    ifstream in;
0022    ///in.open("./EGSnrc/Al_1033keV_EGSnrc.ascii");
0023    in.open("./EGSnrc/Ta_1000keV_EGSnrc.ascii");
0024
0025    TMarker *pt;
0026    Double_t x, y;
0027    // First indicate number of data
0028    int nbdata = 0;
0029    in >> nbdata;
0030    for ( int i = 0 ; i < nbdata ; i++ ) {
0031       in >> x >> y ;
0032       if (!in.good()) break;
0033       pt = new TMarker(x,y,22); // 22 for triangle TMatker
0034       pt->SetMarkerColor(kRed);
0035       pt->Draw();
0036    }
0037    in.close();
0038
0039    // Print the histograms legend
0040    TLegend* legend = new TLegend(0.6,0.5,0.8,0.70);
0041    legend->AddEntry(h1,"ref10 ","l");
0042    legend->AddEntry(pt,"EGSnrc","P");
0043    legend->Draw();
0044 }