Geant4  10.02.p03
Functions | Variables
comparison_ascii.C File Reference
#include <vector>
Include dependency graph for comparison_ascii.C:

Go to the source code of this file.

Functions

cout<< "Reading file \ "<< doseFileExp<< "\ ... ";Long64_t nlines=ntupleExperimental -> ReadFile (doseFileExp, "depthExp:EdepExp")
 
 if (nlines<=0)
 
 printf ("%d Experimental points found\, nlines)
 
ntupleExperimental SetBranchAddress ("EdepExp", &EdepExp)
 
ntupleExperimental SetBranchAddress ("depthExp", &depthExp)
 
ntupleExperimental GetEntry (0)
 
 for (Int_t l=0;l< nentries;l++)
 
ntupleExperimental Reset ()
 
in open (doseFileSim)
 
 if (!in.is_open())
 
TNtupleSim Fill (f1, f2, f3, f4)
 
 while (in.good())
 
in close ()
 
TNtupleSim SetBranchAddress ("dose", &dose)
 
TNtupleSim SetBranchAddress ("iX", &iX)
 
vec_iX clear ()
 
 printf ("%d Simulated points found\, vec_iX.size())
 
TNtupleSim SetMarkerStyle (26)
 
TNtupleSim SetMarkerSize (0.8)
 
ntupleExperimental SetMarkerStyle (4)
 
ntupleExperimental SetMarkerColor (2)
 
ntupleExperimental Draw ("EdepExp:depthExp")
 
TNtupleSim Draw ("dose:iX","","same")
 
leg SetTextSize (0.035)
 
leg SetFillColor (0)
 
leg AddEntry (ntupleExperimental, "Experiment", "P")
 
leg AddEntry (TNtupleSim, "Simulation", "P")
 
leg Draw ()
 

Variables

ifstream in
 
TFile * file = new TFile("Dose.root","RECREATE")
 
TNtuple * ntupleExperimental = new TNtuple("ntupleExperimental","Protons, exp. data", "depthExp:EdepExp")
 
vector< Float_t > vec_dose
 
vector< Float_t > vec_iX
 
TString doseFileExp = "../../../experimentalData/proton/BraggPeak/62MeVInWater.out"
 
Float_t depthExp
 
Float_t EdepExp
 
Int_t nentries = (Int_t)ntupleExperimental -> GetEntries()
 
Float_t normFactor = EdepExp
 
TString doseFileSim = "../../../SimulationOutputs/proton/BraggPeak/Dose.out"
 
TNtuple * TNtupleSim = new TNtuple("SimTree","dose from ascii file", "iX:jY:kZ:dose")
 
Char_t n [5]
 
Float_t f1
 
Float_t f2
 
Float_t f3
 
Float_t f4
 
 nlines = 0
 
 do
 
Float_t iX
 
Float_t dose
 
Float_t sumDose = 0.
 
Float_t norm = 0.
 
Int_t oldX = iX
 
TCanvas * c1 = new TCanvas ("c1","c1",200,10,600,400)
 
 leg = new TLegend(0.50,0.60,0.20,0.70)
 

Function Documentation

◆ AddEntry() [1/2]

leg AddEntry ( ntupleExperimental  ,
"Experiment"  ,
"P"   
)

◆ AddEntry() [2/2]

leg AddEntry ( TNtupleSim  ,
"Simulation"  ,
"P"   
)

◆ clear()

vec_dose clear ( )
Here is the caller graph for this function:

◆ close()

in close ( )
Here is the caller graph for this function:

◆ Draw() [1/3]

ntupleExperimental Draw ( "EdepExp:depthExp"  )

◆ Draw() [2/3]

TNtupleSim Draw ( "dose:iX"  ,
""  ,
"same"   
)

◆ Draw() [3/3]

leg Draw ( )

◆ Fill()

TNtupleSim Fill ( f1  ,
f2  ,
f3  ,
f4   
)
Here is the caller graph for this function:

◆ for()

for ( )

Definition at line 32 of file comparison_ascii.C.

33  {
34  ntupleExperimental -> GetEntry(l);
35  vec_dose.push_back(EdepExp);
36  vec_iX.push_back(depthExp);
37  }
EdepExp
Float_t EdepExp
Definition: comparison_ascii.C:23
vec_dose
vector< Float_t > vec_dose
Definition: comparison_ascii.C:13
ntupleExperimental
TNtuple * ntupleExperimental
Definition: comparison_ascii.C:11
vec_iX
vector< Float_t > vec_iX
Definition: comparison_ascii.C:13
depthExp
Float_t depthExp
Definition: comparison_ascii.C:23
GetEntry
ntupleExperimental GetEntry(0)
Here is the call graph for this function:

◆ GetEntry()

TNtupleSim GetEntry ( 0  )
Here is the caller graph for this function:

◆ if() [1/2]

if ( nlines<=  0)

Definition at line 19 of file comparison_ascii.C.

19 {cout << "Error: Check file \"" << doseFileExp << "\"\n"; return;}
doseFileExp
TString doseFileExp
Definition: comparison_ascii.C:15
Here is the call graph for this function:
Here is the caller graph for this function:

◆ if() [2/2]

if ( !in.  is_open())

Definition at line 54 of file comparison_ascii.C.

54 {cout << "Error: Check file \"" << doseFileSim << "\"\n"; return;}
doseFileSim
TString doseFileSim
Definition: comparison_ascii.C:50

◆ open()

in open ( doseFileSim  )
Here is the caller graph for this function:

◆ printf() [1/2]

printf ( "%d Experimental points found\  ,
nlines   
)
Here is the caller graph for this function:

◆ printf() [2/2]

printf ( "%d Simulated points found\  ,
vec_iX.  size() 
)

◆ ReadFile()

cout<< "Reading file \" " << doseFileExp << "\" ... "; Long64_t nlines = ntupleExperimental -> ReadFile ( doseFileExp  ,
"depthExp:EdepExp"   
)

◆ Reset()

TNtupleSim Reset ( )
Here is the caller graph for this function:

◆ SetBranchAddress() [1/4]

ntupleExperimental SetBranchAddress ( "EdepExp"  ,
&  EdepExp 
)

◆ SetBranchAddress() [2/4]

ntupleExperimental SetBranchAddress ( "depthExp"  ,
&  depthExp 
)

◆ SetBranchAddress() [3/4]

TNtupleSim SetBranchAddress ( "dose"  ,
&  dose 
)

◆ SetBranchAddress() [4/4]

TNtupleSim SetBranchAddress ( "iX"  ,
&  iX 
)

◆ SetFillColor()

leg SetFillColor ( 0  )

◆ SetMarkerColor()

ntuple2 SetMarkerColor ( 2  )
Here is the caller graph for this function:

◆ SetMarkerSize()

ntupleExperimental SetMarkerSize ( 0.  8)

◆ SetMarkerStyle() [1/2]

TNtupleSim SetMarkerStyle ( 26  )

◆ SetMarkerStyle() [2/2]

ntupleExperimental SetMarkerStyle ( 4  )

◆ SetTextSize()

leg SetTextSize ( 0.  035)

◆ while()

while ( in.  good())
Here is the caller graph for this function:

Variable Documentation

◆ c1

TCanvas* c1 = new TCanvas ("c1","c1",200,10,600,400)

Definition at line 113 of file comparison_ascii.C.

◆ depthExp

Float_t depthExp

Definition at line 23 of file comparison_ascii.C.

◆ do

do
Initial value:
{
in >> f1 >> f2 >> f3 >> f4
f4
Float_t f4
Definition: comparison_ascii.C:56
f1
Float_t f1
Definition: comparison_ascii.C:56
f2
Float_t f2
Definition: comparison_ascii.C:56
in
ifstream in
Definition: comparison_ascii.C:1
f3
Float_t f3
Definition: comparison_ascii.C:56

Definition at line 61 of file comparison_ascii.C.

◆ dose

Float_t dose

Definition at line 71 of file comparison_ascii.C.

◆ doseFileExp

TString doseFileExp = "../../../experimentalData/proton/BraggPeak/62MeVInWater.out"

Definition at line 15 of file comparison_ascii.C.

◆ doseFileSim

TString doseFileSim = "../../../SimulationOutputs/proton/BraggPeak/Dose.out"

Definition at line 50 of file comparison_ascii.C.

◆ EdepExp

Float_t EdepExp

Definition at line 23 of file comparison_ascii.C.

◆ f1

Float_t f1

Definition at line 56 of file comparison_ascii.C.

◆ f2

Float_t f2

Definition at line 56 of file comparison_ascii.C.

◆ f3

Float_t f3

Definition at line 56 of file comparison_ascii.C.

◆ f4

Float_t f4

Definition at line 56 of file comparison_ascii.C.

◆ file

TFile* file = new TFile("Dose.root","RECREATE")

Definition at line 5 of file comparison_ascii.C.

◆ in

ifstream in
Initial value:
{
gROOT->Reset()

Definition at line 1 of file comparison_ascii.C.

◆ iX

Float_t iX

Definition at line 71 of file comparison_ascii.C.

◆ leg

leg = new TLegend(0.50,0.60,0.20,0.70)

Definition at line 131 of file comparison_ascii.C.

◆ n

cout<< "Reading file \" " << doseFileSim << "\" ... "; in > n n n n

Definition at line 55 of file comparison_ascii.C.

◆ nentries

nentries = (Int_t)ntupleExperimental -> GetEntries()

Definition at line 29 of file comparison_ascii.C.

◆ nlines

nlines = 0

Definition at line 57 of file comparison_ascii.C.

◆ norm

norm = 0.

Definition at line 71 of file comparison_ascii.C.

◆ normFactor

Float_t normFactor = EdepExp

Definition at line 31 of file comparison_ascii.C.

◆ ntupleExperimental

TNtuple* ntupleExperimental = new TNtuple("ntupleExperimental","Protons, exp. data", "depthExp:EdepExp")

Definition at line 11 of file comparison_ascii.C.

◆ oldX

Int_t oldX = iX

Definition at line 79 of file comparison_ascii.C.

◆ sumDose

Float_t sumDose = 0.

Definition at line 71 of file comparison_ascii.C.

◆ TNtupleSim

TNtuple* TNtupleSim = new TNtuple("SimTree","dose from ascii file", "iX:jY:kZ:dose")

Definition at line 51 of file comparison_ascii.C.

◆ vec_dose

vector<Float_t> vec_dose

Definition at line 13 of file comparison_ascii.C.

◆ vec_iX

vector<Float_t> vec_iX

Definition at line 13 of file comparison_ascii.C.