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)
	DecayAvatar (G4INCL::Particle *aParticle, G4INCL::Particle *bParticle, 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 ()
virtual G4INCL::IChannel *	getChannel ()=0
FinalState *	getFinalState ()
void	fillFinalState (FinalState *fs)
virtual void	preInteraction ()=0
virtual void	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 Member Functions inherited from G4INCL::InteractionAvatar
static void	deleteBackupParticles ()
	Release the memory allocated for the backup particles.
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 ()=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 ()
void	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
Nucleus *	theNucleus
Particle *	particle1
Particle *	particle2
ThreeVector	boostVector
G4double	oldTotalEnergy
G4double	oldXSec
G4bool	isPiN
G4double	weight
ParticleList	modified
ParticleList	created
ParticleList	modifiedAndCreated
ParticleList	Destroyed
ParticleList	ModifiedAndDestroyed
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 55 of file G4INCLDecayAvatar.hh.

Constructor & Destructor Documentation

DecayAvatar() [1/2]

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

Definition at line 52 of file G4INCLDecayAvatar.cc.

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

DecayAvatar() [2/2]

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

Definition at line 59 of file G4INCLDecayAvatar.cc.

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

~DecayAvatar()

G4INCL::DecayAvatar::~DecayAvatar ( )

virtual

Definition at line 66 of file G4INCLDecayAvatar.cc.

                            {
 
  }

Member Function Documentation

dump()

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

virtual

Implements G4INCL::IAvatar.

Definition at line 207 of file G4INCLDecayAvatar.cc.

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

fillFinalState()

void G4INCL::DecayAvatar::fillFinalState

(

FinalState *

fs

)

getChannel()

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

virtual

Implements G4INCL::InteractionAvatar.

Definition at line 70 of file G4INCLDecayAvatar.cc.

                                          {
    if(!particle2){
       if(particle1->isDelta()) {
         INCL_DEBUG("DeltaDecayChannel chosen." << '\n');
         return new DeltaDecayChannel(particle1, incidentDirection);
       }
       else if(particle1->isEta() || particle1->isOmega()) {
         INCL_DEBUG("PionResonanceDecayChannel chosen." << '\n');
         return new PionResonanceDecayChannel(particle1, incidentDirection);
       }
       else if(particle1->getType() == SigmaZero) {
         INCL_DEBUG("SigmaZeroDecayChannel chosen." << '\n');
         return new SigmaZeroDecayChannel(particle1, incidentDirection);
       }
       else if(particle1->getType() == KZero || particle1->getType() == KZeroBar) {
         INCL_DEBUG("NeutralKaonDecayChannel chosen." << '\n');
         return new NeutralKaonDecayChannel(particle1);
       }
    }
    else if(((particle1->isAntiKaon() || particle1->isSigma()) && particle2->isNucleon()) || ((particle2->isAntiKaon() || particle2->isSigma()) && particle1->isNucleon())){
      INCL_DEBUG("StrangeAbsorbtion." << '\n');
      return new StrangeAbsorbtionChannel(particle1, particle2);
    }
    return NULL;
  }

getParticles()

ParticleList G4INCL::DecayAvatar::getParticles ( ) const

inlinevirtual

Implements G4INCL::IAvatar.

Definition at line 67 of file G4INCLDecayAvatar.hh.

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

postInteraction()

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

virtual

Implements G4INCL::IAvatar.

Definition at line 100 of file G4INCLDecayAvatar.cc.

                                                  {
    // Make sure we have at least two particles in the final state
    // Removed because of neutral kaon decay
    
// assert((fs->getModifiedParticles().size() + fs->getCreatedParticles().size() - fs->getDestroyedParticles().size() >= 2) || ((*fs->getModifiedParticles().begin())->getType() == KShort || (*fs->getModifiedParticles().begin())->getType() == KLong ));
    //assert((fs->getModifiedParticles().size() + fs->getCreatedParticles().size() - fs->getDestroyedParticles().size() >= 1));
    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
      modified = fs->getModifiedParticles();
      created = fs->getCreatedParticles();
      Destroyed = fs->getDestroyedParticles();
      modifiedAndCreated = modified;
      modifiedAndCreated.insert(modifiedAndCreated.end(), created.begin(), created.end());
      ModifiedAndDestroyed = modified;
      ModifiedAndDestroyed.insert(ModifiedAndDestroyed.end(), Destroyed.begin(), Destroyed.end());
      
      std::vector<G4int> newBiasCollisionVector;
      newBiasCollisionVector = ModifiedAndDestroyed.getParticleListBiasVector();
      for(ParticleIter i=modifiedAndCreated.begin(), e=modifiedAndCreated.end(); i!=e; ++i ) {
        (*i)->setBiasCollisionVector(newBiasCollisionVector);
      }
      // 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;
  }

preInteraction()

void G4INCL::DecayAvatar::preInteraction ( )

virtual

Implements G4INCL::IAvatar.

Definition at line 96 of file G4INCLDecayAvatar.cc.

                                   {
    InteractionAvatar::preInteraction();
  }

---

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

• G4INCLDecayAvatar.hh
• G4INCLDecayAvatar.cc