00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032 #ifndef G4StatMFChannel_h
00033 #define G4StatMFChannel_h 1
00034
00035 #include <deque>
00036
00037 #include "G4StatMFParameters.hh"
00038 #include "G4StatMFFragment.hh"
00039
00040
00041 class G4StatMFChannel {
00042
00043 public:
00044
00045 G4StatMFChannel();
00046
00047
00048 ~G4StatMFChannel();
00049
00050 private:
00051
00052
00053 G4StatMFChannel(const G4StatMFChannel & right);
00054
00055
00056 G4StatMFChannel & operator=(const G4StatMFChannel & right);
00057
00058 G4bool operator==(const G4StatMFChannel & right) const;
00059 G4bool operator!=(const G4StatMFChannel & right) const;
00060
00061 public:
00062
00063 void CreateFragment(G4int A, G4int Z);
00064
00065 inline size_t GetMultiplicity(void) { return _theFragments.size();}
00066
00067
00068 G4bool CheckFragments(void);
00069
00070 G4double GetFragmentsCoulombEnergy(void);
00071
00072 G4double GetFragmentsEnergy(G4double T) const;
00073
00074 G4FragmentVector * GetFragments(G4int anA, G4int anZ, G4double T);
00075
00076 private:
00077
00078
00079 void CoulombImpulse(G4int anA, G4int anZ, G4double T);
00080
00081 void PlaceFragments(G4int anA);
00082
00083 void SolveEqOfMotion(G4int anA, G4int anZ, G4double T);
00084
00085
00086
00087
00088 void FragmentsMomenta(G4int NF, G4int idx, G4double T);
00089
00090
00091
00092
00093 G4ThreeVector IsotropicVector(G4double Magnitude = 1.0);
00094
00095
00096
00097 G4ThreeVector RotateMomentum(G4ThreeVector Pa, G4ThreeVector V,
00098 G4ThreeVector P);
00099
00100 private:
00101
00102 std::deque<G4StatMFFragment*> _theFragments;
00103
00104 G4int _NumOfNeutralFragments;
00105
00106 G4int _NumOfChargedFragments;
00107
00108 struct DeleteFragment
00109 {
00110 template<typename T>
00111 void operator()(const T* ptr) const
00112 {
00113 delete ptr;
00114 }
00115 };
00116
00117 };
00118
00119 #endif
00120
00121
00122
00123
00124
00125
00126