G4AngularDistributionPP Class Reference

#include <G4AngularDistributionPP.hh>

Inheritance diagram for G4AngularDistributionPP:

Public Member Functions
	G4AngularDistributionPP ()
virtual	~G4AngularDistributionPP ()
virtual G4double	CosTheta (G4double s, G4double m1, G4double m2) const
virtual G4double	Phi () const
Public Member Functions inherited from G4VAngularDistribution
	G4VAngularDistribution ()
virtual	~G4VAngularDistribution ()
virtual G4double	CosTheta (G4double s, G4double m1, G4double m2) const =0
virtual G4double	Phi () const

Detailed Description

Definition at line 40 of file G4AngularDistributionPP.hh.

Constructor & Destructor Documentation

G4AngularDistributionPP()

G4AngularDistributionPP::G4AngularDistributionPP ( )

inline

Definition at line 49 of file G4AngularDistributionPP.hh.

{ }

~G4AngularDistributionPP()

virtual G4AngularDistributionPP::~G4AngularDistributionPP ( )

inlinevirtual

Definition at line 51 of file G4AngularDistributionPP.hh.

{ }

Member Function Documentation

CosTheta()

G4double G4AngularDistributionPP::CosTheta ( G4double  s, G4double  m1, G4double  m2  ) const

virtual

Implements G4VAngularDistribution.

Definition at line 37 of file G4AngularDistributionPP.cc.

{
    G4int verboseLevel=1;
 
    G4double ek= ((S - sqr(m_1) -sqr(m_2) )/(2*m_1) - m_1 )/GeV   ;    // kinetic energy in GeV
 
    // Find energy bin
 
    G4int je1 = 0;
    G4int je2 = NENERGY - 1;
     G4int iterationsLeft=2*NENERGY +1;
    do {
      G4int midBin = (je1 + je2)/2;
      if (ek < elab[midBin])
        je2 = midBin;
      else
        je1 = midBin;
    } while ( (je2 - je1) > 1 && --iterationsLeft > 0 );    /* Loop checking, 30-Oct-2015, G.Folger */
     if ( iterationsLeft <= 0 ) {
         G4Exception("G4AngularDistributionPP", "im_r_matrix020", FatalException,
             "Problem with energy bin (elab) data");
 
     }
    //    G4int j;
    //std::abs(ek-elab[je1]) < std::abs(ek-elab[je2]) ? j = je1 : j = je2;
    G4double delab = elab[je2] - elab[je1];
 
    // Sample the angle
 
    G4double sample = G4UniformRand();
    G4int ke1 = 0;
    G4int ke2 = NANGLE - 1;
    G4double dsig = sig[je2][0] - sig[je1][0];
    G4double rc = dsig/delab;
    G4double b = sig[je1][0] - rc*elab[je1];
    G4double sigint1 = rc*ek + b;
    G4double sigint2 = 0.;
 
    if (verboseLevel > 1) G4cout << "sample=" << sample << G4endl
                                 << ek << " " << ke1 << " " << ke2 << " "
                                 << sigint1 << " " << sigint2 << G4endl;
 
    iterationsLeft= 2*NANGLE +1;
     do {
      G4int midBin = (ke1 + ke2)/2;
      dsig = sig[je2][midBin] - sig[je1][midBin];
      rc = dsig/delab;
      b = sig[je1][midBin] - rc*elab[je1];
      G4double sigint = rc*ek + b;
      if (sample < sigint) {
        ke2 = midBin;
        sigint2 = sigint;
      }
      else {
        ke1 = midBin;
        sigint1 = sigint;
      }
      if (verboseLevel > 1)G4cout << ke1 << " " << ke2 << " "
                                  << sigint1 << " " << sigint2 << G4endl;
    } while ( (ke2 - ke1) > 1 && --iterationsLeft > 0);    /* Loop checking, 30-Oct-2015, G.Folger */
     if ( iterationsLeft <= 0 ) {
         G4Exception("G4AngularDistributionPP", "im_r_matrix021", FatalException,
             "Problem with angular distribution (sig) data");
     }
 
    // sigint1 and sigint2 should be recoverable from above loop
 
    //    G4double dsig = sig[je2][ke1] - sig[je1][ke1];
    //    G4double rc = dsig/delab;
    //    G4double b = sig[je1][ke1] - rc*elab[je1];
    //    G4double sigint1 = rc*ek + b;
 
    //    G4double dsig = sig[je2][ke2] - sig[je1][ke2];
    //    G4double rc = dsig/delab;
    //    G4double b = sig[je1][ke2] - rc*elab[je1];
    //    G4double sigint2 = rc*ek + b;
 
    dsig = sigint2 - sigint1;
    rc = 1./dsig;
    b = ke1 - rc*sigint1;
    G4double kint = rc*sample + b;
    G4double theta = (0.5 + kint)*pi/180.;
 
    //    G4int k;
    //std::abs(sample-sig[j][ke1]) < std::abs(sample-sig[j][ke2]) ? k = ke1 : k = ke2;
    //    G4double theta = (0.5 + k)*pi/180.;
 
    if (verboseLevel > 1) {
      G4cout << "   energy bin " << je1 << " energy=" << elab[je1] << G4endl;
      G4cout << "   angle bin " << kint << " angle=" << theta/degree << G4endl;
    }
    G4double costh= std::cos(theta);
//    G4cout << "scattering angle hpw "<<costh<<" "<<s<<" "<<m1<<" "<<m2<<G4endl;
    return costh;
}

Phi()

G4double G4AngularDistributionPP::Phi ( ) const

virtual

Reimplemented from G4VAngularDistribution.

Definition at line 133 of file G4AngularDistributionPP.cc.

{
    return twopi * G4UniformRand();
}

---

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

• G4AngularDistributionPP.hh
• G4AngularDistributionPPData.hh
• G4AngularDistributionPP.cc