G4FTFModel Class Reference

#include <G4FTFModel.hh>

Inheritance diagram for G4FTFModel:

Public Member Functions
	G4FTFModel (const G4String &modelName="FTF")
	~G4FTFModel ()
void	Init (const G4Nucleus &aNucleus, const G4DynamicParticle &aProjectile)
G4ExcitedStringVector *	GetStrings ()
G4V3DNucleus *	GetWoundedNucleus () const
virtual void	ModelDescription (std::ostream &) const
Public Member Functions inherited from G4VPartonStringModel
	G4VPartonStringModel (const G4String &modelName="Parton String Model")
virtual	~G4VPartonStringModel ()
void	SetFragmentationModel (G4VStringFragmentation *aModel)
G4KineticTrackVector *	Scatter (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)
virtual G4V3DNucleus *	GetWoundedNucleus () const =0
virtual void	ModelDescription (std::ostream &outFile) const
Public Member Functions inherited from G4VHighEnergyGenerator
	G4VHighEnergyGenerator (const G4String &modelName="High Energy Generator")
virtual	~G4VHighEnergyGenerator ()
virtual G4V3DNucleus *	GetWoundedNucleus () const =0
virtual G4KineticTrackVector *	Scatter (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)=0
std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double AbsoluteLevel)
virtual void	ModelDescription (std::ostream &) const
virtual G4String	GetModelName () const

Additional Inherited Members
Protected Member Functions inherited from G4VPartonStringModel
virtual void	Init (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)=0
virtual G4ExcitedStringVector *	GetStrings ()=0
void	SetThisPointer (G4VPartonStringModel *aPointer)
G4bool	EnergyAndMomentumCorrector (G4KineticTrackVector *Output, G4LorentzVector &TotalCollisionMomentum)

Detailed Description

Definition at line 63 of file G4FTFModel.hh.

Constructor & Destructor Documentation

G4FTFModel()

G4FTFModel::G4FTFModel

(

const G4String &

modelName = "FTF"

)

Definition at line 56 of file G4FTFModel.cc.

                                               :G4VPartonStringModel(modelName),
                         theExcitation(new G4DiffractiveExcitation()),
                         theElastic(new G4ElasticHNScattering()),
                         theAnnihilation(new G4FTFAnnihilation())
{
    G4VPartonStringModel::SetThisPointer(this);
        theParameters=0;
    NumberOfInvolvedNucleon=0;
        NumberOfInvolvedNucleonOfProjectile=0;
    SetEnergyMomentumCheckLevels(2*perCent, 150*MeV);
}

~G4FTFModel()

G4FTFModel::~G4FTFModel

(

)

Definition at line 70 of file G4FTFModel.cc.

{
// Because FTF model can be called for various particles
// theParameters must be erased at the end of each call.
// Thus the delete is also in G4FTFModel::GetStrings() method
   if( theParameters   != 0 ) delete theParameters; 
   if( theExcitation   != 0 ) delete theExcitation;
   if( theElastic      != 0 ) delete theElastic; 
   if( theAnnihilation != 0 ) delete theAnnihilation;
 
// Erasing of strings created at annihilation
   if(theAdditionalString.size() != 0)
   {
    std::for_each(theAdditionalString.begin(), theAdditionalString.end(), 
                  DeleteVSplitableHadron());
   }
   theAdditionalString.clear();
 
// Erasing of target involved nucleons
   if( NumberOfInvolvedNucleon != 0)
   {
    for(G4int i=0; i < NumberOfInvolvedNucleon; i++)
    {
     G4VSplitableHadron * aNucleon = TheInvolvedNucleon[i]->GetSplitableHadron();
     if(aNucleon) delete aNucleon;
    }
   }
 
// Erasing of projectile involved nucleons
/*
   if( NumberOfInvolvedNucleonOfProjectile != 0)
   {
    for(G4int i=0; i < NumberOfInvolvedNucleonOfProjectile; i++)
    {
     G4VSplitableHadron * aNucleon = TheInvolvedNucleonOfProjectile[i]->GetSplitableHadron();
     if(aNucleon) delete aNucleon;
    }
   }
*/
}

Member Function Documentation

GetStrings()

G4ExcitedStringVector * G4FTFModel::GetStrings ( )

virtual

Implements G4VPartonStringModel.

Definition at line 206 of file G4FTFModel.cc.

{ 
        G4ExcitedStringVector * theStrings(0);
 
    theParticipants.GetList(theProjectile,theParameters);
//        StoreInvolvedNucleon();
//G4cout<<" GetList theProjectile.GetMomentum()  GetBaryonNumber() "<<theProjectile.GetMomentum()<<" "<<theProjectile.GetDefinition()->GetBaryonNumber()<<G4endl;
        G4bool Success(true);
 
        if((std::abs(theProjectile.GetDefinition()->GetBaryonNumber()) <= 1) &&
                    (theProjectile.GetDefinition()->GetBaryonNumber() != -1)   )
        { // Standard variant of FTF for projectile hadron/nucleon
//G4cout<<"Standard variant of FTF for projectile hadron/nucleon"<<G4endl;
         ReggeonCascade(); 
//G4cout<<"Success after Reggeon "<<Success<<" PutOnMasShell"<<G4endl;
         Success=PutOnMassShell(); 
//G4cout<<"Success after PutOn "<<Success<<" GetResid"<<G4endl;
         GetResidualNucleus();
        } 
//G4cout<<"Success after GetN "<<Success<<G4endl;
//G4int Uzhi; G4cin>>Uzhi;
        if(theProjectile.GetDefinition()->GetBaryonNumber() > 1)
        { // Variant of FTF for projectile nuclei
//G4cout<<"Variant of FTF for projectile nuclei"<<G4endl;
         StoreInvolvedNucleon();
         ReggeonCascade(); 
         Success=PutOnMassShell(); 
         GetResidualNucleus();
        } 
 
//        G4bool LowE_Anti_Ion(false);
        if(theProjectile.GetDefinition()->GetBaryonNumber() <= -1) 
        { // Projectile is Anti-baryon or Anti-Nucleus
//G4cout<<"Projectile is Anti-baryon or Anti-Nucleus "<<G4endl;
//G4cout<<"Be4 Store"<<G4endl;
         StoreInvolvedNucleon();
         if(theProjectile.GetTotalMomentum()/
            std::abs(theProjectile.GetDefinition()->GetBaryonNumber()) > 5000.*MeV)
         {// High energy interaction
//G4cout<<"High energy interaction "<<G4endl;
//G4cout<<"Regeon "<<G4endl;
          ReggeonCascade(); 
//G4cout<<"Put on mass "<<G4endl;
          Success=PutOnMassShell(); 
//G4cout<<"Residual "<<G4endl;
          GetResidualNucleus();
         }
         else
         {
//G4cout<<"Low energy interaction "<<G4endl;
//          LowE_Anti_Ion=true;
          Success=true;
         }
        }
//G4cout<<"Before Excite Success "<<Success<<G4endl;
        Success=Success && ExciteParticipants();
//G4cout<<"Success ExciteParticipants()? "<<Success<<G4endl;
//        if(LowE_Anti_Ion) Success=Success && GetResidualNucleusAfterAnnihilation();
 
        if( Success )
        {       
      theStrings = BuildStrings();
//G4cout<<"BuildStrings OK"<<G4endl;
          if( theParameters != 0 )
          {
           delete theParameters;
           theParameters=0;
          }
         }
/*
        if( Success )
        {
         if( ExciteParticipants() )
         {
//G4cout<<"Excite partic OK"<<G4endl;
      theStrings = BuildStrings();
//G4cout<<"Build String OK"<<G4endl;
          if(LowE_Anti_Ion) GetResidualNucleusAfterAnnihilation();
 
          if( theParameters != 0 )
          {
           delete theParameters;
           theParameters=0;
          }
         } else                      // if( ExciteParticipants() )
         {     Success=false;}
        } else                       // if( Success )
        {      Success=false;}
*/
        if(!Success)   
        {
           // -------------- Erase the projectile ----------------
//G4cout<<"Erase Proj"<<G4endl;
     std::vector<G4VSplitableHadron *> primaries;
 
     theParticipants.StartLoop();    // restart a loop 
         while ( theParticipants.Next() ) 
     {
        const G4InteractionContent & interaction=theParticipants.GetInteraction();
 
                     //  do not allow for duplicates ...
        if ( primaries.end() == std::find(primaries.begin(), primaries.end(),
                                                   interaction.GetProjectile()) )
            primaries.push_back(interaction.GetProjectile());                
         }
         std::for_each(primaries.begin(), primaries.end(), DeleteVSplitableHadron());
         primaries.clear();
        }
// -------------- Cleaning of the memory --------------
        G4VSplitableHadron * aNucleon = 0;
// -------------- Erase the projectile nucleon --------
/*
        G4VSplitableHadron * aNucleon = 0;
        for(G4int i=0; i < NumberOfInvolvedNucleonOfProjectile; i++)
        {
           aNucleon = TheInvolvedNucleonOfProjectile[i]->GetSplitableHadron();
           if(aNucleon) delete aNucleon;
        } 
 
        NumberOfInvolvedNucleonOfProjectile=0;
*/  // Maybe it will be needed latter------------------
 
// -------------- Erase the target nucleons -----------
//G4cout<<"Erase Target Ninv "<<NumberOfInvolvedNucleon<<G4endl;
        for(G4int i=0; i < NumberOfInvolvedNucleon; i++)
        {
           aNucleon = TheInvolvedNucleon[i]->GetSplitableHadron();
           if(aNucleon) delete aNucleon;
        } 
 
        NumberOfInvolvedNucleon=0;
//G4cout<<"Go to fragmentation"<<G4endl;
        return theStrings;
 
}

GetWoundedNucleus()

G4V3DNucleus * G4FTFModel::GetWoundedNucleus ( ) const

inlinevirtual

Implements G4VPartonStringModel.

Definition at line 121 of file G4FTFModel.hh.

{
    return theParticipants.GetWoundedNucleus();
}

Init()

void G4FTFModel::Init ( const G4Nucleus &  aNucleus, const G4DynamicParticle &  aProjectile  )

virtual

Implements G4VPartonStringModel.

Definition at line 112 of file G4FTFModel.cc.

{
    theProjectile = aProjectile;  
 
        G4double PlabPerParticle(0.);  // Laboratory momentum Pz per particle/nucleon
 
/*
G4cout<<"FTF init Pro Name "<<theProjectile.GetDefinition()->GetParticleName()<<G4endl;
G4cout<<"FTF init Pro Mass "<<theProjectile.GetMass()<<" "<<theProjectile.GetMomentum()<<G4endl;
G4cout<<"FTF init Pro B Q  "<<theProjectile.GetDefinition()->GetBaryonNumber()<<" "<<(G4int) theProjectile.GetDefinition()->GetPDGCharge()<<G4endl; 
G4cout<<"FTF init A Z "<<aNucleus.GetA_asInt()<<" "<<aNucleus.GetZ_asInt()<<G4endl;
G4cout<<"             "<<aNucleus.GetN()<<" "<<aNucleus.GetZ()<<G4endl;
//G4int Uzhi; G4cin>>Uzhi;
*/
 
        theParticipants.SetProjectileNucleus(0);
        theParticipants.Init(aNucleus.GetA_asInt(),aNucleus.GetZ_asInt());
 
        if(std::abs(theProjectile.GetDefinition()->GetBaryonNumber()) <= 1) 
        { // Projectile is a hadron
         PlabPerParticle=theProjectile.GetMomentum().z();
 
//         S = sqr( theProjectile.GetMass() ) + sqr( ProtonMass ) +
//                 2*ProtonMass*theProjectile.GetTotalEnergy();
        }
 
 
        if(theProjectile.GetDefinition()->GetBaryonNumber() > 1) 
        { // Projectile is a nucleus
         theParticipants.InitProjectileNucleus(
                      theProjectile.GetDefinition()->GetBaryonNumber(),
              (G4int) theProjectile.GetDefinition()->GetPDGCharge()    );
 
         G4ThreeVector BoostVector=theProjectile.GetMomentum()/theProjectile.GetTotalEnergy();
         theParticipants.theProjectileNucleus->DoLorentzBoost(BoostVector);
 
         PlabPerParticle=theProjectile.GetMomentum().z()/
                         theProjectile.GetDefinition()->GetBaryonNumber();
 
//         S =         2.*sqr( ProtonMass ) + 2*ProtonMass*
//             theProjectile.GetTotalEnergy()/theProjectile.GetDefinition()->GetBaryonNumber();
        }
 
        if(theProjectile.GetDefinition()->GetBaryonNumber() < -1) 
        { // Projectile is an anti-nucleus
         theParticipants.InitProjectileNucleus(
             std::abs(        theProjectile.GetDefinition()->GetBaryonNumber()),
             std::abs((G4int) theProjectile.GetDefinition()->GetPDGCharge())    );
 
         G4ThreeVector BoostVector=theProjectile.GetMomentum()/theProjectile.GetTotalEnergy();
 
         theParticipants.theProjectileNucleus->StartLoop();
         G4Nucleon * aNucleon;
         while ( (aNucleon = theParticipants.theProjectileNucleus->GetNextNucleon()) )
         {
          if(aNucleon->GetDefinition()->GetPDGEncoding() == 2212) // Proton
          {aNucleon->SetParticleType(G4AntiProton::AntiProton());} 
 
          if(aNucleon->GetDefinition()->GetPDGEncoding() == 2112) // Neutron
          {aNucleon->SetParticleType(G4AntiNeutron::AntiNeutron());} 
         }   // end of while (theParticipant.theProjectileNucleus->GetNextNucleon())
 
         theParticipants.theProjectileNucleus->DoLorentzBoost(BoostVector);
 
         PlabPerParticle=         theProjectile.GetMomentum().z()/
                         std::abs(theProjectile.GetDefinition()->GetBaryonNumber());
 
//         S =        2.*sqr( ProtonMass ) + 2*ProtonMass*
//                      theProjectile.GetTotalEnergy()/
//             std::abs(theProjectile.GetDefinition()->GetBaryonNumber());
        }
 
// ------------------------------------------------------------------------
      if( theParameters != 0 ) delete theParameters;
      theParameters = new G4FTFParameters(theProjectile.GetDefinition(),
                                          aNucleus.GetA_asInt(),aNucleus.GetZ_asInt(),
                                          PlabPerParticle);
//G4cout<<" End Init "<<theProjectile.GetMomentum()<<G4endl;
// To turn on/off (1/0) elastic scattering close/open ...
//theParameters->SetProbabilityOfElasticScatt(0.); 
//G4cout<<" etProbabilityOfElasticScatt "<<theParameters->GetProbabilityOfElasticScatt()<<G4endl;
//G4cout<<" INIT ";
//G4int Uzhi; G4cin>>Uzhi;
 
   if(theAdditionalString.size() != 0)
   {
    std::for_each(theAdditionalString.begin(), theAdditionalString.end(), 
                  DeleteVSplitableHadron());
   }
   theAdditionalString.clear();
//G4cout<<" End Init theProjectile.GetMomentum()"<<theProjectile.GetMomentum()<<G4endl;
}

ModelDescription()

void G4FTFModel::ModelDescription ( std::ostream &  desc ) const

virtual

Reimplemented from G4VPartonStringModel.

Definition at line 1745 of file G4FTFModel.cc.

{
    desc << "please add description here" << G4endl;
}

---

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

• G4FTFModel.hh
• G4FTFModel.cc