Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4INCL::OmegaNElasticChannel Class Reference

#include <G4INCLOmegaNElasticChannel.hh>

+ Inheritance diagram for G4INCL::OmegaNElasticChannel:

Public Member Functions

 OmegaNElasticChannel (Particle *, Particle *)
 
virtual ~OmegaNElasticChannel ()
 
void fillFinalState (FinalState *fs)
 
- Public Member Functions inherited from G4INCL::IChannel
 IChannel ()
 
virtual ~IChannel ()
 
FinalStategetFinalState ()
 
virtual void fillFinalState (FinalState *fs)=0
 

Detailed Description

Definition at line 47 of file G4INCLOmegaNElasticChannel.hh.

Constructor & Destructor Documentation

◆ OmegaNElasticChannel()

G4INCL::OmegaNElasticChannel::OmegaNElasticChannel ( Particle p1,
Particle p2 
)

Definition at line 47 of file G4INCLOmegaNElasticChannel.cc.

48 : particle1(p1), particle2(p2)
49 {
50
51 }

◆ ~OmegaNElasticChannel()

G4INCL::OmegaNElasticChannel::~OmegaNElasticChannel ( )
virtual

Definition at line 53 of file G4INCLOmegaNElasticChannel.cc.

53 {
54
55 }

Member Function Documentation

◆ fillFinalState()

void G4INCL::OmegaNElasticChannel::fillFinalState ( FinalState fs)
virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLOmegaNElasticChannel.cc.

57 {
58 Particle * nucleon;
59 Particle * omega;
60 if(particle1->isNucleon()) {
61 nucleon = particle1;
62 omega = particle2;
63 } else {
64 nucleon = particle2;
65 omega = particle1;
66 }
67
68 G4double sh=nucleon->getEnergy()+omega->getEnergy();
69 G4double mn=nucleon->getMass();
70 G4double me=omega->getMass();
71 G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
72 nucleon->setEnergy(en);
73 G4double ee=std::sqrt(en*en-mn*mn+me*me);
74 omega->setEnergy(ee);
75 G4double pn=std::sqrt(en*en-mn*mn);
76
77 ThreeVector mom_nucleon;
78
79// Isotropy
80 mom_nucleon = Random::normVector(pn);
81
82
83 nucleon->setMomentum(mom_nucleon);
84 omega->setMomentum(-mom_nucleon);
85
86 fs->addModifiedParticle(nucleon);
87 fs->addModifiedParticle(omega);
88
89 }
double G4double
Definition: G4Types.hh:83
G4bool isNucleon() const
ThreeVector normVector(G4double norm=1.)
G4bool nucleon(G4int ityp)

The documentation for this class was generated from the following files: