LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
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().