00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // $Id$ 00027 // 00028 // 00029 //J.M. Quesada (August2008). Based on: 00030 // 00031 // Hadronic Process: Nuclear De-excitations 00032 // by V. Lara (Oct 1998) 00033 // 00034 // 17-11-2010 V.Ivanchenko in constructor replace G4VEmissionProbability by 00035 // G4EvaporationProbability and do not new and delete probability 00036 // object at each call; use G4Pow 00037 00038 #ifndef G4EvaporationChannel_h 00039 #define G4EvaporationChannel_h 1 00040 00041 #include "G4VEvaporationChannel.hh" 00042 #include "G4EvaporationProbability.hh" 00043 #include "G4VCoulombBarrier.hh" 00044 00045 class G4EvaporationLevelDensityParameter; 00046 00047 class G4EvaporationChannel : public G4VEvaporationChannel 00048 { 00049 public: 00050 // constructor 00051 G4EvaporationChannel(G4int theA, G4int theZ, const G4String & aName, 00052 G4EvaporationProbability * aEmissionStrategy, 00053 G4VCoulombBarrier * aCoulombBarrier); 00054 00055 // destructor 00056 virtual ~G4EvaporationChannel(); 00057 00058 inline void SetEmissionStrategy(G4EvaporationProbability * aEmissionStrategy) 00059 {theEvaporationProbabilityPtr = aEmissionStrategy;} 00060 00061 inline void SetCoulombBarrierStrategy(G4VCoulombBarrier * aCoulombBarrier) 00062 {theCoulombBarrierPtr = aCoulombBarrier;} 00063 00064 protected: 00065 // default constructor 00066 G4EvaporationChannel(); 00067 00068 public: 00069 00070 // virtual void Initialize(const G4Fragment & fragment); 00071 00072 virtual G4double GetEmissionProbability(G4Fragment* fragment); 00073 00074 G4FragmentVector * BreakUp(const G4Fragment & theNucleus); 00075 00076 inline G4double GetMaximalKineticEnergy(void) const 00077 { return MaximalKineticEnergy; } 00078 00079 private: 00080 00081 // Calculate Binding Energy for separate fragment from nucleus 00082 G4double CalcBindingEnergy(G4int anA, G4int aZ); 00083 00084 // Calculate maximal kinetic energy that can be carried by fragment (in MeV) 00085 G4double CalcMaximalKineticEnergy(G4double U); 00086 00087 // Samples fragment kinetic energy. 00088 G4double GetKineticEnergy(const G4Fragment & aFragment); 00089 00090 // This has to be removed and put in Random Generator 00091 G4ThreeVector IsotropicVector(G4double Magnitude = 1.0); 00092 00093 G4EvaporationChannel(const G4EvaporationChannel & right); 00094 const G4EvaporationChannel & operator=(const G4EvaporationChannel & right); 00095 G4bool operator==(const G4EvaporationChannel & right) const; 00096 G4bool operator!=(const G4EvaporationChannel & right) const; 00097 00098 // Data Members 00099 // ************ 00100 private: 00101 00102 // This data member define the channel. 00103 // They are intializated at object creation (constructor) time. 00104 00105 // Atomic Number of ejectile 00106 G4int theA; 00107 00108 // Charge of ejectile 00109 G4int theZ; 00110 00111 G4double EvaporatedMass; 00112 G4double ResidualMass; 00113 00114 // For evaporation probability calcualation 00115 G4EvaporationProbability * theEvaporationProbabilityPtr; 00116 00117 // For Level Density calculation 00118 // G4bool MyOwnLevelDensity; 00119 G4VLevelDensityParameter * theLevelDensityPtr; 00120 00121 // For Coulomb Barrier calculation 00122 G4VCoulombBarrier * theCoulombBarrierPtr; 00123 G4double CoulombBarrier; 00124 00125 //--------------------------------------------------- 00126 00127 // These values depend on the nucleus that is being evaporated. 00128 // They are calculated through the Initialize method which takes as parameters 00129 // the atomic number, charge and excitation energy of nucleus. 00130 00131 // Residual Mass Number 00132 G4int ResidualA; 00133 00134 // Residual Charge 00135 G4int ResidualZ; 00136 00137 // Emission Probability 00138 G4double EmissionProbability; 00139 00140 // Maximal Kinetic Energy that can be carried by fragment 00141 G4double MaximalKineticEnergy; 00142 00143 }; 00144 00145 00146 #endif