#include <G4InuclParamMomDst.hh>

Inheritance diagram for G4InuclParamMomDst:

Collaboration 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::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.

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

G4VMultiBodyMomDst::G4VMultiBodyMomDst

G4VMultiBodyMomDst(const G4String &name, G4int verbose=0)

Definition: G4VMultiBodyMomDst.cc:35

G4InuclParamMomDst::coeffPR

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

Definition: G4InuclParamMomDst.hh:56

G4InuclParamMomDst::coeffPS

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

Definition: G4InuclParamMomDst.hh:57

◆ ~G4InuclParamMomDst()

virtual G4InuclParamMomDst::~G4InuclParamMomDst ( )

inlinevirtual

Definition at line 51 of file G4InuclParamMomDst.hh.

51 {;}

Here is the call graph for this function:

Member Function Documentation

◆ GetMomentum()

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

virtual

Implements G4VMultiBodyMomDst.

Definition at line 49 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);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ coeffPR

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

protected

Definition at line 56 of file G4InuclParamMomDst.hh.

◆ coeffPS

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