Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4DNASecondOrderReaction Class Reference

#include <G4DNASecondOrderReaction.hh>

+ Inheritance diagram for G4DNASecondOrderReaction:

Classes

struct  SecondOrderReactionState
 

Public Member Functions

 G4DNASecondOrderReaction (const G4String &aName="DNASecondOrderReaction", G4ProcessType type=fDecay)
 
virtual ~G4DNASecondOrderReaction ()
 
 G4DNASecondOrderReaction (const G4DNASecondOrderReaction &)
 
G4DNASecondOrderReactionoperator= (const G4DNASecondOrderReaction &)
 
void StartTracking (G4Track *)
 
void SetReaction (const G4MolecularConfiguration *, const G4Material *, double)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4VParticleChangePostStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChangeAtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4VParticleChangeAlongStepDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VITProcess
 G4VITProcess (const G4String &name, G4ProcessType type=fNotDefined)
 
virtual ~G4VITProcess ()
 
 G4VITProcess (const G4VITProcess &other)
 
G4VITProcessoperator= (const G4VITProcess &other)
 
G4bool operator== (const G4VITProcess &right) const
 
G4bool operator!= (const G4VITProcess &right) const
 
size_t GetProcessID () const
 
G4shared_ptr< G4ProcessState_LockGetProcessState ()
 
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 ()
 
G4VProcessoperator= (const G4VProcess &)=delete
 
G4bool operator== (const G4VProcess &right) const
 
G4bool operator!= (const G4VProcess &right) const
 
virtual G4VParticleChangePostStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChangeAlongStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChangeAtRestDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0
 
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 G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4StringGetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4StringGetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManagerGetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
virtual void ProcessDescription (std::ostream &outfile) const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcessGetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Protected Attributes

G4bool fIsInitialized
 
G4double fReturnedValue
 
const std::vector< double > * fpMoleculeDensity
 
G4double fReactionRate
 
G4double fConcentration
 
G4double fMolarMassOfMaterial
 
G4ParticleChange fParticleChange
 
const G4MolecularConfigurationfpMolecularConfiguration
 
const G4MaterialfpMaterial
 
- Protected Attributes inherited from G4VITProcess
G4shared_ptr< G4ProcessStatefpState
 
G4bool fProposesTimeStep
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManageraProcessManager = nullptr
 
G4VParticleChangepParticleChange = nullptr
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft = -1.0
 
G4double currentInteractionLength = -1.0
 
G4double theInitialNumberOfInteractionLength = -1.0
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType = fNotDefined
 
G4int theProcessSubType = -1
 
G4double thePILfactor = 1.0
 
G4int verboseLevel = 0
 
G4bool enableAtRestDoIt = true
 
G4bool enableAlongStepDoIt = true
 
G4bool enablePostStepDoIt = true
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VITProcess
static const size_t & GetMaxProcessIndex ()
 
- Static Public Member Functions inherited from G4VProcess
static const G4StringGetProcessTypeName (G4ProcessType)
 
- 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 prevStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Detailed Description

Definition at line 54 of file G4DNASecondOrderReaction.hh.

Constructor & Destructor Documentation

◆ G4DNASecondOrderReaction() [1/2]

G4DNASecondOrderReaction::G4DNASecondOrderReaction ( const G4String aName = "DNASecondOrderReaction",
G4ProcessType  type = fDecay 
)

Definition at line 67 of file G4DNASecondOrderReaction.cc.

67 :
68 G4VITProcess(aName,type)
69{
70 Create();
71}

◆ ~G4DNASecondOrderReaction()

G4DNASecondOrderReaction::~G4DNASecondOrderReaction ( )
virtual

Definition at line 79 of file G4DNASecondOrderReaction.cc.

80{
81 ;
82}

◆ G4DNASecondOrderReaction() [2/2]

G4DNASecondOrderReaction::G4DNASecondOrderReaction ( const G4DNASecondOrderReaction rhs)

Definition at line 73 of file G4DNASecondOrderReaction.cc.

73 :
74 G4VITProcess(rhs)
75{
76 Create();
77}

Member Function Documentation

◆ AlongStepDoIt()

virtual G4VParticleChange * G4DNASecondOrderReaction::AlongStepDoIt ( const G4Track ,
const G4Step  
)
inlinevirtual

Implements G4VProcess.

Definition at line 102 of file G4DNASecondOrderReaction.hh.

105 {return 0;}

◆ AlongStepGetPhysicalInteractionLength()

virtual G4double G4DNASecondOrderReaction::AlongStepGetPhysicalInteractionLength ( const G4Track ,
G4double  ,
G4double  ,
G4double ,
G4GPILSelection  
)
inlinevirtual

Implements G4VProcess.

Definition at line 93 of file G4DNASecondOrderReaction.hh.

99 { return -1.0; }

◆ AtRestDoIt()

virtual G4VParticleChange * G4DNASecondOrderReaction::AtRestDoIt ( const G4Track ,
const G4Step  
)
inlinevirtual

Implements G4VProcess.

Definition at line 87 of file G4DNASecondOrderReaction.hh.

90 {return 0;}

◆ AtRestGetPhysicalInteractionLength()

virtual G4double G4DNASecondOrderReaction::AtRestGetPhysicalInteractionLength ( const G4Track ,
G4ForceCondition  
)
inlinevirtual

Implements G4VProcess.

Definition at line 82 of file G4DNASecondOrderReaction.hh.

85 { return -1.0; }

◆ BuildPhysicsTable()

void G4DNASecondOrderReaction::BuildPhysicsTable ( const G4ParticleDefinition )
virtual

Reimplemented from G4VITProcess.

Definition at line 97 of file G4DNASecondOrderReaction.cc.

98{
100 fMolarMassOfMaterial = fpMaterial->GetMassOfMolecule()*CLHEP::Avogadro*1e3;
101 fIsInitialized = true;
102}
const std::vector< G4double > * GetNumMolPerVolTableFor(const G4Material *) const
Retrieve a table of molecular densities (number of molecules per unit volume) in the G4 unit system f...
static G4DNAMolecularMaterial * Instance()
const std::vector< double > * fpMoleculeDensity
G4double GetMassOfMolecule() const
Definition: G4Material.hh:239

◆ operator=()

G4DNASecondOrderReaction & G4DNASecondOrderReaction::operator= ( const G4DNASecondOrderReaction rhs)

Definition at line 83 of file G4DNASecondOrderReaction.cc.

84{
85 if (this == &rhs) return *this; // handle self assignment
86
87 //assignment operator
88 return *this;
89}

◆ PostStepDoIt()

G4VParticleChange * G4DNASecondOrderReaction::PostStepDoIt ( const G4Track track,
const G4Step  
)
virtual

Implements G4VProcess.

Definition at line 244 of file G4DNASecondOrderReaction.cc.

245{
246 G4Molecule* molecule = GetMolecule(track);
247#ifdef G4VERBOSE
248 if(verboseLevel > 1)
249 {
250 G4cout << "___________" << G4endl;
251 G4cout << ">>> Beginning of G4DNASecondOrderReaction verbose" << G4endl;
252 G4cout << ">>> Returned value : " << G4BestUnit(fReturnedValue,"Time") << G4endl;
253 G4cout << ">>> Time Step : " << G4BestUnit(G4VScheduler::Instance()->GetTimeStep(),"Time") << G4endl;
254 G4cout << ">>> Reaction : " << molecule->GetName() << " + " << fpMaterial->GetName() << G4endl;
255 G4cout << ">>> End of G4DNASecondOrderReaction verbose <<<" << G4endl;
256 }
257#endif
262 State(fPreviousTimeAtPreStepPoint) = -1;
263 return &fParticleChange;
264}
#define State(X)
G4Molecule * GetMolecule(const G4Track &track)
Definition: G4Molecule.cc:76
#define G4BestUnit(a, b)
@ fStopAndKill
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
virtual void AddIndirectDamage(const G4String &baseName, const G4Molecule *molecule, const G4ThreeVector &position, double time)
Definition: G4DNADamage.cc:95
static G4DNADamage * Instance()
Definition: G4DNADamage.cc:56
const G4String & GetName() const
Definition: G4Material.hh:175
const G4String & GetName() const
Definition: G4Molecule.cc:338
virtual void Initialize(const G4Track &)
const G4ThreeVector & GetPosition() const
G4double GetGlobalTime() const
void ProposeTrackStatus(G4TrackStatus status)
G4int verboseLevel
Definition: G4VProcess.hh:356
static G4VScheduler * Instance()
Definition: G4VScheduler.cc:48
#define DBL_MAX
Definition: templates.hh:62

◆ PostStepGetPhysicalInteractionLength()

G4double G4DNASecondOrderReaction::PostStepGetPhysicalInteractionLength ( const G4Track track,
G4double  previousStepSize,
G4ForceCondition condition 
)
virtual

Implements G4VProcess.

Definition at line 129 of file G4DNASecondOrderReaction.cc.

132{
133 // G4cout << "G4DNASecondOrderReaction::PostStepGetPhysicalInteractionLength" << G4endl;
134 // G4cout << "For reaction : " << fpMaterial->GetName() << " + " << fpMolecularConfiguration->GetName() << G4endl;
135
136 //_______________________________________________________________________
137 // Check whether the track is in the good material (maybe composite material)
138 const G4Material* material = track.GetMaterial();
139
140 G4Molecule* mol = GetMolecule(track);
141 if(!mol) return DBL_MAX;
143 {
144 // G4cout <<"mol->GetMolecularConfiguration() != fpMolecularConfiguration" << G4endl;
145 return DBL_MAX;
146 }
147
148 G4double molDensity = (*fpMoleculeDensity)[material->GetIndex()];
149
150 if(molDensity == 0.0) // ie : not found
151 {
152 if(State(fIsInGoodMaterial))
153 {
155 //State(fPreviousTimeAtPreStepPoint) = -1;
156 State(fIsInGoodMaterial) = false;
157 }
158
159 // G4cout << " Material " << fpMaterial->GetName() << " not found "
160 // <<" | name of current material : " << material->GetName()
161 // << G4endl;
162
163 return DBL_MAX; // Becareful return here !!
164 }
165
166 // G4cout << " Va calculer le temps d'interaction " << G4endl;
167
168 State(fIsInGoodMaterial) = true;
169
170 // fConcentration = molDensity/fMolarMassOfMaterial;
171 fConcentration = molDensity/CLHEP::Avogadro;
172 // G4cout << "Concentration : " << fConcentration / (g/mole)<< G4endl;
173
174 //_______________________________________________________________________
175 // Either initialize the lapse of time left
176 // meaning
177 // => the track enters for the first time in the material
178 // or substract the previous time step to the previously calculated lapse
179 // of time left
180 // meaning
181 // => the track has not left this material since the previous call
182 G4double previousTimeStep(-1.);
183
184 if(State(fPreviousTimeAtPreStepPoint) != -1)
185 {
186 previousTimeStep = track.GetGlobalTime() -
187 State(fPreviousTimeAtPreStepPoint) ;
188 }
189
190 State(fPreviousTimeAtPreStepPoint) = track.GetGlobalTime();
191
192 // condition is set to "Not Forced"
193 *pForceCond = NotForced;
194
195 if (
196 (previousTimeStep < 0.0) ||
197 (fpState->theNumberOfInteractionLengthLeft<=0.0)) {
198 // beggining of tracking (or just after DoIt of this process)
200 } else if ( previousTimeStep > 0.0) {
201 // get mean free path
202 // subtract NumberOfInteractionLengthLeft
203 SubtractNumberOfInteractionLengthLeft( previousTimeStep );
204 } else {
205 // zero time step
206 // Force trigerring the process
207 //*pForceCond = Forced;
208 }
209
210 fpState->currentInteractionLength = 1/(fReactionRate*fConcentration);
211
212 // G4cout << "fpState->currentInteractionLength = "
213 // << fpState->currentInteractionLength << G4endl;
214
215 G4double value;
216 if (fpState->currentInteractionLength <DBL_MAX) {
217 value = fpState->theNumberOfInteractionLengthLeft
218 * (fpState->currentInteractionLength);
219 } else {
220 value = DBL_MAX;
221 }
222#ifdef G4VERBOSE
223 if (verboseLevel>2) {
224 G4cout << "G4VITRestDiscreteProcess::PostStepGetPhysicalInteractionLength ";
225 G4cout << "[ " << GetProcessName() << "]" <<G4endl;
226 track.GetDynamicParticle()->DumpInfo();
227 G4cout << " in Material " << track.GetMaterial()->GetName() <<G4endl;
228 G4cout << "InteractionLength= " << value/cm <<"[cm] " <<G4endl;
229 }
230#endif
231
232// G4cout << "currentInteractionLength : " << fpState->currentInteractionLength << G4endl;
233// G4cout << "Material : " << fpMaterial->GetName()
234// << "ID: " << track.GetTrackID()
235// << " Returned time : " << G4BestUnit(value,"Time") << G4endl;
236
237 if(value < fReturnedValue)
238 fReturnedValue = value;
239
240 return value*-1;
241 // multiple by -1 to indicate to the tracking system that we are returning a time
242}
@ NotForced
double G4double
Definition: G4Types.hh:83
const G4MolecularConfiguration * fpMolecularConfiguration
void DumpInfo(G4int mode=0) const
size_t GetIndex() const
Definition: G4Material.hh:258
const G4MolecularConfiguration * GetMolecularConfiguration() const
Definition: G4Molecule.cc:532
G4Material * GetMaterial() const
const G4DynamicParticle * GetDynamicParticle() const
virtual void SubtractNumberOfInteractionLengthLeft(G4double previousStepSize)
G4shared_ptr< G4ProcessState > fpState
virtual void ResetNumberOfInteractionLengthLeft()
const G4String & GetProcessName() const
Definition: G4VProcess.hh:382

◆ SetReaction()

void G4DNASecondOrderReaction::SetReaction ( const G4MolecularConfiguration molConf,
const G4Material mat,
double  reactionRate 
)

Definition at line 113 of file G4DNASecondOrderReaction.cc.

115{
117 {
118 G4ExceptionDescription exceptionDescription ;
119 exceptionDescription << "G4DNASecondOrderReaction was already initialised. ";
120 exceptionDescription << "You cannot set a reaction after initialisation.";
121 G4Exception("G4DNASecondOrderReaction::SetReaction","G4DNASecondOrderReaction001",
122 FatalErrorInArgument,exceptionDescription);
123 }
124 fpMolecularConfiguration = molConf;
125 fpMaterial = mat;
126 fReactionRate = reactionRate;
127}
@ FatalErrorInArgument
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40

◆ StartTracking()

void G4DNASecondOrderReaction::StartTracking ( G4Track track)
virtual

Reimplemented from G4VITProcess.

Definition at line 105 of file G4DNASecondOrderReaction.cc.

106{
108 G4VITProcess::fpState.reset(new SecondOrderReactionState());
110}
virtual void StartTracking(G4Track *)
Definition: G4VITProcess.cc:85
virtual void StartTracking(G4Track *)
Definition: G4VProcess.cc:87

Member Data Documentation

◆ fConcentration

G4double G4DNASecondOrderReaction::fConcentration
protected

Definition at line 126 of file G4DNASecondOrderReaction.hh.

Referenced by PostStepGetPhysicalInteractionLength().

◆ fIsInitialized

G4bool G4DNASecondOrderReaction::fIsInitialized
protected

Definition at line 120 of file G4DNASecondOrderReaction.hh.

Referenced by BuildPhysicsTable(), and SetReaction().

◆ fMolarMassOfMaterial

G4double G4DNASecondOrderReaction::fMolarMassOfMaterial
protected

Definition at line 127 of file G4DNASecondOrderReaction.hh.

Referenced by BuildPhysicsTable().

◆ fParticleChange

G4ParticleChange G4DNASecondOrderReaction::fParticleChange
protected

Definition at line 128 of file G4DNASecondOrderReaction.hh.

Referenced by PostStepDoIt().

◆ fpMaterial

const G4Material* G4DNASecondOrderReaction::fpMaterial
protected

Definition at line 131 of file G4DNASecondOrderReaction.hh.

Referenced by BuildPhysicsTable(), PostStepDoIt(), and SetReaction().

◆ fpMolecularConfiguration

const G4MolecularConfiguration* G4DNASecondOrderReaction::fpMolecularConfiguration
protected

◆ fpMoleculeDensity

const std::vector<double>* G4DNASecondOrderReaction::fpMoleculeDensity
protected

Definition at line 124 of file G4DNASecondOrderReaction.hh.

Referenced by BuildPhysicsTable().

◆ fReactionRate

G4double G4DNASecondOrderReaction::fReactionRate
protected

◆ fReturnedValue

G4double G4DNASecondOrderReaction::fReturnedValue
protected

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