#include <G4INCLPiNToEtaChannel.hh>
Definition at line 47 of file G4INCLPiNToEtaChannel.hh.
G4INCL::PiNToEtaChannel::~PiNToEtaChannel |
( |
| ) |
|
|
virtual |
Implements G4INCL::IChannel.
Definition at line 57 of file G4INCLPiNToEtaChannel.cc.
78 G4double sh=nucleon->getEnergy()+pion->getEnergy();
81 G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
82 nucleon->setEnergy(en);
83 G4double ee=std::sqrt(en*en-mn*mn+me*me);
102 G4double f1= -0.0000288627*ECM*ECM+0.09155289*ECM-72.25436;
103 G4double b1=(f1-(f1/(1.5-0.5*std::pow((ECM-1580.)/95.,2))))/2.;
117 fteta=(a1*x1*x1+b1*x1+c1)/interg1;
119 if (u1*f1/interg1 < fteta) {
134 G4double interg2=0.1716182902205207;
145 fteta=((a2*x1*x1+b2*x1+c2)*(0.5+(std::atan(10*(x1+dev)))/
pi) + vert)/interg2;
156 ThreeVector mom_nucleon(
157 pn*std::sin(teta)*std::cos(fi),
158 pn*std::sin(teta)*std::sin(fi),
163 nucleon->setMomentum(-mom_nucleon);
164 pion->setMomentum(mom_nucleon);
166 fs->addModifiedParticle(nucleon);
167 fs->addModifiedParticle(pion);
G4bool nucleon(G4int ityp)
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
static constexpr double pi
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
The documentation for this class was generated from the following files: