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 () const
G4INCL::FinalState *	getFinalState () const
virtual void	preInteraction ()
virtual FinalState *	postInteraction (FinalState *)
ParticleList	getParticles () const
std::string	dump () const

---

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 184 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

(

)

const [virtual]

Implements G4INCL::InteractionAvatar.

Definition at line 58 of file G4INCLDecayAvatar.cc.

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

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

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

(

)

const

ParticleList G4INCL::DecayAvatar::getParticles

(

)

const [inline, virtual]

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 *

)

[virtual]

Reimplemented from G4INCL::InteractionAvatar.

Definition at line 72 of file G4INCLDecayAvatar.cc.

References G4INCL::FinalState::addModifiedParticle(), DEBUG, G4INCL::InteractionAvatar::enforceEnergyConservation(), G4INCL::Store::getBook(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Nucleus::getStore(), G4INCL::FinalState::getValidity(), 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) {
        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(); i != created.end(); ++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.
          DEBUG("No nucleus, continuing anyway." << std::endl);
        }
      } else {
        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(); i!=modified.end(); ++i)
          emissionFS->addModifiedParticle(*i);

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

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

          // Restore the state of both particles
          restoreParticles();

          // Delete newly created particles
          for( ParticleIter i = created.begin(); i != created.end(); ++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.
        }
        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]

Reimplemented from G4INCL::InteractionAvatar.

Definition at line 68 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

---