LArSoft  v09_93_00
Liquid Argon Software toolkit - https://larsoft.org/
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
plotDend.C File Reference

Go to the source code of this file.

Functions

gStyle SetOptStat (0000)
 
c1 Divide (2, 1)
 
 while (1)
 
 fclose (fp)
 
cout<< " Max and Min Energy deposits (keV) == "<< edepMax<< " ; "<< edepMin<< endl;cout<< " Max and Min Dose deposits (Gy) == "<< doseMax<< " ; "<< doseMin<< endl;cout<< " Maximum Basal Distance (um) == "<< distMaxB<< " "<< endl;cout<< " Maximum Apical Distance (um) == "<< distMaxA<< " "<< endl;c1-> cd (1)
 
h1 Draw ("P")
 
h1 SetMarkerSize (2)
 
h1 SetMarkerColor (4)
 
h1 SetMarkerStyle (27)
 
h1 SetFillStyle (3005)
 
h1 GetXaxis () -> SetTitle("Distance from Soma (um)")
 
h1 GetYaxis () -> SetRangeUser(edepMin, edepMax+3.)
 
 if (distMaxB > 0.)
 
 if (distMaxA > 0.)
 
c1 cd (2)
 
h2 SetMarkerColor (kRed)
 

Variables

 c1 = new TCanvas ("c1","",20,20,1200,600)
 
Int_t ncols =0
 
Int_t nlines = 0
 
FILE * fp = fopen("Dend3DEdep.out","r")
 
Float_t posX
 
Float_t posY
 
Float_t posZ
 
Float_t distB
 
Float_t distA
 
Float_t EdepR
 
Float_t DoseR
 
Float_t distMaxA = -1e-9
 
Float_t distMaxB = -1e-9
 
Float_t edepMax = -1e-9
 
Float_t doseMax = -1e-9
 
Float_t edepMin = 1e9
 
Float_t doseMin = 1e9
 
 h1 = new TProfile("Energy", "Energy deposits (keV) in dendritic compartments",1000,-1000,1000,0.001,1000)
 
 h2 = new TProfile("Dose", "Dose deposits (Gy) in dendritic compartments",1000,-1000,1000,0.001,1000)
 

Function Documentation

cout<< " Max and Min Energy deposits (keV) == " << edepMax << " ; "<< edepMin<<endl;cout << " Max and Min Dose deposits (Gy) == " << doseMax << " ; "<< doseMin<<endl;cout << " Maximum Basal Distance (um) == " << distMaxB << " "<<endl;cout << " Maximum Apical Distance (um) == " << distMaxA << " "<<endl;c1-> cd ( )
c1 cd ( )
c1 Divide ( ,
 
)
h1 Draw ( "P"  )
fclose ( fp  )

Referenced by while().

h1 GetXaxis ( ) -> SetTitle("Distance from Soma (um)")
h1 GetYaxis ( ) -> SetRangeUser(edepMin, edepMax+3.)
if ( distMaxB  ,
0.   
)

Definition at line 64 of file plotDend.C.

65 {
66  TLatex text(-distMaxB,edepMax-1.,"Basal");
67  text.DrawClone();
68 }
Float_t distMaxB
Definition: plotDend.C:23
Float_t edepMax
Definition: plotDend.C:24
if ( distMaxA  ,
0.   
)

Definition at line 69 of file plotDend.C.

References c1, h1, and h2.

70 {
71  TLatex text(distMaxA/3.,edepMax-2.,"Apical");
72  text.DrawClone();
73 }
Float_t distMaxA
Definition: plotDend.C:22
Float_t edepMax
Definition: plotDend.C:24
h2 SetFillStyle ( 3005  )
h1 SetMarkerColor ( )
h2 SetMarkerColor ( kRed  )
h1 SetMarkerSize ( )
h2 SetMarkerStyle ( 27  )
gStyle SetOptStat ( 0000  )
while ( )

Definition at line 31 of file plotDend.C.

References c1, distA, distB, DoseR, EdepR, fclose(), h1, and h2.

32  {
33  ncols = fscanf(fp," %f %f %f %f %f %f %f",&posX, &posY, &posZ, &distA, &distB, &EdepR, &DoseR);
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 dendrite
43  h1->Fill(distA, EdepR); // Apical dendrite
44  h2->Fill(-distB, DoseR);
45  h2->Fill(distA, DoseR);
46  nlines++;
47  }
Float_t posZ
Definition: plotDend.C:20
Float_t EdepR
Definition: plotDend.C:21
h1
Definition: plotDend.C:29
Float_t posY
Definition: plotDend.C:20
h2
Definition: plotDend.C:30
Int_t nlines
Definition: plotDend.C:17
Float_t edepMin
Definition: plotDend.C:26
Float_t doseMin
Definition: plotDend.C:27
Float_t DoseR
Definition: plotDend.C:21
Float_t distB
Definition: plotDend.C:21
FILE * fp
Definition: plotDend.C:19
Float_t distMaxA
Definition: plotDend.C:22
Int_t ncols
Definition: plotDend.C:16
Float_t distMaxB
Definition: plotDend.C:23
Float_t posX
Definition: plotDend.C:20
Float_t doseMax
Definition: plotDend.C:25
Float_t distA
Definition: plotDend.C:21
Float_t edepMax
Definition: plotDend.C:24

Variable Documentation

c1 = new TCanvas ("c1","",20,20,1200,600)

Definition at line 13 of file plotDend.C.

Referenced by if(), and while().

Float_t distA

Definition at line 21 of file plotDend.C.

Referenced by while().

Float_t distB

Definition at line 21 of file plotDend.C.

Referenced by while().

Float_t distMaxA = -1e-9

Definition at line 22 of file plotDend.C.

Float_t distMaxB = -1e-9

Definition at line 23 of file plotDend.C.

Float_t doseMax = -1e-9

Definition at line 25 of file plotDend.C.

Float_t doseMin = 1e9

Definition at line 27 of file plotDend.C.

Float_t DoseR

Definition at line 21 of file plotDend.C.

Referenced by while().

Float_t edepMax = -1e-9

Definition at line 24 of file plotDend.C.

Float_t edepMin = 1e9

Definition at line 26 of file plotDend.C.

Float_t EdepR

Definition at line 21 of file plotDend.C.

Referenced by while().

FILE* fp = fopen("Dend3DEdep.out","r")

Definition at line 19 of file plotDend.C.

h1 = new TProfile("Energy", "Energy deposits (keV) in dendritic compartments",1000,-1000,1000,0.001,1000)

Definition at line 29 of file plotDend.C.

Referenced by if(), and while().

h2 = new TProfile("Dose", "Dose deposits (Gy) in dendritic compartments",1000,-1000,1000,0.001,1000)

Definition at line 30 of file plotDend.C.

Referenced by if(), and while().

Int_t ncols =0

Definition at line 16 of file plotDend.C.

Int_t nlines = 0

Definition at line 17 of file plotDend.C.

Float_t posX

Definition at line 20 of file plotDend.C.

Referenced by trkf::TrackTrajectoryAlg::TrackTrajectory().