G4Pythia6Decayer Class Reference

#include <G4Pythia6Decayer.hh>

Inheritance 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)
{
/// Standard constructor

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

G4Pythia6Decayer::~G4Pythia6Decayer ( )

virtual

Destructor

Definition at line 71 of file G4Pythia6Decayer.cc.

{
/// Destructor

  delete fDecayProductsArray;
}

Member Function Documentation

void G4Pythia6Decayer::ForceDecayType

(

EDecayType

decayType

)

Force a given decay type

Definition at line 618 of file G4Pythia6Decayer.cc.

Referenced by G4Pythia6DecayerMessenger::SetNewValue().

{ 
/// Force a given decay type

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

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

virtual

Import decay products

Implements G4VExtDecayer.

Definition at line 539 of file G4Pythia6Decayer.cc.

References Pythia6Particle::fKF, Pythia6Particle::fKS, G4cout, G4endl, G4Track::GetDefinition(), G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4Track::GetMomentum(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGEncoding(), G4DynamicParticle::GetTotalEnergy(), python.hepunit::GeV, G4DecayProducts::PushProducts(), status, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

{
/// Import decay products

  // 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;
}

void G4Pythia6Decayer::SetVerboseLevel ( G4int  verboseLevel )

inline

Definition at line 63 of file G4Pythia6Decayer.hh.

Referenced by P6DExtDecayerPhysics::ConstructProcess(), and G4Pythia6DecayerMessenger::SetNewValue().

{ fVerboseLevel =  verboseLevel; }

---

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

• G4Pythia6Decayer.hh
• G4Pythia6Decayer.cc