include/Geant4/G4DNASecondOrderReaction.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // Author: Mathieu Karamitros
0028
0029 // The code is developed in the framework of the ESA AO7146
0030 //
0031 // We would be very happy hearing from you, send us your feedback! :)
0032 //
0033 // In order for Geant4-DNA to be maintained and still open-source,
0034 // article citations are crucial.
0035 // If you use Geant4-DNA chemistry and you publish papers about your software,
0036 // in addition to the general paper on Geant4-DNA:
0037 //
0038 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
0039 //
0040 // we would be very happy if you could please also cite the following
0041 // reference papers on chemistry:
0042 //
0043 // J. Comput. Phys. 274 (2014) 841-882
0044 // Prog. Nucl. Sci. Tec. 2 (2011) 503-508
0045
0046
0047 #ifndef G4DNASECONDORDERREACTION_HH
0048 #define G4DNASECONDORDERREACTION_HH
0049
0050 #include "G4VITProcess.hh"
0051
0052 class G4MolecularConfiguration;
0053
0054 class G4DNASecondOrderReaction : public G4VITProcess
0055 {
0056 public:
0057     G4DNASecondOrderReaction(const G4String& aName =  "DNASecondOrderReaction",
0058                              G4ProcessType type = fDecay);
0059     ~G4DNASecondOrderReaction() override;
0060
0061     G4IT_ADD_CLONE(G4VITProcess,G4DNASecondOrderReaction)
0062
0063     G4DNASecondOrderReaction(const G4DNASecondOrderReaction&);
0064     G4DNASecondOrderReaction& operator=(const G4DNASecondOrderReaction&);
0065     void StartTracking(G4Track*) override;
0066
0067     void SetReaction(const G4MolecularConfiguration*, const G4Material*, double /*reactionRate*/);
0068
0069 public :
0070    void BuildPhysicsTable(const G4ParticleDefinition&) override;
0071    G4double PostStepGetPhysicalInteractionLength(
0072                            const G4Track& track,
0073                G4double   previousStepSize,
0074                G4ForceCondition* condition
0075               ) override;
0076
0077    G4VParticleChange* PostStepDoIt(
0078                const G4Track& ,
0079                const G4Step&
0080               ) override;
0081
0082    G4double AtRestGetPhysicalInteractionLength(
0083                            const G4Track& ,
0084                G4ForceCondition*
0085               ) override{ return -1.0; }
0086
0087    G4VParticleChange* AtRestDoIt(
0088                const G4Track& ,
0089                const G4Step&
0090               ) override{return nullptr;}
0091
0092    //  no operation in  AlongStepDoIt
0093    G4double AlongStepGetPhysicalInteractionLength(
0094                            const G4Track&,
0095                G4double  ,
0096                G4double  ,
0097                G4double& ,
0098                            G4GPILSelection*
0099                           ) override{ return -1.0; }
0100
0101    //  no operation in  AlongStepDoIt
0102    G4VParticleChange* AlongStepDoIt(
0103                const G4Track& ,
0104                const G4Step&
0105                           ) override {return nullptr;}
0106
0107 protected:
0108     struct SecondOrderReactionState : public G4ProcessState
0109     {
0110         SecondOrderReactionState();
0111         ~SecondOrderReactionState() override= default;
0112         G4double fPreviousTimeAtPreStepPoint;
0113         G4bool fIsInGoodMaterial;
0114     };
0115
0116 private :
0117     void Create();
0118
0119 protected:
0120     G4bool fIsInitialized;
0121
0122     G4double fReturnedValue;
0123
0124     const std::vector<double>* fpMoleculeDensity;
0125     G4double fReactionRate;
0126     G4double fConcentration;
0127     G4double fMolarMassOfMaterial;
0128     G4ParticleChange fParticleChange;
0129
0130     const G4MolecularConfiguration* fpMolecularConfiguration;
0131     const G4Material* fpMaterial;
0132 };
0133
0134 #endif // G4DNASECONDORDERREACTION_HH