34 #define INCLXX_IN_GEANT4_MODE 1
48 : particle1(p1), particle2(p2)
93 G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
95 G4double ee=std::sqrt(en*en-mn*mn+me*me);
121 G4double f1= -0.0000288627*ECM*ECM+0.09155289*ECM-72.25436;
122 G4double b1=(f1-(f1/(1.5-0.5*std::pow((ECM-1580.)/95.,2))))/2.;
136 fteta=(a1*x1*x1+b1*x1+c1)/interg1;
138 if (u1*f1/interg1 < fteta) {
153 G4double interg2=0.1716182902205207;
164 fteta=((a2*x1*x1+b2*x1+c2)*(0.5+(std::atan(10*(x1+dev)))/
pi) + vert)/interg2;
176 pn*std::sin(teta)*std::cos(fi),
177 pn*std::sin(teta)*std::sin(fi),
G4double getMass() const
Get the cached particle mass.
G4double getEnergy() const
virtual ~EtaNToPiNChannel()
void setEnergy(G4double energy)
void fillFinalState(FinalState *fs)
G4bool nucleon(G4int ityp)
G4INCL::ParticleType getType() const
void setType(ParticleType t)
static constexpr double pi
EtaNToPiNChannel(Particle *, Particle *)
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
void addModifiedParticle(Particle *p)
virtual void setMomentum(const G4INCL::ThreeVector &momentum)