34 #define INCLXX_IN_GEANT4_MODE 1
52 :theParticle(p), incidentDirection(dir)
82 if(std::abs(*ctet_par) > 1.0) (*ctet_par) =
Math::sign(*ctet_par);
83 (*stet_par) = std::sqrt(1.-(*ctet_par)*(*ctet_par));
104 else if (drnd < 0.7265) {
110 else if (drnd < 0.9575) {
130 else if (drnd < 0.9845) {
167 q1=(b1*t1+b2*t2*cfi)/beta;
168 q2=(b2*t1-b1*t2*cfi)/beta;
169 q3=(b3*t1/beta-t2*sfi);
185 Particle *createdParticle =
new Particle(createdType, createdMomentum, createdPosition);
193 else if (nbpart == 3) {
201 list.push_back(Pion1);
202 list.push_back(Pion2);
void addModifiedParticle(Particle *p)
const G4INCL::ThreeVector & getPosition() const
G4double getEnergy() const
static G4double computeDecayTime(Particle *p)
G4double getMass() const
Get the cached particle mass.
PionResonanceDecayChannel(Particle *, ThreeVector const &)
void generate(const G4double sqrtS, ParticleList &particles)
Generate an event in the CM system.
void fillFinalState(FinalState *fs)
static const G4double angularSlope
G4double momentumInCM(Particle const *const p1, Particle const *const p2)
gives the momentum in the CM frame of two particles.
static constexpr double m
void sampleAngles(G4double *, G4double *, G4double *)
G4double getINCLMass(const G4int A, const G4int Z)
Get INCL nuclear mass (in MeV/c^2)
void setType(ParticleType t)
virtual void setMomentum(const G4INCL::ThreeVector &momentum)
ThreeVector const incidentDirection
G4INCL::ParticleType getType() const
G4double adjustEnergyFromMomentum()
Recompute the energy to match the momentum.
void addCreatedParticle(Particle *p)
virtual ~PionResonanceDecayChannel()
const G4double hc
[MeV*fm]