14 gStyle->SetPalette(1);
15 gROOT->SetStyle(
"Plain");
16 gStyle->SetOptStat(000000);
18 c1 =
new TCanvas (
"c1",
"",60,60,800,800);
21 FILE *
fp = fopen(
"s.txt",
"r");
27 TNtuple *
ntuple =
new TNtuple(
"ntuple",
"s",
"radius:thick:E:sc:ssc:sn:ssn");
30 ncols = fscanf(fp,
"%f %f %f %f %f %f %f",&radius,&thick,&E,&sc,&ssc,&sn,&ssn);
32 ntuple->Fill(radius,thick,E,sc,ssc,sn,ssn);
41 TH2F *
h2 =
new TH2F (
"h2",
"",2,99.999,1e4,2,1e0,1e3);
44 ntuple->SetMarkerStyle(24);
45 ntuple->SetMarkerSize(1.);
46 ntuple->Draw(
"sc:E",
"",
"LPsame");
48 ntuple->SetMarkerStyle(20);
49 ntuple->SetMarkerSize(1.);
50 ntuple->Draw(
"sn:E",
"",
"LPsame");
52 h2->GetXaxis()->SetLabelSize(0.025);
53 h2->GetYaxis()->SetLabelSize(0.025);
54 h2->GetXaxis()->SetTitleSize(0.035);
55 h2->GetYaxis()->SetTitleSize(0.035);
56 h2->GetXaxis()->SetTitleOffset(1.4);
57 h2->GetYaxis()->SetTitleOffset(1.4);
58 h2->GetXaxis()->SetTitle(
"E (eV)");
59 h2->GetYaxis()->SetTitle(
"S (Gy/Bq.s)");
61 TText *
pt1 =
new TText(200,200,
"Cytoplasm");
64 TText *
pt2 =
new TText(400,15,
"Nucleus");