11 gStyle->SetOptStat(0000);
13 c1 =
new TCanvas (
"c1",
"",20,20,1200,600);
19 FILE *
fp = fopen(
"Dend3DEdep.out",
"r");
29 h1 =
new TProfile(
"Energy",
"Energy deposits (keV) in dendritic compartments",1000,-1000,1000,0.001,1000);
30 h2 =
new TProfile(
"Dose",
"Dose deposits (Gy) in dendritic compartments",1000,-1000,1000,0.001,1000);
33 ncols = fscanf(fp,
" %f %f %f %f %f %f %f",&posX, &posY, &posZ, &distA, &distB, &EdepR, &DoseR);
35 if (distMaxA < distA ) distMaxA =
distA ;
36 if (distMaxB < distB ) distMaxB =
distB ;
37 if (edepMax < EdepR ) edepMax =
EdepR ;
38 if (doseMax < DoseR ) doseMax =
DoseR ;
39 if (edepMin > EdepR ) edepMin =
EdepR ;
40 if (doseMin > DoseR ) doseMin =
DoseR ;
42 h1->Fill(-distB, EdepR);
43 h1->Fill(distA, EdepR);
44 h2->Fill(-distB, DoseR);
45 h2->Fill(distA, DoseR);
49 cout <<
" Max and Min Energy deposits (keV) == " << edepMax <<
" ; "<< edepMin<<endl;
50 cout <<
" Max and Min Dose deposits (Gy) == " << doseMax <<
" ; "<< doseMin<<endl;
51 cout <<
" Maximum Basal Distance (um) == " << distMaxB <<
" "<<endl;
52 cout <<
" Maximum Apical Distance (um) == " << distMaxA <<
" "<<endl;
58 h1->SetMarkerColor(4);
59 h1->SetMarkerStyle(27);
60 h1->SetFillStyle(3005);
62 h1->GetXaxis()->SetTitle(
"Distance from Soma (um)");
63 h1->GetYaxis()->SetRangeUser(edepMin, edepMax+3.);
66 TLatex text(-distMaxB,edepMax-1.,
"Basal");
71 TLatex text(distMaxA/3.,edepMax-2.,
"Apical");
74 h1->GetXaxis()->SetRangeUser(-distMaxB-10., distMaxA+10.);
81 h2->SetMarkerColor(kRed);
82 h2->SetMarkerStyle(27);
83 h2->SetFillStyle(3005);
85 h2->GetXaxis()->SetTitle(
"Distance from Soma (um)");
86 h2->GetYaxis()->SetRangeUser(doseMin, doseMax+0.3);
89 TLatex text(-distMaxB,doseMax-0.1,
"Basal");
94 TLatex text(distMaxA/3.,doseMax-0.2,
"Apical");
97 h2->GetXaxis()->SetRangeUser(-distMaxB-10., distMaxA+10.);