G4DeltaAngle Class Reference

#include <G4DeltaAngle.hh>

Inheritance diagram for G4DeltaAngle:

Public Member Functions
	G4DeltaAngle (const G4String &name="")
virtual	~G4DeltaAngle ()
virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, const G4Material *mat=0)
void	PrintGeneratorInformation () const
Public Member Functions inherited from G4VEmAngularDistribution
	G4VEmAngularDistribution (const G4String &name)
virtual	~G4VEmAngularDistribution ()
virtual G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
const G4String &	GetName () const

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection

Detailed Description

Definition at line 57 of file G4DeltaAngle.hh.

Constructor & Destructor Documentation

G4DeltaAngle::G4DeltaAngle

(

const G4String &

name = ""

)

Definition at line 60 of file G4DeltaAngle.cc.

References G4Electron::Electron().

  : G4VEmAngularDistribution("deltaVI")
{
  fElectron = G4Electron::Electron();
  nprob = 26;
  prob.resize(nprob,0.0);
}    

G4DeltaAngle::~G4DeltaAngle ( )

virtual

Definition at line 68 of file G4DeltaAngle.cc.

{}

Member Function Documentation

void G4DeltaAngle::PrintGeneratorInformation

(

)

const

Definition at line 154 of file G4DeltaAngle.cc.

{} 

G4ThreeVector & G4DeltaAngle::SampleDirection ( const G4DynamicParticle *  dp, G4double  kinEnergyFinal, G4int  Z, const G4Material *  mat = 0  )

virtual

Implements G4VEmAngularDistribution.

Definition at line 72 of file G4DeltaAngle.cc.

References CLHEP::HepLorentzVector::beta(), G4InuclSpecialFunctions::bindingEnergy(), CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), python.hepunit::electron_mass_c2, G4VEmAngularDistribution::fLocalDirection, G4Log(), G4UniformRand, G4InuclParticleNames::gam, CLHEP::HepLorentzVector::gamma(), G4AtomicShells::GetBindingEnergy(), G4DynamicParticle::GetMomentumDirection(), G4AtomicShells::GetNumberOfElectrons(), G4AtomicShells::GetNumberOfShells(), G4DynamicParticle::GetParticleDefinition(), G4ParticleDefinition::GetPDGMass(), G4DynamicParticle::GetTotalEnergy(), G4DynamicParticle::GetTotalMomentum(), CLHEP::HepLorentzVector::mag(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), python.hepunit::twopi, CLHEP::Hep3Vector::unit(), test::x, CLHEP::HepLorentzVector::x(), CLHEP::HepLorentzVector::y(), and CLHEP::HepLorentzVector::z().

{
  G4int nShells = G4AtomicShells::GetNumberOfShells(Z);
  if(nShells> nprob) {
    nprob = nShells;
    prob.resize(nprob,0.0);
  }
  G4int idx;
  G4double sum = 0.0;
  for(idx=0; idx<nShells; ++idx) {
    sum += G4AtomicShells::GetNumberOfElectrons(Z, idx)
      /G4AtomicShells::GetBindingEnergy(Z, idx);
    prob[idx] = sum;
  }
  sum *= G4UniformRand();
  for(idx=0; idx<nShells; ++idx) {
    if(sum <= prob[idx]) { break; }
  }
  G4double bindingEnergy = G4AtomicShells::GetBindingEnergy(Z, idx);
  G4double mass = dp->GetParticleDefinition()->GetPDGMass();

  G4ThreeVector bst(0.0,0.0,0.0);
  G4double cost, en, mom;

  do {
  
    // the atomic electron
    G4double x = -G4Log(G4UniformRand());
    G4double eKinEnergy = bindingEnergy*x;
    G4double ePotEnergy = bindingEnergy*(1.0 + x);
    G4double e = kinEnergyFinal + ePotEnergy + electron_mass_c2;

    G4double totEnergy = dp->GetTotalEnergy();
    G4double totMomentum = dp->GetTotalMomentum();
    if(dp->GetParticleDefinition() == fElectron) {
      totEnergy += ePotEnergy;
      totMomentum = sqrt((totEnergy + electron_mass_c2)
             *(totEnergy - electron_mass_c2));
    }
 
    G4double eTotMomentum = sqrt(eKinEnergy*(eKinEnergy + 2*electron_mass_c2));
    G4double phi = G4UniformRand()*twopi;
    G4double costet = 2*G4UniformRand() - 1;
    G4double sintet = sqrt((1 - costet)*(1 + costet));
 
    G4LorentzVector lv0(eTotMomentum*sintet*cos(phi),
            eTotMomentum*sintet*sin(phi),
            eTotMomentum*costet + totMomentum,
            eKinEnergy + electron_mass_c2 + totEnergy);
    bst = lv0.boostVector();

    G4double m0  = lv0.mag();
    G4double bet = lv0.beta();
    G4double gam = lv0.gamma();

    en = 0.5*(m0*m0 - mass*mass + electron_mass_c2*electron_mass_c2)/m0;
    mom = sqrt((en + electron_mass_c2)*(en - electron_mass_c2));

    cost= (e/gam - en)/(mom*bet);

    //G4cout << "e= " << e << " gam= " << gam << " en= " << en 
    //     << " mom= " << mom << "  beta= " << bet << " cost= " << cost 
    //     << G4endl; 

  } while(std::fabs(cost) > 1.0);

  G4double sint = sqrt((1 - cost)*(1 + cost));
  G4double phi  = twopi*G4UniformRand(); 

  G4LorentzVector lv1(sint*std::cos(phi)*mom, sint*std::sin(phi)*mom,
              mom*cost, en);
  lv1.boost(bst);

  fLocalDirection.set(lv1.x(), lv1.y(), lv1.z());
  fLocalDirection = fLocalDirection.unit();
  fLocalDirection.rotateUz(dp->GetMomentumDirection());

  return fLocalDirection;
}

---

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

• G4DeltaAngle.hh
• G4DeltaAngle.cc