G4AngularDistributionPP Class Reference

#include <G4AngularDistributionPP.hh>

Inheritance diagram for G4AngularDistributionPP:

Collaboration 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 ()

Private Types
enum	{ NENERGY =40, NENERGYC =22, NANGLE =180 }

Static Private Attributes
static const G4float	sig [NENERGY][NANGLE]
static const G4float	pcm [NENERGY]
static const G4float	elab [NENERGY]
static const G4float	dsigmax [NENERGY]
static const G4float	sigtot [NENERGY]

Detailed Description

Definition at line 40 of file G4AngularDistributionPP.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
NENERGY
NENERGYC
NANGLE

Definition at line 45 of file G4AngularDistributionPP.hh.

{ NENERGY=40, NENERGYC=22, NANGLE=180 };

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.

{ }

Here is the call graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

Phi()

G4double G4AngularDistributionPP::Phi ( ) const

virtual

Reimplemented from G4VAngularDistribution.

Definition at line 133 of file G4AngularDistributionPP.cc.

{
    return twopi * G4UniformRand();
}

Here is the caller graph for this function:

Member Data Documentation

dsigmax

const G4float G4AngularDistributionPP::dsigmax[NENERGY]

staticprivate

Definition at line 61 of file G4AngularDistributionPP.hh.

elab

const G4float G4AngularDistributionPP::elab[NENERGY]

staticprivate

Definition at line 61 of file G4AngularDistributionPP.hh.

pcm

const G4float G4AngularDistributionPP::pcm[NENERGY]

staticprivate

Definition at line 61 of file G4AngularDistributionPP.hh.

sig

const G4float G4AngularDistributionPP::sig[NENERGY][NANGLE]

staticprivate

Definition at line 60 of file G4AngularDistributionPP.hh.

sigtot

const G4float G4AngularDistributionPP::sigtot[NENERGY]

staticprivate

Definition at line 61 of file G4AngularDistributionPP.hh.

---

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

• G4AngularDistributionPP.hh
• G4AngularDistributionPP.cc