Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4INCL::OmegaNElasticChannel Class Reference

#include <G4INCLOmegaNElasticChannel.hh>

Inheritance diagram for G4INCL::OmegaNElasticChannel:
Collaboration 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 ()
 

Detailed Description

Definition at line 47 of file G4INCLOmegaNElasticChannel.hh.

Constructor & Destructor Documentation

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

Definition at line 47 of file G4INCLOmegaNElasticChannel.cc.

48  : particle1(p1), particle2(p2)
49  {
50 
51  }
G4INCL::OmegaNElasticChannel::~OmegaNElasticChannel ( )
virtual

Definition at line 53 of file G4INCLOmegaNElasticChannel.cc.

53  {
54 
55  }

Member Function Documentation

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  }
G4double getEnergy() const
ThreeVector normVector(G4double norm=1.)
G4bool nucleon(G4int ityp)
G4bool isNucleon() const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:


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