G4INCL::NDeltaOmegaProductionChannel Class Reference

#include <G4INCLNDeltaOmegaProductionChannel.hh>

Inheritance diagram for G4INCL::NDeltaOmegaProductionChannel:

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

Detailed Description

Definition at line 48 of file G4INCLNDeltaOmegaProductionChannel.hh.

Constructor & Destructor Documentation

NDeltaOmegaProductionChannel()

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

Definition at line 51 of file G4INCLNDeltaOmegaProductionChannel.cc.

    : particle1(p1), particle2(p2)
  {}

~NDeltaOmegaProductionChannel()

G4INCL::NDeltaOmegaProductionChannel::~NDeltaOmegaProductionChannel ( )

virtual

Definition at line 55 of file G4INCLNDeltaOmegaProductionChannel.cc.

{}

Member Function Documentation

fillFinalState()

void G4INCL::NDeltaOmegaProductionChannel::fillFinalState ( FinalState *  fs )

virtual

Unlike NN -> NDelta, NN -> NDeltaOmega is drawn from a phase-space generator

Implements G4INCL::IChannel.

Definition at line 97 of file G4INCLNDeltaOmegaProductionChannel.cc.

                                                                 {
  
/**
*
* Unlike NN -> NDelta, NN -> NDeltaOmega is drawn from a phase-space generator
*
**/
 
  G4int is1=ParticleTable::getIsospin(particle1->getType());
  G4int is2=ParticleTable::getIsospin(particle2->getType());
  
  ParticleList list;
  list.push_back(particle1);
  list.push_back(particle2);
  
//  isospin Repartition of N and Delta;
  G4double ecm = KinematicsUtils::totalEnergyInCM(particle1, particle2);  
  const G4int isospin = is1+is2;
  
  G4double rndm = 0.0;
  G4double xmdel = sampleDeltaMass(ecm);
 
  G4int index2=0;
  if (isospin == 0) { // pn case
   rndm = Random::shoot();
   if (rndm < 0.5) index2=1;
  }
 
  if (isospin == 0) {
   if(index2 == 1) {
    G4int isi=is1;
    is1=is2;
    is2=isi;
   }
//   particle1->setHelicity(0.0);
  } else {
   rndm = Random::shoot();
   if (rndm >= 0.25) {
    is1=3*is1;
    is2=-is2;
   }
//   particle1->setHelicity(ctet*ctet);
  }
  
  if(is1 == ParticleTable::getIsospin(DeltaMinus)) {
   particle1->setType(DeltaMinus);
  } else if(is1 == ParticleTable::getIsospin(DeltaZero)) {
   particle1->setType(DeltaZero);
  } else if(is1 == ParticleTable::getIsospin(DeltaPlus)) {
   particle1->setType(DeltaPlus);
  } else if(is1 == ParticleTable::getIsospin(DeltaPlusPlus)) {
   particle1->setType(DeltaPlusPlus);
  }
  
  if(is2 == ParticleTable::getIsospin(Proton)) {
   particle2->setType(Proton);
  } else if(is2 == ParticleTable::getIsospin(Neutron)) {
   particle2->setType(Neutron);
  }
  
  if(particle1->isDelta()) particle1->setMass(xmdel);
  if(particle2->isDelta()) particle2->setMass(xmdel);
  
  const ThreeVector &rcolnucleon1 = particle1->getPosition();
  const ThreeVector &rcolnucleon2 = particle2->getPosition();
  const ThreeVector rcol = (rcolnucleon1+rcolnucleon2)*0.5;
  const ThreeVector zero;
  Particle *omega = new Particle(Omega,zero,rcol);
  list.push_back(omega);
  fs->addCreatedParticle(omega);
  
  const G4double sqrtS = KinematicsUtils::totalEnergyInCM(particle1, particle2);
  G4int biasIndex = ((Random::shoot()<0.5) ? 0 : 1);
  PhaseSpaceGenerator::generateBiased(sqrtS, list, biasIndex, angularSlope);
  
  const ThreeVector vz(0.0,0.0,1.0);
  G4double ctet=(particle1->getMomentum().dot(vz))/particle1->getMomentum().mag();
  if (isospin == 0)
   particle1->setHelicity(0.0);
  else
   particle1->setHelicity(ctet*ctet);
  
  fs->addModifiedParticle(particle1);
  fs->addModifiedParticle(particle2);
  
 }

---

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

• G4INCLNDeltaOmegaProductionChannel.hh
• G4INCLNDeltaOmegaProductionChannel.cc