Functions
gStyle	SetPalette (1)

gROOT	SetStyle ("Plain")

gStyle	SetOptStat (00000)

	if (fp)

c1	cd (1)

ntuple	Draw ("y>>hfy","","")

hfy	Scale (1./(hfy->GetEntries() *hfy->GetBinWidth(1)))

hfy	SetTitle ("f(y) (um/keV)")

hfy	GetXaxis () -> SetTitle("y (keV/um)")

hfy	SetFillColor (1)

hfy	Draw ("")

	for (Int_t j=0;j< hfy->GetNbinsX();j++) norm =normLogfy+(log(10)/B)*hlogyfy->GetBinContent(j)

cout<< endl;cout<< "** I - Results from lin-lin histograms **"<< endl;cout<< endl;cout<< "---> sum of	f (hdy->SetTitle("d() (um/keV)") y)

axis	Set (bins, new_bins)

ntuple	SetBranchAddress ("radius",&radius)

ntuple	SetBranchAddress ("y",&y)

ntuple	SetBranchAddress ("z",&z)

	for (Int_t i=0;i< nentries;i++)

hlogyfy	Scale (1./(hlogyfy->GetEntries()))

hlogy2fy	Scale (1./(hlogy2fy->GetEntries()))

hlogydy	Scale (1./(hlogydy->GetEntries()))

hlogy3fy	Scale (1./(hlogy3fy->GetEntries()))

gPad	SetLogx ()

hlogyfy	SetLineColor (2)

hlogyfy	SetFillColor (2)

hlogyfy	SetTitle ("y*f(y)")

cout<< "---> sum of Log	f (y)

Variables
Int_t	B =1000

Int_t	min =-2

Int_t	max =2

FILE *	fp = fopen("yz.root","r")

	else

cout<< "* it will be created from merged ROOT files "<< endl;system("rm -rf yz.root");system("hadd yz.root yz_.root");}c1=new TCanvas("c1","", 60, 60, 800, 600);c1.Divide(3, 2);TFile f("yz.root");TNtuple ntuple;ntuple=(TNtuple ) f->	Get ("yz")

Double_t	norm =0

delete	new_bins

Double_t	radius

Double_t	y

Double_t	z

Int_t	nentries = (Int_t)ntuple->GetEntries()

Double_t	normLogfy =0

Double_t	logyF =0

Function Documentation

◆ cd()

c1 cd ( 1 )

◆ Draw() [1/2]

ntuple Draw	(	"y>>hfy"	,
		""	,
		""
	)

◆ Draw() [2/2]

hlogyfy Draw ( "" )

◆ f() [1/2]

cout<< endl;cout << "**** I - Results from lin-lin histograms ****" << endl;cout << endl;cout << "---> sum of f ( hdy->SetTitle("d() (um/keV)") y )

Definition at line 91 of file extended/medical/dna/microyz/plot.C.

 {
   new_bins[i] = TMath::Power(10, from + i * width);
 } 

◆ f() [2/2]

cout<< "---> sum of Log f ( y )

Definition at line 286 of file extended/medical/dna/microyz/plot.C.

 {
  hlogy2fy->SetBinContent(i,
 //(log(10)/B)*
 hlogy2fy->GetBinCenter(i)
 *hlogy2fy->GetBinCenter(i)
 *hlogy2fy->GetBinContent(i)
 /(TMath::Power(10, from + (i+1) * width)-TMath::Power(10, from + i * width))
 );
 }

◆ for() [1/2]

for ( Int_t j = 0;j< hfy->GetNbinsX();j++ ) =normLogfy+(log(10)/B)*hlogyfy->GetBinContent(j)

Definition at line 263 of file extended/medical/dna/microyz/plot.C.

 {
  hlogyfy->SetBinContent(i,
 //(log(10)/B)*
 hlogyfy->GetBinCenter(i)
 *hlogyfy->GetBinContent(i)
 /(TMath::Power(10, from + (i+1) * width)-TMath::Power(10, from + i * width))
 );
 }

◆ for() [2/2]

for ( )

Definition at line 241 of file extended/medical/dna/microyz/plot.C.

 {
  ntuple->GetEntry(i);
       hlogyfy->Fill(y);
       hlogy2fy->Fill(y);
       hlogydy->Fill(y);
       hlogy3fy->Fill(y);
 
 }

◆ GetXaxis()

hfy GetXaxis ( ) -> SetTitle("y (keV/um)")

◆ if()

if ( fp )

Definition at line 31 of file extended/medical/dna/microyz/plot.C.

          {
   // exists
   cout << "*** Notice: the output file yz.root exists ***"<< endl;
   fclose(fp);
 } else {

Here is the call graph for this function:

◆ Scale() [1/5]

hfy Scale ( 1./ hfy->GetEntries() *hfy->GetBinWidth(1) )

◆ Scale() [2/5]

hlogyfy Scale ( 1./ hlogyfy->GetEntries() )

◆ Scale() [3/5]

hlogy2fy Scale ( 1./ hlogy2fy->GetEntries() )

◆ Scale() [4/5]

hlogydy Scale ( 1./ hlogydy->GetEntries() )

◆ Scale() [5/5]

hlogy3fy Scale ( 1./ hlogy3fy->GetEntries() )

◆ Set()

axis4 Set	(	bins	,
		new_bins
	)

◆ SetBranchAddress() [1/3]

ntuple SetBranchAddress	(	"radius"	,
		&	radius
	)

◆ SetBranchAddress() [2/3]

ntuple SetBranchAddress	(	"y"	,
		&	y
	)

◆ SetBranchAddress() [3/3]

ntuple SetBranchAddress	(	"z"	,
		&	z
	)

◆ SetFillColor() [1/2]

hfy SetFillColor ( 1 )

◆ SetFillColor() [2/2]

hlogyfy SetFillColor ( 2 )

◆ SetLineColor()

hlogyfy SetLineColor ( 2 )

◆ SetLogx()

gPad SetLogx ( )

◆ SetOptStat()

gStyle SetOptStat ( 00000 )

◆ SetPalette()

gStyle SetPalette ( 1 )

◆ SetStyle()

gROOT SetStyle ( "Plain" )

◆ SetTitle() [1/2]

hfy SetTitle ( "f(y) (um/keV)" )

◆ SetTitle() [2/2]

hlogyfy SetTitle ( "y*f(y)" )

Variable Documentation

◆ B

Int_t B =1000

Definition at line 25 of file extended/medical/dna/microyz/plot.C.

◆ else

else

Initial value:

{

cout << "*** Notice: the output file yz.root does not exist ***"<< endl

Definition at line 35 of file extended/medical/dna/microyz/plot.C.

◆ fp

FILE* fp = fopen("yz.root","r")

Definition at line 30 of file extended/medical/dna/microyz/plot.C.

◆ Get

cout<< "*** it will be created from merged ROOT files ***"<< endl; system ("rm -rf yz.root"); system ("hadd yz.root yz_*.root");}c1 = new TCanvas ("c1","",60,60,800,600);c1.Divide(3,2);TFile f("yz.root"); TNtuple* ntuple;ntuple = (TNtuple*)f-> Get("yz")

Definition at line 66 of file extended/medical/dna/microyz/plot.C.