G4Pythia6Decayer Class Reference

#include <G4Pythia6Decayer.hh>

Inheritance diagram for G4Pythia6Decayer:

Collaboration diagram for G4Pythia6Decayer:

Public Member Functions
	G4Pythia6Decayer ()
virtual	~G4Pythia6Decayer ()
virtual G4DecayProducts *	ImportDecayProducts (const G4Track &track)
void	ForceDecayType (EDecayType decayType)
void	SetVerboseLevel (G4int verboseLevel)
Public Member Functions inherited from G4VExtDecayer
	G4VExtDecayer (const G4String &name="")
virtual	~G4VExtDecayer ()
const G4String &	GetName () const

Additional Inherited Members
Protected Attributes inherited from G4VExtDecayer
G4String	decayerName

Detailed Description

Pythia6 decayer

Implements the G4VExtDecayer abstract class using the Pythia6 interface. According to TPythia6Decayer class in Root: http://root.cern.ch/ see http://root.cern.ch/root/License.html

Definition at line 53 of file G4Pythia6Decayer.hh.

Constructor & Destructor Documentation

G4Pythia6Decayer::G4Pythia6Decayer

(

)

Standard constructor

Definition at line 55 of file G4Pythia6Decayer.cc.

  : G4VExtDecayer("G4Pythia6Decayer"),
    fMessenger(this),
    fVerboseLevel(0),
    fDecayType(fgkDefaultDecayType),
    fDecayProductsArray(0)
{

  fDecayProductsArray = new ParticleVector();
  
  ForceDecay(fDecayType);
}

G4Pythia6Decayer::~G4Pythia6Decayer ( )

virtual

Destructor

Definition at line 71 of file G4Pythia6Decayer.cc.

{

  delete fDecayProductsArray;
}

Member Function Documentation

void G4Pythia6Decayer::ForceDecayType

(

EDecayType

decayType

)

Force a given decay type

Definition at line 617 of file G4Pythia6Decayer.cc.

{ 

  // Do nothing if the decay type is not different from current one
  if ( decayType == fDecayType ) return;
  
  fDecayType =  decayType; 
  ForceDecay(fDecayType);
}

Here is the caller graph for this function:

G4DecayProducts * G4Pythia6Decayer::ImportDecayProducts ( const G4Track &  track )

virtual

Import decay products

Implements G4VExtDecayer.

Definition at line 538 of file G4Pythia6Decayer.cc.

{

  // get particle momentum
  G4ThreeVector momentum = track.GetMomentum(); 
  G4double etot = track.GetDynamicParticle()->GetTotalEnergy();;  
  CLHEP::HepLorentzVector p;    
  p[0] = momentum.x() / GeV;
  p[1] = momentum.y() / GeV;
  p[2] = momentum.z() / GeV;
  p[3] = etot         / GeV;
  
  // get particle PDG
  // ask G4Pythia6Decayer to get PDG encoding 
  // (in order to get PDG from extended TDatabasePDG
  // in case the standard PDG code is not defined)
  G4ParticleDefinition* particleDef = track.GetDefinition();
  G4int pdgEncoding = particleDef->GetPDGEncoding();

  // let Pythia6Decayer decay the particle
  // and import the decay products
  Decay(pdgEncoding, p);
  G4int nofParticles = ImportParticles(fDecayProductsArray);
  
  if ( fVerboseLevel > 0 ) {
    G4cout << "nofParticles: " <<  nofParticles << G4endl;
  }  

  // convert decay products Pythia6Particle type 
  // to G4DecayProducts  
  G4DecayProducts* decayProducts
    = new G4DecayProducts(*(track.GetDynamicParticle()));

  G4int counter = 0;
  for (G4int i=0; i<nofParticles; i++) {

    // get particle from ParticleVector
    Pythia6Particle* particle = (*fDecayProductsArray)[i];
      
    G4int status = particle->fKS;
    G4int pdg = particle->fKF;
    if ( status>0 && status<11 && 
         std::abs(pdg)!=12 && std::abs(pdg)!=14 && std::abs(pdg)!=16 ) {
      // pass to tracking final particles only;
      // skip neutrinos

      if ( fVerboseLevel > 0 ) {
        G4cout << "  " << i << "th particle PDG: " << pdg << "   ";
      }  
            
      // create G4DynamicParticle 
      G4DynamicParticle* dynamicParticle 
        = CreateDynamicParticle(particle);

      if (dynamicParticle) {

        if ( fVerboseLevel > 0 ) {
          G4cout << "  G4 particle name: " 
                 << dynamicParticle->GetDefinition()->GetParticleName()
                 << G4endl;
        }         

        // add dynamicParticle to decayProducts
        decayProducts->PushProducts(dynamicParticle);
        
        counter++;
      }
    }       
  }                             
  if ( fVerboseLevel > 0 ) {
    G4cout << "nofParticles for tracking: " <<  counter << G4endl;
  }  
     
  return decayProducts;
}

Here is the call graph for this function:

void G4Pythia6Decayer::SetVerboseLevel ( G4int  verboseLevel )

inline

Definition at line 63 of file G4Pythia6Decayer.hh.

{ fVerboseLevel =  verboseLevel; }

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/examples/extended/eventgenerator/pythia/decayer6/include/G4Pythia6Decayer.hh
• source/geant4.10.03.p02/examples/extended/eventgenerator/pythia/decayer6/src/G4Pythia6Decayer.cc