#include <G4INCLInteractionAvatar.hh>

Inheritance diagram for G4INCL::InteractionAvatar:

Data Structures
class	ViolationEEnergyFunctor
	RootFunctor-derived object for enforcing energy conservation in delta production. More...

class	ViolationEMomentumFunctor
	RootFunctor-derived object for enforcing energy conservation in N-N. More...

Public Member Functions
virtual std::string	dump () const =0

void	fillFinalState (FinalState *fs)

FinalState *	getFinalState ()

long	getID () const

virtual ParticleList	getParticles () const =0

G4double	getTime () const

AvatarType	getType () const

	InteractionAvatar (G4double, G4INCL::Nucleus , G4INCL::Particle )

	InteractionAvatar (G4double, G4INCL::Nucleus , G4INCL::Particle , G4INCL::Particle *)

G4bool	isACollision () const

G4bool	isADecay () const

void	setType (AvatarType t)

std::string	toString ()

virtual	~InteractionAvatar ()

Static Public Member Functions
static void	deleteBackupParticles ()
	Release the memory allocated for the backup particles. More...

Static Public Attributes
static const G4double	locEAccuracy = 1.E-4
	Target accuracy in the determination of the local-energy Q-value. More...

static const G4int	maxIterLocE = 50
	Max number of iterations for the determination of the local-energy Q-value. More...

Protected Member Functions
G4bool	bringParticleInside (Particle *const p)

G4bool	enforceEnergyConservation (FinalState *const fs)
	Enforce energy conservation. More...

virtual G4INCL::IChannel *	getChannel ()=0

void	postInteraction (FinalState *)

void	preInteraction ()

void	preInteractionBlocking ()
	Store the state of the particles before the interaction. More...

void	preInteractionLocalEnergy (Particle *const p)
	Apply local-energy transformation, if appropriate. More...

void	restoreParticles () const
	Restore the state of both particles. More...

G4bool	shouldUseLocalEnergy () const
	true if the given avatar should use local energy More...

Protected Attributes
ThreeVector	boostVector

ParticleList	created

ParticleList	Destroyed

G4bool	isPiN

ParticleList	modified

ParticleList	modifiedAndCreated

ParticleList	ModifiedAndDestroyed

G4double	oldTotalEnergy

G4double	oldXSec

Particle *	particle1

Particle *	particle2

Nucleus *	theNucleus

G4double	theTime

G4double	weight

Static Protected Attributes
static G4ThreadLocal Particle *	backupParticle1 = NULL

static G4ThreadLocal Particle *	backupParticle2 = NULL

Private Attributes
long	ID

AvatarType	type

RootFunctor *	violationEFunctor

Static Private Attributes
static G4ThreadLocal long	nextID = 1

Detailed Description

Definition at line 61 of file G4INCLInteractionAvatar.hh.

Constructor & Destructor Documentation

◆ InteractionAvatar() [1/2]

G4INCL::InteractionAvatar::InteractionAvatar	(	G4double	time,
		G4INCL::Nucleus *	n,
		G4INCL::Particle *	p1
	)

Definition at line 65 of file G4INCLInteractionAvatar.cc.

    : IAvatar(time), theNucleus(n),
    particle1(p1), particle2(NULL),
    isPiN(false),
    weight(1.),
    violationEFunctor(NULL)
  {
  }

◆ InteractionAvatar() [2/2]

G4INCL::InteractionAvatar::InteractionAvatar	(	G4double	time,
		G4INCL::Nucleus *	n,
		G4INCL::Particle *	p1,
		G4INCL::Particle *	p2
	)

Definition at line 74 of file G4INCLInteractionAvatar.cc.

    : IAvatar(time), theNucleus(n),
    particle1(p1), particle2(p2),
    isPiN((p1->isPion() && p2->isNucleon()) || (p2->isPion() && p1->isNucleon())),
    weight(1.),
    violationEFunctor(NULL)
  {
  }

◆ ~InteractionAvatar()

G4INCL::InteractionAvatar::~InteractionAvatar ( )

virtual

Definition at line 84 of file G4INCLInteractionAvatar.cc.

84 {

85 }

Member Function Documentation

◆ bringParticleInside()

G4bool G4INCL::InteractionAvatar::bringParticleInside ( Particle *const p )

protected

Definition at line 137 of file G4INCLInteractionAvatar.cc.

                                                                  {
    if(!theNucleus)
      return false;
 
    ThreeVector pos = p->getPosition();
    p->rpCorrelate();
    G4double pos2 = pos.mag2();
    const G4double r = theNucleus->getSurfaceRadius(p);
    short iterations=0;
    const short maxIterations=50;
 
    if(pos2 < r*r) return true;
 
    while( pos2 >= r*r && iterations<maxIterations ) /* Loop checking, 10.07.2015, D.Mancusi */
    {
      pos *= std::sqrt(r*r*0.9801/pos2); // 0.9801 == 0.99*0.99
      pos2 = pos.mag2();
      iterations++;
    }
    if( iterations < maxIterations)
    {
      INCL_DEBUG("Particle position vector length was : " << p->getPosition().mag() << ", rescaled to: " << pos.mag() << '\n');
      p->setPosition(pos);
      return true;
    }
    else
      return false;
  }

References G4INCL::Particle::getPosition(), G4INCL::Nucleus::getSurfaceRadius(), INCL_DEBUG, G4INCL::ThreeVector::mag(), G4INCL::RootFinder::anonymous_namespace{G4INCLRootFinder.cc}::maxIterations, pos, G4INCL::Particle::rpCorrelate(), G4INCL::Particle::setPosition(), and theNucleus.

Referenced by postInteraction().

◆ deleteBackupParticles()

void G4INCL::InteractionAvatar::deleteBackupParticles ( )

static

Release the memory allocated for the backup particles.

Definition at line 87 of file G4INCLInteractionAvatar.cc.

                                                {
    delete backupParticle1;
    if(backupParticle2)
      delete backupParticle2;
    backupParticle1 = NULL;
    backupParticle2 = NULL;
  }

References backupParticle1, and backupParticle2.

Referenced by G4INCL::INCL::~INCL().

◆ dump()

virtual std::string G4INCL::IAvatar::dump ( ) const

pure virtualinherited

Implemented in G4INCL::BinaryCollisionAvatar, G4INCL::DecayAvatar, G4INCL::ParticleEntryAvatar, and G4INCL::SurfaceAvatar.

Referenced by G4INCL::CascadeAction::beforeAvatarDefaultAction().

◆ enforceEnergyConservation()

G4bool G4INCL::InteractionAvatar::enforceEnergyConservation ( FinalState *const fs )

protected

Enforce energy conservation.

Final states generated by the channels might violate energy conservation because of different reasons (energy-dependent potentials, local energy...). This conservation law must therefore be enforced by hand. We do so by rescaling the momenta of the final-state particles in the CM frame. If this turns out to be impossible, this method returns false.

Returns: true if the algorithm succeeded

Definition at line 359 of file G4INCLInteractionAvatar.cc.

                                                                           {
    // Set up the violationE calculation
    const G4bool manyBodyFinalState = (modifiedAndCreated.size() > 1);
    
    if(manyBodyFinalState)
      violationEFunctor = new ViolationEMomentumFunctor(theNucleus, modifiedAndCreated, fs->getTotalEnergyBeforeInteraction(), boostVector, shouldUseLocalEnergy());
    else {
      Particle * const p = modified.front();
      // The following condition is necessary for the functor to work
      // correctly. A similar condition exists in INCL4.6.
      if(p->getMass() < ParticleTable::minDeltaMass)
        return false;
      violationEFunctor = new ViolationEEnergyFunctor(theNucleus, p, fs->getTotalEnergyBeforeInteraction(), shouldUseLocalEnergy());
    }
 
    // Apply the root-finding algorithm
    const RootFinder::Solution theSolution = RootFinder::solve(violationEFunctor, 1.0);
    if(theSolution.success) { // Apply the solution
      (*violationEFunctor)(theSolution.x);
    } else if(theNucleus){
      INCL_DEBUG("Couldn't enforce energy conservation after an interaction, root-finding algorithm failed." << '\n');
      theNucleus->getStore()->getBook().incrementEnergyViolationInteraction();
    }
    delete violationEFunctor;
    violationEFunctor = NULL;
    return theSolution.success;
  }

References boostVector, G4INCL::Store::getBook(), G4INCL::Particle::getMass(), G4INCL::Nucleus::getStore(), G4INCL::FinalState::getTotalEnergyBeforeInteraction(), INCL_DEBUG, G4INCL::Book::incrementEnergyViolationInteraction(), G4INCL::ParticleTable::minDeltaMass, modified, modifiedAndCreated, shouldUseLocalEnergy(), G4INCL::RootFinder::solve(), G4INCL::RootFinder::Solution::success, theNucleus, violationEFunctor, and G4INCL::RootFinder::Solution::x.

Referenced by postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ fillFinalState()

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

inherited

Definition at line 96 of file G4INCLIAvatar.cc.

                                             {
    INCL_DEBUG("Random seeds before preInteraction: " << Random::getSeeds() << '\n');
    preInteraction();
    INCL_DEBUG("Random seeds before getChannel: " << Random::getSeeds() << '\n');
    IChannel *c = getChannel();
    if( !c )
      return;
    INCL_DEBUG("Random seeds before getFinalState: " << Random::getSeeds() << '\n');
    c->fillFinalState(fs);
    INCL_DEBUG("Random seeds before postInteraction: " << Random::getSeeds() << '\n');
    postInteraction(fs);
    delete c;
  }

References G4INCL::IChannel::fillFinalState(), G4INCL::IAvatar::getChannel(), G4INCL::Random::getSeeds(), INCL_DEBUG, G4INCL::IAvatar::postInteraction(), and G4INCL::IAvatar::preInteraction().

Referenced by G4INCL::INCL::cascade(), and G4INCL::IAvatar::getFinalState().

◆ getChannel()

virtual G4INCL::IChannel * G4INCL::InteractionAvatar::getChannel ( )

protectedpure virtual

Implements G4INCL::IAvatar.

Implemented in G4INCL::BinaryCollisionAvatar, and G4INCL::DecayAvatar.

◆ getFinalState()

FinalState * G4INCL::IAvatar::getFinalState ( )

inherited

Definition at line 90 of file G4INCLIAvatar.cc.

                                     {
    FinalState *fs = new FinalState;
    fillFinalState(fs);
    return fs;
  }

References G4INCL::IAvatar::fillFinalState().

Referenced by G4INCL::INCL::makeCompoundNucleus().

◆ getID()

long G4INCL::IAvatar::getID ( ) const

inlineinherited

Definition at line 94 of file G4INCLIAvatar.hh.

94{ return ID; };

G4INCL::IAvatar::ID

long ID

Definition: G4INCLIAvatar.hh:98

Referenced by G4INCL::CascadeAction::afterAvatarDefaultAction(), G4INCL::CascadeAction::beforeAvatarDefaultAction(), and G4INCL::IAvatar::toString().

◆ getParticles()

virtual ParticleList G4INCL::IAvatar::getParticles ( ) const

pure virtualinherited

Implemented in G4INCL::BinaryCollisionAvatar, G4INCL::DecayAvatar, G4INCL::ParticleEntryAvatar, and G4INCL::SurfaceAvatar.

Referenced by G4INCL::Store::add(), G4INCL::Store::addParticleEntryAvatar(), G4INCL::Store::removeAvatar(), and G4INCL::IAvatar::toString().

◆ getTime()

G4double G4INCL::IAvatar::getTime ( ) const

inlineinherited

Definition at line 84 of file G4INCLIAvatar.hh.

84{ return theTime; };

G4INCL::IAvatar::theTime

G4double theTime

Definition: G4INCLIAvatar.hh:102

Referenced by G4INCL::AvatarDumpAction::afterAvatarUserAction(), G4INCL::Store::avatarComparisonPredicate(), G4INCL::StandardPropagationModel::propagate(), and G4INCL::IAvatar::toString().

◆ getType()

AvatarType G4INCL::IAvatar::getType ( ) const

inlineinherited

Definition at line 90 of file G4INCLIAvatar.hh.

90{ return type; };

G4INCL::IAvatar::type

AvatarType type

Definition: G4INCLIAvatar.hh:99

Referenced by G4INCL::AvatarDumpAction::afterAvatarUserAction(), G4INCL::CascadeAction::beforeAvatarDefaultAction(), shouldUseLocalEnergy(), and G4INCL::IAvatar::toString().

◆ isACollision()

G4bool G4INCL::IAvatar::isACollision ( ) const

inlineinherited

Definition at line 91 of file G4INCLIAvatar.hh.

91{ return (type==CollisionAvatarType); };

G4INCL::CollisionAvatarType

@ CollisionAvatarType

Definition: G4INCLIAvatar.hh:67

References G4INCL::CollisionAvatarType.

Referenced by G4INCL::CascadeAction::afterAvatarDefaultAction().

◆ isADecay()

G4bool G4INCL::IAvatar::isADecay ( ) const

inlineinherited

Definition at line 92 of file G4INCLIAvatar.hh.

92{ return (type==DecayAvatarType); };

G4INCL::DecayAvatarType

@ DecayAvatarType

Definition: G4INCLIAvatar.hh:68

References G4INCL::DecayAvatarType.

Referenced by G4INCL::CascadeAction::afterAvatarDefaultAction().

◆ postInteraction()

void G4INCL::InteractionAvatar::postInteraction ( FinalState * fs )

protectedvirtual

Implements G4INCL::IAvatar.

Definition at line 166 of file G4INCLInteractionAvatar.cc.

                                                        {
    INCL_DEBUG("postInteraction: final state: " << '\n' << fs->print() << '\n');
    modified = fs->getModifiedParticles();
    created = fs->getCreatedParticles();
    Destroyed = fs->getDestroyedParticles();
    modifiedAndCreated = modified;
    modifiedAndCreated.insert(modifiedAndCreated.end(), created.begin(), created.end());
    ModifiedAndDestroyed = modified;
    ModifiedAndDestroyed.insert(ModifiedAndDestroyed.end(), Destroyed.begin(), Destroyed.end());
 
    // Boost back to lab
    modifiedAndCreated.boost(-boostVector);
 
    // If there is no Nucleus, just return
    if(!theNucleus) return;
 
    // Mark pions and kaons that have been created outside their well (we will force them
    // to be emitted later).
    for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
      if(((*i)->isPion() || (*i)->isKaon() || (*i)->isAntiKaon()) && (*i)->getPosition().mag() > theNucleus->getSurfaceRadius(*i)) {
        (*i)->makeParticipant();
        (*i)->setOutOfWell();
        fs->addOutgoingParticle(*i);
        INCL_DEBUG("Pion was created outside its potential well." << '\n'
            << (*i)->print());
      }
 
    // Try to enforce energy conservation
    fs->setTotalEnergyBeforeInteraction(oldTotalEnergy);
    G4bool success = enforceEnergyConservation(fs);
    if(!success) {
      INCL_DEBUG("Enforcing energy conservation: failed!" << '\n');
 
      // Restore the state of the initial particles
      restoreParticles();
 
      // Delete newly created particles
      for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
        delete *i;
 
      fs->reset();
      fs->makeNoEnergyConservation();
      fs->setTotalEnergyBeforeInteraction(0.0);
 
      return; // Interaction is blocked. Return an empty final state.
    }
    INCL_DEBUG("Enforcing energy conservation: success!" << '\n');
 
    INCL_DEBUG("postInteraction after energy conservation: final state: " << '\n' << fs->print() << '\n');
 
    // Check that outgoing delta resonances can decay to pi-N
    for(ParticleIter i=modified.begin(), e=modified.end(); i!=e; ++i )
      if((*i)->isDelta() &&
          (*i)->getMass() < ParticleTable::minDeltaMass) {
        INCL_DEBUG("Mass of the produced delta below decay threshold; forbidding collision. deltaMass=" <<
            (*i)->getMass() << '\n');
 
        // Restore the state of the initial particles
        restoreParticles();
 
        // Delete newly created particles
        for(ParticleIter j=created.begin(), end=created.end(); j!=end; ++j )
          delete *j;
 
        fs->reset();
        fs->makeNoEnergyConservation();
        fs->setTotalEnergyBeforeInteraction(0.0);
 
        return; // Interaction is blocked. Return an empty final state.
      }
 
    INCL_DEBUG("Random seeds before Pauli blocking: " << Random::getSeeds() << '\n');
    // Test Pauli blocking
    G4bool isBlocked = Pauli::isBlocked(modifiedAndCreated, theNucleus);
 
    if(isBlocked) {
      INCL_DEBUG("Pauli: Blocked!" << '\n');
 
      // Restore the state of the initial particles
      restoreParticles();
 
      // Delete newly created particles
      for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
        delete *i;
 
      fs->reset();
      fs->makePauliBlocked();
      fs->setTotalEnergyBeforeInteraction(0.0);
 
      return; // Interaction is blocked. Return an empty final state.
    }
    INCL_DEBUG("Pauli: Allowed!" << '\n');
 
    // Test CDPP blocking
    G4bool isCDPPBlocked = Pauli::isCDPPBlocked(created, theNucleus);
 
    if(isCDPPBlocked) {
      INCL_DEBUG("CDPP: Blocked!" << '\n');
 
      // Restore the state of the initial particles
      restoreParticles();
 
      // Delete newly created particles
      for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
        delete *i;
 
      fs->reset();
      fs->makePauliBlocked();
      fs->setTotalEnergyBeforeInteraction(0.0);
 
      return; // Interaction is blocked. Return an empty final state.
    }
    INCL_DEBUG("CDPP: Allowed!" << '\n');
 
    // If all went well, try to bring particles inside the nucleus...
    for(ParticleIter i=modifiedAndCreated.begin(), e=modifiedAndCreated.end(); i!=e; ++i )
    {
      // ...except for pions beyond their surface radius.
      if((*i)->isOutOfWell()) continue;
 
      const G4bool successBringParticlesInside = bringParticleInside(*i);
      if( !successBringParticlesInside ) {
        INCL_ERROR("Failed to bring particle inside the nucleus!" << '\n');
      }
    }
 
    // Collision accepted!
    // Biasing of the final state
    std::vector<G4int> newBiasCollisionVector;
    newBiasCollisionVector = ModifiedAndDestroyed.getParticleListBiasVector();
    if(std::fabs(weight-1.) > 1E-6){
      newBiasCollisionVector.push_back(Particle::nextBiasedCollisionID);
      Particle::FillINCLBiasVector(1./weight);
      weight = 1.; // useless?
    }
    for(ParticleIter i=modifiedAndCreated.begin(), e=modifiedAndCreated.end(); i!=e; ++i ) {
      (*i)->setBiasCollisionVector(newBiasCollisionVector);
      if(!(*i)->isOutOfWell()) {
        // Decide if the particle should be made into a spectator
        // (Back to spectator)
        G4bool goesBackToSpectator = false;
        if((*i)->isNucleon() && theNucleus->getStore()->getConfig()->getBackToSpectator()) {
          G4double threshold = (*i)->getPotentialEnergy();
          if((*i)->getType()==Proton)
            threshold += Math::twoThirds*theNucleus->getTransmissionBarrier(*i);
          if((*i)->getKineticEnergy() < threshold)
            goesBackToSpectator = true;
        }
        // Thaw the particle propagation
        (*i)->thawPropagation();
 
        // Increment or decrement the participant counters
        if(goesBackToSpectator) {
          INCL_DEBUG("The following particle goes back to spectator:" << '\n'
              << (*i)->print() << '\n');
          if(!(*i)->isTargetSpectator()) {
            theNucleus->getStore()->getBook().decrementCascading();
          }
          (*i)->makeTargetSpectator();
        } else {
          if((*i)->isTargetSpectator()) {
            theNucleus->getStore()->getBook().incrementCascading();
          }
          (*i)->makeParticipant();
        }
      }
    }
    ParticleList destroyed = fs->getDestroyedParticles();
    for(ParticleIter i=destroyed.begin(), e=destroyed.end(); i!=e; ++i )
      if(!(*i)->isTargetSpectator())
        theNucleus->getStore()->getBook().decrementCascading();
    return;
  }

Referenced by G4INCL::BinaryCollisionAvatar::postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ preInteraction()

void G4INCL::InteractionAvatar::preInteraction ( )

protectedvirtual

Implements G4INCL::IAvatar.

Definition at line 122 of file G4INCLInteractionAvatar.cc.

                                         {
    preInteractionBlocking();
 
    preInteractionLocalEnergy(particle1);
 
    if(particle2) {
      preInteractionLocalEnergy(particle2);
      boostVector = KinematicsUtils::makeBoostVector(particle1, particle2);
      particle2->boost(boostVector);
    } else {
      boostVector = particle1->getMomentum()/particle1->getEnergy();
    }
    particle1->boost(boostVector);
  }

References G4INCL::Particle::boost(), boostVector, G4INCL::Particle::getEnergy(), G4INCL::Particle::getMomentum(), G4INCL::KinematicsUtils::makeBoostVector(), particle1, particle2, preInteractionBlocking(), and preInteractionLocalEnergy().

Referenced by G4INCL::BinaryCollisionAvatar::preInteraction(), and G4INCL::DecayAvatar::preInteraction().

◆ preInteractionBlocking()

void G4INCL::InteractionAvatar::preInteractionBlocking ( )

protected

Store the state of the particles before the interaction.

If the interaction cannot be realised for any reason, we will need to restore the particle state as it was before. This is done by calling the restoreParticles() method.

Definition at line 95 of file G4INCLInteractionAvatar.cc.

                                                 {
    if(backupParticle1)
      (*backupParticle1) = (*particle1);
    else
      backupParticle1 = new Particle(*particle1);
 
    if(particle2) {
      if(backupParticle2)
        (*backupParticle2) = (*particle2);
      else
        backupParticle2 = new Particle(*particle2);
 
      oldTotalEnergy = particle1->getEnergy() + particle2->getEnergy()
        - particle1->getPotentialEnergy() - particle2->getPotentialEnergy();
      oldXSec = CrossSections::total(particle1, particle2);
    } else {
      oldTotalEnergy = particle1->getEnergy() - particle1->getPotentialEnergy();
    }
  }

References backupParticle1, backupParticle2, G4INCL::Particle::getEnergy(), G4INCL::Particle::getPotentialEnergy(), oldTotalEnergy, oldXSec, particle1, particle2, and G4INCL::CrossSections::total().

Referenced by preInteraction().

◆ preInteractionLocalEnergy()

void G4INCL::InteractionAvatar::preInteractionLocalEnergy ( Particle *const p )

protected

Apply local-energy transformation, if appropriate.

Parameters

p	particle to apply the transformation to

Definition at line 115 of file G4INCLInteractionAvatar.cc.

                                                                      {
    if(!theNucleus || p->isMeson()) return; // Local energy does not make any sense without a nucleus
 
    if(shouldUseLocalEnergy())
      KinematicsUtils::transformToLocalEnergyFrame(theNucleus, p);
  }

References G4INCL::Particle::isMeson(), shouldUseLocalEnergy(), theNucleus, and G4INCL::KinematicsUtils::transformToLocalEnergyFrame().

Referenced by preInteraction().

◆ restoreParticles()

void G4INCL::InteractionAvatar::restoreParticles ( ) const

protected

Restore the state of both particles.

The state must first be stored by calling preInteractionBlocking().

Definition at line 340 of file G4INCLInteractionAvatar.cc.

                                                 {
    (*particle1) = (*backupParticle1);
    if(particle2)
      (*particle2) = (*backupParticle2);
  }

References particle2.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ setType()

void G4INCL::IAvatar::setType ( AvatarType t )

inlineinherited

Definition at line 93 of file G4INCLIAvatar.hh.

93{ type = t; };

Referenced by G4INCL::BinaryCollisionAvatar::BinaryCollisionAvatar(), G4INCL::DecayAvatar::DecayAvatar(), G4INCL::ParticleEntryAvatar::ParticleEntryAvatar(), and G4INCL::SurfaceAvatar::SurfaceAvatar().

◆ shouldUseLocalEnergy()

G4bool G4INCL::InteractionAvatar::shouldUseLocalEnergy ( ) const

protected

true if the given avatar should use local energy

Definition at line 346 of file G4INCLInteractionAvatar.cc.

                                                       {
    if(!theNucleus) return false;
    LocalEnergyType theLocalEnergyType;
    if(getType()==DecayAvatarType || isPiN)
      theLocalEnergyType = theNucleus->getStore()->getConfig()->getLocalEnergyPiType();
    else
      theLocalEnergyType = theNucleus->getStore()->getConfig()->getLocalEnergyBBType();
 
    const G4bool firstAvatar = (theNucleus->getStore()->getBook().getAcceptedCollisions() == 0);
    return ((theLocalEnergyType == FirstCollisionLocalEnergy && firstAvatar) ||
            theLocalEnergyType == AlwaysLocalEnergy);
  }

References G4INCL::AlwaysLocalEnergy, G4INCL::DecayAvatarType, G4INCL::FirstCollisionLocalEnergy, G4INCL::Book::getAcceptedCollisions(), G4INCL::Store::getBook(), G4INCL::Store::getConfig(), G4INCL::Config::getLocalEnergyBBType(), G4INCL::Config::getLocalEnergyPiType(), G4INCL::Nucleus::getStore(), G4INCL::IAvatar::getType(), isPiN, and theNucleus.

Referenced by enforceEnergyConservation(), preInteractionLocalEnergy(), G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), and G4INCL::InteractionAvatar::ViolationEEnergyFunctor::setParticleEnergy().

◆ toString()

std::string G4INCL::IAvatar::toString ( )

inherited

Definition at line 73 of file G4INCLIAvatar.cc.

                              {
    std::stringstream entry;
    std::stringstream particleString;
    ParticleList const &pl = getParticles();
    G4int numberOfParticles = 0;
    for(ParticleIter i=pl.begin(), e=pl.end(); i!=e; ++i) {
      numberOfParticles++;
      particleString << (*i)->getID() << " ";
    }
    if(numberOfParticles == 1) particleString << "-1";
    entry << getID() << " "
      << getType() << " "
      << getTime() << " "
      << particleString.str();
    return entry.str();
  }

References G4INCL::IAvatar::getID(), G4INCL::IAvatar::getParticles(), G4INCL::IAvatar::getTime(), G4INCL::IAvatar::getType(), and numberOfParticles.

Field Documentation

◆ backupParticle1

G4ThreadLocal Particle * G4INCL::InteractionAvatar::backupParticle1 = NULL

staticprotected

Definition at line 108 of file G4INCLInteractionAvatar.hh.

Referenced by deleteBackupParticles(), G4INCL::BinaryCollisionAvatar::postInteraction(), and preInteractionBlocking().

◆ backupParticle2

G4ThreadLocal Particle * G4INCL::InteractionAvatar::backupParticle2 = NULL

staticprotected

Definition at line 108 of file G4INCLInteractionAvatar.hh.

Referenced by deleteBackupParticles(), G4INCL::BinaryCollisionAvatar::postInteraction(), and preInteractionBlocking().

◆ boostVector

ThreeVector G4INCL::InteractionAvatar::boostVector

protected

Definition at line 109 of file G4INCLInteractionAvatar.hh.

Referenced by enforceEnergyConservation(), G4INCL::BinaryCollisionAvatar::getChannel(), postInteraction(), preInteraction(), and G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta().

◆ created

ParticleList G4INCL::InteractionAvatar::created

protected

Definition at line 223 of file G4INCLInteractionAvatar.hh.

Referenced by postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ Destroyed

ParticleList G4INCL::InteractionAvatar::Destroyed

protected

Definition at line 223 of file G4INCLInteractionAvatar.hh.

Referenced by postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ ID

long G4INCL::IAvatar::ID

privateinherited

Definition at line 98 of file G4INCLIAvatar.hh.

Referenced by G4INCL::IAvatar::IAvatar().

◆ isPiN

G4bool G4INCL::InteractionAvatar::isPiN

protected

Definition at line 111 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), and shouldUseLocalEnergy().

◆ locEAccuracy

const G4double G4INCL::InteractionAvatar::locEAccuracy = 1.E-4

static

Target accuracy in the determination of the local-energy Q-value.

Definition at line 68 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), and G4INCL::InteractionAvatar::ViolationEEnergyFunctor::setParticleEnergy().

◆ maxIterLocE

const G4int G4INCL::InteractionAvatar::maxIterLocE = 50

static

Max number of iterations for the determination of the local-energy Q-value.

Definition at line 70 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), and G4INCL::InteractionAvatar::ViolationEEnergyFunctor::setParticleEnergy().

◆ modified

ParticleList G4INCL::InteractionAvatar::modified

protected

Definition at line 223 of file G4INCLInteractionAvatar.hh.

Referenced by enforceEnergyConservation(), postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ modifiedAndCreated

ParticleList G4INCL::InteractionAvatar::modifiedAndCreated

protected

Definition at line 223 of file G4INCLInteractionAvatar.hh.

Referenced by enforceEnergyConservation(), postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ ModifiedAndDestroyed

ParticleList G4INCL::InteractionAvatar::ModifiedAndDestroyed

protected

Definition at line 223 of file G4INCLInteractionAvatar.hh.

Referenced by postInteraction(), and G4INCL::DecayAvatar::postInteraction().

◆ nextID

G4ThreadLocal long G4INCL::IAvatar::nextID = 1

staticprivateinherited

Definition at line 100 of file G4INCLIAvatar.hh.

Referenced by G4INCL::IAvatar::IAvatar().

◆ oldTotalEnergy

G4double G4INCL::InteractionAvatar::oldTotalEnergy

protected

Definition at line 110 of file G4INCLInteractionAvatar.hh.

Referenced by postInteraction(), G4INCL::DecayAvatar::postInteraction(), and preInteractionBlocking().

◆ oldXSec

G4double G4INCL::InteractionAvatar::oldXSec

protected

Definition at line 110 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::postInteraction(), and preInteractionBlocking().

◆ particle1

Particle* G4INCL::InteractionAvatar::particle1

protected

Definition at line 107 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::dump(), G4INCL::DecayAvatar::dump(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::BinaryCollisionAvatar::getParticles(), G4INCL::DecayAvatar::getParticles(), G4INCL::DecayAvatar::postInteraction(), G4INCL::BinaryCollisionAvatar::preInteraction(), preInteraction(), and preInteractionBlocking().

◆ particle2

Particle * G4INCL::InteractionAvatar::particle2

protected

Definition at line 107 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::dump(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::BinaryCollisionAvatar::getParticles(), G4INCL::BinaryCollisionAvatar::preInteraction(), preInteraction(), preInteractionBlocking(), and restoreParticles().

◆ theNucleus

Nucleus* G4INCL::InteractionAvatar::theNucleus

protected

Definition at line 106 of file G4INCLInteractionAvatar.hh.

Referenced by bringParticleInside(), enforceEnergyConservation(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::BinaryCollisionAvatar::postInteraction(), postInteraction(), G4INCL::DecayAvatar::postInteraction(), preInteractionLocalEnergy(), G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), G4INCL::InteractionAvatar::ViolationEEnergyFunctor::setParticleEnergy(), and shouldUseLocalEnergy().

◆ theTime

G4double G4INCL::IAvatar::theTime

protectedinherited

Definition at line 102 of file G4INCLIAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::dump(), G4INCL::DecayAvatar::dump(), G4INCL::ParticleEntryAvatar::dump(), and G4INCL::SurfaceAvatar::dump().

◆ type

AvatarType G4INCL::IAvatar::type

privateinherited

Definition at line 99 of file G4INCLIAvatar.hh.

◆ violationEFunctor

RootFunctor* G4INCL::InteractionAvatar::violationEFunctor

private

Definition at line 208 of file G4INCLInteractionAvatar.hh.

Referenced by enforceEnergyConservation().

◆ weight

G4double G4INCL::InteractionAvatar::weight

protected

Definition at line 112 of file G4INCLInteractionAvatar.hh.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), and postInteraction().

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

source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLInteractionAvatar.hh
source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLInteractionAvatar.cc

Data Structures

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ InteractionAvatar() [1/2]

◆ InteractionAvatar() [2/2]

◆ ~InteractionAvatar()

Member Function Documentation

◆ bringParticleInside()

◆ deleteBackupParticles()

◆ dump()

◆ enforceEnergyConservation()

◆ fillFinalState()

◆ getChannel()

◆ getFinalState()

◆ getID()

◆ getParticles()

◆ getTime()

◆ getType()

◆ isACollision()

◆ isADecay()

◆ postInteraction()

◆ preInteraction()

◆ preInteractionBlocking()

◆ preInteractionLocalEnergy()

◆ restoreParticles()

◆ setType()

◆ shouldUseLocalEnergy()

◆ toString()

Field Documentation

◆ backupParticle1

◆ backupParticle2

◆ boostVector

◆ created

◆ Destroyed

◆ ID

◆ isPiN

◆ locEAccuracy

◆ maxIterLocE

◆ modified

◆ modifiedAndCreated

◆ ModifiedAndDestroyed

◆ nextID

◆ oldTotalEnergy

◆ oldXSec

◆ particle1

◆ particle2

◆ theNucleus

◆ theTime

◆ type

◆ violationEFunctor

◆ weight