G4DNAChemistryManager Class Reference

#include <G4DNAChemistryManager.hh>

Public Member Functions
void	WriteInto (const G4String &, std::ios_base::openmode mode=std::ios_base::out)
void	CloseFile ()
G4bool	IsChemistryActived ()
void	SetChemistryActivation (G4bool)
void	CreateWaterMolecule (ElectronicModification, G4int, const G4Track *)
void	CreateSolvatedElectron (const G4Track *, G4ThreeVector *finalPosition=0)
void	PushMolecule (G4Molecule *&molecule, double time, const G4ThreeVector &position, int parentID)
void	PushMoleculeAtParentTimeAndPlace (G4Molecule *&molecule, const G4Track *)

Static Public Member Functions
static G4DNAChemistryManager *	Instance ()
static void	DeleteInstance ()

Protected Member Functions
G4DNAWaterExcitationStructure *	GetExcitationLevel ()
G4DNAWaterIonisationStructure *	GetIonisationLevel ()

Friends
class	std::auto_ptr< G4DNAChemistryManager >

Detailed Description

G4DNAChemistryManager is called from the physics models. It creates the water molecules and the solvated electrons and and send them to G4ITStepManager to be treated in the chemistry stage. For this, the fActiveChemistry flag needs to be on. It is also possible to give already molecule's pointers already built. G4DNAChemistryManager will then be in charge of creating the track and loading it to the IT system. The user can also ask to create a file containing a information about the creation of water molecules and solvated electrons.

Definition at line 71 of file G4DNAChemistryManager.hh.

Member Function Documentation

CloseFile()

void G4DNAChemistryManager::CloseFile

(

)

Close the file specified with WriteInto

Definition at line 117 of file G4DNAChemistryManager.cc.

{
    if (fOutput.is_open())
    {
        fOutput.close();
    }
    fWriteFile = false;
}

CreateSolvatedElectron()

void G4DNAChemistryManager::CreateSolvatedElectron	(	const G4Track *	theIncomingTrack,
		G4ThreeVector *	finalPosition = 0
	)

On the same idea as the previous method but for solvated electron. This method should be used by the physics model of the ElectronSolvatation process.

Definition at line 209 of file G4DNAChemistryManager.cc.

{
    if(fWriteFile)
    {
        fOutput << setw(11)<< theIncomingTrack->GetTrackID()
                << setw(10)<< "e_aq"
                << setw(14)<< -1
                << std::setprecision(2) << std::fixed
                << setw(13)<< theIncomingTrack->GetKineticEnergy()/eV
                << std::setprecision(6) << std::scientific
                << setw(22)<< (theIncomingTrack->GetPosition().x())/nanometer
                << setw(22)<< (theIncomingTrack->GetPosition().y())/nanometer
                << setw(22)<< (theIncomingTrack->GetPosition().z())/nanometer  ;
 
        if(finalPosition != 0)
        {
            fOutput<< setw(14)<< (finalPosition->x())/nanometer
                   << setw(14)<< (finalPosition->y())/nanometer
                   << setw(14)<< (finalPosition->z())/nanometer ;
        }
 
        fOutput << G4endl;
    }
 
    if(fActiveChemistry)
    {
        G4Molecule* e_aq = new G4Molecule(G4Electron_aq::Definition());
        G4Track * e_aqTrack(0);
        if(finalPosition)
        {
            e_aqTrack  = e_aq->BuildTrack(picosecond,*finalPosition);
        }
        else
        {
            e_aqTrack  = e_aq->BuildTrack(picosecond,theIncomingTrack->GetPosition());
        }
        e_aqTrack -> SetTrackStatus(fAlive);
        e_aqTrack -> SetParentID(theIncomingTrack->GetTrackID());
        G4ITTrackHolder::Instance()->PushTrack(e_aqTrack);
        G4ITManager<G4Molecule>::Instance()->Push(e_aqTrack);
    }
}

Referenced by G4DNASancheSolvatationModel::SampleSecondaries(), and G4DNATransformElectronModel::SampleSecondaries().

CreateWaterMolecule()

void G4DNAChemistryManager::CreateWaterMolecule	(	ElectronicModification	modification,
		G4int	electronicLevel,
		const G4Track *	theIncomingTrack
	)

Method used by DNA physics model to create a water molecule. The ElectronicModification is a flag telling wheter the molecule is ionized or excited, the electronic level is calculated by the model and the IncomingTrack is the track responsible for the creation of this molecule, for instance an electron.

Definition at line 144 of file G4DNAChemistryManager.cc.

{
    if(fWriteFile)
    {
        G4double energy = -1.;
 
        switch (modification)
        {
        case eDissociativeAttachment:
            energy = -1;
            break;
        case eExcitedMolecule :
            energy = GetExcitationLevel()->ExcitationEnergy(electronicLevel);
            break;
        case eIonizedMolecule :
            energy = GetIonisationLevel()->IonisationEnergy(electronicLevel);
            break;
        }
 
        fOutput << setw(11) << left << theIncomingTrack->GetTrackID()
                << setw(10) << "H2O"
                << left << modification
                << internal <<":"
                << right <<electronicLevel
                << left
                << setw(11) << ""
                << std::setprecision(2) << std::fixed
                << setw(13) << energy/eV
                << std::setprecision(6) << std::scientific
                << setw(22) << (theIncomingTrack->GetPosition().x())/nanometer
                << setw(22) << (theIncomingTrack->GetPosition().y())/nanometer
                << setw(22) << (theIncomingTrack->GetPosition().z())/nanometer
                << G4endl;
    }
 
    if(fActiveChemistry)
    {
        G4Molecule * H2O = new G4Molecule (G4H2O::Definition());
 
        switch (modification)
        {
        case eDissociativeAttachment:
            H2O -> AddElectron(5,1);
            break;
        case eExcitedMolecule :
            H2O -> ExciteMolecule(electronicLevel);
            break;
        case eIonizedMolecule :
            H2O -> IonizeMolecule(electronicLevel);
            break;
        }
 
        G4Track * H2OTrack = H2O->BuildTrack(1*picosecond,
                                             theIncomingTrack->GetPosition());
 
        H2OTrack -> SetParentID(theIncomingTrack->GetTrackID());
        H2OTrack -> SetTrackStatus(fStopButAlive);
        H2OTrack -> SetKineticEnergy(0.);
 
        G4ITTrackHolder::Instance()->PushTrack(H2OTrack);
    }
}

Referenced by G4DNABornExcitationModel::SampleSecondaries(), G4DNABornIonisationModel::SampleSecondaries(), G4DNAEmfietzoglouExcitationModel::SampleSecondaries(), G4DNAMeltonAttachmentModel::SampleSecondaries(), G4DNAMillerGreenExcitationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), and G4DNARuddIonisationModel::SampleSecondaries().

DeleteInstance()

void G4DNAChemistryManager::DeleteInstance ( )

static

You should rather use DeleteInstance than the destructor of this class

Definition at line 85 of file G4DNAChemistryManager.cc.

{
    if(fInstance.get())
        fInstance.reset();
}

Referenced by G4EmDNAPhysicsChemistry::~G4EmDNAPhysicsChemistry().

GetExcitationLevel()

G4DNAWaterExcitationStructure * G4DNAChemistryManager::GetExcitationLevel ( )

protected

Definition at line 126 of file G4DNAChemistryManager.cc.

{
    if(!fExcitationLevel)
    {
        fExcitationLevel = new G4DNAWaterExcitationStructure;
    }
    return fExcitationLevel;
}

Referenced by CreateWaterMolecule().

GetIonisationLevel()

G4DNAWaterIonisationStructure * G4DNAChemistryManager::GetIonisationLevel ( )

protected

Definition at line 135 of file G4DNAChemistryManager.cc.

{
    if(!fIonisationLevel)
    {
        fIonisationLevel = new G4DNAWaterIonisationStructure;
    }
    return fIonisationLevel;
}

Referenced by CreateWaterMolecule().

Instance()

G4DNAChemistryManager * G4DNAChemistryManager::Instance ( )

static

Definition at line 64 of file G4DNAChemistryManager.cc.

{
    if(!fInstance.get()) fInstance = auto_ptr<G4DNAChemistryManager>(new G4DNAChemistryManager());
    return fInstance.get();
}

Referenced by G4EmDNAPhysicsChemistry::ConstructProcess(), G4DNABornExcitationModel::SampleSecondaries(), G4DNABornIonisationModel::SampleSecondaries(), G4DNAEmfietzoglouExcitationModel::SampleSecondaries(), G4DNAMeltonAttachmentModel::SampleSecondaries(), G4DNAMillerGreenExcitationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4DNASancheSolvatationModel::SampleSecondaries(), and G4DNATransformElectronModel::SampleSecondaries().

IsChemistryActived()

G4bool G4DNAChemistryManager::IsChemistryActived ( )

inline

Definition at line 154 of file G4DNAChemistryManager.hh.

{
    return fActiveChemistry;
}

PushMolecule()

void G4DNAChemistryManager::PushMolecule	(	G4Molecule *&	molecule,
		double	time,
		const G4ThreeVector &	position,
		int	parentID
	)

WARNING : In case chemistry is not activated, PushMolecule will take care of deleting the transfered molecule. Before calling this method, it is also possible to check if the chemistry is activated through IsChemistryActived(). This method will create the track corresponding to the transfered molecule and will be in charge of loading the new track to the system.

Definition at line 255 of file G4DNAChemistryManager.cc.

{
    if(fWriteFile)
    {
        fOutput << setw(11)<< parentID
                << setw(10)<< molecule->GetName()
                << setw(14)<< -1
                << std::setprecision(2) << std::fixed
                << setw(13)<< -1
                << std::setprecision(6) << std::scientific
                << setw(22)<< (position.x())/nanometer
                << setw(22)<< (position.y())/nanometer
                << setw(22)<< (position.z())/nanometer;
        fOutput << G4endl;
    }
 
    if(fActiveChemistry)
    {
        G4Track* track = molecule->BuildTrack(time,position);
        track -> SetTrackStatus(fAlive);
        track -> SetParentID(parentID);
        G4ITTrackHolder::Instance()->PushTrack(track);
        G4ITManager<G4Molecule>::Instance()->Push(track);
    }
    else
    {
        delete molecule;
        molecule = 0;
    }
}

PushMoleculeAtParentTimeAndPlace()

void G4DNAChemistryManager::PushMoleculeAtParentTimeAndPlace	(	G4Molecule *&	molecule,
		const G4Track *	theIncomingTrack
	)

WARNING : In case chemistry is not activated, PushMoleculeAtParentTimeAndPlace will take care of deleting the transfered molecule. Before calling this method, it is also possible to check if the chemistry is activated through IsChemistryActived(). This method will create the track corresponding to the transfered molecule and will be in charge of loading the new track to the system.

Definition at line 287 of file G4DNAChemistryManager.cc.

{
    if(fWriteFile)
    {
        fOutput << setw(11)<< theIncomingTrack->GetTrackID()
                << setw(10)<< molecule->GetName()
                << setw(14)<< -1
                << std::setprecision(2) << std::fixed
                << setw(13)<< theIncomingTrack->GetKineticEnergy()/eV
                << std::setprecision(6) << std::scientific
                << setw(22)<< (theIncomingTrack->GetPosition().x())/nanometer
                << setw(22)<< (theIncomingTrack->GetPosition().y())/nanometer
                << setw(22)<< (theIncomingTrack->GetPosition().z())/nanometer  ;
        fOutput << G4endl;
    }
 
    if(fActiveChemistry)
    {
        G4Track* track = molecule->BuildTrack(theIncomingTrack->GetGlobalTime(),theIncomingTrack->GetPosition());
        track -> SetTrackStatus(fAlive);
        track -> SetParentID(theIncomingTrack->GetTrackID());
        G4ITTrackHolder::Instance()->PushTrack(track);
        G4ITManager<G4Molecule>::Instance()->Push(track);
    }
    else
    {
        delete molecule;
        molecule = 0;
    }
}

SetChemistryActivation()

void G4DNAChemistryManager::SetChemistryActivation ( G4bool  flag )

inline

Definition at line 159 of file G4DNAChemistryManager.hh.

{
    fActiveChemistry = flag;
}

Referenced by G4EmDNAPhysicsChemistry::ConstructProcess().

WriteInto()

void G4DNAChemistryManager::WriteInto	(	const G4String &	,
		std::ios_base::openmode	mode = std::ios_base::out
	)

Tells the chemMan to write into a file the position and electronic state of the water molecule and the position thermalized or not of the solvated electron

Definition at line 91 of file G4DNAChemistryManager.cc.

{
    fOutput.open(output.data(), mode);
    fOutput << std::setprecision(6) << std::scientific;
 
    fOutput << setw(11) << left << "#Parent ID"
            << setw(10) << "Molecule"
            << setw(14) << "Elec Modif"
            << setw(13) << "Energy (eV)"
            << setw(22) << "X pos of parent [nm]"
            << setw(22) << "Y pos of parent [nm]"
            << setw(22) << "Z pos of parent [nm]"
            << setw(14) << "X pos [nm]"
            << setw(14) << "Y pos [nm]"
            << setw(14) << "Z pos [nm]"
            << G4endl
            << setw(21) << "#"
            << setw(13) << "1)io/ex=0/1"
            << G4endl
            << setw(21) << "#"
            << setw(13) << "2)level=0...5"
            << G4endl;
    fWriteFile = true;
}

Friends And Related Function Documentation

std::auto_ptr< G4DNAChemistryManager >

friend class std::auto_ptr< G4DNAChemistryManager >

friend

Definition at line 159 of file G4DNAChemistryManager.hh.

Referenced by Instance().

---

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

• G4DNAChemistryManager.hh
• G4DNAChemistryManager.cc