Changeset 1340 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4HETCTriton.hh
- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4HETCTriton.hh
r819 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4HETCTriton.hh,v 1.3 2010/08/28 15:16:55 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-03-ref-09 $ 28 // 26 29 // by V. Lara 30 // 31 // Modified: 32 // 23.08.2010 V.Ivanchenko general cleanup, move constructor and destructor 33 // the source, use G4Pow 27 34 28 35 #ifndef G4HETCTriton_h … … 31 38 #include "G4HETCChargedFragment.hh" 32 39 #include "G4ReactionProduct.hh" 33 #include "G4Triton.hh"34 35 40 #include "G4TritonCoulombBarrier.hh" 36 37 41 38 42 class G4HETCTriton : public G4HETCChargedFragment 39 43 { 40 44 public: 41 // default constructor42 G4HETCTriton():G4HETCChargedFragment(3,1,&theTritonCoulombBarrier,"Triton") {}43 45 44 // copy constructor 45 G4HETCTriton(const G4HETCTriton &right): G4HETCChargedFragment(right) {} 46 G4HETCTriton(); 46 47 47 // destructor 48 ~G4HETCTriton() {} 48 ~G4HETCTriton(); 49 50 protected: 51 52 virtual G4double GetAlpha(); 53 54 virtual G4double GetBeta(); 55 56 virtual G4double GetSpinFactor(); 57 58 virtual G4double K(const G4Fragment & aFragment); 59 60 private: 49 61 50 62 // operators 51 const G4HETCTriton & operator=(const G4HETCTriton &right) 52 { 53 if (&right != this) this->G4HETCChargedFragment::operator=(right); 54 return *this; 55 } 63 G4HETCTriton(const G4HETCTriton &right); 64 const G4HETCTriton & operator=(const G4HETCTriton &right); 65 G4bool operator==(const G4HETCTriton &right) const; 66 G4bool operator!=(const G4HETCTriton &right) const; 56 67 57 G4bool operator==(const G4HETCTriton &right) const58 {59 return G4HETCChargedFragment::operator==(right);60 }61 62 63 G4bool operator!=(const G4HETCTriton &right) const64 {65 return G4HETCChargedFragment::operator!=(right);66 }67 68 69 G4ReactionProduct * GetReactionProduct() const70 {71 G4ReactionProduct * theReactionProduct =72 new G4ReactionProduct(G4Triton::TritonDefinition());73 theReactionProduct->SetMomentum(GetMomentum().vect());74 theReactionProduct->SetTotalEnergy(GetMomentum().e());75 #ifdef PRECOMPOUND_TEST76 theReactionProduct->SetCreatorModel("G4PrecompoundModel");77 #endif78 return theReactionProduct;79 }80 81 private:82 virtual G4double GetAlpha()83 {84 G4double C = 0.0;85 G4double aZ = GetZ() + GetRestZ();86 if (aZ >= 70)87 {88 C = 0.10;89 }90 else91 {92 C = ((((0.15417e-06*aZ) - 0.29875e-04)*aZ + 0.21071e-02)*aZ - 0.66612e-01)*aZ + 0.98375;93 }94 95 return 1.0 + C/3.0;96 }97 98 virtual G4double GetBeta()99 {100 return -GetCoulombBarrier();101 }102 103 104 105 virtual G4double GetSpinFactor()106 {107 // 2s+1108 return 2.0;109 }110 111 virtual G4double K(const G4Fragment& aFragment);112 113 private:114 115 68 G4TritonCoulombBarrier theTritonCoulombBarrier; 116 117 69 }; 118 70
Note: See TracChangeset
for help on using the changeset viewer.