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 #ifndef G4DiffractiveExcitation_h
00030 #define G4DiffractiveExcitation_h 1
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041 #include "globals.hh"
00042 class G4VSplitableHadron;
00043 class G4ExcitedString;
00044 #include "G4FTFParameters.hh"
00045 #include "G4ElasticHNScattering.hh"
00046 #include "G4ThreeVector.hh"
00047
00048 class G4DiffractiveExcitation
00049 {
00050
00051 public:
00052
00053 G4DiffractiveExcitation();
00054 virtual ~G4DiffractiveExcitation();
00055
00056 virtual G4bool ExciteParticipants (G4VSplitableHadron *aPartner,
00057 G4VSplitableHadron * bPartner,
00058 G4FTFParameters *theParameters,
00059 G4ElasticHNScattering *theElastic) const;
00060
00061 virtual void CreateStrings (G4VSplitableHadron * aHadron,
00062 G4bool isProjectile,
00063 G4ExcitedString * &FirstString,
00064 G4ExcitedString * &SecondString,
00065 G4FTFParameters *theParameters) const;
00066 private:
00067
00068 G4DiffractiveExcitation(const G4DiffractiveExcitation &right);
00069
00070 G4ThreeVector GaussianPt(G4double AveragePt2, G4double maxPtSquare) const;
00071 G4double ChooseP(G4double Pmin, G4double Pmax) const;
00072 G4double GetQuarkFractionOfKink(G4double zmin, G4double zmax) const;
00073 void UnpackMeson( G4int IdPDG, G4int &Q1, G4int &Q2) const;
00074 void UnpackBaryon(G4int IdPDG, G4int &Q1, G4int &Q2, G4int &Q3) const;
00075 G4int NewNucleonId(G4int Q1, G4int Q2, G4int Q3) const;
00076
00077 const G4DiffractiveExcitation & operator=(const G4DiffractiveExcitation &right);
00078 int operator==(const G4DiffractiveExcitation &right) const;
00079 int operator!=(const G4DiffractiveExcitation &right) const;
00080
00081 };
00082
00083 #endif