G4VLongitudinalStringDecay Class Reference

#include <G4VLongitudinalStringDecay.hh>

Inheritance diagram for G4VLongitudinalStringDecay:

Collaboration diagram for G4VLongitudinalStringDecay:

Public Member Functions
	G4VLongitudinalStringDecay ()
virtual	~G4VLongitudinalStringDecay ()
virtual G4KineticTrackVector *	FragmentString (const G4ExcitedString &theString)=0
G4int	SampleQuarkFlavor (void)
G4ThreeVector	SampleQuarkPt (G4double ptMax=-1.)
G4KineticTrackVector *	DecayResonans (G4KineticTrackVector *aHadrons)
void	SetSigmaTransverseMomentum (G4double aQT)
void	SetStrangenessSuppression (G4double aValue)
void	SetDiquarkSuppression (G4double aValue)
void	SetDiquarkBreakProbability (G4double aValue)
void	SetVectorMesonProbability (G4double aValue)
void	SetSpinThreeHalfBarionProbability (G4double aValue)
void	SetScalarMesonMixings (std::vector< G4double > aVector)
void	SetVectorMesonMixings (std::vector< G4double > aVector)
void	SetStringTensionParameter (G4double aValue)

Protected Types
typedef std::pair < G4ParticleDefinition *, G4ParticleDefinition * >	pDefPair
typedef G4ParticleDefinition *(G4HadronBuilder::*	Pcreate )(G4ParticleDefinition *, G4ParticleDefinition *)

Protected Member Functions
virtual void	SetMassCut (G4double aValue)
G4KineticTrackVector *	LightFragmentationTest (const G4ExcitedString *const theString)
G4double	FragmentationMass (const G4FragmentingString *const string, Pcreate build=0, pDefPair *pdefs=0)
G4ParticleDefinition *	FindParticle (G4int Encoding)
virtual void	Sample4Momentum (G4LorentzVector *Mom, G4double Mass, G4LorentzVector *AntiMom, G4double AntiMass, G4double InitialMass)=0
virtual G4bool	StopFragmenting (const G4FragmentingString *const string)=0
virtual G4bool	IsFragmentable (const G4FragmentingString *const string)=0
virtual G4bool	SplitLast (G4FragmentingString *string, G4KineticTrackVector *LeftVector, G4KineticTrackVector *RightVector)=0
G4ExcitedString *	CPExcited (const G4ExcitedString &string)
G4KineticTrack *	Splitup (G4FragmentingString *string, G4FragmentingString *&newString)
G4ParticleDefinition *	QuarkSplitup (G4ParticleDefinition *decay, G4ParticleDefinition *&created)
virtual G4ParticleDefinition *	DiQuarkSplitup (G4ParticleDefinition *decay, G4ParticleDefinition *&created)=0
pDefPair	CreatePartonPair (G4int NeedParticle, G4bool AllowDiquarks=true)
virtual G4LorentzVector *	SplitEandP (G4ParticleDefinition *pHadron, G4FragmentingString *string, G4FragmentingString *newString)=0
virtual G4double	GetLightConeZ (G4double zmin, G4double zmax, G4int PartonEncoding, G4ParticleDefinition *pHadron, G4double Px, G4double Py)=0
void	CalculateHadronTimePosition (G4double theInitialStringMass, G4KineticTrackVector *)
void	ConstructParticle ()
G4ParticleDefinition *	CreateHadron (G4int id1, G4int id2, G4bool theGivenSpin, G4int theSpin)
G4double	GetDiquarkSuppress ()
G4double	GetDiquarkBreakProb ()
G4double	GetStrangeSuppress ()
G4double	GetClusterMass ()
G4int	GetClusterLoopInterrupt ()
G4double	GetStringTensionParameter ()

Protected Attributes
G4double	MassCut
G4double	ClusterMass
G4double	SigmaQT
G4double	DiquarkSuppress
G4double	DiquarkBreakProb
G4double	SmoothParam
G4double	StrangeSuppress
G4int	StringLoopInterrupt
G4int	ClusterLoopInterrupt
G4HadronBuilder *	hadronizer
G4double	pspin_meson
G4double	pspin_barion
std::vector< G4double >	vectorMesonMix
std::vector< G4double >	scalarMesonMix
G4bool	PastInitPhase
G4double	Kappa

Detailed Description

Definition at line 46 of file G4VLongitudinalStringDecay.hh.

Member Typedef Documentation

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

protected

Definition at line 75 of file G4VLongitudinalStringDecay.hh.

typedef std::pair<G4ParticleDefinition*, G4ParticleDefinition*> G4VLongitudinalStringDecay::pDefPair

protected

Definition at line 71 of file G4VLongitudinalStringDecay.hh.

Constructor & Destructor Documentation

G4VLongitudinalStringDecay::G4VLongitudinalStringDecay

(

)

Definition at line 65 of file G4VLongitudinalStringDecay.cc.

{
  MassCut  = 0.35*GeV; 
  ClusterMass = 0.15*GeV;

  SmoothParam      = 0.9; 
  StringLoopInterrupt    = 1000;
  ClusterLoopInterrupt   =  500;

  // Changable Parameters below.
  SigmaQT = 0.5 * GeV;  // 0.5 0.1
   
  StrangeSuppress  = 0.44;    //  27 % strange quarks produced, ie. u:d:s=1:1:0.27
  DiquarkSuppress  = 0.07;
  DiquarkBreakProb = 0.1;
   
  //... pspin_meson is probability to create pseudo-scalar meson 
  pspin_meson = 0.5;

  //... pspin_barion is probability to create 1/2 barion 
  pspin_barion = 0.5;

  //... vectorMesonMix[] is quark mixing parameters for vector mesons (Variable spin = 3)
  vectorMesonMix.resize(6);
  vectorMesonMix[0] = 0.0;  //AR-20Oct2014 : it was 0.5
  vectorMesonMix[1] = 0.0;
  vectorMesonMix[2] = 0.0;  //AR-20Oct2014 : it was 0.5
  vectorMesonMix[3] = 0.0;
  vectorMesonMix[4] = 1.0;
  vectorMesonMix[5] = 1.0; 

  //... scalarMesonMix[] is quark mixing parameters for scalar mesons (Variable spin=1)
  scalarMesonMix.resize(6);
  scalarMesonMix[0] = 0.5; 
  scalarMesonMix[1] = 0.25; 
  scalarMesonMix[2] = 0.5; 
  scalarMesonMix[3] = 0.25; 
  scalarMesonMix[4] = 1.0; 
  scalarMesonMix[5] = 0.5; 

  // Parameters may be changed until the first fragmentation starts
  PastInitPhase=false;
  hadronizer = new G4HadronBuilder(pspin_meson,pspin_barion,scalarMesonMix,vectorMesonMix);
  Kappa = 1.0 * GeV/fermi;
}

G4VLongitudinalStringDecay::~G4VLongitudinalStringDecay ( )

virtual

Definition at line 112 of file G4VLongitudinalStringDecay.cc.

{
  delete hadronizer;
}

Member Function Documentation

void G4VLongitudinalStringDecay::CalculateHadronTimePosition ( G4double  theInitialStringMass, G4KineticTrackVector *  Hadrons  )

protected

Definition at line 487 of file G4VLongitudinalStringDecay.cc.

{
  // 'yo-yo' formation time
  //const G4double kappa = 1.0 * GeV/fermi/4.;      
  G4double kappa = GetStringTensionParameter();
  for(size_t c1 = 0; c1 < Hadrons->size(); c1++)
  {
    G4double SumPz = 0; 
    G4double SumE  = 0;
    for(size_t c2 = 0; c2 < c1; c2++)
    {
      SumPz += Hadrons->operator[](c2)->Get4Momentum().pz();
      SumE  += Hadrons->operator[](c2)->Get4Momentum().e();   
    } 
    G4double HadronE  = Hadrons->operator[](c1)->Get4Momentum().e();
    G4double HadronPz = Hadrons->operator[](c1)->Get4Momentum().pz();
    Hadrons->operator[](c1)->SetFormationTime( 
        (theInitialStringMass - 2.*SumPz + HadronE - HadronPz)/(2.*kappa)/c_light ); 

    G4ThreeVector aPosition(0, 0, (theInitialStringMass - 2.*SumE  - HadronE + HadronPz)/(2.*kappa));
    Hadrons->operator[](c1)->SetPosition(aPosition);

  } 
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::ConstructParticle ( )

protected

G4ExcitedString * G4VLongitudinalStringDecay::CPExcited ( const G4ExcitedString &  string )

protected

Definition at line 289 of file G4VLongitudinalStringDecay.cc.

{
  G4Parton *Left=new G4Parton(*in.GetLeftParton());
  G4Parton *Right=new G4Parton(*in.GetRightParton());
  return new G4ExcitedString(Left,Right,in.GetDirection());
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

protected

G4VLongitudinalStringDecay::pDefPair G4VLongitudinalStringDecay::CreatePartonPair ( G4int  NeedParticle, G4bool  AllowDiquarks = true  )

protected

Definition at line 446 of file G4VLongitudinalStringDecay.cc.

{
  //  NeedParticle = +1 for Particle, -1 for Antiparticle

  if ( AllowDiquarks && G4UniformRand() < DiquarkSuppress )
  {
    // Create a Diquark - AntiDiquark pair , first in pair is anti to IsParticle
    G4int q1  = SampleQuarkFlavor();
    G4int q2  = SampleQuarkFlavor();
    G4int spin = (q1 != q2 && G4UniformRand() <= 0.5)? 1 : 3;
    // Convention: quark with higher PDG number is first
    G4int PDGcode = (std::max(q1,q2) * 1000 + std::min(q1,q2) * 100 + spin) * NeedParticle;
    return pDefPair (FindParticle(-PDGcode),FindParticle(PDGcode));
  } else {
    // Create a Quark - AntiQuark pair, first in pair  IsParticle
    G4int PDGcode=SampleQuarkFlavor()*NeedParticle;
    return pDefPair (FindParticle(PDGcode),FindParticle(-PDGcode));
  }

}

Here is the call graph for this function:

Here is the caller graph for this function:

G4KineticTrackVector* G4VLongitudinalStringDecay::DecayResonans

(

G4KineticTrackVector *

aHadrons

)

virtual G4ParticleDefinition* G4VLongitudinalStringDecay::DiQuarkSplitup ( G4ParticleDefinition *  decay, G4ParticleDefinition *&  created  )

protectedpure virtual

Here is the caller graph for this function:

G4ParticleDefinition * G4VLongitudinalStringDecay::FindParticle ( G4int  Encoding )

protected

Definition at line 262 of file G4VLongitudinalStringDecay.cc.

{
  G4ParticleDefinition* ptr = G4ParticleTable::GetParticleTable()->FindParticle(Encoding);
  if (ptr == NULL)
  {
    G4cout << "Particle with encoding "<<Encoding<<" does not exist!!!"<<G4endl;
    throw G4HadronicException(__FILE__, __LINE__, "Check your particle table");
  }
  return ptr;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4VLongitudinalStringDecay::FragmentationMass ( const G4FragmentingString *const  string, Pcreate  build = 0, pDefPair *  pdefs = 0  )

protected

Definition at line 204 of file G4VLongitudinalStringDecay.cc.

{
  G4double mass;
  static G4ThreadLocal G4bool NeedInit(true);
  static G4ThreadLocal std::vector<double> *nomix_G4MT_TLS_ = 0 ; 
  if (!nomix_G4MT_TLS_) nomix_G4MT_TLS_ = new  std::vector<double>;
  std::vector<double> &nomix = *nomix_G4MT_TLS_;
  static G4ThreadLocal G4HadronBuilder * minMassHadronizer;
  if ( NeedInit ) 
  {
    NeedInit = false;
    nomix.resize(6);
    for ( G4int i=0; i<6 ; i++ ) nomix[i]=0;

    //minMassHadronizer=new G4HadronBuilder(pspin_meson,pspin_barion,nomix,nomix);
    minMassHadronizer=hadronizer;
  }

  if ( build==0 ) build=&G4HadronBuilder::BuildLowSpin;

  G4ParticleDefinition *Hadron1, *Hadron2=0;

  if (!string->FourQuarkString() )
  {
    // spin 0 meson or spin 1/2 barion will be built

    Hadron1 = (minMassHadronizer->*build)(string->GetLeftParton(), string->GetRightParton());

    #ifdef debug_VStringDecay
    G4cout<<"Quarks at the string ends "<<string->GetLeftParton()->GetParticleName()
          <<" "<<string->GetRightParton()->GetParticleName()<<G4endl;
    G4cout<<"(G4VLongitudinalStringDecay) Hadron "<<Hadron1->GetParticleName()<<" "<<Hadron1->GetPDGMass()<<G4endl;
    #endif

    mass= (Hadron1)->GetPDGMass();
  } else {
    //... string is qq--qqbar: Build two stable hadrons,
    //...    with extra uubar or ddbar quark pair
    G4int iflc = (G4UniformRand() < 0.5)? 1 : 2;
    if (string->GetLeftParton()->GetPDGEncoding() < 0) iflc = -iflc;

    //... theSpin = 4; spin 3/2 baryons will be built
    Hadron1 = (minMassHadronizer->*build)(string->GetLeftParton(), FindParticle(iflc));
    Hadron2 = (minMassHadronizer->*build)(string->GetRightParton(), FindParticle(-iflc));
    mass = (Hadron1)->GetPDGMass() + (Hadron2)->GetPDGMass();
  }
    
  if ( pdefs != 0 ) 
  { // need to return hadrons as well....
    pdefs->first  = Hadron1;
    pdefs->second = Hadron2;
  }
       
  return mass;
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4KineticTrackVector* G4VLongitudinalStringDecay::FragmentString ( const G4ExcitedString &  theString )

pure virtual

Implemented in G4LundStringFragmentation, and G4QGSMFragmentation.

Here is the caller graph for this function:

G4int G4VLongitudinalStringDecay::GetClusterLoopInterrupt ( )

inlineprotected

Definition at line 164 of file G4VLongitudinalStringDecay.hh.

{return ClusterLoopInterrupt;};

G4double G4VLongitudinalStringDecay::GetClusterMass ( )

inlineprotected

Definition at line 163 of file G4VLongitudinalStringDecay.hh.

{return ClusterMass;};

G4double G4VLongitudinalStringDecay::GetDiquarkBreakProb ( )

inlineprotected

Definition at line 161 of file G4VLongitudinalStringDecay.hh.

{return DiquarkBreakProb;};

G4double G4VLongitudinalStringDecay::GetDiquarkSuppress ( )

inlineprotected

Definition at line 160 of file G4VLongitudinalStringDecay.hh.

{return DiquarkSuppress;};

virtual G4double G4VLongitudinalStringDecay::GetLightConeZ ( G4double  zmin, G4double  zmax, G4int  PartonEncoding, G4ParticleDefinition *  pHadron, G4double  Px, G4double  Py  )

protectedpure virtual

G4double G4VLongitudinalStringDecay::GetStrangeSuppress ( )

inlineprotected

Definition at line 162 of file G4VLongitudinalStringDecay.hh.

{return StrangeSuppress;};

G4double G4VLongitudinalStringDecay::GetStringTensionParameter ( )

inlineprotected

Definition at line 166 of file G4VLongitudinalStringDecay.hh.

{return Kappa;};

Here is the caller graph for this function:

virtual G4bool G4VLongitudinalStringDecay::IsFragmentable ( const G4FragmentingString *const  string )

protectedpure virtual

G4KineticTrackVector * G4VLongitudinalStringDecay::LightFragmentationTest ( const G4ExcitedString *const  theString )

protected

Definition at line 147 of file G4VLongitudinalStringDecay.cc.

{
  // Check string decay threshold   
  G4KineticTrackVector * result=0;  // return 0 when string exceeds the mass cut
    
  pDefPair hadrons((G4ParticleDefinition *)0,(G4ParticleDefinition *)0);

  G4FragmentingString aString(*string);
  if ( sqr(FragmentationMass(&aString,0,&hadrons)+MassCut) < aString.Mass2()) {
    return 0;
  }

  // The string mass is very low ---------------------------
    
  result=new G4KineticTrackVector;
        
  if ( hadrons.second ==0 ) {
    // Substitute string by light hadron, Note that Energy is not conserved here!

    #ifdef debug_VStringDecay
    G4cout << "VlongSF Warning replacing string by single hadron (G4VLongitudinalStringDecay)" <<G4endl;
    G4cout << hadrons.first->GetParticleName()<<G4endl
       << "string .. " << string->Get4Momentum() << " " 
       << string->Get4Momentum().m() << G4endl;
    #endif      

    G4ThreeVector   Mom3 = string->Get4Momentum().vect();
    G4LorentzVector Mom(Mom3, std::sqrt(Mom3.mag2() + sqr(hadrons.first->GetPDGMass())));
    result->push_back(new G4KineticTrack(hadrons.first, 0, string->GetPosition(), Mom));
  } else {
    //... string was qq--qqbar type: Build two stable hadrons,

    #ifdef debug_VStringDecay
    G4cout << "VlongSF Warning replacing qq-qqbar string by TWO hadrons (G4VLongitudinalStringDecay)" 
       << hadrons.first->GetParticleName() << " / " 
       << hadrons.second->GetParticleName()
       << "string .. " << string->Get4Momentum() << " " 
       << string->Get4Momentum().m() << G4endl;
    #endif            

    G4LorentzVector  Mom1, Mom2;
    Sample4Momentum(&Mom1, hadrons.first->GetPDGMass(), &Mom2,hadrons.second->GetPDGMass(),
                string->Get4Momentum().mag());

    result->push_back(new G4KineticTrack(hadrons.first, 0, string->GetPosition(), Mom1));
    result->push_back(new G4KineticTrack(hadrons.second, 0, string->GetPosition(), Mom2));

    G4ThreeVector Velocity = string->Get4Momentum().boostVector();
    result->Boost(Velocity);          
  }

  return result;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4ParticleDefinition * G4VLongitudinalStringDecay::QuarkSplitup ( G4ParticleDefinition *  decay, G4ParticleDefinition *&  created  )

protected

Definition at line 376 of file G4VLongitudinalStringDecay.cc.

{
  // if we have a quark, we need antiquark (or diquark)
  G4int IsParticle=(decay->GetPDGEncoding()>0) ? -1 : +1;  
                                                            
  pDefPair QuarkPair = CreatePartonPair(IsParticle);
  created = QuarkPair.second;
  return hadronizer->Build(QuarkPair.first, decay);
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual void G4VLongitudinalStringDecay::Sample4Momentum ( G4LorentzVector *  Mom, G4double  Mass, G4LorentzVector *  AntiMom, G4double  AntiMass, G4double  InitialMass  )

protectedpure virtual

Here is the caller graph for this function:

G4int G4VLongitudinalStringDecay::SampleQuarkFlavor

(

void

)

Definition at line 438 of file G4VLongitudinalStringDecay.cc.

{
  return (1 + (int)(G4UniformRand()/StrangeSuppress));
}

Here is the caller graph for this function:

G4ThreeVector G4VLongitudinalStringDecay::SampleQuarkPt

(

G4double

ptMax = -1.

)

Definition at line 469 of file G4VLongitudinalStringDecay.cc.

{
  G4double Pt;
  if ( ptMax < 0 ) {
    // sample full gaussian
    Pt = -G4Log(G4UniformRand());
  } else {
    // sample in limited range
    Pt = -G4Log(G4RandFlat::shoot(G4Exp(-sqr(ptMax)/sqr(SigmaQT)), 1.));
  }
  Pt = SigmaQT * std::sqrt(Pt);
  G4double phi = 2.*pi*G4UniformRand();
  return G4ThreeVector(Pt * std::cos(phi),Pt * std::sin(phi),0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetDiquarkBreakProbability

(

G4double

aValue

)

Definition at line 550 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "4VLongitudinalStringDecay::SetDiquarkBreakProbability after FragmentString() not allowed");
  } else {
    DiquarkBreakProb = aValue;
  }
}

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetDiquarkSuppression

(

G4double

aValue

)

Definition at line 538 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "4VLongitudinalStringDecay::SetDiquarkSuppression after FragmentString() not allowed");
  } else {
    DiquarkSuppress = aValue;
  }
}

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetMassCut ( G4double  aValue )

protectedvirtual

Definition at line 140 of file G4VLongitudinalStringDecay.cc.

{MassCut=aValue;}

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetScalarMesonMixings

(

std::vector< G4double >

aVector

)

Definition at line 590 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "G4VLongitudinalStringDecay::SetScalarMesonMixings after FragmentString() not allowed");
  } else {
    if ( aVector.size() < 6 )  
      throw G4HadronicException(__FILE__, __LINE__, 
                                "G4VLongitudinalStringDecay::SetScalarMesonMixings( argument Vector too small");
    scalarMesonMix[0] = aVector[0];
    scalarMesonMix[1] = aVector[1];
    scalarMesonMix[2] = aVector[2];
    scalarMesonMix[3] = aVector[3];
    scalarMesonMix[4] = aVector[4];
    scalarMesonMix[5] = aVector[5];
    delete hadronizer;
    hadronizer = new G4HadronBuilder(pspin_meson,pspin_barion,scalarMesonMix,vectorMesonMix);
  }
}

void G4VLongitudinalStringDecay::SetSigmaTransverseMomentum

(

G4double

aQT

)

Definition at line 514 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "4VLongitudinalStringDecay::SetSigmaTransverseMomentum after FragmentString() not allowed");
  } else {
    SigmaQT = aValue;
  }
}

void G4VLongitudinalStringDecay::SetSpinThreeHalfBarionProbability

(

G4double

aValue

)

Definition at line 576 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "G4VLongitudinalStringDecay::SetSpinThreeHalfBarionProbability after FragmentString() not allowed");
  } else {
    pspin_barion = aValue;
    delete hadronizer;
    hadronizer = new G4HadronBuilder(pspin_meson,pspin_barion,scalarMesonMix,vectorMesonMix);
  }
}

void G4VLongitudinalStringDecay::SetStrangenessSuppression

(

G4double

aValue

)

Definition at line 526 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "4VLongitudinalStringDecay::SetStrangenessSuppression after FragmentString() not allowed");
  } else {
    StrangeSuppress = aValue;
  }
}

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetStringTensionParameter

(

G4double

aValue

)

Definition at line 633 of file G4VLongitudinalStringDecay.cc.

{
  Kappa = aValue * GeV/fermi;
}   

Here is the caller graph for this function:

void G4VLongitudinalStringDecay::SetVectorMesonMixings

(

std::vector< G4double >

aVector

)

Definition at line 612 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "G4VLongitudinalStringDecay::SetVectorMesonMixings after FragmentString() not allowed");
  } else {
    if ( aVector.size() < 6 ) 
      throw G4HadronicException(__FILE__, __LINE__, 
                                "G4VLongitudinalStringDecay::SetVectorMesonMixings( argument Vector too small");
    vectorMesonMix[0] = aVector[0];
    vectorMesonMix[1] = aVector[1];
    vectorMesonMix[2] = aVector[2];
    vectorMesonMix[3] = aVector[3];
    vectorMesonMix[4] = aVector[4];
    vectorMesonMix[5] = aVector[5];
    delete hadronizer;
    hadronizer = new G4HadronBuilder(pspin_meson,pspin_barion,scalarMesonMix,vectorMesonMix);
  }
}

void G4VLongitudinalStringDecay::SetVectorMesonProbability

(

G4double

aValue

)

Definition at line 562 of file G4VLongitudinalStringDecay.cc.

{
  if ( PastInitPhase ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                              "G4VLongitudinalStringDecay::SetVectorMesonProbability after FragmentString() not allowed");
  } else {
    pspin_meson = aValue;
    delete hadronizer;
    hadronizer = new G4HadronBuilder(pspin_meson,pspin_barion,scalarMesonMix,vectorMesonMix);
  }
}

virtual G4LorentzVector* G4VLongitudinalStringDecay::SplitEandP ( G4ParticleDefinition *  pHadron, G4FragmentingString *  string, G4FragmentingString *  newString  )

protectedpure virtual

Here is the caller graph for this function:

virtual G4bool G4VLongitudinalStringDecay::SplitLast ( G4FragmentingString *  string, G4KineticTrackVector *  LeftVector, G4KineticTrackVector *  RightVector  )

protectedpure virtual

G4KineticTrack * G4VLongitudinalStringDecay::Splitup ( G4FragmentingString *  string, G4FragmentingString *&  newString  )

protected

Definition at line 299 of file G4VLongitudinalStringDecay.cc.

{

  #ifdef debug_VStringDecay
  G4cout<<G4endl;
  G4cout<<"Start SplitUP (G4VLongitudinalStringDecay) ========================="<<G4endl;
  G4cout<<"String partons: " <<string->GetLeftParton()->GetPDGEncoding()<<" "
                             <<string->GetRightParton()->GetPDGEncoding()<<" "
        <<"Direction "       <<string->GetDecayDirection()<<G4endl;
  #endif

  //... random choice of string end to use for creating the hadron (decay)   
  G4int SideOfDecay = (G4UniformRand() < 0.5)? 1: -1;
  if (SideOfDecay < 0)
  {
    string->SetLeftPartonStable();
  } else {
    string->SetRightPartonStable();
  }

  G4ParticleDefinition *newStringEnd;
  G4ParticleDefinition * HadronDefinition;
  if (string->DecayIsQuark())
  {
    HadronDefinition= QuarkSplitup(string->GetDecayParton(), newStringEnd);
  } else {
    HadronDefinition= DiQuarkSplitup(string->GetDecayParton(), newStringEnd);
  }      

  #ifdef debug_VStringDecay
  G4cout<<"The parton "<<string->GetDecayParton()->GetPDGEncoding()<<" "
        <<" produces hadron "<<HadronDefinition->GetParticleName()
        <<" and is transformed to "<<newStringEnd->GetPDGEncoding()<<G4endl;
  G4cout<<"The side of the string decay Left/Right (1/-1) "<<SideOfDecay<<G4endl;
  #endif
  // create new String from old, ie. keep Left and Right order, but replace decay

  newString=new G4FragmentingString(*string,newStringEnd); // To store possible quark containt of new string

  #ifdef debug_VStringDecay
  G4cout<<"An attempt to determine its energy (SplitEandP)"<<G4endl;
  #endif

  G4LorentzVector* HadronMomentum=SplitEandP(HadronDefinition, string, newString);

  delete newString; newString=0;
    
  G4KineticTrack * Hadron =0;
  if ( HadronMomentum != 0 ) {

    #ifdef debug_VStringDecay                     
    G4cout<<"The attempt was successful"<<G4endl;
    #endif

    G4ThreeVector Pos;
    Hadron = new G4KineticTrack(HadronDefinition, 0,Pos, *HadronMomentum);
   
    newString=new G4FragmentingString(*string,newStringEnd,HadronMomentum);

    delete HadronMomentum;
  } else {

    #ifdef debug_VStringDecay
    G4cout<<"The attempt was not successful !!!"<<G4endl;
    #endif
  }

  #ifdef debug_VStringDecay
  G4cout<<"End SplitUP (G4VLongitudinalStringDecay) ====================="<<G4endl;
  #endif

  return Hadron;
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4bool G4VLongitudinalStringDecay::StopFragmenting ( const G4FragmentingString *const  string )

protectedpure virtual

Member Data Documentation

G4int G4VLongitudinalStringDecay::ClusterLoopInterrupt

protected

Definition at line 177 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::ClusterMass

protected

Definition at line 170 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::DiquarkBreakProb

protected

Definition at line 173 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::DiquarkSuppress

protected

Definition at line 172 of file G4VLongitudinalStringDecay.hh.

G4HadronBuilder* G4VLongitudinalStringDecay::hadronizer

protected

Definition at line 178 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::Kappa

protected

Definition at line 187 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::MassCut

protected

Definition at line 166 of file G4VLongitudinalStringDecay.hh.

G4bool G4VLongitudinalStringDecay::PastInitPhase

protected

Definition at line 185 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::pspin_barion

protected

Definition at line 181 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::pspin_meson

protected

Definition at line 180 of file G4VLongitudinalStringDecay.hh.

std::vector<G4double> G4VLongitudinalStringDecay::scalarMesonMix

protected

Definition at line 183 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::SigmaQT

protected

Definition at line 171 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::SmoothParam

protected

Definition at line 174 of file G4VLongitudinalStringDecay.hh.

G4double G4VLongitudinalStringDecay::StrangeSuppress

protected

Definition at line 175 of file G4VLongitudinalStringDecay.hh.

G4int G4VLongitudinalStringDecay::StringLoopInterrupt

protected

Definition at line 176 of file G4VLongitudinalStringDecay.hh.

std::vector<G4double> G4VLongitudinalStringDecay::vectorMesonMix

protected

Definition at line 182 of file G4VLongitudinalStringDecay.hh.

---

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

• geant4.10.03.p01/source/processes/hadronic/models/parton_string/hadronization/include/G4VLongitudinalStringDecay.hh
• geant4.10.03.p01/source/processes/hadronic/models/parton_string/hadronization/src/G4VLongitudinalStringDecay.cc