G4INCL::DecayAvatar Class Reference

#include <G4INCLDecayAvatar.hh>

Inheritance diagram for G4INCL::DecayAvatar:

Collaboration diagram for G4INCL::DecayAvatar:

Public Member Functions
	DecayAvatar (G4INCL::Particle *aParticle, G4double time, G4INCL::Nucleus *aNucleus, G4bool force=false)
virtual	~DecayAvatar ()
IChannel *	getChannel ()
void	fillFinalState (FinalState *fs)
virtual void	preInteraction ()
virtual void	postInteraction (FinalState *fs)
ParticleList	getParticles () const
std::string	dump () const
Public Member Functions inherited from G4INCL::InteractionAvatar
	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 ()
FinalState *	getFinalState ()
void	fillFinalState (FinalState *fs)
G4double	getTime () const
AvatarType	getType () const
G4bool	isACollision () const
G4bool	isADecay () const
void	setType (AvatarType t)
long	getID () const
std::string	toString ()

Private Member Functions
	INCL_DECLARE_ALLOCATION_POOL (DecayAvatar)

Private Attributes
G4bool	forced
ThreeVector const	incidentDirection

Additional Inherited Members
Static Public Member Functions inherited from G4INCL::InteractionAvatar
static void	deleteBackupParticles ()
	Release the memory allocated for the backup particles. More...
Static Public Attributes inherited from G4INCL::InteractionAvatar
static const G4double	locEAccuracy
	Target accuracy in the determination of the local-energy Q-value. More...
static const G4int	maxIterLocE
	Max number of iterations for the determination of the local-energy Q-value. More...
Protected Member Functions inherited from G4INCL::InteractionAvatar
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 ()
void	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...
	INCL_DECLARE_ALLOCATION_POOL (InteractionAvatar)
Protected Member Functions inherited from G4INCL::IAvatar
	INCL_DECLARE_ALLOCATION_POOL (IAvatar)
Protected Attributes inherited from G4INCL::InteractionAvatar
Nucleus *	theNucleus
Particle *	particle1
Particle *	particle2
ThreeVector	boostVector
G4double	oldTotalEnergy
G4double	oldXSec
G4bool	isPiN
ParticleList	modified
ParticleList	created
ParticleList	modifiedAndCreated
Protected Attributes inherited from G4INCL::IAvatar
G4double	theTime
Static Protected Attributes inherited from G4INCL::InteractionAvatar
static G4ThreadLocal Particle *	backupParticle1
static G4ThreadLocal Particle *	backupParticle2

Detailed Description

Decay avatar

The reflection avatar is created when a particle reaches the boundary of the nucleus. At this point it can either be reflected from the boundary or exit the nucleus.

Definition at line 55 of file G4INCLDecayAvatar.hh.

Constructor & Destructor Documentation

DecayAvatar()

G4INCL::DecayAvatar::DecayAvatar	(	G4INCL::Particle *	aParticle,
		G4double	time,
		G4INCL::Nucleus *	aNucleus,
		G4bool	force = false
	)

Definition at line 48 of file G4INCLDecayAvatar.cc.

    : InteractionAvatar(time, n, aParticle), forced(force),
      incidentDirection(aParticle->getMomentum())
  {
    setType(DecayAvatarType);
  }

Here is the call graph for this function:

~DecayAvatar()

G4INCL::DecayAvatar::~DecayAvatar ( )

virtual

Definition at line 55 of file G4INCLDecayAvatar.cc.

                            {

  }

Member Function Documentation

dump()

std::string G4INCL::DecayAvatar::dump ( ) const

virtual

Implements G4INCL::IAvatar.

Definition at line 169 of file G4INCLDecayAvatar.cc.

                                    {
    std::stringstream ss;
    ss << "(avatar " << theTime << " 'decay" << '\n'
      << "(list " << '\n'
      << particle1->dump()
      << "))" << '\n';
    return ss.str();
  }

Here is the call graph for this function:

Here is the caller graph for this function:

fillFinalState()

void G4INCL::DecayAvatar::fillFinalState

(

FinalState *

fs

)

getChannel()

G4INCL::IChannel * G4INCL::DecayAvatar::getChannel ( )

virtual

Implements G4INCL::InteractionAvatar.

Definition at line 59 of file G4INCLDecayAvatar.cc.

                                          {
    if(particle1->isDelta()) {
      INCL_DEBUG("DeltaDecayChannel chosen." << '\n');
      return new DeltaDecayChannel(particle1, incidentDirection);
    }
    else
      return NULL;
  }

Here is the call graph for this function:

getParticles()

ParticleList G4INCL::DecayAvatar::getParticles ( ) const

inlinevirtual

Implements G4INCL::IAvatar.

Definition at line 66 of file G4INCLDecayAvatar.hh.

                                      {
      ParticleList theParticleList;
      theParticleList.push_back(particle1);
      return theParticleList;
    }

Here is the call graph for this function:

INCL_DECLARE_ALLOCATION_POOL()

G4INCL::DecayAvatar::INCL_DECLARE_ALLOCATION_POOL ( DecayAvatar  )

private

postInteraction()

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

virtual

Implements G4INCL::IAvatar.

Definition at line 72 of file G4INCLDecayAvatar.cc.

                                                  {
    // Make sure we have at least two particles in the final state
// assert(fs->getModifiedParticles().size() + fs->getCreatedParticles().size() - fs->getDestroyedParticles().size() >= 2);

    if(!forced) { // Normal decay

      // Call the postInteraction method of the parent class
      // (provides Pauli blocking and enforces energy conservation)
      InteractionAvatar::postInteraction(fs);

      if(fs->getValidity() == PauliBlockedFS)
        /* If the decay was Pauli-blocked, make sure the propagation model
         * generates a new decay avatar on the next call to propagate().
         *
         * \bug{Note that we don't generate new decay avatars for deltas that
         * could not satisfy energy conservation. This is in keeping with
         * INCL4.6, but doesn't seem to make much sense to me (DM), as energy
         * conservation can be impossible to satisfy due to weird local-energy
         * conditions, for example, that evolve with time.}
         */
        fs->addModifiedParticle(particle1);

    } else { // Forced decay
      created = fs->getCreatedParticles();
      modified = fs->getModifiedParticles();

      // Try to enforce energy conservation
      fs->setTotalEnergyBeforeInteraction(oldTotalEnergy);
      const G4bool success = enforceEnergyConservation(fs);
      if(!success) {
        INCL_DEBUG("Enforcing energy conservation: failed!" << '\n');

        if(theNucleus) {
          // 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.
        } else {
          // If there is no nucleus we have to continue anyway, even if energy
          // conservation failed. We cannot afford producing unphysical
          // remnants.
          INCL_DEBUG("No nucleus, continuing anyway." << '\n');
        }
      } else {
        INCL_DEBUG("Enforcing energy conservation: success!" << '\n');
      }

      if(theNucleus) {
        // Test CDPP blocking
        G4bool isCDPPBlocked = Pauli::isCDPPBlocked(created, theNucleus);

        if(isCDPPBlocked) {
          INCL_DEBUG("CDPP: Blocked!" << '\n');

          // Restore the state of both 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 there is a nucleus, increment the counters
    if(theNucleus) {
      switch(fs->getValidity()) {
        case PauliBlockedFS:
          theNucleus->getStore()->getBook().incrementBlockedDecays();
          break;
        case NoEnergyConservationFS:
        case ParticleBelowFermiFS:
        case ParticleBelowZeroFS:
          break;
        case ValidFS:
          theNucleus->getStore()->getBook().incrementAcceptedDecays();
      }
    }
    return;
  }

Here is the call graph for this function:

preInteraction()

void G4INCL::DecayAvatar::preInteraction ( )

virtual

Implements G4INCL::IAvatar.

Definition at line 68 of file G4INCLDecayAvatar.cc.

                                   {
    InteractionAvatar::preInteraction();
  }

Here is the call graph for this function:

Member Data Documentation

forced

G4bool G4INCL::DecayAvatar::forced

private

Definition at line 74 of file G4INCLDecayAvatar.hh.

incidentDirection

ThreeVector const G4INCL::DecayAvatar::incidentDirection

private

Definition at line 75 of file G4INCLDecayAvatar.hh.

---

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

• G4INCLDecayAvatar.hh
• G4INCLDecayAvatar.cc