#include <G4DNAMolecularDissociation.hh>

Inheritance diagram for G4DNAMolecularDissociation:

Public Member Functions
	G4DNAMolecularDissociation (const G4String &processName="DNAMolecularDecay", G4ProcessType type=fDecay)

virtual	~G4DNAMolecularDissociation ()

virtual G4bool	IsApplicable (const G4ParticleDefinition &)

G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)

G4VParticleChange *	AtRestDoIt (const G4Track &track, const G4Step &step)

void	SetVerbose (G4int)

void	SetDecayDisplacer (const G4ParticleDefinition , G4VMolecularDecayDisplacer )

G4VMolecularDecayDisplacer *	GetDecayDisplacer (const G4ParticleDefinition *)

Public Member Functions inherited from G4VITRestProcess
	G4VITRestProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VITRestProcess (const G4VITRestProcess &)

virtual	~G4VITRestProcess ()

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &, G4double, G4ForceCondition *)

virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)

Public Member Functions inherited from G4VITProcess
	G4VITProcess (const G4String &name, G4ProcessType type=fNotDefined)

virtual	~G4VITProcess ()

	G4VITProcess (const G4VITProcess &other)

G4VITProcess &	operator= (const G4VITProcess &other)

G4int	operator== (const G4VITProcess &right) const

G4int	operator!= (const G4VITProcess &right) const

size_t	GetProcessID () const

G4ProcessState_Lock *	GetProcessState ()

void	SetProcessState (G4ProcessState_Lock *aProcInfo)

virtual void	StartTracking (G4Track *)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

G4double	GetInteractionTimeLeft ()

virtual void	ResetNumberOfInteractionLengthLeft ()

G4bool	ProposesTimeStep () const

Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)

	G4VProcess (const G4VProcess &right)

virtual	~G4VProcess ()

G4int	operator== (const G4VProcess &right) const

G4int	operator!= (const G4VProcess &right) const

G4double	GetCurrentInteractionLength () const

void	SetPILfactor (G4double value)

G4double	GetPILfactor () const

G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)

G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)

G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual void	PreparePhysicsTable (const G4ParticleDefinition &)

virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)

const G4String &	GetProcessName () const

G4ProcessType	GetProcessType () const

void	SetProcessType (G4ProcessType)

G4int	GetProcessSubType () const

void	SetProcessSubType (G4int)

virtual void	EndTracking ()

virtual void	SetProcessManager (const G4ProcessManager *)

virtual const G4ProcessManager *	GetProcessManager ()

G4double	GetNumberOfInteractionLengthLeft () const

G4double	GetTotalNumberOfInteractionLengthTraversed () const

G4bool	isAtRestDoItIsEnabled () const

G4bool	isAlongStepDoItIsEnabled () const

G4bool	isPostStepDoItIsEnabled () const

virtual void	DumpInfo () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

virtual void	SetMasterProcess (G4VProcess *masterP)

const G4VProcess *	GetMasterProcess () const

virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)

virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Member Functions
virtual G4VParticleChange *	DecayIt (const G4Track &, const G4Step &)

virtual G4double	GetMeanLifeTime (const G4Track &, G4ForceCondition *)

Protected Member Functions inherited from G4VITRestProcess
	G4VITRestProcess ()

G4VITRestProcess &	operator= (const G4VITRestProcess &right)

Protected Member Functions inherited from G4VITProcess
void	RetrieveProcessInfo ()

void	CreateInfo ()

virtual void	ClearInteractionTimeLeft ()

virtual void	ClearNumberOfInteractionLengthLeft ()

void	SetInstantiateProcessState (G4bool flag)

G4bool	InstantiateProcessState ()

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Additional Inherited Members
Static Public Member Functions inherited from G4VITProcess
static const size_t &	GetMaxProcessIndex ()

Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Attributes inherited from G4VITProcess
G4ProcessState *	fpState

G4bool	fProposesTimeStep

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

Detailed Description

G4DNAMolecularDissociation should be called only for molecules. It will dissociate the molecules using the decay associated to this molecule and if a displacement scheme has been registered, it will place the products to the expected position.

Definition at line 56 of file G4DNAMolecularDissociation.hh.

Constructor & Destructor Documentation

G4DNAMolecularDissociation::G4DNAMolecularDissociation	(	const G4String &	processName = `"DNAMolecularDecay"`,
		G4ProcessType	type = `fDecay`
	)

Definition at line 48 of file G4DNAMolecularDissociation.cc.

References G4VProcess::aParticleChange, G4VProcess::enableAlongStepDoIt, G4VProcess::enableAtRestDoIt, G4VProcess::enablePostStepDoIt, G4cout, G4DNAMolecularDissociation(), G4endl, G4VProcess::pParticleChange, G4VProcess::SetProcessSubType(), and G4VProcess::verboseLevel.

Referenced by G4DNAMolecularDissociation().

                                                              : G4VITRestProcess(processName, type)
 {
     // set Process Sub Type
     SetProcessSubType(59); // DNA sub-type
     enableAlongStepDoIt = false;
     enablePostStepDoIt = false;
     enableAtRestDoIt=true;
 
     fVerbose = 0 ;
 
 #ifdef G4VERBOSE
     if (verboseLevel>1)
     {
         G4cout << "G4MolecularDecayProcess constructor " << "  Name:" << processName << G4endl;
     }
 #endif
 
     pParticleChange = &aParticleChange;
 
     fDecayAtFixedTime = true ;
 }

G4DNAMolecularDissociation::~G4DNAMolecularDissociation ( )

virtual

Definition at line 71 of file G4DNAMolecularDissociation.cc.

 {
     DecayDisplacementMap::iterator it = fDecayDisplacementMap.begin();
 
     for( ; it != fDecayDisplacementMap.end() ; it++)
     {
         if(it->second)
         {
             delete it->second;
             it->second = 0;
         }
     }
     fDecayDisplacementMap.clear();
 }

Member Function Documentation

G4VParticleChange * G4DNAMolecularDissociation::AtRestDoIt	(	const G4Track &	track,
		const G4Step &	step
	)

inlinevirtual

Reimplemented from G4VITRestProcess.

Definition at line 120 of file G4DNAMolecularDissociation.hh.

References G4VITProcess::ClearInteractionTimeLeft(), G4VITProcess::ClearNumberOfInteractionLengthLeft(), and DecayIt().

 {
     ClearNumberOfInteractionLengthLeft();
     ClearInteractionTimeLeft();
     return DecayIt(track, step);
 }

G4double G4DNAMolecularDissociation::AtRestGetPhysicalInteractionLength	(	const G4Track &	track,
		G4ForceCondition *	condition
	)

inlinevirtual

Reimplemented from G4VITRestProcess.

Definition at line 108 of file G4DNAMolecularDissociation.hh.

References G4VITRestProcess::AtRestGetPhysicalInteractionLength(), and GetMeanLifeTime().

 {
     if(fDecayAtFixedTime)
     {
         return GetMeanLifeTime(track, condition);
     }
 
     return G4VITRestProcess::AtRestGetPhysicalInteractionLength(track, condition);
 }

G4VParticleChange * G4DNAMolecularDissociation::DecayIt	(	const G4Track &	track,
		const G4Step &
	)

protectedvirtual

Definition at line 121 of file G4DNAMolecularDissociation.cc.

References G4VProcess::aParticleChange, G4Molecule::BuildTrack(), G4ElectronOccupancy::DumpInfo(), fAlive, FatalErrorInArgument, FatalException, fStopAndKill, G4cout, G4endl, G4Exception(), G4UniformRand, G4Molecule::GetDefinition(), G4Molecule::GetElectronOccupancy(), G4MolecularDecayChannel::GetEnergy(), GetMolecule(), G4MolecularDecayChannel::GetName(), G4Molecule::GetName(), G4MolecularDecayChannel::GetNbProducts(), G4Track::GetPosition(), G4MolecularDecayChannel::GetProbability(), G4MolecularDecayChannel::GetProduct(), G4VMolecularDecayDisplacer::GetProductsDisplacement(), G4Track::GetTrackID(), G4ParticleChange::Initialize(), G4ITManager< T >::Instance(), python.hepunit::picosecond, G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeTrackStatus(), and G4VParticleChange::SetNumberOfSecondaries().

Referenced by AtRestDoIt().

 {
     // DEBUG
     //    G4cout << "Is calling G4MolecularDecayProcess::DecayIt" << G4endl;
 
     aParticleChange.Initialize(track);
     const G4Molecule * theMotherMolecule = GetMolecule(track);
     const G4MoleculeDefinition* moleculeDefinition = theMotherMolecule->GetDefinition();
 
     //    DEBUG
     //        G4cout <<"Calling G4MolecularDecayProcess::DecayIt"<<G4endl;
     //        G4cout << "The mother molecule state : " << G4endl;
     //        theMotherMolecule -> PrintState();
 
     if(moleculeDefinition-> GetDecayTable())
     {
         const vector<const G4MolecularDecayChannel*>* DecayVector =
                 (theMotherMolecule -> GetDecayChannel());
 
         if(DecayVector == 0)
         {
             G4ExceptionDescription exceptionDescription;
             theMotherMolecule->GetElectronOccupancy()->DumpInfo();
             exceptionDescription << "No decay channel was found for the molecule : " << theMotherMolecule-> GetName() << G4endl;
             G4Exception("G4DNAMolecularDissociation::DecayIt", "G4DNAMolecularDissociation::NoDecayChannel",FatalException,exceptionDescription);
             return &aParticleChange;
         }
 
         G4int DecayVectorSize = DecayVector-> size();
         //        DEBUG
         //            G4cout<< "Number of decay channels : " << DecayVectorSize<<G4endl;
         G4double RdmValue = G4UniformRand();
 
         const G4MolecularDecayChannel* decayChannel(0);
         G4int i=0;
         do
         {
             decayChannel = (*DecayVector)[i];
             if(RdmValue < decayChannel->GetProbability()) break;
             RdmValue -= decayChannel->GetProbability();
             i++;
         }
         while(i< DecayVectorSize);
 
         //        DEBUG
         //            G4cout<< "Selected Decay channel : " << decayChannel->GetName()<<G4endl;
 
         G4double decayEnergy = decayChannel->GetEnergy();
         G4int nbProducts = decayChannel->GetNbProducts();
 
         if(decayEnergy)
         {
             //            DEBUG
             //                G4cout<<"Deposit energy :" <<decayChannel->GetEnergy()/eV << " eV"<<G4endl;
 
             aParticleChange.ProposeLocalEnergyDeposit(decayChannel->GetEnergy());
         }
 
         if(nbProducts)
         {
 
             //            DEBUG
             //                G4cout<<"Number of products :" <<nbProducts<<G4endl;
 
             vector<G4ThreeVector> ProductsDisplacement(nbProducts);
             G4ThreeVector theMotherMoleculeDisplacement;
 
             DecayDisplacementMap::iterator it = fDecayDisplacementMap.find(moleculeDefinition);
 
             if(it!=fDecayDisplacementMap.end())
             {
                 G4VMolecularDecayDisplacer* displacer = it->second;
                 ProductsDisplacement = displacer->GetProductsDisplacement(decayChannel);
                 theMotherMoleculeDisplacement = displacer-> GetMotherMoleculeDisplacement(decayChannel);
             }
             else
             {
                 G4ExceptionDescription errMsg;
                 errMsg << "No G4MolecularDecayProcess::theDecayDisplacementMap["
                        << theMotherMolecule->GetName() +"]" ;
                 G4Exception("G4MolecularDecayProcess::DecayIt","DNAMolecularDecay001",FatalErrorInArgument, errMsg);
             }
 
             aParticleChange.SetNumberOfSecondaries(nbProducts);
 
 #ifdef G4VERBOSE
             if(fVerbose)
             {
                 G4cout<<"Decay Process : "
                      << theMotherMolecule->GetName()
                      << " (trackID :" << track.GetTrackID() << ") "
                      << decayChannel->GetName()
                      << G4endl;
             }
 #endif
 
             for (G4int j=0; j<nbProducts ; j++)
             {
                 G4Molecule* product = new G4Molecule(*decayChannel->GetProduct(j));
 
                 // create a new track object
                 // Be carefull as this processes is dedicated to be used in the DNA module
                 // The molecular decay will happen one picosecond after the start of the simulation.
                 // This may be seen as a bug and will be hopefully improve in the next releases
                 G4Track* secondary = product->BuildTrack(picosecond,track.GetPosition()
                                                          + theMotherMoleculeDisplacement + ProductsDisplacement[j]);
 
                 secondary-> SetTrackStatus(fAlive);
 #ifdef G4VERBOSE
                 if(fVerbose)
                 {
                     G4cout<<"Product : "<< product->GetName()<<G4endl;
                 }
 #endif
                 // add the secondary track in the List
                 aParticleChange.G4VParticleChange::AddSecondary(secondary);
                 G4ITManager<G4Molecule>::Instance()->Push(secondary);
             }
 #ifdef G4VERBOSE
             if(fVerbose)
                 G4cout<<"-------------"<<G4endl;
 #endif
         }
         //        DEBUG
         //        else if(decayEnergy)
         //        {
         //            G4cout << "No products for this channel" << G4endl ;
         //            G4cout<<"-------------"<<G4endl;
         //        }
         else if(!decayEnergy && !nbProducts)
         {
             G4ExceptionDescription errMsg;
             errMsg << "There is no products and no energy specified in the molecular decay channel";
             G4Exception("G4MolecularDecayProcess::DecayIt","DNAMolecularDecay002",FatalErrorInArgument, errMsg);
         }
     }
 
     aParticleChange.ProposeTrackStatus(fStopAndKill);
 
     return &aParticleChange;
 }

G4VMolecularDecayDisplacer * G4DNAMolecularDissociation::GetDecayDisplacer ( const G4ParticleDefinition * molDef )

Definition at line 270 of file G4DNAMolecularDissociation.cc.

 {
     return fDecayDisplacementMap[molDef] ;
 }

G4double G4DNAMolecularDissociation::GetMeanLifeTime	(	const G4Track &	track,
		G4ForceCondition *
	)

protectedvirtual

Implements G4VITRestProcess.

Definition at line 114 of file G4DNAMolecularDissociation.cc.

References GetMolecule(), and G4Track::GetProperTime().

Referenced by AtRestGetPhysicalInteractionLength().

 {
     G4double output = GetMolecule(track)-> GetDecayTime() - track.GetProperTime() ;
     return (output > 0 ? output : 0 );
 }

G4bool G4DNAMolecularDissociation::IsApplicable ( const G4ParticleDefinition & aParticleType )

virtual

Reimplemented from G4VProcess.

Definition at line 94 of file G4DNAMolecularDissociation.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), and G4ParticleDefinition::GetParticleType().

 {
     if(aParticleType.GetParticleType()=="Molecule")
     {
 #ifdef G4VERBOSE
         if(fVerbose>1)
         {
             G4cout<<"G4MolecularDecay::IsApplicable(";
             G4cout<<aParticleType.GetParticleName()<<",";
             G4cout<<aParticleType.GetParticleType()<<")"<<G4endl;
         }
 #endif
         return(true);
     }
     else
     {
         return false;
     }
 }

void G4DNAMolecularDissociation::SetDecayDisplacer	(	const G4ParticleDefinition *	molDef,
		G4VMolecularDecayDisplacer *	aDisplacer
	)

Definition at line 265 of file G4DNAMolecularDissociation.cc.

 {
     fDecayDisplacementMap[molDef] = aDisplacer;
 }

void G4DNAMolecularDissociation::SetVerbose ( G4int verbose )

inline

Definition at line 103 of file G4DNAMolecularDissociation.hh.

 {
     fVerbose = verbose ;
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation