G4GammaParticipants Class Reference

#include <G4GammaParticipants.hh>

Inheritance diagram for G4GammaParticipants:

Collaboration diagram for G4GammaParticipants:

Public Member Functions
virtual	~G4GammaParticipants ()
Public Member Functions inherited from G4QGSParticipants
	G4QGSParticipants ()
	G4QGSParticipants (const G4QGSParticipants &right)
const G4QGSParticipants &	operator= (const G4QGSParticipants &right)
virtual	~G4QGSParticipants ()
int	operator== (const G4QGSParticipants &right) const
int	operator!= (const G4QGSParticipants &right) const
virtual void	DoLorentzBoost (G4ThreeVector aBoost)
G4PartonPair *	GetNextPartonPair ()
void	BuildInteractions (const G4ReactionProduct &thePrimary)
void	StartPartonPairLoop ()
Public Member Functions inherited from G4VParticipants
	G4VParticipants ()
	G4VParticipants (const G4VParticipants &right)
virtual	~G4VParticipants ()
const G4VParticipants &	operator= (const G4VParticipants &right)
int	operator== (const G4VParticipants &right) const
int	operator!= (const G4VParticipants &right) const
virtual void	Init (G4int theZ, G4int theA)
virtual void	SetNucleus (G4V3DNucleus *aNucleus)
virtual G4V3DNucleus *	GetWoundedNucleus () const
virtual void	InitProjectileNucleus (G4int theZ, G4int theA)
virtual void	SetProjectileNucleus (G4V3DNucleus *aNucleus)
virtual G4V3DNucleus *	GetProjectileNucleus () const

Private Member Functions
virtual G4VSplitableHadron *	SelectInteractions (const G4ReactionProduct &thePrimary)

Additional Inherited Members
Public Attributes inherited from G4VParticipants
G4V3DNucleus *	theNucleus
G4V3DNucleus *	theProjectileNucleus
Protected Types inherited from G4QGSParticipants
enum	{ SOFT, DIFFRACTIVE }
Protected Member Functions inherited from G4QGSParticipants
void	SplitHadrons ()
void	PerformSoftCollisions ()
void	PerformDiffractiveCollisions ()
G4bool	IsSingleDiffractive ()
Protected Attributes inherited from G4QGSParticipants
std::vector< G4InteractionContent * >	theInteractions
std::vector< G4VSplitableHadron * >	theTargets
std::vector< G4PartonPair * >	thePartonPairs
G4SingleDiffractiveExcitation	theSingleDiffExcitation
G4QGSDiffractiveExcitation	theDiffExcitaton
G4int	ModelMode
G4ThreeVector	theBoost
const G4int	nCutMax
const G4double	ThresholdParameter
const G4double	QGSMThreshold
const G4double	theNucleonRadius

Detailed Description

Definition at line 37 of file G4GammaParticipants.hh.

Constructor & Destructor Documentation

~G4GammaParticipants()

virtual G4GammaParticipants::~G4GammaParticipants ( )

inlinevirtual

Definition at line 40 of file G4GammaParticipants.hh.

{}

Here is the call graph for this function:

Member Function Documentation

SelectInteractions()

G4VSplitableHadron * G4GammaParticipants::SelectInteractions ( const G4ReactionProduct &  thePrimary )

privatevirtual

Reimplemented from G4QGSParticipants.

Definition at line 40 of file G4GammaParticipants.cc.

{
    // Check reaction threshold  - goes to CheckThreshold
    G4VSplitableHadron* aProjectile = new G4QGSMSplitableHadron(thePrimary, TRUE); // @@@ check the TRUE

    const std::vector<G4Nucleon>& theTargetNuc = theNucleus->GetNucleons();
    G4LorentzVector aPrimaryMomentum(thePrimary.GetMomentum(), thePrimary.GetTotalEnergy());
    if((!(aPrimaryMomentum.e()>-1)) && (!(aPrimaryMomentum.e()<1)) )
    {
        throw G4HadronicException(__FILE__, __LINE__,
                "G4GammaParticipants::SelectInteractions: primary nan energy.");
    }
    G4double S = (aPrimaryMomentum + theTargetNuc[0].Get4Momentum()).mag2();
    G4double ThresholdMass = thePrimary.GetMass() + theTargetNuc[0].GetDefinition()->GetPDGMass();
    ModelMode = SOFT;
    if (sqr(ThresholdMass + ThresholdParameter) > S)
    {
        ModelMode = DIFFRACTIVE;
        //throw G4HadronicException(__FILE__, __LINE__, "Initial energy is too low. The 4-vectors of the input are inconsistant with the particle masses.");
    }
    if (sqr(ThresholdMass + QGSMThreshold) > S) // thus only diffractive in cascade!
    {
        ModelMode = DIFFRACTIVE;
    }

    // first find the collisions HPW
    std::for_each(theInteractions.begin(), theInteractions.end(), DeleteInteractionContent());
    theInteractions.clear();
    G4int totalCuts = 0;

    #ifdef debug_G4GammaParticipants
        G4double eK = thePrimary.GetKineticEnergy()/GeV;
        G4int nucleonCount = theTargetNuc.size(); // debug
    #endif

    G4int theCurrent = static_cast<G4int> (theTargetNuc.size()*G4UniformRand());
    const G4Nucleon& pNucleon = theTargetNuc[theCurrent];
    G4QGSMSplitableHadron* aTarget = new G4QGSMSplitableHadron(pNucleon);
    theTargets.push_back(aTarget);
    const_cast<G4Nucleon&>(pNucleon).Hit(aTarget);
    if ( (0.06 > G4UniformRand() &&(ModelMode==SOFT)) || (ModelMode==DIFFRACTIVE ) )
    {
        // diffractive interaction occurs
        if(IsSingleDiffractive())
        {
            theSingleDiffExcitation.ExciteParticipants(aProjectile, aTarget);
        }
        else
        {
            theDiffExcitaton.ExciteParticipants(aProjectile, aTarget);
        }
        G4InteractionContent * aInteraction = new G4InteractionContent(aProjectile);
        aInteraction->SetTarget(aTarget);
        theInteractions.push_back(aInteraction);
        aInteraction->SetNumberOfDiffractiveCollisions(1);
        totalCuts += 1;
    }
    else
    {
        // nondiffractive soft interaction occurs
        aTarget->IncrementCollisionCount(1);
        aProjectile->IncrementCollisionCount(1);
        G4InteractionContent * aInteraction = new G4InteractionContent(aProjectile);
        aInteraction->SetTarget(aTarget);
        aInteraction->SetNumberOfSoftCollisions(1);
        theInteractions.push_back(aInteraction);
        totalCuts += 1;
    }
    return aProjectile;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4GammaParticipants.hh
• G4GammaParticipants.cc