G4INCL::AntinucleiAtrestEntryChannel Class Reference

#include <G4INCLAntinucleiAtrestEntryChannel.hh>

Inheritance diagram for G4INCL::AntinucleiAtrestEntryChannel:

Public Member Functions
	AntinucleiAtrestEntryChannel (Nucleus *n, Cluster *ac, ThreeVector pos1, ThreeVector pos2)
	AntinucleiAtrestEntryChannel (Nucleus *n, Particle *p)
virtual	~AntinucleiAtrestEntryChannel ()
void	fillFinalState (FinalState *fs)
ThreeVector	getAnnihilationPosition (ThreeVector nbarPos, ThreeVector pbarPos)
ParticleList	makeMesonStar ()
IAvatarList	bringMesonStar (ParticleList const &pL, Nucleus *const n)
Public Member Functions inherited from G4INCL::IChannel
	IChannel ()
virtual	~IChannel ()
FinalState *	getFinalState ()

Detailed Description

Definition at line 53 of file G4INCLAntinucleiAtrestEntryChannel.hh.

Constructor & Destructor Documentation

AntinucleiAtrestEntryChannel() [1/2]

G4INCL::AntinucleiAtrestEntryChannel::AntinucleiAtrestEntryChannel	(	Nucleus *	n,
		Cluster *	ac,
		ThreeVector	pos1,
		ThreeVector	pos2 )

Definition at line 52 of file G4INCLAntinucleiAtrestEntryChannel.cc.

      :theNucleus(n), theantiComposite(ac), Posnbar(pos1), Pospbar(pos2){}

AntinucleiAtrestEntryChannel() [2/2]

G4INCL::AntinucleiAtrestEntryChannel::AntinucleiAtrestEntryChannel	(	Nucleus *	n,
		Particle *	p )

Definition at line 55 of file G4INCLAntinucleiAtrestEntryChannel.cc.

:theNucleus(n),Meson(p){}

~AntinucleiAtrestEntryChannel()

G4INCL::AntinucleiAtrestEntryChannel::~AntinucleiAtrestEntryChannel ( )

virtual

Definition at line 57 of file G4INCLAntinucleiAtrestEntryChannel.cc.

{}

Member Function Documentation

bringMesonStar()

IAvatarList G4INCL::AntinucleiAtrestEntryChannel::bringMesonStar	(	ParticleList const &	pL,
		Nucleus *const	n )

Definition at line 115 of file G4INCLAntinucleiAtrestEntryChannel.cc.

                                                                                                     {
        ThreeVector ann_position = getAnnihilationPosition(Posnbar,Pospbar);
        IAvatarList theAvatarList;
        G4int cnt=1;
        if (ann_position.getX() == 999. && ann_position.getY() == 999. && ann_position.getZ() == 999.){
            INCL_DEBUG("Particle are close to each other : 2 sources of annihilation "<< '\n');
            for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p){
                if(cnt <= pbarListSize){
                    (*p)->setPosition(Pospbar);
                }
                else
                    (*p)->setPosition(Posnbar);
                theAvatarList.push_back(new ParticleEntryAvatar(0.0, n, *p, ADAR));
                cnt++;
            }
        } 
        else{ 
            for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p){
                (*p)->setPosition(ann_position);
                theAvatarList.push_back(new ParticleEntryAvatar(0.0, n, *p, ADAR));
            }
        }
        return theAvatarList;
    }

Referenced by G4INCL::StandardPropagationModel::shootCompositeAtrest().

fillFinalState()

void G4INCL::AntinucleiAtrestEntryChannel::fillFinalState ( FinalState * fs )

virtual

Implements G4INCL::IChannel.

Definition at line 140 of file G4INCLAntinucleiAtrestEntryChannel.cc.

                                                                   {
        const G4double energyBefore = Meson->getEnergy();
        fs->addEnteringParticle(Meson);
        INCL_DEBUG("Entering antiComposite annihilation product added " << '\n');
        fs->setTotalEnergyBeforeInteraction(energyBefore);
    }

getAnnihilationPosition()

ThreeVector G4INCL::AntinucleiAtrestEntryChannel::getAnnihilationPosition	(	ThreeVector	nbarPos,
		ThreeVector	pbarPos )

Definition at line 59 of file G4INCLAntinucleiAtrestEntryChannel.cc.

                                                                                                             { //Choose between the pbar or nbar annihilation position
        if((nbarPos - pbarPos).mag2() <= ParticleTable::getLargestNuclearRadius(-theantiComposite->getA(),- theantiComposite->getZ())){
        //If the annihilions positions are close (the radius of a deuteron) then we have 2 sources for the meson star
            return ThreeVector(999.,999.,999.);
        }
        if(nbarPos.mag2() <= pbarPos.mag2())
            return nbarPos;
        else
            return pbarPos;
    }

Referenced by bringMesonStar().

makeMesonStar()

ParticleList G4INCL::AntinucleiAtrestEntryChannel::makeMesonStar

(

)

Definition at line 70 of file G4INCLAntinucleiAtrestEntryChannel.cc.

                                                            {
        ParticleList Antiparticles = theantiComposite->getParticles();
        Particle *nbar=nullptr;
        Particle *pbar=nullptr;
        for(ParticleIter p =Antiparticles.begin(), e=Antiparticles.end(); p!=e; ++p){
            if((*p)->getType()==antiProton)
                pbar = *p;
            else if((*p)->getType()==antiNeutron)
                nbar = *p;
            else
                INCL_ERROR("ERROR : something else than antiNeutron or antiProton in antiComposite");
        }
        PbarAtrestEntryChannel *pbarChannel = new PbarAtrestEntryChannel(theNucleus, pbar);
        NbarAtrestEntryChannel *nbarChannel = new NbarAtrestEntryChannel(theNucleus, nbar);
        
        ParticleList TotalStarList = pbarChannel->makeMesonStar(); //pbar in first because polarisation of the dbar (Coulomb)
        ParticleList nbarMesonStar = nbarChannel->makeMesonStar();
        pbarListSize = (G4int)TotalStarList.size();
 
        for(ParticleIter p=nbarMesonStar.begin(), e=nbarMesonStar.end(); p!=e; ++p){
            TotalStarList.push_back(*p);
        }
        Pospbar = pbarChannel->getAnnihilationPosition();
        Posnbar = nbarChannel->getAnnihilationPosition();
        
        G4double EnergyofFinalMesonStar = 0;
        G4int a=theNucleus->getA();
        G4int z=theNucleus->getZ();
        G4int stra=theNucleus->getS();
        if(theNucleus->getAnnihilationType()==DNbarPPbarPType){
            EnergyofFinalMesonStar = theantiComposite->getMass() + (ParticleTable::getTableMass(a+2,z+2,stra)- ParticleTable::getTableMass(a,z+1,stra));
        }
        else if(theNucleus->getAnnihilationType()==DNbarPPbarNType || theNucleus->getAnnihilationType()==DNbarNPbarPType){
            //Correction for all but they all cancel out !  
            EnergyofFinalMesonStar = theantiComposite->getMass() + (ParticleTable::getTableMass(a+2,z+1,stra)- ParticleTable::getTableMass(a,z,stra));
 
        }
        else if(theNucleus->getAnnihilationType()==DNbarNPbarNType){
            EnergyofFinalMesonStar = theantiComposite->getMass() + (ParticleTable::getTableMass(a+2,z,stra) - ParticleTable::getTableMass(a,z-1,stra));
        }
        PhaseSpaceGenerator::generate(EnergyofFinalMesonStar, TotalStarList);
        
        return TotalStarList; 
    }

Referenced by G4INCL::StandardPropagationModel::shootCompositeAtrest().

---

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

• G4INCLAntinucleiAtrestEntryChannel.hh
• G4INCLAntinucleiAtrestEntryChannel.cc