34 #define INCLXX_IN_GEANT4_MODE 1
47 :theParticle(p), incidentDirection(dir)
56 if(m > 1500.0) gg = 200.0;
59 const G4double psf = std::pow(qqq, 3)/(std::pow(qqq, 3) + 5832000.0);
66 unsigned long loopCounter = 0;
67 const unsigned long maxLoopCounter = 10000000;
70 if(std::abs(*ctet_par) > 1.0) (*ctet_par) =
Math::sign(*ctet_par);
72 }
while(loopCounter<maxLoopCounter &&
Random::shoot() > ((1.0 + 3.0 * hel * (*ctet_par) * (*ctet_par)) /(1.0 + 3.0 * hel)));
73 (*stet_par) = std::sqrt(1.-(*ctet_par)*(*ctet_par));
107 sampleAngles(&ctet, &stet, &fi);
116 sal = incidentDirection.
perp()/beta;
122 G4double t1 = ctet+cal*stet*sfi/sal;
124 q1=(b1*t1+b2*t2*cfi)/beta;
125 q2=(b2*t1-b1*t2*cfi)/beta;
126 q3=(b3*t1/beta-t2*sfi);
135 switch(theParticle->
getType()) {
void fillFinalState(FinalState *fs)
void addCreatedParticle(Particle *p)
G4double getMass() const
Get the cached particle mass.
const G4double hc
[MeV*fm]
G4double adjustEnergyFromMomentum()
Recompute the energy to match the momentum.
DeltaDecayChannel(Particle *, ThreeVector const &)
virtual ~DeltaDecayChannel()
G4double getEnergy() const
G4double momentumInCM(Particle const *const p1, Particle const *const p2)
gives the momentum in the CM frame of two particles.
static G4double computeDecayTime(Particle *p)
static constexpr double m
const G4double effectivePionMass
G4double getINCLMass(const G4int A, const G4int Z)
Get INCL nuclear mass (in MeV/c^2)
G4INCL::ParticleType getType() const
const G4INCL::ThreeVector & getPosition() const
void setType(ParticleType t)
const G4double effectiveNucleonMass
void addModifiedParticle(Particle *p)
void setHelicity(G4double h)
virtual void setMomentum(const G4INCL::ThreeVector &momentum)