#include <G4DalitzDecayChannel.hh>

Inheritance diagram for G4DalitzDecayChannel:

Public Member Functions
	G4DalitzDecayChannel (const G4String &theParentName, G4double theBR, const G4String &theLeptonName, const G4String &theAntiLeptonName)

virtual	~G4DalitzDecayChannel ()

virtual G4DecayProducts *	DecayIt (G4double)

Public Member Functions inherited from G4VDecayChannel
	G4VDecayChannel (const G4String &aName, G4int Verbose=1)

	G4VDecayChannel (const G4String &aName, const G4String &theParentName, G4double theBR, G4int theNumberOfDaughters, const G4String &theDaughterName1, const G4String &theDaughterName2="", const G4String &theDaughterName3="", const G4String &theDaughterName4="")

virtual	~G4VDecayChannel ()

G4int	operator== (const G4VDecayChannel &right) const

G4int	operator!= (const G4VDecayChannel &right) const

G4int	operator< (const G4VDecayChannel &right) const

const G4String &	GetKinematicsName () const

G4double	GetBR () const

G4int	GetNumberOfDaughters () const

G4ParticleDefinition *	GetParent ()

G4ParticleDefinition *	GetDaughter (G4int anIndex)

G4int	GetAngularMomentum ()

const G4String &	GetParentName () const

const G4String &	GetDaughterName (G4int anIndex) const

G4double	GetParentMass () const

G4double	GetDaughterMass (G4int anIndex) const

void	SetParent (const G4ParticleDefinition *particle_type)

void	SetParent (const G4String &particle_name)

void	SetBR (G4double value)

void	SetNumberOfDaughters (G4int value)

void	SetDaughter (G4int anIndex, const G4ParticleDefinition *particle_type)

void	SetDaughter (G4int anIndex, const G4String &particle_name)

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

void	DumpInfo ()

Protected Member Functions
	G4DalitzDecayChannel (const G4DalitzDecayChannel &)

G4DalitzDecayChannel &	operator= (const G4DalitzDecayChannel &)

Protected Member Functions inherited from G4VDecayChannel
void	ClearDaughtersName ()

void	FillDaughters ()

void	FillParent ()

	G4VDecayChannel ()

	G4VDecayChannel (const G4VDecayChannel &)

G4VDecayChannel &	operator= (const G4VDecayChannel &)

Additional Inherited Members
Protected Attributes inherited from G4VDecayChannel
G4String	kinematics_name

G4double	rbranch

G4int	numberOfDaughters

G4String *	parent_name

G4String **	daughters_name

G4ParticleTable *	particletable

G4int	verboseLevel

G4ParticleDefinition *	G4MT_parent

G4ParticleDefinition **	G4MT_daughters

G4double	G4MT_parent_mass

G4double *	G4MT_daughters_mass

Static Protected Attributes inherited from G4VDecayChannel
static const G4String	noName = " "

Detailed Description

Definition at line 43 of file G4DalitzDecayChannel.hh.

Constructor & Destructor Documentation

G4DalitzDecayChannel::G4DalitzDecayChannel	(	const G4String &	theParentName,
		G4double	theBR,
		const G4String &	theLeptonName,
		const G4String &	theAntiLeptonName
	)

Definition at line 53 of file G4DalitzDecayChannel.cc.

References G4DalitzDecayChannel(), G4VDecayChannel::SetBR(), G4VDecayChannel::SetDaughter(), G4VDecayChannel::SetNumberOfDaughters(), and G4VDecayChannel::SetParent().

Referenced by G4DalitzDecayChannel().

                    :G4VDecayChannel("Dalitz Decay",1)
 {
   // set names for daughter particles
   SetParent(theParentName);
   SetBR(theBR);
   SetNumberOfDaughters(3);
   G4String gammaName = "gamma";
   SetDaughter(idGamma, gammaName);
   SetDaughter(idLepton, theLeptonName);
   SetDaughter(idAntiLepton, theAntiLeptonName);
 }

G4DalitzDecayChannel::~G4DalitzDecayChannel ( )

virtual

Definition at line 70 of file G4DalitzDecayChannel.cc.

71 {

72 }

G4DalitzDecayChannel::G4DalitzDecayChannel ( const G4DalitzDecayChannel & right )

protected

Definition at line 74 of file G4DalitzDecayChannel.cc.

References G4DalitzDecayChannel().

                                                                            :
   G4VDecayChannel(right)
 {
 }

Member Function Documentation

G4DecayProducts * G4DalitzDecayChannel::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 106 of file G4DalitzDecayChannel.cc.

References CLHEP::HepLorentzVector::boost(), G4DecayProducts::DumpInfo(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::FillParent(), G4cout, G4endl, G4VDecayChannel::G4MT_daughters, G4VDecayChannel::G4MT_parent, G4UniformRand, G4DynamicParticle::Get4Momentum(), G4ParticleDefinition::GetPDGMass(), G4VDecayChannel::GetVerboseLevel(), G4PhaseSpaceDecayChannel::Pmx(), G4DecayProducts::PushProducts(), python.hepunit::rad, G4DynamicParticle::Set4Momentum(), python.hepunit::twopi, test::x, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

 {
 #ifdef G4VERBOSE
   if (GetVerboseLevel()>1) G4cout << "G4DalitzDecayChannel::DecayIt ";
 #endif 
   if (G4MT_parent == 0) FillParent();  
   if (G4MT_daughters == 0) FillDaughters();
 
   // parent mass
   G4double parentmass = G4MT_parent->GetPDGMass();
  
  //create parent G4DynamicParticle at rest
   G4ThreeVector dummy;
   G4DynamicParticle * parentparticle = new G4DynamicParticle( G4MT_parent, dummy, 0.0);
  
   //daughters'mass
   G4double leptonmass = G4MT_daughters[idLepton]->GetPDGMass(); 
 
  // Generate t ( = std::exp(x):mass Square of (l+ l-) system) 
   G4double xmin  = 2.0*std::log(2.0*leptonmass);
   G4double xmax  = 2.0*std::log(parentmass);
   G4double wmax = 1.5;
   G4double x, w, ww, w1, w2, w3, t;
   do {
     x = G4UniformRand()*(xmax-xmin) + xmin;
     w = G4UniformRand()*wmax;
     t = std::exp(x);
     w1 = (1.0-4.0*leptonmass*leptonmass/t);
     if ( w1 > 0.0) {
       w2 = ( 1.0 + 2.0*leptonmass*leptonmass/t );
       w3 = ( 1.0 - t/parentmass/parentmass );
       w3 = w3 * w3 * w3;
       ww = w3 * w2 * std::sqrt(w1);
     } else {
       ww = 0.0;
     }
   } while (w > ww);    
  
   // calculate gamma momentum
   G4double Pgamma = 
       G4PhaseSpaceDecayChannel::Pmx(parentmass, 0.0, std::sqrt(t)); 
   G4double costheta = 2.*G4UniformRand()-1.0;
   G4double sintheta = std::sqrt((1.0 - costheta)*(1.0 + costheta));
   G4double phi  = twopi*G4UniformRand()*rad;
   G4ThreeVector gdirection(sintheta*std::cos(phi),sintheta*std::sin(phi),costheta);
 
   //create G4DynamicParticle for gamma 
   G4DynamicParticle * gammaparticle
       = new G4DynamicParticle(G4MT_daughters[idGamma] , gdirection, Pgamma);
 
   // calcurate beta of (l+ l-)system
   G4double beta = Pgamma/(parentmass-Pgamma);
 
   // calculate lepton momentum in the rest frame of (l+ l-)system
   G4double Plepton = 
       G4PhaseSpaceDecayChannel::Pmx(std::sqrt(t),leptonmass, leptonmass); 
   G4double Elepton = std::sqrt(Plepton*Plepton + leptonmass*leptonmass );
   costheta = 2.*G4UniformRand()-1.0;
   sintheta = std::sqrt((1.0 - costheta)*(1.0 + costheta));
   phi  = twopi*G4UniformRand()*rad;
   G4ThreeVector ldirection(sintheta*std::cos(phi),sintheta*std::sin(phi),costheta);
   //create G4DynamicParticle for leptons  in the rest frame of (l+ l-)system
   G4DynamicParticle * leptonparticle 
     = new G4DynamicParticle(G4MT_daughters[idLepton] , 
                 ldirection, Elepton-leptonmass );
   G4DynamicParticle * antileptonparticle 
     = new G4DynamicParticle(G4MT_daughters[idAntiLepton] , 
                 -1.0*ldirection, Elepton-leptonmass );
   //boost leptons in the rest frame of the parent 
   G4LorentzVector p4 = leptonparticle->Get4Momentum();
   p4.boost( -1.0*gdirection.x()*beta, -1.0*gdirection.y()*beta, -1.0*gdirection.z()*beta);
   leptonparticle->Set4Momentum(p4);
   p4 = antileptonparticle->Get4Momentum();
   p4.boost( -1.0*gdirection.x()*beta, -1.0*gdirection.y()*beta, -1.0*gdirection.z()*beta);
   antileptonparticle->Set4Momentum(p4);
 
   //create G4Decayproducts
   G4DecayProducts *products = new G4DecayProducts(*parentparticle);
   delete parentparticle;
   products->PushProducts(gammaparticle);
   products->PushProducts(leptonparticle);
   products->PushProducts(antileptonparticle);
 
 #ifdef G4VERBOSE
   if (GetVerboseLevel()>1) {
      G4cout << "G4DalitzDecayChannel::DecayIt ";
      G4cout << "  create decay products in rest frame " <<G4endl;
      products->DumpInfo();
   }
 #endif
   return products;
 }

G4DalitzDecayChannel & G4DalitzDecayChannel::operator= ( const G4DalitzDecayChannel & right )

protected

Definition at line 79 of file G4DalitzDecayChannel.cc.

References G4VDecayChannel::ClearDaughtersName(), G4VDecayChannel::daughters_name, G4VDecayChannel::kinematics_name, G4VDecayChannel::numberOfDaughters, G4VDecayChannel::parent_name, G4VDecayChannel::rbranch, and G4VDecayChannel::verboseLevel.

 {
   if (this != &right) { 
     kinematics_name = right.kinematics_name;
     verboseLevel = right.verboseLevel;
     rbranch = right.rbranch;
 
     // copy parent name
     parent_name = new G4String(*right.parent_name);
 
     // clear daughters_name array
     ClearDaughtersName();
 
     // recreate array
     numberOfDaughters = right.numberOfDaughters;
     if ( numberOfDaughters >0 ) {
       if (daughters_name !=0) ClearDaughtersName();
       daughters_name = new G4String*[numberOfDaughters];
       //copy daughters name
       for (G4int index=0; index < numberOfDaughters; index++) {
           daughters_name[index] = new G4String(*right.daughters_name[index]);
       }
     }
   }
   return *this;
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation