Gabriele Cosmo
41 lines
983 B
C++
41 lines
983 B
C++
|
|
{
|
|
gROOT->Reset();
|
|
|
|
// Draw histogram fill by Geant4 TestBruce simulation
|
|
TFile f("./Au3.local.root");
|
|
TH1D* h1d = (TH1D*) f.Get("4");
|
|
h1d->SetTitle("Fluence distribution of 13 MeV e- in Au (93.7 mg/cm2)");
|
|
h1d->GetXaxis()->SetTitle("r (mm)");
|
|
h1d->GetYaxis()->SetTitle("Fluence");
|
|
h1d->SetStats(kFALSE); // Eliminate statistics box
|
|
h1d->Draw("HIST");
|
|
|
|
/* data
|
|
* Bruce et al.
|
|
*/
|
|
|
|
ifstream in;
|
|
in.open("../data/Au3.13MeV.ascii");
|
|
|
|
TMarker *pt;
|
|
Double_t x, y;
|
|
// First indicate number of data
|
|
int nbdata = 0;
|
|
in >> nbdata;
|
|
for ( int i = 0 ; i < nbdata ; i++ ) {
|
|
in >> x >> y ;
|
|
if (!in.good()) break;
|
|
pt = new TMarker(x,y,22); // 22 for triangle TMatker
|
|
pt->SetMarkerColor(kRed);
|
|
pt->Draw();
|
|
}
|
|
in.close();
|
|
|
|
// Print the histograms legend
|
|
TLegend* legend = new TLegend(0.6,0.55,0.8,0.68);
|
|
legend->AddEntry(h1d,"msc95","l");
|
|
legend->AddEntry(pt,"Faddegon data","P");
|
|
legend->Draw();
|
|
}
|