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
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 ()
virtual G4INCL::IChannel *	getChannel () const =0
G4INCL::FinalState *	getFinalState ()
virtual void	preInteraction ()=0
virtual FinalState *	postInteraction (FinalState *)=0
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 ()

Additional Inherited Members
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.
static const G4int	maxIterLocE = 50
	Max number of iterations for the determination of the local-energy Q-value.
Protected Member Functions inherited from G4INCL::InteractionAvatar
virtual G4INCL::IChannel *	getChannel () const =0
G4bool	bringParticleInside (Particle *const p)
void	preInteractionLocalEnergy (Particle *const p)
	Apply local-energy transformation, if appropriate.
void	preInteractionBlocking ()
	Store the state of the particles before the interaction.
void	preInteraction ()
FinalState *	postInteraction (FinalState *)
void	restoreParticles () const
	Restore the state of both particles.
G4bool	shouldUseLocalEnergy () const
	true if the given avatar should use local energy
G4bool	enforceEnergyConservation (FinalState *const fs)
	Enforce energy conservation.
Protected Attributes inherited from G4INCL::InteractionAvatar
G4INCL::Nucleus *	theNucleus
G4INCL::Particle *	particle1
G4INCL::Particle *	particle2
ThreeVector	boostVector
ParticleType	oldParticle1Type
ParticleType	oldParticle2Type
G4double	oldParticle1Energy
G4double	oldParticle2Energy
G4double	oldTotalEnergy
G4double	oldXSec
G4double	oldParticle1Potential
G4double	oldParticle2Potential
G4double	oldParticle1Mass
G4double	oldParticle2Mass
G4double	oldParticle1Helicity
G4double	oldParticle2Helicity
ThreeVector	oldParticle1Momentum
ThreeVector	oldParticle2Momentum
ThreeVector	oldParticle1Position
ThreeVector	oldParticle2Position
G4bool	isPiN
Protected Attributes inherited from G4INCL::IAvatar
G4double	theTime

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 49 of file G4INCLDecayAvatar.cc.

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

~DecayAvatar()

G4INCL::DecayAvatar::~DecayAvatar ( )

virtual

Definition at line 56 of file G4INCLDecayAvatar.cc.

                            {
 
  }

Member Function Documentation

dump()

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

virtual

Implements G4INCL::IAvatar.

Definition at line 186 of file G4INCLDecayAvatar.cc.

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

getChannel()

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

virtual

Implements G4INCL::InteractionAvatar.

Definition at line 60 of file G4INCLDecayAvatar.cc.

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

getFinalState()

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

(

)

const

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

postInteraction()

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

virtual

Implements G4INCL::IAvatar.

Definition at line 74 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)
      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;
  }

preInteraction()

void G4INCL::DecayAvatar::preInteraction ( )

virtual

Implements G4INCL::IAvatar.

Definition at line 70 of file G4INCLDecayAvatar.cc.

                                   {
    InteractionAvatar::preInteraction();
  }

---

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

• G4INCLDecayAvatar.hh
• G4INCLDecayAvatar.cc