IORTAnalysisManager Class Reference

#include <IORTAnalysisManager.hh>

Public Member Functions
	~IORTAnalysisManager ()
void	book ()
void	SetAnalysisFileName (G4String)
void	FillEnergyDeposit (G4int voxelXId, G4int voxelYId, G4int voxelZId, G4double energyDeposit)
void	BraggPeak (G4int, G4double)
	Fill 1D histogram with the Bragg peak in the phantom. More...
void	SecondaryProtonEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary protons. More...
void	SecondaryNeutronEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary neutrons. More...
void	SecondaryAlphaEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary alpha particles. More...
void	SecondaryGammaEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary gamma. More...
void	SecondaryElectronEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary electrons. More...
void	SecondaryTritonEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary tritons. More...
void	SecondaryDeuteronEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary deuterons. More...
void	SecondaryPionEnergyDeposit (G4int slice, G4double energy)
	Fill 1D histogram with the energy deposit of secondary pions. More...
void	electronEnergyDistribution (G4double secondaryParticleKineticEnergy)
	Energy distribution of secondary electrons originated in the phantom. More...
void	gammaEnergyDistribution (G4double secondaryParticleKineticEnergy)
	Energy distribution of secondary gamma originated in the phantom. More...
void	deuteronEnergyDistribution (G4double secondaryParticleKineticEnergy)
	Energy distribution of secondary deuterons originated in the phantom. More...
void	tritonEnergyDistribution (G4double secondaryParticleKineticEnergy)
	Energy distribution of secondary tritons originated in the phantom. More...
void	alphaEnergyDistribution (G4double secondaryParticleKineticEnergy)
	Energy distribution of secondary alpha originated in the phantom. More...
void	heliumEnergy (G4double secondaryParticleKineticEnergy)
	Energy distribution of the helium (He3 and alpha) particles after the phantom. More...
void	hydrogenEnergy (G4double secondaryParticleKineticEnergy)
	Energy distribution of the hydrogen (proton, d, t) particles after the phantom. More...
void	FillKineticFragmentTuple (G4int i, G4int j, G4int k, G4int A, G4double Z, G4double kinEnergy)
void	FillKineticEnergyPrimaryNTuple (G4int i, G4int j, G4int k, G4double kinEnergy)
	Energy by voxel, mass number A and atomic number Z. More...
void	FillVoxelFragmentTuple (G4int i, G4int j, G4int k, G4int A, G4double Z, G4double energy, G4double fluence)
void	FillFragmentTuple (G4int A, G4double Z, G4double energy, G4double posX, G4double posY, G4double posZ)
	Energy ntuple. More...
void	FillLetFragmentTuple (G4int i, G4int j, G4int k, G4int A, G4double Z, G4double letT, G4double letD)
	let ntuple More...
void	genericIonInformation (G4int, G4double, G4int, G4double)
void	ThintargetBeamDisp (G4double, G4double)
void	startNewEvent ()
	Tell the analysis manager that a new event is starting. More...
void	setGeometryMetaData (G4double, G4double, G4double)
	from the detector construction information about the geometry can be written as metadata More...
void	setBeamMetaData (G4double, G4double)
	metadata about the beam can be written this way More...
void	flush ()
	Close the .hbk file with the histograms and the ntuples. More...

Static Public Member Functions
static IORTAnalysisManager *	GetInstance ()

Detailed Description

A class for connecting the simulation to an analysis package.

Definition at line 56 of file IORTAnalysisManager.hh.

Constructor & Destructor Documentation

IORTAnalysisManager::~IORTAnalysisManager

(

)

Definition at line 58 of file IORTAnalysisManager.cc.

{
  if (fMess)
    delete fMess; 
  delete G4AnalysisManager::Instance(); 
}

Member Function Documentation

void IORTAnalysisManager::alphaEnergyDistribution

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of secondary alpha originated in the phantom.

Definition at line 311 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(14,energy);
}

void IORTAnalysisManager::book

(

)

Book the histograms and ntuples.

Definition at line 82 of file IORTAnalysisManager.cc.

{
  // Create analysis manager
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->SetVerboseLevel(1);
  man->SetFirstHistoId(1);
  man->SetFirstNtupleId(1);

  man->OpenFile(analysisFileName);
  
  // Create the histograms with the energy deposit along the X axis
  //ID=1 <different waterthicknesses are accoutned for in ROOT-analysis stage>
  man->CreateH1("braggPeak","slice, energy", 400, 0., 80); //
  //ID=2
  man->CreateH1("h20","Secondary protons - slice, energy", 400, 0., 400.);
  //ID=3
  man->CreateH1("h30","Secondary neutrons - slice, energy", 400, 0., 400.);
  //ID=4
  man->CreateH1("h40","Secondary alpha - slice, energy", 400, 0., 400.);
  //ID=5
  man->CreateH1("h50","Secondary gamma - slice, energy", 400, 0., 400.);
  //ID=6
  man->CreateH1("h60","Secondary electron - slice, energy", 400, 0., 400.);
  //ID=7
  man->CreateH1("h70","Secondary triton - slice, energy", 400, 0., 400.);
  //ID=8
  man->CreateH1("h80","Secondary deuteron - slice, energy", 400, 0., 400.);
  //ID=9
  man->CreateH1("h90","Secondary pion - slice, energy", 400, 0., 400.);
  //ID=10
  man->CreateH1("h100","Energy distribution of secondary electrons", 70, 0., 70.);
  //ID=11
  man->CreateH1("h110","Energy distribution of secondary photons", 70, 0., 70.);
  //ID=12
  man->CreateH1("h120","Energy distribution of secondary deuterons", 70, 0., 70.);
  //ID = 13
  man->CreateH1("h130","Energy distribution of secondary tritons", 70, 0., 70.);
  //ID = 14
  man->CreateH1("h140","Energy distribution of secondary alpha particles", 70, 0., 70.);
  //ID = 15
  man->CreateH1("heliumEnergyAfterPhantom",
        "Energy distribution of secondary helium fragments after the phantom",
        70, 0., 500.);
  //ID= 16
  man->CreateH1("hydrogenEnergyAfterPhantom",
        "Energy distribution of secondary helium fragments after the phantom",
        70, 0., 500.);

  //Now the ntuples
  //ID = 1
  man->CreateNtuple("kinFragNtuple", 
            "Kinetic energy by voxel & fragment");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleIColumn("A");
  man->CreateNtupleDColumn("Z");
  man->CreateNtupleDColumn("kineticEnergy");
  man->FinishNtuple();

  //ID = 2
  man->CreateNtuple("kineticEnergyPrimaryNtuple", 
            "Kinetic energy by voxel of primary");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleDColumn("kineticEnergy");
  man->FinishNtuple();

  //ID = 3
  man->CreateNtuple("doseFragNtuple",
            "Energy deposit by voxel & fragment");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleIColumn("A");
  man->CreateNtupleDColumn("Z");
  man->CreateNtupleDColumn("energy");
  man->FinishNtuple();

  // ID =4
  man->CreateNtuple("fluenceFragNtuple", 
            "Fluence by voxel & fragment");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleIColumn("A");
  man->CreateNtupleDColumn("Z");
  man->CreateNtupleDColumn("fluence");
  man->FinishNtuple();

  // ID=5
  man->CreateNtuple("letFragNtuple", 
            "Let by voxel & fragment");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleIColumn("A");
  man->CreateNtupleDColumn("Z");
  man->CreateNtupleDColumn("letT");
  man->CreateNtupleDColumn("letD");
  man->FinishNtuple();

  //ID=6
  man->CreateNtuple("theROOTNtuple", 
            "Energy deposit by slice");
  man->CreateNtupleIColumn("i");
  man->CreateNtupleIColumn("j");
  man->CreateNtupleIColumn("k");
  man->CreateNtupleDColumn("energy");
  man->FinishNtuple();

  //ID=7
  man->CreateNtuple("theROOTIonTuple",
            "Generic ion information");
  man->CreateNtupleIColumn("a");
  man->CreateNtupleDColumn("z");
  man->CreateNtupleIColumn("occupancy");
  man->CreateNtupleDColumn("energy");
  man->FinishNtuple();

  //ID=8
  man->CreateNtuple("fragmentNtuple",
            "Fragments");
  man->CreateNtupleIColumn("A");
  man->CreateNtupleDColumn("Z");
  man->CreateNtupleDColumn("energy");
  man->CreateNtupleDColumn("posX");
  man->CreateNtupleDColumn("posY");
  man->CreateNtupleDColumn("posZ");
  man->FinishNtuple();

}

Here is the call graph for this function:

void IORTAnalysisManager::BraggPeak	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the Bragg peak in the phantom.

Definition at line 231 of file IORTAnalysisManager.cc.

{
  //FIXME
  G4AnalysisManager::Instance()->FillH1(1,slice,energy);
  //histo1->SetBinContent(slice, energy); //This uses setbincontent instead of fill to get labels correct
}

void IORTAnalysisManager::deuteronEnergyDistribution

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of secondary deuterons originated in the phantom.

Definition at line 299 of file IORTAnalysisManager.cc.

{
 G4AnalysisManager::Instance()->FillH1(12,energy);
}

void IORTAnalysisManager::electronEnergyDistribution

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of secondary electrons originated in the phantom.

Definition at line 287 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(10,energy);
}

void IORTAnalysisManager::FillEnergyDeposit	(	G4int	voxelXId,
		G4int	voxelYId,
		G4int	voxelZId,
		G4double	energyDeposit
	)

Fill the ntuple with the energy deposit in the phantom

Definition at line 217 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(6,0,i);
  man->FillNtupleIColumn(6,1,j);
  man->FillNtupleIColumn(6,2,k);
  man->FillNtupleDColumn(6,3,energy);
  man->AddNtupleRow(6);  
}

Here is the call graph for this function:

void IORTAnalysisManager::FillFragmentTuple	(	G4int	A,
		G4double	Z,
		G4double	energy,
		G4double	posX,
		G4double	posY,
		G4double	posZ
	)

Energy ntuple.

Definition at line 401 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(8,0,A);
  man->FillNtupleDColumn(8,1,Z);
  man->FillNtupleDColumn(8,2,energy);
  man->FillNtupleDColumn(8,3,posX);
  man->FillNtupleDColumn(8,4,posY);
  man->FillNtupleDColumn(8,5,posZ);
  man->AddNtupleRow(8); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTAnalysisManager::FillKineticEnergyPrimaryNTuple	(	G4int	i,
		G4int	j,
		G4int	k,
		G4double	kinEnergy
	)

Energy by voxel, mass number A and atomic number Z.

Definition at line 346 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(2,0,i);
  man->FillNtupleIColumn(2,1,j);
  man->FillNtupleIColumn(2,2,k);
  man->FillNtupleDColumn(2,3,kinEnergy);
  man->AddNtupleRow(2);  
}

Here is the call graph for this function:

void IORTAnalysisManager::FillKineticFragmentTuple	(	G4int	i,
		G4int	j,
		G4int	k,
		G4int	A,
		G4double	Z,
		G4double	kinEnergy
	)

Definition at line 331 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(1,0,i);
  man->FillNtupleIColumn(1,1,j);
  man->FillNtupleIColumn(1,2,k);
  man->FillNtupleIColumn(1,3,A);
  man->FillNtupleDColumn(1,4,Z);
  man->FillNtupleDColumn(1,5,kinEnergy);
  man->AddNtupleRow(1);  
}

Here is the call graph for this function:

void IORTAnalysisManager::FillLetFragmentTuple	(	G4int	i,
		G4int	j,
		G4int	k,
		G4int	A,
		G4double	Z,
		G4double	letT,
		G4double	letD
	)

let ntuple

Definition at line 386 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(5,0,i);
  man->FillNtupleIColumn(5,1,j);
  man->FillNtupleIColumn(5,2,k);
  man->FillNtupleIColumn(5,3,A);
  man->FillNtupleDColumn(5,4,Z);
  man->FillNtupleDColumn(5,5,letT);
  man->FillNtupleDColumn(5,6,letD);
  man->AddNtupleRow(5);  
}

Here is the call graph for this function:

void IORTAnalysisManager::FillVoxelFragmentTuple	(	G4int	i,
		G4int	j,
		G4int	k,
		G4int	A,
		G4double	Z,
		G4double	energy,
		G4double	fluence
	)

Definition at line 361 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(3,0,i);
  man->FillNtupleIColumn(3,1,j);
  man->FillNtupleIColumn(3,2,k);
  man->FillNtupleIColumn(3,3,A);
  man->FillNtupleDColumn(3,4,Z);
  man->FillNtupleDColumn(3,5,energy);
  man->AddNtupleRow(3);
  
  
  // Fill the ntuple containing the voxel, mass and atomic number and the fluence
  if (i==1 && Z==1) {
    man->FillNtupleIColumn(4,0,i);
    man->FillNtupleIColumn(4,1,j);
    man->FillNtupleIColumn(4,2,k);
    man->FillNtupleIColumn(4,3,A);
    man->FillNtupleDColumn(4,4,Z);
    man->FillNtupleDColumn(4,5,fluence);
    man->AddNtupleRow(4);    
  }
}

Here is the call graph for this function:

void IORTAnalysisManager::flush

(

)

Close the .hbk file with the histograms and the ntuples.

Definition at line 450 of file IORTAnalysisManager.cc.

{
  // Save histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->Write();
  man->CloseFile();
  eventCounter = 0;
}

Here is the call graph for this function:

void IORTAnalysisManager::gammaEnergyDistribution

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of secondary gamma originated in the phantom.

Definition at line 293 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(11,energy);
}

void IORTAnalysisManager::genericIonInformation	(	G4int	a,
		G4double	z,
		G4int	electronOccupancy,
		G4double	energy
	)

Definition at line 415 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleIColumn(7,0,a);
  man->FillNtupleDColumn(7,1,z);
  man->FillNtupleIColumn(7,2,electronOccupancy);
  man->FillNtupleDColumn(7,3,energy);
  man->AddNtupleRow(7);  
}

Here is the call graph for this function:

IORTAnalysisManager * IORTAnalysisManager::GetInstance ( )

static

Get the pointer to the analysis manager.

Definition at line 67 of file IORTAnalysisManager.cc.

{
  if (instance == 0) instance = new IORTAnalysisManager;
  return instance;
}

Here is the caller graph for this function:

void IORTAnalysisManager::heliumEnergy

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of the helium (He3 and alpha) particles after the phantom.

Definition at line 317 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(15,secondaryParticleKineticEnergy);
}

Here is the caller graph for this function:

void IORTAnalysisManager::hydrogenEnergy

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of the hydrogen (proton, d, t) particles after the phantom.

Definition at line 323 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(16,secondaryParticleKineticEnergy);
}

Here is the caller graph for this function:

void IORTAnalysisManager::SecondaryAlphaEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary alpha particles.

Definition at line 251 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(4,slice,energy);
}

void IORTAnalysisManager::SecondaryDeuteronEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary deuterons.

Definition at line 275 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(8,slice,energy);
}

void IORTAnalysisManager::SecondaryElectronEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary electrons.

Definition at line 263 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(6,slice,energy);
}

void IORTAnalysisManager::SecondaryGammaEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary gamma.

Definition at line 257 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(5,slice,energy);
}

void IORTAnalysisManager::SecondaryNeutronEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary neutrons.

Definition at line 245 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(3,slice,energy);
}

void IORTAnalysisManager::SecondaryPionEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary pions.

Definition at line 281 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(9,slice,energy);
}

void IORTAnalysisManager::SecondaryProtonEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary protons.

Definition at line 239 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(2,slice,energy);
}

void IORTAnalysisManager::SecondaryTritonEnergyDeposit	(	G4int	slice,
		G4double	energy
	)

Fill 1D histogram with the energy deposit of secondary tritons.

Definition at line 269 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(7,slice,energy);
}

void IORTAnalysisManager::SetAnalysisFileName

(

G4String

aFileName

)

Set name for the analysis file .root (used by macro)

Definition at line 75 of file IORTAnalysisManager.cc.

{
  analysisFileName = aFileName;
}

void IORTAnalysisManager::setBeamMetaData	(	G4double	meanKineticEnergy,
		G4double	sigmaEnergy
	)

metadata about the beam can be written this way

Definition at line 443 of file IORTAnalysisManager.cc.

{
  beamEnergy = meanKineticEnergy;
  energyError = sigmaEnergy;
}

void IORTAnalysisManager::setGeometryMetaData	(	G4double	endDetectorPosition,
		G4double	waterThickness,
		G4double	phantomCenter
	)

from the detector construction information about the geometry can be written as metadata

Definition at line 434 of file IORTAnalysisManager.cc.

{
  detectorDistance = endDetectorPosition;
  phantomDepth = waterThickness;
  phantomCenterDistance = phantomCenter;
}

void IORTAnalysisManager::startNewEvent

(

)

Tell the analysis manager that a new event is starting.

Definition at line 429 of file IORTAnalysisManager.cc.

{
  eventCounter++;
}

Here is the caller graph for this function:

void IORTAnalysisManager::ThintargetBeamDisp	(	G4double	,
		G4double
	)

void IORTAnalysisManager::tritonEnergyDistribution

(

G4double

secondaryParticleKineticEnergy

)

Energy distribution of secondary tritons originated in the phantom.

Definition at line 305 of file IORTAnalysisManager.cc.

{
  G4AnalysisManager::Instance()->FillH1(13,energy);
}

---

The documentation for this class was generated from the following files:

• source/geant4.10.03.p03/examples/advanced/iort_therapy/include/IORTAnalysisManager.hh
• source/geant4.10.03.p03/examples/advanced/iort_therapy/src/IORTAnalysisManager.cc