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 #include "globals.hh"
00030 #include "G4KineticTrack.hh"
00031 #include "G4VCrossSectionSource.hh"
00032 #include "G4Proton.hh"
00033 #include "G4Neutron.hh"
00034 #include "G4XAqmElastic.hh"
00035 #include "G4AngularDistribution.hh"
00036 #include "G4ThreeVector.hh"
00037 #include "G4LorentzVector.hh"
00038 #include "G4LorentzRotation.hh"
00039 #include "G4KineticTrackVector.hh"
00040 #include "G4XResonance.hh"
00041 #include "G4ConcreteNNTwoBodyResonance.hh"
00042
00043 G4ConcreteNNTwoBodyResonance::G4ConcreteNNTwoBodyResonance(const G4ParticleDefinition* aPrimary,
00044 const G4ParticleDefinition* bPrimary,
00045 const G4ParticleDefinition* aSecondary,
00046 const G4ParticleDefinition* bSecondary,
00047 const G4VXResonanceTable& sigmaTable)
00048 : thePrimary1(aPrimary), thePrimary2(bPrimary)
00049 {
00050 theOutGoing.push_back(aSecondary);
00051 theOutGoing.push_back(bSecondary);
00052
00053 crossSectionSource = new G4XResonance(aPrimary, bPrimary,
00054 aSecondary->GetPDGiIsospin(),
00055 aSecondary->GetPDGiSpin(),
00056 aSecondary->GetPDGMass(),
00057 bSecondary->GetPDGiIsospin(),
00058 bSecondary->GetPDGiSpin(),
00059 bSecondary->GetPDGMass(),
00060 aSecondary->GetParticleName(),
00061 bSecondary->GetParticleName(),
00062 sigmaTable);
00063 }
00064
00065 G4ConcreteNNTwoBodyResonance::~G4ConcreteNNTwoBodyResonance()
00066 {
00067 delete crossSectionSource;
00068 }
00069
00070 G4bool G4ConcreteNNTwoBodyResonance::IsInCharge(const G4KineticTrack& trk1,
00071 const G4KineticTrack& trk2) const
00072 {
00073 if (trk1.GetDefinition()==thePrimary1 && trk2.GetDefinition()==thePrimary2) return true;
00074 if (trk1.GetDefinition()==thePrimary2 && trk2.GetDefinition()==thePrimary1) return true;
00075 return false;
00076 }
00077