#include <G4DNAMolecularDissociation.hh>

Inheritance diagram for G4DNAMolecularDissociation:

Collaboration diagram for G4DNAMolecularDissociation:

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

virtual	~G4DNAMolecularDissociation ()

virtual G4bool	IsApplicable (const G4ParticleDefinition &)

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

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

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

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

void	SetVerbose (G4int)

void	SetDecayDisplacer (const G4ParticleDefinition , G4VMolecularDecayDisplacer )

G4VMolecularDecayDisplacer *	GetDecayDisplacer (const G4ParticleDefinition *)

void	SetDisplacer (const G4ParticleDefinition , G4VMolecularDecayDisplacer )

G4VMolecularDecayDisplacer *	GetDisplacer (const G4ParticleDefinition *)

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

	G4VITRestDiscreteProcess (const G4VITRestDiscreteProcess &)

virtual	~G4VITRestDiscreteProcess ()

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

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

G4shared_ptr< G4ProcessState_Lock >	GetProcessState ()

void	SetProcessState (G4shared_ptr< G4ProcessState_Lock > aProcInfo)

void	ResetProcessState ()

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 *)

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

Protected Member Functions inherited from G4VITRestDiscreteProcess
	G4VITRestDiscreteProcess ()

G4VITRestDiscreteProcess &	operator= (const G4VITRestDiscreteProcess &right)

Protected Member Functions inherited from G4VITProcess
void	RetrieveProcessInfo ()

void	CreateInfo ()

template<typename T >
T *	GetState ()

virtual void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)

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
G4shared_ptr< 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 66 of file G4DNAMolecularDissociation.hh.

Constructor & Destructor Documentation

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

Definition at line 52 of file G4DNAMolecularDissociation.cc.

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

Here is the call graph for this function:

G4DNAMolecularDissociation::~G4DNAMolecularDissociation ( )

virtual

Definition at line 80 of file G4DNAMolecularDissociation.cc.

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

Member Function Documentation

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

inlinevirtual

Reimplemented from G4VITRestDiscreteProcess.

Definition at line 155 of file G4DNAMolecularDissociation.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

inlinevirtual

Reimplemented from G4VITRestDiscreteProcess.

Definition at line 142 of file G4DNAMolecularDissociation.hh.

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

Here is the call graph for this function:

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

protectedvirtual

Definition at line 141 of file G4DNAMolecularDissociation.cc.

 {
   // 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 G4MolecularDissociationChannel*>* DecayVector =
         (theMotherMolecule->GetDecayChannel());
 
     if(DecayVector == 0)
     {
       G4ExceptionDescription exceptionDescription;
       theMotherMolecule->PrintState();
       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 G4MolecularDissociationChannel* 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;
 
       DisplacementMap::iterator it =
           fDisplacementMap.find(moleculeDefinition);
 
       if(it != fDisplacementMap.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);
 
 //      G4cout << " nbProducts = " << nbProducts << G4endl;
 //      theMotherMolecule->PrintState();
 
 #ifdef G4VERBOSE
       if(fVerbose)
       {
         G4cout << "Decay Process : " << theMotherMolecule->GetName()
                << " (trackID :" << track.GetTrackID() << ") "
                << decayChannel->GetName() << G4endl;
       }
 #endif
 
       G4ITNavigator* navigator =
           G4ITTransportationManager::GetTransportationManager()
               ->GetNavigatorForTracking();
 
       for(G4int j = 0; j < nbProducts; j++)
       {
         G4Molecule* product = new G4Molecule(decayChannel->GetProduct(j));
 
         G4ThreeVector displacement = theMotherMoleculeDisplacement
             + ProductsDisplacement[j];
         double mag_displacement = displacement.mag();
         G4ThreeVector displacement_direction =
           displacement/(mag_displacement+1e-30);
         
         double prNewSafety = DBL_MAX;
 
         //double step =
         navigator->CheckNextStep(track.GetPosition(),
                                  displacement_direction,
                                  mag_displacement,
                                  prNewSafety);
 
         //if(prNewSafety < mag_displacement || step < mag_displacement)
         mag_displacement = min(prNewSafety*0.8, mag_displacement);
 
         G4ThreeVector product_pos = track.GetPosition()
                        + displacement_direction * mag_displacement;
         
         const G4AffineTransform& transform = navigator
             ->GetGlobalToLocalTransform();
 
         G4ThreeVector localPoint =
           transform.TransformPoint(product_pos); //track.GetPosition());
 
         if(track.GetTouchable()->GetSolid()->Inside(localPoint) !=
            EInside::kInside)
         {
           G4cout << "Mother volume: "
                  << track.GetTouchable()->GetVolume()->GetName()
                  << G4endl;
           G4Exception("G4DNAMolecularDissociation::DecayIt",
                       "OUTSIDE_OF_MOTHER_VOLUME",
                       FatalException,
                       "Product has been placed outside of the volume "
                       "containing the mother molecule");
         }
 
         G4Track* secondary =
             product->BuildTrack(track.GetGlobalTime(),
                                 product_pos);
 
         secondary->SetTrackStatus(fAlive);
 #ifdef G4VERBOSE
         if(fVerbose)
         {
           G4cout << "Product : " << product->GetName() << G4endl;
         }
 #endif
         // add the secondary track in the List
         aParticleChange.G4VParticleChange::AddSecondary(secondary);
       }
 #ifdef G4VERBOSE
       if(fVerbose) G4cout << "-------------" << G4endl;
 #endif
     }
     //DEBUG
     else if(fVerbose && 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 "
      "dissociation channel";
      G4Exception("G4MolecularDissociationProcess::DecayIt",
      "DNAMolecularDissociation002",
      FatalErrorInArgument,
      errMsg);
      }
      */
   }
 
   aParticleChange.ProposeTrackStatus(fStopAndKill);
 
   return &aParticleChange;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 354 of file G4DNAMolecularDissociation.cc.

 {
   return fDisplacementMap[molDef];
 }

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

Definition at line 372 of file G4DNAMolecularDissociation.cc.

 {
   return fDisplacementMap[molDef];
 }

virtual G4double G4DNAMolecularDissociation::GetMeanFreePath	(	const G4Track &	,
		G4double	,
		G4ForceCondition *
	)

inlineprotectedvirtual

Implements G4VITRestDiscreteProcess.

Definition at line 116 of file G4DNAMolecularDissociation.hh.

     {
       return 0;
     }

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

protectedvirtual

Implements G4VITRestDiscreteProcess.

Definition at line 132 of file G4DNAMolecularDissociation.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4DNAMolecularDissociation::IsApplicable ( const G4ParticleDefinition & aParticleType )

virtual

Reimplemented from G4VProcess.

Definition at line 109 of file G4DNAMolecularDissociation.cc.

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

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from G4VITRestDiscreteProcess.

Definition at line 93 of file G4DNAMolecularDissociation.hh.

     {
       return AtRestDoIt(track, step);
     }

Here is the call graph for this function:

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

virtual

Reimplemented from G4VITRestDiscreteProcess.

Definition at line 380 of file G4DNAMolecularDissociation.cc.

 {
   return 0; //-1*(-log(1-G4UniformRand())*100*1e-15*s);
 }

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

Definition at line 344 of file G4DNAMolecularDissociation.cc.

 {
   fDisplacementMap[molDef] = aDisplacer;
 }

Here is the caller graph for this function:

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

Definition at line 362 of file G4DNAMolecularDissociation.cc.

 {
   fDisplacementMap[molDef] = aDisplacer;
 }

void G4DNAMolecularDissociation::SetVerbose ( G4int verbose )

inline

Definition at line 137 of file G4DNAMolecularDissociation.hh.

 {
     fVerbose = verbose ;
 }

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

geant4.10.03.p01/source/processes/electromagnetic/dna/processes/include/G4DNAMolecularDissociation.hh
geant4.10.03.p01/source/processes/electromagnetic/dna/processes/src/G4DNAMolecularDissociation.cc

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation