#include <G4LundStringFragmentation.hh>

Inheritance diagram for G4LundStringFragmentation:

Public Member Functions
	G4LundStringFragmentation ()

virtual	~G4LundStringFragmentation ()

virtual G4KineticTrackVector *	FragmentString (const G4ExcitedString &theString)

Public Member Functions inherited from G4VLongitudinalStringDecay
	G4VLongitudinalStringDecay (const G4String &name="StringDecay")

virtual	~G4VLongitudinalStringDecay ()

G4HadFinalState *	ApplyYourself (const G4HadProjectile &, G4Nucleus &) final

void	AddNewParticles ()

void	EraseNewParticles ()

void	SetMinMasses ()

void	SetMinimalStringMass (const G4FragmentingString *const string)

void	SetMinimalStringMass2 (const G4double aValue)

G4int	SampleQuarkFlavor (void)

G4ThreeVector	SampleQuarkPt (G4double ptMax=-1.)

void	SetSigmaTransverseMomentum (G4double aQT)

void	SetStrangenessSuppression (G4double aValue)

void	SetDiquarkSuppression (G4double aValue)

void	SetDiquarkBreakProbability (G4double aValue)

void	SetSpinThreeHalfBarionProbability (G4double aValue)

void	SetScalarMesonMixings (std::vector< G4double > aVector)

void	SetVectorMesonMixings (std::vector< G4double > aVector)

void	SetStringTensionParameter (G4double aValue)

void	SetProbCCbar (G4double aValue)

void	SetProbEta_c (G4double aValue)

void	SetProbBBbar (G4double aValue)

void	SetProbEta_b (G4double aValue)

Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")

virtual	~G4HadronicInteraction ()

virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)

virtual G4bool	IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)

G4double	GetMinEnergy () const

G4double	GetMinEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMinEnergy (G4double anEnergy)

void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)

G4double	GetMaxEnergy () const

G4double	GetMaxEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMaxEnergy (const G4double anEnergy)

void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int value)

const G4String &	GetModelName () const

void	DeActivateFor (const G4Material *aMaterial)

void	ActivateFor (const G4Material *aMaterial)

void	DeActivateFor (const G4Element *anElement)

void	ActivateFor (const G4Element *anElement)

G4bool	IsBlocked (const G4Material *aMaterial) const

G4bool	IsBlocked (const G4Element *anElement) const

void	SetRecoilEnergyThreshold (G4double val)

G4double	GetRecoilEnergyThreshold () const

virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const

virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const

void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)

virtual void	ModelDescription (std::ostream &outFile) const

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	InitialiseModel ()

	G4HadronicInteraction (const G4HadronicInteraction &right)=delete

const G4HadronicInteraction &	operator= (const G4HadronicInteraction &right)=delete

G4bool	operator== (const G4HadronicInteraction &right) const =delete

G4bool	operator!= (const G4HadronicInteraction &right) const =delete

Additional Inherited Members
Public Attributes inherited from G4VLongitudinalStringDecay
G4double	Mass_of_light_quark

G4double	Mass_of_s_quark

G4double	Mass_of_c_quark

G4double	Mass_of_b_quark

G4double	Mass_of_string_junction

G4double	minMassQQbarStr [5][5]

G4double	minMassQDiQStr [5][5][5]

G4double	MinimalStringMass

G4double	MinimalStringMass2

G4int	Qcharge [5]

G4int	Meson [5][5][7]

G4double	MesonWeight [5][5][7]

G4int	Baryon [5][5][5][4]

G4double	BaryonWeight [5][5][5][4]

G4double	Prob_QQbar [5]

G4int	DecayQuark

G4int	NewQuark

G4ParticleDefinition *	FS_LeftHadron [350]

G4ParticleDefinition *	FS_RightHadron [350]

G4double	FS_Weight [350]

G4int	NumberOf_FS

Protected Types inherited from G4VLongitudinalStringDecay
typedef std::pair< G4ParticleDefinition , G4ParticleDefinition >	pDefPair

typedef G4ParticleDefinition (G4HadronBuilder::)	Pcreate(G4ParticleDefinition , G4ParticleDefinition )

Protected Member Functions inherited from G4VLongitudinalStringDecay
virtual void	SetMassCut (G4double aValue)

G4double	GetMassCut ()

G4KineticTrackVector *	ProduceOneHadron (const G4ExcitedString *const theString)

G4double	PossibleHadronMass (const G4FragmentingString const string, Pcreate build=0, pDefPair pdefs=0)

G4ParticleDefinition *	FindParticle (G4int Encoding)

G4ExcitedString *	CopyExcited (const G4ExcitedString &string)

virtual G4ParticleDefinition *	QuarkSplitup (G4ParticleDefinition decay, G4ParticleDefinition &created)

pDefPair	CreatePartonPair (G4int NeedParticle, G4bool AllowDiquarks=true)

void	CalculateHadronTimePosition (G4double theInitialStringMass, G4KineticTrackVector *)

void	ConstructParticle ()

G4ParticleDefinition *	CreateHadron (G4int id1, G4int id2, G4bool theGivenSpin, G4int theSpin)

G4double	GetDiquarkSuppress ()

G4double	GetDiquarkBreakProb ()

G4double	GetStrangeSuppress ()

G4int	GetClusterLoopInterrupt ()

G4double	GetProbCCbar ()

G4double	GetProbEta_c ()

G4double	GetProbBBbar ()

G4double	GetProbEta_b ()

G4double	GetStringTensionParameter ()

Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)

G4bool	IsBlocked () const

void	Block ()

Protected Attributes inherited from G4VLongitudinalStringDecay
G4double	MassCut

G4double	SigmaQT

G4double	DiquarkSuppress

G4double	DiquarkBreakProb

G4double	StrangeSuppress

G4int	StringLoopInterrupt

G4int	ClusterLoopInterrupt

G4HadronBuilder *	hadronizer

std::vector< G4double >	pspin_meson

G4double	pspin_barion

std::vector< G4double >	vectorMesonMix

std::vector< G4double >	scalarMesonMix

G4double	ProbCCbar

G4double	ProbEta_c

G4double	ProbBBbar

G4double	ProbEta_b

G4double	ProbCB

G4double	MaxMass

G4bool	PastInitPhase

G4double	Kappa

std::vector< G4ParticleDefinition * >	NewParticles

Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange

G4int	verboseLevel

G4double	theMinEnergy

G4double	theMaxEnergy

G4bool	isBlocked

Detailed Description

Definition at line 41 of file G4LundStringFragmentation.hh.

Constructor & Destructor Documentation

◆ G4LundStringFragmentation()

G4LundStringFragmentation::G4LundStringFragmentation ( )

Definition at line 49 of file G4LundStringFragmentation.cc.

  : G4VLongitudinalStringDecay("LundStringFragmentation")
{
    SetMassCut(210.*MeV);   //  Mpi + Delta
                            // For ProduceOneHadron it is required
                            // that no one pi-meson can be produced.
    SigmaQT = 0.435 * GeV;
    Tmt = 190.0 * MeV;       
 
    SetStringTensionParameter(1.*GeV/fermi);
    SetDiquarkBreakProbability(0.3);
 
    SetStrangenessSuppression((1.0 - 0.12)/2.0);
    SetDiquarkSuppression(0.07);
 
    // Check if charmed and bottom hadrons are enabled: if this is the case, then
    // set the non-zero probabilities for c-cbar and b-bbar creation from the vacuum,
    // else set them to 0.0. If these probabilities are/aren't zero then charmed or bottom
    // hadrons can't/can be created during the string fragmentation of ordinary
    // (i.e. not heavy) projectile hadron nuclear reactions.
    if ( G4HadronicParameters::Instance()->EnableBCParticles() ) {
      SetProbCCbar(0.0002);  // According to O.I. Piskunova Yad. Fiz. 56 (1993) 1094; tuned by Uzhi Oct. 2022
      SetProbBBbar(5.0e-5);  // According to O.I. Piskunova Yad. Fiz. 56 (1993) 1094
    } else {
      SetProbCCbar(0.0);
      SetProbBBbar(0.0);
    }
    
    SetMinMasses();  // For treating of small string decays
}

◆ ~G4LundStringFragmentation()

G4LundStringFragmentation::~G4LundStringFragmentation ( )

virtual

Definition at line 1364 of file G4LundStringFragmentation.cc.

1365{}

Member Function Documentation

◆ FragmentString()

G4KineticTrackVector * G4LundStringFragmentation::FragmentString ( const G4ExcitedString & theString )

virtual

Implements G4VLongitudinalStringDecay.

Definition at line 82 of file G4LundStringFragmentation.cc.

{
    // Can no longer modify Parameters for Fragmentation.
 
    PastInitPhase=true;
 
    G4FragmentingString  aString(theString);
    SetMinimalStringMass(&aString);
 
        #ifdef debug_LUNDfragmentation
        G4cout<<G4endl<<"LUND StringFragmentation ------------------------------------"<<G4endl;
        G4cout<<G4endl<<"LUND StringFragm: String Mass "
                      <<theString.Get4Momentum().mag()<<G4endl
                      <<"4Mom "<<theString.Get4Momentum()<<G4endl
                      <<"------------------------------------"<<G4endl;
        G4cout<<"String ends Direct "<<theString.GetLeftParton()->GetPDGcode()<<" "
                                     <<theString.GetRightParton()->GetPDGcode()<<" "
                                     <<theString.GetDirection()<< G4endl;
        G4cout<<"Left  mom "<<theString.GetLeftParton()->Get4Momentum()<<G4endl;
        G4cout<<"Right mom "<<theString.GetRightParton()->Get4Momentum()<<G4endl<<G4endl;
        G4cout<<"Check for Fragmentation "<<G4endl;
        #endif
 
    G4KineticTrackVector * LeftVector(0);
 
    if (!aString.IsAFourQuarkString() && !IsItFragmentable(&aString))
    {
                #ifdef debug_LUNDfragmentation
                G4cout<<"Non fragmentable - the string is converted to one hadron "<<G4endl;
                #endif
                // SetMassCut(210.*MeV);  // For ProduceOneHadron it is required
                                          // that no one pi-meson can be produced.
 
        G4double Mcut = GetMassCut();
        SetMassCut(10000.*MeV);
        LeftVector=ProduceOneHadron(&theString);
        SetMassCut(Mcut);
 
        if ( LeftVector )
        {
          if ( LeftVector->size() > 0)
                  {
                LeftVector->operator[](0)->SetFormationTime(theString.GetTimeOfCreation());
                LeftVector->operator[](0)->SetPosition(theString.GetPosition());
                  }
          if (LeftVector->size() > 1)
                  {
                // 2 hadrons created from qq-qqbar are stored
            LeftVector->operator[](1)->SetFormationTime(theString.GetTimeOfCreation());
            LeftVector->operator[](1)->SetPosition(theString.GetPosition());
          }
        }       
        return LeftVector;
    }
 
        #ifdef debug_LUNDfragmentation
        G4cout<<"The string will be fragmented. "<<G4endl;
        #endif
 
    // The string can fragment. At least two particles can be produced.
                   LeftVector =new G4KineticTrackVector;
    G4KineticTrackVector * RightVector=new G4KineticTrackVector;
 
        G4bool success = Loop_toFragmentString(theString, LeftVector, RightVector);
 
    if ( ! success )
    {
        std::for_each(LeftVector->begin(), LeftVector->end(), DeleteKineticTrack());
        LeftVector->clear();
        std::for_each(RightVector->begin(), RightVector->end(), DeleteKineticTrack());
        delete RightVector;
        return LeftVector;
    }
 
    // Join Left- and RightVector into LeftVector in correct order.
    while (!RightVector->empty())
    {
        LeftVector->push_back(RightVector->back());
        RightVector->erase(RightVector->end()-1);
    }
    delete RightVector;
 
    return LeftVector;
}

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