#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 1360 of file G4LundStringFragmentation.cc.

1361{}

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: