Geant4-11
Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
G4INCL::EtaNToPiPiNChannel Class Reference

#include <G4INCLEtaNToPiPiNChannel.hh>

Inheritance diagram for G4INCL::EtaNToPiPiNChannel:
G4INCL::IChannel

Public Member Functions

 EtaNToPiPiNChannel (Particle *, Particle *)
 
void fillFinalState (FinalState *fs)
 
FinalState * getFinalState ()
 
virtual ~EtaNToPiPiNChannel ()
 

Private Member Functions

 INCL_DECLARE_ALLOCATION_POOL (EtaNToPiPiNChannel)
 

Private Attributes

G4int ind2
 
Particle * particle1
 
Particle * particle2
 

Static Private Attributes

static const G4double angularSlope = 15.
 

Detailed Description

Definition at line 47 of file G4INCLEtaNToPiPiNChannel.hh.

Constructor & Destructor Documentation

◆ EtaNToPiPiNChannel()

G4INCL::EtaNToPiPiNChannel::EtaNToPiPiNChannel ( Particle *  p1,
Particle *  p2 
)

Definition at line 51 of file G4INCLEtaNToPiPiNChannel.cc.

52 : ind2(0),
53 particle1(p1),
54 particle2(p2)
55 {}

◆ ~EtaNToPiPiNChannel()

G4INCL::EtaNToPiPiNChannel::~EtaNToPiPiNChannel ( )
virtual

Definition at line 57 of file G4INCLEtaNToPiPiNChannel.cc.

57 {
58
59 }

Member Function Documentation

◆ fillFinalState()

void G4INCL::EtaNToPiPiNChannel::fillFinalState ( FinalState *  fs)
virtual

Implements G4INCL::IChannel.

Definition at line 61 of file G4INCLEtaNToPiPiNChannel.cc.

61 {
62
63 Particle * nucleon;
64 Particle * eta;
65 if(particle1->isNucleon()) {
66 nucleon = particle1;
67 eta = particle2;
68 } else {
69 nucleon = particle2;
70 eta = particle1;
71 }
72
73 const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, eta);
74
75 ind2=ParticleTable::getIsospin(nucleon->getType());
76
77 eta->setType(PiZero);
78
79 ParticleType pionType;
80 G4double rs1=Random::shoot();
81 if (ind2 == 1) {
82 if (rs1*6. > 5.) pionType=PiZero;
83 else if (rs1*6. > 3.) {
84 pionType=PiPlus;
85 ind2=-ind2;
86 }
87 else {
88 pionType=PiPlus;
89 eta->setType(PiMinus);
90 }
91 }
92 else {
93 if (rs1*6. > 5.) pionType=PiZero;
94 else if (rs1*6. > 3.) {
95 pionType=PiMinus;
96 ind2=-ind2;
97 }
98 else {
99 pionType=PiPlus;
100 eta->setType(PiMinus);
101 }
102 }
103
104
105 const ParticleType tn=ParticleTable::getNucleonType(ind2);
106 nucleon->setType(tn);
107 ParticleList list;
108 list.push_back(nucleon);
109 list.push_back(eta);
110 const ThreeVector &rcolpion = eta->getPosition();
111 const ThreeVector zero;
112 Particle *newPion = new Particle(pionType,zero,rcolpion);
113 list.push_back(newPion);
114 fs->addModifiedParticle(nucleon);
115 fs->addModifiedParticle(eta);
116 fs->addCreatedParticle(newPion);
117
118 PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
119
120 }
G4double
double G4double
Definition: G4Types.hh:83
G4INCL::Particle::isNucleon
G4bool isNucleon() const
Definition: G4INCLParticle.hh:303
G4INCL::KinematicsUtils::totalEnergyInCM
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
Definition: G4INCLKinematicsUtils.cc:94
G4INCL::ParticleTable::getIsospin
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
Definition: G4INCLParticleTable.cc:478
G4INCL::ParticleTable::getNucleonType
ParticleType getNucleonType(const G4int isosp)
Get the type of nucleon.
Definition: G4INCLParticleTable.cc:1253
G4INCL::PhaseSpaceGenerator::generateBiased
void generateBiased(const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
Generate a biased event in the CM system.
Definition: G4INCLPhaseSpaceGenerator.cc:98
G4INCL::Random::shoot
G4double shoot()
Definition: G4INCLRandom.cc:93
G4InuclParticleNames::nucleon
G4bool nucleon(G4int ityp)
Definition: G4InuclParticleNames.hh:94
anonymous_namespace{G4CascadeDeexciteBase.cc}::zero
static const G4LorentzVector zero(0., 0., 0., 0.)

References G4INCL::FinalState::addCreatedParticle(), G4INCL::FinalState::addModifiedParticle(), angularSlope, G4INCL::PhaseSpaceGenerator::generateBiased(), G4INCL::ParticleTable::getIsospin(), G4INCL::ParticleTable::getNucleonType(), G4INCL::Particle::getPosition(), ind2, G4INCL::Particle::isNucleon(), G4InuclParticleNames::nucleon(), particle1, particle2, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Particle::setType(), G4INCL::Random::shoot(), G4INCL::KinematicsUtils::totalEnergyInCM(), and anonymous_namespace{G4CascadeDeexciteBase.cc}::zero.

◆ getFinalState()

FinalState * G4INCL::IChannel::getFinalState ( )
inherited

Definition at line 50 of file G4INCLIChannel.cc.

50 {
51 FinalState *fs = new FinalState;
52 fillFinalState(fs);
53 return fs;
54 }
G4INCL::IChannel::fillFinalState
virtual void fillFinalState(FinalState *fs)=0

References G4INCL::IChannel::fillFinalState().

◆ INCL_DECLARE_ALLOCATION_POOL()

G4INCL::EtaNToPiPiNChannel::INCL_DECLARE_ALLOCATION_POOL ( EtaNToPiPiNChannel  )
private

Field Documentation

◆ angularSlope

const G4double G4INCL::EtaNToPiPiNChannel::angularSlope = 15.
staticprivate

Definition at line 58 of file G4INCLEtaNToPiPiNChannel.hh.

Referenced by fillFinalState().

◆ ind2

G4int G4INCL::EtaNToPiPiNChannel::ind2
private

Definition at line 55 of file G4INCLEtaNToPiPiNChannel.hh.

Referenced by fillFinalState().

◆ particle1

Particle* G4INCL::EtaNToPiPiNChannel::particle1
private

Definition at line 56 of file G4INCLEtaNToPiPiNChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::EtaNToPiPiNChannel::particle2
private

Definition at line 56 of file G4INCLEtaNToPiPiNChannel.hh.

Referenced by fillFinalState().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3