G4InuclParamMomDst Class Reference

#include <G4InuclParamMomDst.hh>

Inheritance diagram for G4InuclParamMomDst:

Public Member Functions
	G4InuclParamMomDst (const G4String &name, const G4double(&pqprC)[2][4][4], const G4double(&psC)[2][3], G4int verbose=0)
virtual	~G4InuclParamMomDst ()
virtual G4double	GetMomentum (G4int ptype, const G4double &ekin) const
Public Member Functions inherited from G4VMultiBodyMomDst
	G4VMultiBodyMomDst (const G4String &name, G4int verbose=0)
virtual	~G4VMultiBodyMomDst ()
virtual void	setVerboseLevel (G4int verbose=0)
virtual const G4String &	GetName () const

Protected Attributes
const G4double(&	coeffPR )[2][4][4]
const G4double(&	coeffPS )[2][3]
Protected Attributes inherited from G4VMultiBodyMomDst
G4String	theName
G4int	verboseLevel

Detailed Description

Definition at line 42 of file G4InuclParamMomDst.hh.

Constructor & Destructor Documentation

G4InuclParamMomDst::G4InuclParamMomDst ( const G4String &  name, const G4double(&)  pqprC[2][4][4], const G4double(&)  psC[2][3], G4int  verbose = 0  )

inline

Definition at line 45 of file G4InuclParamMomDst.hh.

    : G4VMultiBodyMomDst(name, verbose), coeffPR(pqprC), coeffPS(psC) {;}

virtual G4InuclParamMomDst::~G4InuclParamMomDst ( )

inlinevirtual

Definition at line 51 of file G4InuclParamMomDst.hh.

{;}

Member Function Documentation

G4double G4InuclParamMomDst::GetMomentum ( G4int  ptype, const G4double &  ekin  ) const

virtual

Implements G4VMultiBodyMomDst.

Definition at line 49 of file G4InuclParamMomDst.cc.

References G4cout, G4endl, G4Pow::GetInstance(), G4InuclParticleNames::neu, G4Pow::powN(), G4InuclParticleNames::pro, and G4InuclSpecialFunctions::randomInuclPowers().

                                                                       {
  if (verboseLevel>3) {
    G4cout << theName << "::GetMomentum: ptype " << ptype << " ekin " << ekin
       << G4endl;
  }

  G4int JK = (ptype==pro || ptype==neu) ? 0 : 1;    // nucleon vs. other

  if (verboseLevel > 3) G4cout << " JK " << JK << G4endl;

  G4Pow* theG4Pow = G4Pow::GetInstance();   // For convenience

  G4double Spow = randomInuclPowers(ekin, coeffPR[JK]);

  G4double C=0., PS=0.;
  for(G4int im = 0; im < 3; im++) {
    C = coeffPS[JK][im];
    PS += C * theG4Pow->powN(ekin, im);

    if (verboseLevel >3) {
      G4cout << " im " << im << " : coeffPS[JK][im] " << C
         << " ekin^im " << theG4Pow->powN(ekin, im) << G4endl;
    }
  }
  
  G4double PRA = PS * Spow;

  if (verboseLevel > 3) 
    G4cout << " PS " << PS << " Spow = sqrt(S)*(PR+(1-PQ)*S^4) " << Spow
       << " PRA = PS*Spow " << PRA << G4endl;

  return std::fabs(PRA);
}

Field Documentation

const G4double(& G4InuclParamMomDst::coeffPR)[2][4][4]

protected

Definition at line 56 of file G4InuclParamMomDst.hh.

const G4double(& G4InuclParamMomDst::coeffPS)[2][3]

protected

Definition at line 57 of file G4InuclParamMomDst.hh.

---

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

• G4InuclParamMomDst.hh
• G4InuclParamMomDst.cc