#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 G4double	GetMomentum (G4int ptype, const G4double &ekin) const =0

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 41 of file G4InuclParamMomDst.hh.

Constructor & Destructor Documentation

◆ G4InuclParamMomDst()

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

inline

Definition at line 44 of file G4InuclParamMomDst.hh.

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

G4InuclParamMomDst::coeffPS

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

Definition: G4InuclParamMomDst.hh:56

G4InuclParamMomDst::coeffPR

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

Definition: G4InuclParamMomDst.hh:55

G4VMultiBodyMomDst

Definition: G4VMultiBodyMomDst.hh:40

◆ ~G4InuclParamMomDst()

virtual G4InuclParamMomDst::~G4InuclParamMomDst ( )

inlinevirtual

Definition at line 50 of file G4InuclParamMomDst.hh.

50{;}

Member Function Documentation

◆ GetMomentum()

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

virtual

Implements G4VMultiBodyMomDst.

Definition at line 48 of file G4InuclParamMomDst.cc.

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

Member Data Documentation

◆ coeffPR

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

protected

Definition at line 55 of file G4InuclParamMomDst.hh.

Referenced by GetMomentum().

◆ coeffPS

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