| Geant4 Cross Reference |
1 // ------------------------------------------- 1
2 // -------------------------------------------
3 //
4 // *******************************************
5 // To execute this macro under ROOT,
6 // launch ROOT (usually type 'root' at your ma
7 // This macro needs Dend3DEdep.out file :
8 // *******************************************
9 {
10 gROOT->Reset();
11 gStyle->SetOptStat(0000);
12
13 c1 = new TCanvas ("c1","",20,20,1200,600);
14 c1->Divide(2,1);
15
16 Int_t ncols=0;
17 Int_t nlines = 0;
18
19 FILE * fp = fopen("Dend3DEdep.out","r");
20 Float_t posX, posY, posZ ;
21 Float_t distB, distA, EdepR, DoseR;
22 Float_t distMaxA = -1e-9;
23 Float_t distMaxB = -1e-9;
24 Float_t edepMax = -1e-9;
25 Float_t doseMax = -1e-9;
26 Float_t edepMin = 1e9;
27 Float_t doseMin = 1e9;
28
29 h1 = new TProfile("Energy", "Energy deposits (
30 h2 = new TProfile("Dose", "Dose deposits (Gy)
31 while (1)
32 {
33 ncols = fscanf(fp," %f %f %f %f %f %f %f
34 if (ncols < 0) break;
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 ;
41 // ....
42 h1->Fill(-distB, EdepR); // Basal dendr
43 h1->Fill(distA, EdepR); // Apical dendri
44 h2->Fill(-distB, DoseR);
45 h2->Fill(distA, DoseR);
46 nlines++;
47 }
48 fclose(fp);
49 cout << " Max and Min Energy deposits (keV) ==
50 cout << " Max and Min Dose deposits (Gy) == "
51 cout << " Maximum Basal Distance (um) == " <<
52 cout << " Maximum Apical Distance (um) == " <<
53
54 c1->cd(1);
55 h1->Draw("P");
56 //gPad->SetLogy();
57 h1->SetMarkerSize(2);
58 h1->SetMarkerColor(4);
59 h1->SetMarkerStyle(27);
60 h1->SetFillStyle(3005);
61 //h1->GetYaxis()->SetTitle("Energy deposits in
62 h1->GetXaxis()->SetTitle("Distance from Soma (
63 h1->GetYaxis()->SetRangeUser(edepMin, edepMax+
64 if (distMaxB > 0.)
65 {
66 TLatex text(-distMaxB,edepMax-1.,"Basal");
67 text.DrawClone();
68 }
69 if (distMaxA > 0.)
70 {
71 TLatex text(distMaxA/3.,edepMax-2.,"Apical
72 text.DrawClone();
73 }
74 h1->GetXaxis()->SetRangeUser(-distMaxB-10., di
75 //h1->Fit("gaus");
76
77 c1->cd(2);
78 h2->Draw("P");
79 //gPad->SetLogy();
80 h2->SetMarkerSize(2);
81 h2->SetMarkerColor(kRed);
82 h2->SetMarkerStyle(27);
83 h2->SetFillStyle(3005);
84 //h2->GetYaxis()->SetTitle("Dose deposits in b
85 h2->GetXaxis()->SetTitle("Distance from Soma (
86 h2->GetYaxis()->SetRangeUser(doseMin, doseMax+
87 if (distMaxB > 0.)
88 {
89 TLatex text(-distMaxB,doseMax-0.1,"Basal")
90 text.DrawClone();
91 }
92 if (distMaxA > 0.)
93 {
94 TLatex text(distMaxA/3.,doseMax-0.2,"Apica
95 text.DrawClone();
96 }
97 h2->GetXaxis()->SetRangeUser(-distMaxB-10., di
98
99 }
100
101