#include <G4INCLInteractionAvatar.hh>

Inheritance diagram for G4INCL::InteractionAvatar:

Public Member Functions
	InteractionAvatar (G4double, G4INCL::Nucleus , G4INCL::Particle )

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

virtual	~InteractionAvatar ()

Public Member Functions inherited from G4INCL::IAvatar
	IAvatar ()

	IAvatar (G4double time)

virtual	~IAvatar ()

G4INCL::FinalState *	getFinalState ()

G4double	getTime () const

virtual ParticleList	getParticles () const =0

virtual std::string	dump () const =0

AvatarType	getType () const

G4bool	isACollision () const

G4bool	isADecay () const

void	setType (AvatarType t)

long	getID () const

std::string	toString ()

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
virtual G4INCL::IChannel *	getChannel ()=0

G4bool	bringParticleInside (Particle *const p)

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

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

void	preInteraction ()

FinalState *	postInteraction (FinalState *)

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

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

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

Protected Attributes
Nucleus *	theNucleus

Particle *	particle1

Particle *	particle2

ThreeVector	boostVector

G4double	oldTotalEnergy

G4double	oldXSec

G4bool	isPiN

Protected Attributes inherited from G4INCL::IAvatar
G4double	theTime

Static Protected Attributes
static G4ThreadLocal Particle *	backupParticle1 = NULL

static G4ThreadLocal Particle *	backupParticle2 = NULL

Detailed Description

Definition at line 59 of file G4INCLInteractionAvatar.hh.

Constructor & Destructor Documentation

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

Definition at line 64 of file G4INCLInteractionAvatar.cc.

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

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

Definition at line 72 of file G4INCLInteractionAvatar.cc.

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

G4INCL::InteractionAvatar::~InteractionAvatar ( )

virtual

Definition at line 81 of file G4INCLInteractionAvatar.cc.

81 {

82 }

Member Function Documentation

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

protected

Definition at line 137 of file G4INCLInteractionAvatar.cc.

References G4INCL::Particle::getPosition(), G4INCL::Nucleus::getSurfaceRadius(), INCL_DEBUG, G4INCL::ThreeVector::mag(), G4INCL::ThreeVector::mag2(), G4INCL::Particle::rpCorrelate(), G4INCL::Particle::setPosition(), and theNucleus.

Referenced by postInteraction().

                                                                   {
     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 )
     {
       pos *= std::sqrt(r*r*0.99/pos2);
       pos2 = pos.mag2();
       iterations++;
     }
     if( iterations < maxIterations)
     {
       INCL_DEBUG("Particle position vector length was : " << p->getPosition().mag() << ", rescaled to: " << pos.mag() << std::endl);
       p->setPosition(pos);
       return true;
     }
     else
       return false;
   }

void G4INCL::InteractionAvatar::deleteBackupParticles ( )

static

Release the memory allocated for the backup particles.

Definition at line 84 of file G4INCLInteractionAvatar.cc.

References backupParticle1, and backupParticle2.

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

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

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 345 of file G4INCLInteractionAvatar.cc.

References boostVector, G4INCL::ParticleTable::effectiveDeltaDecayThreshold, G4INCL::Store::getBook(), G4INCL::FinalState::getCreatedParticles(), G4INCL::Particle::getMass(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Nucleus::getStore(), INCL_DEBUG, G4INCL::Book::incrementEnergyViolationInteraction(), shouldUseLocalEnergy(), G4INCL::RootFinder::solve(), G4INCL::RootFinder::Solution::success, theNucleus, and G4INCL::RootFinder::Solution::x.

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

                                                                            {
     // Set up the violationE calculation
     ParticleList modified = fs->getModifiedParticles();
     const G4bool manyBodyFinalState = (modified.size() + fs->getCreatedParticles().size() > 1);
     if(manyBodyFinalState)
       violationEFunctor = new ViolationEMomentumFunctor(theNucleus, fs, &boostVector, shouldUseLocalEnergy());
     else {
       Particle const * 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::effectiveDeltaDecayThreshold)
         return false;
       violationEFunctor = new ViolationEEnergyFunctor(theNucleus, fs, 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." << std::endl);
       theNucleus->getStore()->getBook().incrementEnergyViolationInteraction();
     }
     delete violationEFunctor;
     violationEFunctor = NULL;
     return theSolution.success;
   }

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

protectedpure virtual

Implements G4INCL::IAvatar.

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

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

protectedvirtual

Implements G4INCL::IAvatar.

Definition at line 166 of file G4INCLInteractionAvatar.cc.

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

                                                                {
     INCL_DEBUG("postInteraction: final state: " << std::endl << fs->print() << std::endl);
     ParticleList modified = fs->getModifiedParticles();
     ParticleList modifiedAndCreated = modified;
     ParticleList created = fs->getCreatedParticles();
     modifiedAndCreated.insert(modifiedAndCreated.end(), created.begin(), created.end());
 
     if(!isPiN) {
       // Boost back to lab
       for(ParticleIter i=modifiedAndCreated.begin(), e=modifiedAndCreated.end(); i!=e; ++i )
         (*i)->boost(-boostVector);
     }
 
     // If there is no Nucleus, just return
     if(!theNucleus) return fs;
 
     // Mark pions 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)->getPosition().mag() > theNucleus->getSurfaceRadius(*i)) {
         (*i)->makeParticipant();
         (*i)->setOutOfWell();
         fs->addOutgoingParticle(*i);
         INCL_DEBUG("Pion was created outside its potential well." << std::endl
             << (*i)->print());
       }
 
     // Try to enforce energy conservation
     fs->setTotalEnergyBeforeInteraction(oldTotalEnergy);
     G4bool success = true;
     if(!isPiN || shouldUseLocalEnergy())
       success = enforceEnergyConservation(fs);
     if(!success) {
       INCL_DEBUG("Enforcing energy conservation: failed!" << std::endl);
 
       // 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;
 
       FinalState *fsBlocked = new FinalState;
       delete fs;
       fsBlocked->makeNoEnergyConservation();
       fsBlocked->setTotalEnergyBeforeInteraction(0.0);
 
       return fsBlocked; // Interaction is blocked. Return an empty final state.
     }
     INCL_DEBUG("Enforcing energy conservation: success!" << std::endl);
 
     INCL_DEBUG("postInteraction after energy conservation: final state: " << std::endl << fs->print() << std::endl);
 
     // 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::effectiveDeltaDecayThreshold) {
         INCL_DEBUG("Mass of the produced delta below decay threshold; forbidding collision. deltaMass=" <<
             (*i)->getMass() << std::endl);
 
         // 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;
 
         FinalState *fsBlocked = new FinalState;
         delete fs;
         fsBlocked->makeNoEnergyConservation();
         fsBlocked->setTotalEnergyBeforeInteraction(0.0);
 
         return fsBlocked; // Interaction is blocked. Return an empty final state.
       }
 
     INCL_DEBUG("Random seeds before Pauli blocking: " << Random::getSeeds() << std::endl);
     // Test Pauli blocking
     G4bool isBlocked = Pauli::isBlocked(modifiedAndCreated, theNucleus);
 
     if(isBlocked) {
       INCL_DEBUG("Pauli: Blocked!" << std::endl);
 
       // 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;
 
       FinalState *fsBlocked = new FinalState;
       delete fs;
       fsBlocked->makePauliBlocked();
       fsBlocked->setTotalEnergyBeforeInteraction(0.0);
 
       return fsBlocked; // Interaction is blocked. Return an empty final state.
     }
     INCL_DEBUG("Pauli: Allowed!" << std::endl);
 
     // Test CDPP blocking
     G4bool isCDPPBlocked = Pauli::isCDPPBlocked(created, theNucleus);
 
     if(isCDPPBlocked) {
       INCL_DEBUG("CDPP: Blocked!" << std::endl);
 
       // 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;
 
       FinalState *fsBlocked = new FinalState;
       delete fs;
       fsBlocked->makePauliBlocked();
       fsBlocked->setTotalEnergyBeforeInteraction(0.0);
 
       return fsBlocked; // Interaction is blocked. Return an empty final state.
     }
     INCL_DEBUG("CDPP: Allowed!" << std::endl);
 
     // 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!" << std::endl);
       }
     }
 
     // Collision accepted!
     for(ParticleIter i=modifiedAndCreated.begin(), e=modifiedAndCreated.end(); i!=e; ++i ) {
       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:" << std::endl
               << (*i)->print() << std::endl);
           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 fs;
   }

void G4INCL::InteractionAvatar::preInteraction ( )

protectedvirtual

Implements G4INCL::IAvatar.

Definition at line 119 of file G4INCLInteractionAvatar.cc.

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

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

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

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 92 of file G4INCLInteractionAvatar.cc.

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

Referenced by preInteraction().

                                                  {
     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();
     }
   }

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 112 of file G4INCLInteractionAvatar.cc.

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

Referenced by preInteraction().

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

void G4INCL::InteractionAvatar::restoreParticles ( ) const

protected

Restore the state of both particles.

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

Definition at line 339 of file G4INCLInteractionAvatar.cc.

References particle2.

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

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

G4bool G4INCL::InteractionAvatar::shouldUseLocalEnergy ( ) const

inlineprotected

true if the given avatar should use local energy

Definition at line 102 of file G4INCLInteractionAvatar.hh.

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(), postInteraction(), and preInteractionLocalEnergy().

                                           {
         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);
       }