Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
egs.C
이 파일의 문서화 페이지로 가기
1 
2 {
3  gROOT->Reset();
4 
5  // Draw histograms fill by Geant4 TestEm11 simulation
6  //
8  TFile f1("./Ta_1000keV_opt3.root");
9  TH1D* h1 = (TH1D*) f1.Get("8");
10  h1->SetTitle("Depth dose distribution of 1000 keV e- in Ta");
11  h1->GetXaxis()->SetTitle("Edep (Mev.cm2/g) along x/r0 x/r0");
12  h1->GetYaxis()->SetTitle("MeV*cm2/g");
13  h1->SetStats(kFALSE); // Eliminate statistics box
14  h1->SetLineColor(kBlack);
15  h1->Draw("HIST");
16 
17 /* EGSnrc
18 * Yann Perrot
19 */
20 
21  ifstream in;
23  in.open("./EGSnrc/Ta_1000keV_EGSnrc.ascii");
24 
25  TMarker *pt;
27  // First indicate number of data
28  int nbdata = 0;
29  in >> nbdata;
30  for ( int i = 0 ; i < nbdata ; i++ ) {
31  in >> x >> y ;
32  if (!in.good()) break;
33  pt = new TMarker(x,y,22); // 22 for triangle TMatker
34  pt->SetMarkerColor(kRed);
35  pt->Draw();
36  }
37  in.close();
38 
39  // Print the histograms legend
40  TLegend* legend = new TLegend(0.6,0.5,0.8,0.70);
41  legend->AddEntry(h1,"ref10 ","l");
42  legend->AddEntry(pt,"EGSnrc","P");
43  legend->Draw();
44 }
Float_t x
Definition: compare.C:6
G4double Double_t
int nbdata
Definition: egs.C:28
Float_t y
Definition: compare.C:6
TMarker * pt
Definition: egs.C:25
Float_t f1
TLegend * legend
Definition: egs.C:40
ifstream in
Definition: comparison.C:7
TH1F * h1