34 #define INCLXX_IN_GEANT4_MODE 1
52 : particle1(p1), particle2(p2)
58 const G4double ecm = ecmorigin - 783.437;
65 G4double q2=(y-1.157776E6)*(y-6.4E5)/y/4.0;
66 G4double q3=std::pow(std::sqrt(q2), 3.);
67 const G4double f3max=q3/(q3+5.832E6);
74 INCL_WARN(
"NDeltaOmegaProductionChannel::sampleDeltaMass loop was stopped because maximum number of tries was reached. Minimum delta mass "
87 q2=(y-1.157776E6)*(y-6.4E5)/y/4.0;
88 q3=std::pow(std::sqrt(q2), 3.);
116 const G4int isospin = is1+is2;
124 if (rndm < 0.5) index2=1;
167 const ThreeVector rcol = (rcolnucleon1+rcolnucleon2)*0.5;
170 list.push_back(omega);
void addCreatedParticle(Particle *p)
void setMass(G4double mass)
Set the mass of the particle in MeV/c^2.
static const G4int maxTries
G4double dot(const ThreeVector &v) const
Dot product.
G4ThreadLocal G4double minDeltaMassRndm
virtual ~NDeltaOmegaProductionChannel()
const G4INCL::ThreeVector & getMomentum() const
Get the momentum vector.
G4bool isDelta() const
Is it a Delta?
NDeltaOmegaProductionChannel(Particle *, Particle *)
Final state of an interaction.
void fillFinalState(FinalState *fs)
void generateBiased(const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
Generate a biased event in the CM system.
G4double sampleDeltaMass(G4double ecm)
const G4double effectiveDeltaWidth
G4ThreadLocal G4double minDeltaMass
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
G4INCL::ParticleType getType() const
Get the particle type.
const G4INCL::ThreeVector & getPosition() const
Set the position vector.
void setType(ParticleType t)
G4double shoot()
Generate flat distribution of random numbers.
G4double mag() const
Get the length of the vector.
const G4double effectiveDeltaMass
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
static const G4double angularSlope
const G4double effectiveNucleonMass
void addModifiedParticle(Particle *p)
void setHelicity(G4double h)