G4INCL::DecayAvatar Class Reference

#include <G4INCLDecayAvatar.hh>

Inheritance diagram for G4INCL::DecayAvatar:

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

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 = 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 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 ()
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 inherited from G4INCL::InteractionAvatar
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 inherited from G4INCL::InteractionAvatar
static G4ThreadLocal Particle *	backupParticle1 = NULL
static G4ThreadLocal Particle *	backupParticle2 = NULL

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 53 of file G4INCLDecayAvatar.hh.

Constructor & Destructor Documentation

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

Definition at line 47 of file G4INCLDecayAvatar.cc.

References G4INCL::DecayAvatarType, and G4INCL::IAvatar::setType().

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

G4INCL::DecayAvatar::~DecayAvatar ( )

virtual

Definition at line 54 of file G4INCLDecayAvatar.cc.

                            {

  }

Member Function Documentation

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

virtual

Implements G4INCL::IAvatar.

Definition at line 183 of file G4INCLDecayAvatar.cc.

References G4INCL::Particle::dump(), G4INCL::InteractionAvatar::particle1, and G4INCL::IAvatar::theTime.

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

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

virtual

Implements G4INCL::InteractionAvatar.

Definition at line 58 of file G4INCLDecayAvatar.cc.

References INCL_DEBUG, G4INCL::Particle::isDelta(), and G4INCL::InteractionAvatar::particle1.

                                          {
    if(particle1->isDelta()) {
      INCL_DEBUG("DeltaDecayChannel chosen." << std::endl);
      return new DeltaDecayChannel(particle1, incidentDirection);
    }
    else
      return NULL;
  }

G4INCL::FinalState* G4INCL::DecayAvatar::getFinalState

(

)

ParticleList G4INCL::DecayAvatar::getParticles ( ) const

inlinevirtual

Implements G4INCL::IAvatar.

Definition at line 64 of file G4INCLDecayAvatar.hh.

References G4INCL::InteractionAvatar::particle1.

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

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

virtual

Implements G4INCL::IAvatar.

Definition at line 71 of file G4INCLDecayAvatar.cc.

References G4INCL::FinalState::addModifiedParticle(), G4INCL::InteractionAvatar::enforceEnergyConservation(), G4INCL::Store::getBook(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Nucleus::getStore(), G4INCL::FinalState::getValidity(), INCL_DEBUG, G4INCL::Book::incrementAcceptedDecays(), G4INCL::Book::incrementBlockedDecays(), G4INCL::Pauli::isCDPPBlocked(), G4INCL::FinalState::makeNoEnergyConservation(), G4INCL::FinalState::makePauliBlocked(), G4INCL::NoEnergyConservationFS, G4INCL::InteractionAvatar::oldTotalEnergy, G4INCL::InteractionAvatar::particle1, G4INCL::ParticleBelowFermiFS, G4INCL::ParticleBelowZeroFS, G4INCL::PauliBlockedFS, G4INCL::InteractionAvatar::postInteraction(), G4INCL::InteractionAvatar::restoreParticles(), G4INCL::FinalState::setBlockedDelta(), G4INCL::FinalState::setTotalEnergyBeforeInteraction(), G4INCL::InteractionAvatar::theNucleus, and G4INCL::ValidFS.

                                                         {
    // 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)
      fs = 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->setBlockedDelta(particle1);

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

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

        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;

          FinalState *fsBlocked = new FinalState;
          delete fs;
          fsBlocked->makeNoEnergyConservation();
          fsBlocked->setTotalEnergyBeforeInteraction(0.0);

          return fsBlocked; // 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." << std::endl);
        }
      } else {
        INCL_DEBUG("Enforcing energy conservation: success!" << std::endl);
      }

      if(theNucleus) {
        ParticleList modified = fs->getModifiedParticles();

        // Copy the final state, but don't include the pion (as if it had been
        // emitted right away).
        FinalState *emissionFS = new FinalState;
        for(ParticleIter i=modified.begin(), e=modified.end(); i!=e; ++i)
          emissionFS->addModifiedParticle(*i);

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

        if(isCDPPBlocked) {
          INCL_DEBUG("CDPP: Blocked!" << std::endl);

          // Restore the state of both 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;
      delete emissionFS;

          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 (energy conservation enforced and CDPP satisfied),
        // delete the auxiliary final state
        delete emissionFS;

      }
    }

    // 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 fs;
  }

void G4INCL::DecayAvatar::preInteraction ( )

virtual

Implements G4INCL::IAvatar.

Definition at line 67 of file G4INCLDecayAvatar.cc.

References G4INCL::InteractionAvatar::preInteraction().

                                   {
    InteractionAvatar::preInteraction();
  }

---

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

• G4INCLDecayAvatar.hh
• G4INCLDecayAvatar.cc