84 throw G4HadronicException(__FILE__, __LINE__,
"G4NeutronHPDiscreteTwoBody: Unknown ion case 2");
93 if(
theCoeff[i].GetEnergy()>anEnergy)
break;
95 if(it==0||it==nEnergy-1)
97 if(
theCoeff[it].GetRepresentation()==0)
107 else if(
theCoeff[it].GetRepresentation()==12)
122 cosTh = theStore.
Sample();
124 else if(
theCoeff[it].GetRepresentation()==14)
139 cosTh = theStore.
Sample();
143 throw G4HadronicException(__FILE__, __LINE__,
"unknown representation type in Two-body scattering");
150 if(
theCoeff[it].GetRepresentation()==0)
163 else if(
theCoeff[it].GetRepresentation()==12)
205 mu = theBuff1.
GetX(i);
206 y1 = theBuff1.
GetY(i);
207 y2 = theBuff2.
GetY(mu);
213 mu = theBuff2.
GetX(i);
214 y1 = theBuff2.
GetY(i);
215 y2 = theBuff1.
GetY(mu);
219 theStore.
Merge(&theStore1, &theStore2);
220 cosTh = theStore.
Sample();
222 else if(
theCoeff[it].GetRepresentation()==14)
266 mu = theBuff1.
GetX(i);
267 y1 = theBuff1.
GetY(i);
268 y2 = theBuff2.
GetY(mu);
274 mu = theBuff2.
GetX(i);
275 y1 = theBuff2.
GetY(i);
276 y2 = theBuff1.
GetY(mu);
280 theStore.
Merge(&theStore1, &theStore2);
281 cosTh = theStore.
Sample();
285 throw G4HadronicException(__FILE__, __LINE__,
"Two neighbouring distributions with different interpolation");
290 throw G4HadronicException(__FILE__, __LINE__,
"unknown representation type in Two-body scattering, case 2");
306 G4double E1 = (A1+1)*(A1+1)/A1/A1*anEnergy;
314 G4ThreeVector tempVector(mtot*sinth*std::cos(phi), mtot*sinth*std::sin(phi), mtot*std::cos(theta) );
G4ReactionProduct * GetNeutron()
G4double GetY(G4double x)
G4int GetVectorLength() const
void SetInterpolationManager(const G4InterpolationManager &aManager)
G4double GetTotalMomentum() const
CLHEP::Hep3Vector G4ThreeVector
void Init(G4int aScheme, G4int aRange)
void SetKineticEnergy(const G4double en)
void SetMomentum(const G4double x, const G4double y, const G4double z)
void SetCoeff(G4int i, G4int l, G4double coeff)
void Merge(G4NeutronHPVector *active, G4NeutronHPVector *passive)
G4InterpolationManager theManager
G4double GetX(G4int i) const
void SetData(G4int i, G4double x, G4double y)
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4ReactionProduct * Sample(G4double anEnergy, G4double massCode, G4double mass)
G4NeutronHPLegendreTable * theCoeff
void SetY(G4int i, G4double x)
G4InterpolationScheme GetScheme(G4int index) const
static G4Triton * Triton()
G4double SampleDiscreteTwoBody(G4double anEnergy)
G4ReactionProduct * GetTarget()
static G4Proton * Proton()
void SetX(G4int i, G4double e)
static G4Neutron * Neutron()
static const G4double A[nN]
static G4Deuteron * Deuteron()
static G4Positron * Positron()
void SetManager(G4InterpolationManager &aManager)
static G4Electron * Electron()
G4NeutronHPInterpolator theInt