G4INCL::SrcChannel Class Reference

#include <G4INCLSrcChannel.hh>

Inheritance diagram for G4INCL::SrcChannel:

Public Member Functions
	SrcChannel (Particle *p1, Particle *p2, Nucleus *n)
virtual	~SrcChannel ()
void	fillFinalState (FinalState *fs)
void	fillFinalState (FinalState *fs, ParticleType, ParticleType)
Public Member Functions inherited from G4INCL::IChannel
	IChannel ()
virtual	~IChannel ()
FinalState *	getFinalState ()

Detailed Description

Definition at line 57 of file G4INCLSrcChannel.hh.

Constructor & Destructor Documentation

SrcChannel()

G4INCL::SrcChannel::SrcChannel	(	Particle *	p1,
		Particle *	p2,
		Nucleus *	n )

Definition at line 54 of file G4INCLSrcChannel.cc.

    : particle1(p1), particle2(p2), thenucleus(n) {
  fDistSrc = ParticleTable::getsrcPairDistance();
  srcpartner = nullptr;
  ftype1 = UnknownParticle;
  ftype2 = UnknownParticle;
}

~SrcChannel()

G4INCL::SrcChannel::~SrcChannel ( )

virtual

Definition at line 62 of file G4INCLSrcChannel.cc.

{}

Member Function Documentation

fillFinalState() [1/2]

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

virtual

Implements G4INCL::IChannel.

Definition at line 89 of file G4INCLSrcChannel.cc.

                                              {
 
  auto fSource = InteractionAvatar::Instance();
 
  G4double psrcmax = 2.0 * PhysicalConstants::Pf; // Fermi momentum in MeV/c
 
  if (particle1->getSrcPair() > 0) {
 
    if (particle2->getSrcPair() > 0) 
    {
      INCL_ERROR(particle2->print() << " \n");
    }
    srcpartner = findpairpartner(particle1);
 
    if (srcpartner) 
    {
      srcpartner->setSrcPartner();
      fSource->setSrcPartner(srcpartner);
 
      ThreeVector d1 = particle1->getPosition();
      ThreeVector d2 = srcpartner->getPosition();
      auto d = d1 - d2;
      if (d.mag() > fDistSrc) {
        INCL_DEBUG("Distance src > " << fDistSrc << " fm : " << d.mag()
                                     << " \n");
      }
 
      auto x = (fDistSrc - d.mag()) / fDistSrc;
      auto srcp = x * x * psrcmax;
      INCL_DEBUG("Src momentum = " << srcp << " , eventnb: "
                                   << theEventInfo.eventNumber << " \n");
 
      auto pmomentum =
          particle1->getMomentum() / particle1->getMomentum().mag() * srcp;
      particle1->setMomentum(particle1->getMomentum() + pmomentum);
      srcpartner->setMomentum(srcpartner->getMomentum() - pmomentum);
 
      particle1->adjustEnergyFromMomentum();
      srcpartner->adjustEnergyFromMomentum();
      thenucleus->updatePotentialEnergy(srcpartner);
 
      fs->addModifiedParticle(particle1);
      fs->addModifiedParticle(particle2);
      fs->addModifiedParticle(srcpartner);
    }
  } else {
 
    if (particle1->getSrcPair() > 0) 
    {
      INCL_ERROR(particle1->print() << " \n");
    }
    srcpartner = findpairpartner(particle2);
 
    if (srcpartner) 
    {
      srcpartner->setSrcPartner();
      fSource->setSrcPartner(srcpartner);
 
      ThreeVector d1 = particle2->getPosition();
      ThreeVector d2 = srcpartner->getPosition();
      auto d = d1 - d2;
 
      if (d.mag() > fDistSrc) {
        INCL_DEBUG("Distance src > " << fDistSrc << " fm : " << d.mag()
                                     << " \n");
      }
 
      auto x = (fDistSrc - d.mag()) / fDistSrc;
      auto srcp = x * x * psrcmax;
      INCL_DEBUG("Src momentum = " << srcp << " , eventnb: "
                                   << theEventInfo.eventNumber << " \n");
 
      auto pmomentum =
          particle2->getMomentum() / particle2->getMomentum().mag() * srcp;
 
      particle2->setMomentum(particle2->getMomentum() + pmomentum);
      srcpartner->setMomentum(srcpartner->getMomentum() - pmomentum);
 
      particle2->adjustEnergyFromMomentum();
      srcpartner->adjustEnergyFromMomentum();
      thenucleus->updatePotentialEnergy(srcpartner);
 
      fs->addModifiedParticle(particle1);
      fs->addModifiedParticle(particle2);
      fs->addModifiedParticle(srcpartner);
    }
  }
  thenucleus->getStore()->getBook().incrementAcceptedSrcCollisions();
}

Referenced by fillFinalState().

fillFinalState() [2/2]

void G4INCL::SrcChannel::fillFinalState	(	FinalState *	fs,
		ParticleType	type1,
		ParticleType	type2 )

Definition at line 82 of file G4INCLSrcChannel.cc.

                                                    {
  ftype1 = type1;
  ftype2 = type2;
  fillFinalState(fs);
}

---

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

• G4INCLSrcChannel.hh
• G4INCLSrcChannel.cc