Functions
gStyle	SetOptStat (0000)

gStyle	SetOptFit ()

gStyle	SetPalette (1)

gROOT	SetStyle ("Plain")

c1	Divide (4, 3)

	while (1)

	fclose (fp)

c1	cd (1)

h1	Draw ()

h1	GetXaxis () -> SetLabelSize(0.025)

h1	GetYaxis () -> SetLabelSize(0.025)

h1	SetLineColor (3)

h1	SetFillColor (3)

h11	SetLineColor (8)

h11	SetFillColor (8)

h11	Draw ("same")

c1	cd (5)

h2	SetLineColor (2)

h2	SetFillColor (2)

h20	SetLineColor (5)

h20	SetFillColor (5)

c1	cd (7)

ntupleYX	Draw ("Y:X>>hist","vox","colz")

gPad	SetLogz ()

hist	Draw ("colz")

hist	GetZaxis () -> SetLabelSize(0.025)

hist	SetTitle ("Cytoplasm intensity on transverse section")

c1	cd (3)

ntupleYXN	Draw ("Y:X>>hist2","vox","colz")

hist2	SetTitle ("Nucleus intensity on transverse section")

	system ("rm -rf microbeam.root")

	system ("hadd microbeam.root microbeam_*.root")

TFile	f ("microbeam.root")

c1	cd (2)

ntuple3	Project ("h1bis","doseN")

h1bis	Scale (scale)

c1	cd (6)

ntuple3	Project ("h10bis","doseC")

c1	cd (9)

ntuple0	Project ("h2bis","sp")

h2bis	SetTitle ("dE/dx at cell entrance")

h2bis	SetFillColor (4)

h2bis	SetLineColor (4)

h2bis	Fit ("gaus")

gaus	SetLineColor (6)

ntuple2	SetBranchAddress ("x",&x)

ntuple2	SetBranchAddress ("y",&y)

ntuple2	SetBranchAddress ("z",&z)

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

c1	cd (10)

ntupleR	Draw ("X2:Z2","abs(X2)<50","surf3")

c1	cd (11)

ntuple4	Draw ("y0.359060:x0.359060>>histbis","doseV","contz")

histbis	Draw ("contz")

histbis	SetTitle ("Mean energy deposit -transverse- (z axis in eV)")

c1	cd (12)

ntuple4	Draw ("x0.359060:(z+1500/0.162810+21)0.162810>>histter","doseV","contz")

histter	SetTitle ("Mean energy deposit -longitudinal- (z axis in eV)")

c1	cd (4)

ntuple1	Project ("hx","x")

h77	SetTitle ("Beam X position on cell")

c1	cd (8)

ntuple1	Project ("hy","y")

h88	SetTitle ("Beam Y position on cell")

Variables
Double_t	scale = 1/h1bis->Integral()

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

FILE *	fp = fopen("phantom.dat","r")

Float_t	xVox

Float_t	yVox

Float_t	zVox

Float_t	tmp

Float_t	den

Float_t	dose

Float_t	X

Float_t	Y

Float_t	Z

Float_t	vox = 0

Float_t	mat = 0

Float_t	voxelSizeX

Float_t	voxelSizeY

Float_t	voxelSizeZ

TH1F *	h1 = new TH1F("h1","Nucleus marker intensity",100,1,300)

TH1F *	h11 = new TH1F("h11 ","",100,1,300)

TH1F *	h2 = new TH1F("h2","Cytoplasm marker intensity",100,1,300)

TH1F *	h20 = new TH1F("h20 ","",100,1,300)

TNtuple *	ntupleYXN = new TNtuple("NUCLEUS","ntuple","Y:X:vox")

TNtuple *	ntupleZX = new TNtuple("CYTOPLASM","ntuple","Z:X:vox")

TNtuple *	ntupleYX = new TNtuple("CYTOPLASM","ntuple","Y:X:vox")

Int_t	nlines =0

Int_t	ncols =0

TH2F *	hist = new TH2F("hist","hist",50,-20,20,50,-20,20)

TH2F *	hist2 = new TH2F("hist2","hist2",50,-20,20,50,-20,20)

TNtuple *	ntuple0 = (TNtuple*)f.Get("ntuple0")

TNtuple *	ntuple1 = (TNtuple*)f.Get("ntuple1")

TNtuple *	ntuple2 = (TNtuple*)f.Get("ntuple2")

TNtuple *	ntuple3 = (TNtuple*)f.Get("ntuple3")

TNtuple *	ntuple4 = (TNtuple*)f.Get("ntuple4")

TH1F *	h1bis = new TH1F("h1bis","Dose distribution in Nucleus",100,0.001,1.)

TH1F *	h10 = new TH1F("h10bis","Dose distribution in Cytoplasm",100,0.001,.2)

Float_t	d

TH1F *	h2bis = new TH1F("h2bis","Beam stopping power at cell entrance",200,0,300)

Double_t	Xc = -1295.59e3 - 955e3sin(10TMath::Pi()/180)

Double_t	Zc = -1327e3 + 955e3cos(10TMath::Pi()/180)

Double_t	X1

Double_t	Y1

Double_t	Z1

Double_t	X2

Double_t	Y2

Double_t	Z2

TNtuple *	ntupleR = new TNtuple("Rmax","ntuple","Z2:Y2:X2")

Double_t	x

Double_t	y

Double_t	z

Double_t	xx

Double_t	zz

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

TH2F *	histbis = new TH2F("histbis","histbis",50,-20,20,50,-20,20)

TH2F *	histter = new TH2F("histter","histter",50,-20,20,50,-20,20)

TH1F *	h77 = new TH1F("hx","h1",200,-10,10)

TH1F *	h88 = new TH1F("hy","h1",200,-10,10)

Function Documentation

c1 cd ( 1 )

Referenced by cluster::ClusterParamsAlg::GetFinalSlope(), and cluster::ClusterParamsAlg::GetProfileInfo().

c1 cd ( 5 )

c1 cd ( 7 )

c1 cd ( 3 )

c1 cd ( 2 )

c1 cd ( 6 )

c1 cd ( 9 )

c1 cd ( 10 )

c1 cd ( 11 )

c1 cd ( 12 )

c1 cd ( 4 )

c1 cd ( 8 )

c1 Divide	(	4	,
		3
	)

h1 Draw ( )

h11 Draw ( "same" )

ntupleYX Draw	(	"Y:X>>hist"	,
		"vox"	,
		"colz"
	)

hist2 Draw ( "colz" )

ntupleYXN Draw	(	"Y:X>>hist2"	,
		"vox"	,
		"colz"
	)

ntupleR Draw ( "X2:Z2" )

ntuple4 Draw	(	"y0.359060:x0.359060>>histbis"	,
		"doseV"	,
		"contz"
	)

histter Draw ( "contz" )

ntuple4 Draw	(	"x0.359060:(z+1500/0.162810+21)0.162810>>histter"	,
		"doseV"	,
		"contz"
	)

TFile f ( "microbeam.root" )

fclose ( fp )

Referenced by if(), and while().

h88 Fit ( "gaus" )

for ( )

Definition at line 284 of file plot.C.

References c1, x, Xc, y, Y1, z, and Zc.

 {
       ntuple2->GetEntry(i);
       X1=x;
       Y1=y;
       Z1=z;
       xx = X1-Xc;
       zz = Z1-Zc;
       Z2 = zz*cos(10*TMath::Pi()/180)-xx*sin(10*TMath::Pi()/180);
       X2 = zz*sin(10*TMath::Pi()/180)+xx*cos(10*TMath::Pi()/180);
       Y2 = Y1;    
       ntupleR->Fill(Z2,Y2,X2);
 }

h1 GetXaxis ( ) -> SetLabelSize(0.025)

h1 GetYaxis ( ) -> SetLabelSize(0.025)

histter GetZaxis ( ) -> SetLabelSize(0.025)

ntuple3 Project	(	"h1bis"	,
		"doseN"
	)

ntuple3 Project	(	"h10bis"	,
		"doseC"
	)

ntuple0 Project	(	"h2bis"	,
		"sp"
	)

ntuple1 Project	(	"hx"	,
		"x"
	)

ntuple1 Project	(	"hy"	,
		"y"
	)

h88 Scale ( scale )

ntuple2 SetBranchAddress	(	"x"	,
		&	x
	)

ntuple2 SetBranchAddress	(	"y"	,
		&	y
	)

ntuple2 SetBranchAddress	(	"z"	,
		&	z
	)

h1bis SetFillColor ( 3 )

h11 SetFillColor ( 8 )

h10 SetFillColor ( 2 )

h20 SetFillColor ( 5 )

h88 SetFillColor ( 4 )

h1bis SetLineColor ( 3 )

h11 SetLineColor ( 8 )

h2 SetLineColor ( 2 )

h20 SetLineColor ( 5 )

h88 SetLineColor ( 4 )

gaus SetLineColor ( 6 )

gPad SetLogz ( )

gStyle SetOptFit ( )

gStyle SetOptStat ( 0000 )

gStyle SetPalette ( 1 )

gROOT SetStyle ( "Plain" )

hist SetTitle ( "Cytoplasm intensity on transverse section" )

hist2 SetTitle ( "Nucleus intensity on transverse section" )

h2bis SetTitle ( "dE/dx at cell entrance" )

histbis SetTitle ( "Mean energy deposit -transverse- (z axis in eV)" )

histter SetTitle ( "Mean energy deposit -longitudinal- (z axis in eV)" )

h77 SetTitle ( "Beam X position on cell" )

h88 SetTitle ( "Beam Y position on cell" )

system ( "rm -rf microbeam.root" )

Referenced by wc::CellTree::endJob(), and wc::CellTree::initOutput().

system ( "hadd microbeam.root microbeam_*.root" )

while ( 1 )

Definition at line 56 of file plot.C.

References c1, fclose(), mat, voxelSizeX, voxelSizeY, and voxelSizeZ.

    {
       if ( nlines == 0 ) ncols = fscanf(fp,"%f %f %f",&tmp,&tmp,&tmp);
       if ( nlines == 1 ) ncols = fscanf(fp,"%f %f %f",&voxelSizeX,&voxelSizeY,&voxelSizeZ);
       if ( nlines == 2 ) ncols = fscanf(fp,"%f %f %f",&tmp,&tmp,&tmp);
       if ( nlines >= 3 ) ncols = fscanf(fp,"%f %f %f %f %f %f", &X, &Y, &Z, &mat, &den, &vox);
       if (ncols < 0) break;
 
       X= X*voxelSizeX;
       Y= Y*voxelSizeY;
       Z= Z*voxelSizeZ;
   
       if ( mat == 2 )  // noyau
          {
           if (den==1) h1->Fill( vox );
           if (den==2) h11->Fill( vox );
           ntupleYXN->Fill(Y,X,vox);
          }
       if ( mat == 1 ) // cytoplasm
          {
           if (den==1) h2->Fill( vox );
           if (den==2) h20->Fill( vox );
           ntupleZX->Fill(Z,X,vox);
           ntupleYX->Fill(Y,X,vox);
          }
       nlines++;
       
    }