CML2ExpVoxels Class Reference

#include <ML2ExpVoxels.hh>

Public Member Functions
	CML2ExpVoxels (G4bool bHasExperimentalData, G4int saving_in_Selected_Voxels_every_events, G4int seed, G4String FileExperimentalData, G4String FileExperimentalDataOut)
	~CML2ExpVoxels (void)
void	add (G4ThreeVector pos, G4double depEnergy, G4double density)
void	add (const G4Step *aStep)
std::vector< Svoxel >	getVoxels ()
G4int	getMinNumberOfEvents ()
G4int	getMaxNumberOfEvents ()
G4bool	loadData ()
void	setRecycling (int recycling)
void	saveResults (void)
void	resetNEventsInVoxels ()

Detailed Description

Definition at line 51 of file ML2ExpVoxels.hh.

Constructor & Destructor Documentation

CML2ExpVoxels::CML2ExpVoxels	(	G4bool	bHasExperimentalData,
		G4int	saving_in_Selected_Voxels_every_events,
		G4int	seed,
		G4String	FileExperimentalData,
		G4String	FileExperimentalDataOut
	)

Definition at line 48 of file ML2ExpVoxels.cc.

References test::a, and CLHEP::Hep3Vector::set().

                                                                          :startCurve(0),
          stopCurve(0),chi2Factor(0)
{
    char a[10];
    sprintf(a,"%d", seed);
    seedName=(G4String)a;
    saving_in_Selected_Voxels_every_events=saveEvents;
    nRecycling=1;

    
    fullFileOut=FileExperimentalDataOut+seedName+".m";
    fullFileIn=FileExperimentalData;
    nParticle=nTotalEvents=0;

// define the extremes of global-volume containing all experimental voxels
    G4double extr=100000000000.;
    minZone.set(extr, extr, extr);
    maxZone.set(-extr, -extr, -extr);
    bHasExperimentalData=bData;
}

CML2ExpVoxels::~CML2ExpVoxels

(

void

)

Definition at line 70 of file ML2ExpVoxels.cc.

{
    delete [] startCurve;
    delete [] stopCurve;
    delete [] chi2Factor;
        delete [] nVoxelsgeometry;
}

Member Function Documentation

void CML2ExpVoxels::add	(	G4ThreeVector	pos,
		G4double	depEnergy,
		G4double	density
	)

Referenced by CML2Convergence::add().

void CML2ExpVoxels::add

(

const G4Step *

aStep

)

Definition at line 165 of file ML2ExpVoxels.cc.

References density, Sparticle::dir, G4Track::GetDefinition(), G4Material::GetDensity(), G4VPhysicalVolume::GetLogicalVolume(), G4LogicalVolume::GetMaterial(), G4StepPoint::GetMomentumDirection(), G4ParticleDefinition::GetPDGEncoding(), G4StepPoint::GetPhysicalVolume(), G4StepPoint::GetPosition(), G4Step::GetPreStepPoint(), G4Step::GetTotalEnergyDeposit(), G4Step::GetTrack(), CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), CLHEP::Hep3Vector::getZ(), int(), Sparticle::kinEnergy, Sparticle::nPrimaryPart, Sparticle::partPDGE, Sparticle::pos, Sparticle::primaryParticlePDGE, saveResults(), Sparticle::volumeId, and Sparticle::volumeName.

{
    G4ThreeVector pos;
    G4double depEnergy, density, voxelVolume;
    
    pos=aStep->GetPreStepPoint()->GetPosition();
    depEnergy=aStep->GetTotalEnergyDeposit();
    density=aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetDensity();

    G4ThreeVector minPos, maxPos;
    G4bool newEvent=false;
    G4double voxelMass, dose;

// check if the event is inside the global-volume 
    if (minZone.getX()<= pos.getX() && pos.getX()<maxZone.getX() && 
        minZone.getY()<= pos.getY() && pos.getY()<maxZone.getY() && 
        minZone.getZ()<= pos.getZ() && pos.getZ()<maxZone.getZ())
    {
// look for the voxel containing the event
        for (int i=0; i<(int)voxels.size(); i++)
        {
            minPos=voxels[i].pos-voxels[i].halfSize;
            maxPos=voxels[i].pos+voxels[i].halfSize;
            if (minPos.getX()<= pos.getX() && pos.getX()<maxPos.getX() && 
                minPos.getY()<= pos.getY() && pos.getY()<maxPos.getY() && 
                minPos.getZ()<= pos.getZ() && pos.getZ()<maxPos.getZ())
            {
                voxelVolume=voxels[i].halfSize.getX()*voxels[i].halfSize.getY()*voxels[i].halfSize.getZ()*8.;
                voxelMass=density*voxelVolume;
// calculate the dose 
                dose=depEnergy/(voxelMass*nRecycling);
                voxels[i].nEvents++;
                                nVoxelsgeometry[i]++;
                voxels[i].depEnergy+=dose;
                voxels[i].depEnergy2+=dose*dose;
                newEvent=true;

                Sparticle *particle=new Sparticle;
                particle->dir=aStep->GetPreStepPoint()->GetMomentumDirection();
                particle->pos=aStep->GetPreStepPoint()->GetPosition();
                particle->kinEnergy=dose; // I use the same kinEnergy name to store the dose
                particle->nPrimaryPart=-1;
                particle->partPDGE=aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
                particle->primaryParticlePDGE=-1;
                particle->volumeId=i; // voxel index where the dose is accumulating
                particle->volumeName="-1";
            }
        }
        if (newEvent)
        {
// save data
            nTotalEvents++;
            if (nTotalEvents%saving_in_Selected_Voxels_every_events==0 && nTotalEvents>0)
            {
                saveResults();
            }
        }
    }
}

G4int CML2ExpVoxels::getMaxNumberOfEvents

(

)

Definition at line 234 of file ML2ExpVoxels.cc.

References int(), and n.

{
        int n=nVoxelsgeometry[0];
    for (int i=0;i<(int)voxels.size();i++)
    { 
                if (n<nVoxelsgeometry[i]){n = nVoxelsgeometry[i];}
    }
    return n;
}

G4int CML2ExpVoxels::getMinNumberOfEvents

(

)

Definition at line 225 of file ML2ExpVoxels.cc.

References int(), and n.

{
    int n=voxels[0].nEvents;
    for (int i=0;i<(int)voxels.size();i++)
    { 
        if (n>voxels[i].nEvents){n = voxels[i].nEvents;}
    }
    return n;
}

std::vector<Svoxel> CML2ExpVoxels::getVoxels ( )

inline

Definition at line 59 of file ML2ExpVoxels.hh.

{return voxels;}

G4bool CML2ExpVoxels::loadData

(

void

)

Definition at line 77 of file ML2ExpVoxels.cc.

References test::a, Svoxel::depEnergy, Svoxel::depEnergy2, Svoxel::depEnergyNorm, Svoxel::depEnergyNormError, Svoxel::expDose, G4endl, CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), CLHEP::Hep3Vector::getZ(), Svoxel::halfSize, python.hepunit::joule, python.hepunit::kg, Svoxel::nEvents, Svoxel::pos, resetNEventsInVoxels(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), and Svoxel::volumeId.

Referenced by CML2Convergence::CML2Convergence().

{
    bHasExperimentalData=true;
    std::ifstream in;

    Svoxel voxel; voxel.volumeId=0;
    G4ThreeVector pos, halfSize;
    G4double expDose;

    in.open(fullFileIn, std::ios::in);
    if (in)
    {
        G4String appo;
        char a[1000];
        in.getline(a,1000,'\n'); headerText1=(G4String)a;
        in.getline(a,1000,'\n');
        in >> nCurves;
        startCurve=new G4int[nCurves];
        stopCurve=new G4int[nCurves]; 
        chi2Factor=new G4double[nCurves];
        for (int i=0; i< nCurves; i++)
        {
            chi2Factor[i]=0.;
            in >> startCurve[i]; 
            in >> stopCurve[i]; 
            in >> chi2Factor[i];
        }
        in.getline(a,1000,'\n');
        in.getline(a,1000,'\n'); headerText2=(G4String)a;

        while (!in.eof())
        {
            in >> pos;
            in >> halfSize; 
            if (bHasExperimentalData)
            {
                in >> expDose;
                voxel.expDose=expDose/100.*(joule/kg); // input data in cGy
            }
            else
            {
                voxel.expDose=0.;
            }
            voxel.pos=pos;
            voxel.halfSize=halfSize;
            voxel.depEnergy=0.;
            voxel.depEnergy2=0.;
            voxel.nEvents=0;
            voxel.depEnergyNorm=0.;
            voxel.depEnergyNormError=0.;
            voxels.push_back(voxel);

// calculate the actual extremes of the global-volume containing all the experimental data
            if (minZone.getX()>pos.getX()-halfSize.getX())
            {minZone.setX(pos.getX()-halfSize.getX());}
            if (maxZone.getX()<pos.getX()+halfSize.getX())
            {maxZone.setX(pos.getX()+halfSize.getX());}

            if (minZone.getY()>pos.getY()-halfSize.getY())
            {minZone.setY(pos.getY()-halfSize.getY());}
            if (maxZone.getY()<pos.getY()+halfSize.getY())
            {maxZone.setY(pos.getY()+halfSize.getY());}

            if (minZone.getZ()>pos.getZ()-halfSize.getZ())
            {minZone.setZ(pos.getZ()-halfSize.getZ());}
            if (maxZone.getZ()<pos.getZ()+halfSize.getZ())
            {maxZone.setZ(pos.getZ()+halfSize.getZ());}

        }
    }
    else
    {
        std::cout << "ERROR I can't find the experimental data file" << G4endl;
        return false;
    }
    in.close();

        nVoxelsgeometry=new G4int[(G4int) voxels.size()];
        resetNEventsInVoxels();

    return true;
}

void CML2ExpVoxels::resetNEventsInVoxels

(

)

Definition at line 159 of file ML2ExpVoxels.cc.

References int().

Referenced by loadData().

{
    for (int i=0; i<(int) voxels.size(); i++ )
    {nVoxelsgeometry[i]=0;}
}

void CML2ExpVoxels::saveResults

(

void

)

Definition at line 262 of file ML2ExpVoxels.cc.

References testem0::app, G4endl, int(), python.hepunit::joule, python.hepunit::kg, and python.hepunit::mm.

Referenced by add(), and CML2Convergence::saveResults().

{
    if (nTotalEvents>0)
    {
        calculateNormalizedEd(voxels);
        saveHeader();
        std::ofstream out;
        out.open(fullFileOut, std::ios::app);
        out <<"d"<< seedName<<"=["<< G4endl;
        for (int i=0; i<(int)voxels.size(); i++)
        {
            out <<voxels[i].pos.getX()/mm<<'\t'<<voxels[i].pos.getY()/mm<<'\t'<<voxels[i].pos.getZ()/mm<<'\t';
            out <<voxels[i].halfSize.getX()/mm<<'\t'<<voxels[i].halfSize.getY()/mm<<'\t'<<voxels[i].halfSize.getZ()/mm<<'\t';
            out <<voxels[i].expDose/(joule/kg)<<'\t'<<voxels[i].depEnergy/(joule/kg)<<'\t'<<voxels[i].depEnergy2/((joule/kg)*(joule/kg))<<'\t'<<voxels[i].nEvents<< '\t';
            out <<voxels[i].depEnergyNorm/(joule/kg)<<'\t'<<voxels[i].depEnergyNormError/(joule/kg);
            out << G4endl;
        }
        out << "];"<<G4endl;
        out.close();
    }
}

void CML2ExpVoxels::setRecycling ( int  recycling )

inline

Definition at line 66 of file ML2ExpVoxels.hh.

Referenced by CML2Convergence::CML2Convergence().

{nRecycling=recycling;}

---

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

• ML2ExpVoxels.hh
• ML2ExpVoxels.cc