Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 // 28 /// \file TLKModel.hh 29 /// \brief Definition of the TLKModel class 30 31 #ifndef TLKMODEL_HH 32 #define TLKMODEL_HH 33 34 //#include <boost/numeric/odeint.hpp> 35 #include <vector> 36 #include <string> 37 class Damage; 38 class TLKModel 39 { 40 public: 41 42 /// \brief constructor 43 TLKModel(double pLambda1 =-1,double pLambda2=-1, double pBeta1=-1, double pBeta2=-1,double pEta=-1); 44 /// \brief destructor 45 ~TLKModel() = default; 46 47 double GetLambda1(){return fLambda1;}; 48 void SetLambda1(double pVal){fLambda1 = pVal;}; 49 double GetLambda2(){return fLambda2;}; 50 void SetLambda2(double pVal){fLambda2 = pVal;}; 51 double GetBeta1(){return fBeta1;}; 52 void SetBeta1(double pVal){fBeta1 = pVal;}; 53 double GetBeta2(){return fBeta2;}; 54 void SetBeta2(double pVal){fBeta2 = pVal;}; 55 double GetEta(){return fEta;}; 56 void SetEta(double pVal){fEta = pVal;}; 57 58 double GetSingleDSBYield(){return fSingleDSBYield;}; 59 void SetSingleDSBYield(double pVal){fSingleDSBYield = pVal;}; 60 61 double GetComplexDSBYield(){return fComplexDSBYield;}; 62 void SetComplexDSBYield(double pVal){fComplexDSBYield = pVal;}; 63 64 unsigned int GetBpForDSB(){return fBpForDSB;}; 65 void SetBpForDSB(unsigned int pVal){fBpForDSB = pVal;}; 66 67 double GetStartTime(){return fStartTime;}; 68 void SetStartTime(double pVal){fStartTime = pVal;}; 69 double GetStopTime(){return fStopTime;}; 70 void SetStopTime(double pVal){fStopTime = pVal;}; 71 double GetStepTime(){return fStepTime;}; 72 void SetStepTime(double pVal){fStepTime = pVal;}; 73 74 // Compute damage inputs required by TLK model 75 void ComputeAndSetDamageInput(std::vector<Damage>); 76 77 // Compute SF for a given absorbed dose expressed in Gy 78 double ComputeSF(double pDose); 79 80 // Compute a SF Curve 81 void CalculateRepair(double pDoseMax, double pDeltaDose); 82 83 // Write output, dose expressed in Gy 84 void WriteOutput(std::string pFileName); 85 86 void SetDose(double d) {fDose = d;} 87 private: 88 std::vector<double> TLK_odes_system(double t,std::vector<double> y); 89 // TLK parameters 90 // simple DSB repair probability (h-1) 91 double fLambda1{0}; 92 // complex DSB repair probability (h-1) 93 double fLambda2{0}; 94 // simple DSB misrepair probability (h-1) 95 double fBeta1{0}; 96 // complex DSB misrepair probability (h-1) 97 double fBeta2{0}; 98 // binary misrepair probability (h-1) 99 double fEta{0}; 100 101 // DSB Yields in Gy-1 102 double fSingleDSBYield{0}; 103 double fComplexDSBYield{0}; 104 105 // num of bp to consider a DSB for the MakeCluster function 106 // default value is 10 107 unsigned int fBpForDSB{10}; 108 109 // Time for integration 110 double fStartTime{0}; 111 double fStopTime{0}; 112 double fStepTime{0}; 113 114 // SF curve 115 std::vector<std::pair<double,double>> fSFCurve; 116 // store dose deposited in nucleus cell 117 double fDose{0}; 118 }; 119 120 #endif